V2 JOE'S HACKATHON MODEL STACKED.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "import sklearn\n",
    "from sklearn.tree import DecisionTreeRegressor\n",
    "from sklearn.metrics import mean_squared_error, r2_score, mean_absolute_error\n",
    "from sklearn.ensemble import AdaBoostRegressor, BaggingRegressor, RandomForestRegressor, ExtraTreesRegressor, StackingRegressor, GradientBoostingRegressor\n",
    "from sklearn.svm import SVR\n",
    "from sklearn.model_selection import KFold, cross_val_score, RandomizedSearchCV, GridSearchCV, train_test_split, cross_val_predict\n",
    "from sklearn.linear_model import Lasso, Ridge, LinearRegression\n",
    "from sklearn.preprocessing import StandardScaler\n",
    "from sklearn.datasets import make_regression\n",
    "from sklearn.neighbors import KNeighborsRegressor\n",
    "import xgboost as xgb\n",
    "%matplotlib inline\n",
    "pd.set_option('display.max_columns', None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "df1 = pd.read_csv(\"frequency_domain_features_train.csv\")\n",
    "df2 = pd.read_csv(\"heart_rate_non_linear_features_train.csv\")\n",
    "df3 = pd.read_csv(\"time_domain_features_train.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((369289, 12), (369289, 7), (369289, 20))"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df1.shape, df2.shape, df3.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "df4 = df1.merge(df2,how='inner',on='uuid')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(369289, 18)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df4.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "df5 = df4.merge(df3,how='inner',on='uuid')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(369289, 37)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df5.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# RANDOM SAMPLE OF DATAFRAME\n",
    "#df5 = df5.sample(50000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "3"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df5.condition.nunique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7fce7ef41c40>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEGCAYAAABYV4NmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAetElEQVR4nO3dfZRVV5nn8e8voIgaEELpIC+CCdEmaSVDDYkdtWNjA+10S7TBVJZtSMsMJpPYndHpaaMzkpUsZiUdNWO0kx5skJcVkyBJDE5LJyyIyWjCS5FgeInplAmaEiZBoRFfwljkmT/OvnLqcqu4kNr3QvH7rHVWnfOcvffd59669dQ++9xzFRGYmZn1tdOa3QEzM+ufnGDMzCwLJxgzM8vCCcbMzLJwgjEzsywGNrsDJ4oRI0bEuHHjmt0NM7OTyubNm38WES219jnBJOPGjaO9vb3Z3TAzO6lI+nFP+3yKzMzMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MssiUYSWMkPSTpKUnbJf11ig+XtEbSM+nnsFKdayV1SHpa0vRSfLKkrWnfrZKU4oMk3Z3iGySNK9WZkx7jGUlzch2nmZnVlnME0wV8OiJ+D7gAuErSROAzwNqImACsTdukfW3AOcAM4DZJA1JbtwPzgAlpmZHic4F9EXEWcAtwU2prODAfOB+YAswvJzIzM8svW4KJiN0R8XhaPwA8BYwCZgJLU7GlwMVpfSZwV0QcjIjngA5giqSRwJCIeCyKL69ZVlWn0tZKYGoa3UwH1kTE3ojYB6zhcFIyM7MGaMgn+dOpq/OADcCbImI3FElI0htTsVHA+lK1zhT7bVqvjlfqPJ/a6pK0HzijHK9Rp9yveRQjI8aOHXtMxzT5b5YdU3k7dptvvqzZXTCzVyD7JL+k1wP3ANdExC96K1ojFr3Ej7fO4UDEwohojYjWlpaat9IxM7PjlDXBSHoVRXK5IyLuTeEX0mkv0s8XU7wTGFOqPhrYleKja8S71ZE0EBgK7O2lLTMza5CcV5EJWAQ8FRFfKu1aBVSu6poD3F+Kt6Urw8ZTTOZvTKfTDki6ILV5WVWdSluzgHVpnuYBYJqkYWlyf1qKmZlZg+Scg7kQ+BiwVdKWFPsscCOwQtJc4CfAbICI2C5pBbCD4gq0qyLiUKp3JbAEGAysTgsUCWy5pA6KkUtbamuvpBuATanc9RGxN9eBmpnZkbIlmIj4HrXnQgCm9lBnAbCgRrwdOLdG/CVSgqqxbzGwuN7+mplZ3/In+c3MLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8vCCcbMzLJwgjEzsyycYMzMLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsi2wJRtJiSS9K2laK3S1pS1p2Vr5KWdI4Sb8p7fuHUp3JkrZK6pB0qySl+KDUXoekDZLGlerMkfRMWubkOkYzM+tZtq9MBpYAXwWWVQIRcUllXdIXgf2l8j+KiEk12rkdmAesB74DzABWA3OBfRFxlqQ24CbgEknDgflAKxDAZkmrImJfHx6bmZkdRbYRTEQ8AuyttS+NQj4C3NlbG5JGAkMi4rGICIpkdXHaPRNYmtZXAlNTu9OBNRGxNyWVNRRJyczMGqhZczDvAV6IiGdKsfGSnpD0sKT3pNgooLNUpjPFKvueB4iILorR0BnleI063UiaJ6ldUvuePXte6TGZmVlJsxLMpXQfvewGxkbEecCngG9IGgKoRt1IP3va11ud7sGIhRHRGhGtLS0tdXfezMyOruEJRtJA4MPA3ZVYRByMiJ+n9c3Aj4CzKUYfo0vVRwO70nonMKbU5lCKU3K/i9eoY2ZmDdKMEcz7gR9GxO9OfUlqkTQgrb8VmAA8GxG7gQOSLkjzK5cB96dqq4DKFWKzgHVpnuYBYJqkYZKGAdNSzMzMGijbVWSS7gQuAkZI6gTmR8QioI0jJ/ffC1wvqQs4BFwREZULBK6kuCJtMMXVY6tTfBGwXFIHxcilDSAi9kq6AdiUyl1fasvMzBokW4KJiEt7iF9eI3YPcE8P5duBc2vEXwJm91BnMbD4GLprZmZ9zJ/kNzOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8vCCcbMzLJwgjEzsyycYMzMLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyyyJZgJC2W9KKkbaXYdZJ+KmlLWj5Q2netpA5JT0uaXopPlrQ17bs1fXUykgZJujvFN0gaV6ozR9Izaal8rbKZmTVQzhHMEmBGjfgtETEpLd8BkDSR4iuPz0l1bpM0IJW/HZgHTEhLpc25wL6IOAu4BbgptTUcmA+cD0wB5ksa1veHZ2ZmvcmWYCLiEWBvncVnAndFxMGIeA7oAKZIGgkMiYjHIiKAZcDFpTpL0/pKYGoa3UwH1kTE3ojYB6yhdqIzM7OMmjEHc7WkJ9MptMrIYhTwfKlMZ4qNSuvV8W51IqIL2A+c0UtbZmbWQI1OMLcDZwKTgN3AF1NcNcpGL/HjrdONpHmS2iW179mzp7d+m5nZMWpogomIFyLiUES8DHyNYo4EilHGmFLR0cCuFB9dI96tjqSBwFCKU3I9tVWrPwsjojUiWltaWl7JoZmZWZWGJpg0p1LxIaByhdkqoC1dGTaeYjJ/Y0TsBg5IuiDNr1wG3F+qU7lCbBawLs3TPABMkzQsnYKblmJmZtZAA3M1LOlO4CJghKROiiu7LpI0ieKU1U7gEwARsV3SCmAH0AVcFRGHUlNXUlyRNhhYnRaARcBySR0UI5e21NZeSTcAm1K56yOi3osNzMysj2RLMBFxaY3wol7KLwAW1Ii3A+fWiL8EzO6hrcXA4ro7a2Zmfc6f5DczsyycYMzMLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8vCCcbMzLJwgjEzsyycYMzMLAsnGDMzy8IJxszMsnCCMTOzLLIlGEmLJb0oaVspdrOkH0p6UtJ9kt6Q4uMk/UbSlrT8Q6nOZElbJXVIulWSUnyQpLtTfIOkcaU6cyQ9k5Y5uY7RzMx6lnMEswSYURVbA5wbEe8A/gW4trTvRxExKS1XlOK3A/OACWmptDkX2BcRZwG3ADcBSBoOzAfOB6YA8yUN68sDMzOzo8uWYCLiEWBvVezBiOhKm+uB0b21IWkkMCQiHouIAJYBF6fdM4GlaX0lMDWNbqYDayJib0Tso0hq1YnOzMwya+YczMeB1aXt8ZKekPSwpPek2Cigs1SmM8Uq+54HSElrP3BGOV6jTjeS5klql9S+Z8+eV3o8ZmZW0pQEI+lzQBdwRwrtBsZGxHnAp4BvSBoCqEb1qDTTw77e6nQPRiyMiNaIaG1paTmWQzAzs6NoeIJJk+5/Cnw0nfYiIg5GxM/T+mbgR8DZFKOP8mm00cCutN4JjEltDgSGUpyS+128Rh0zM2uQhiYYSTOAvwU+GBG/LsVbJA1I62+lmMx/NiJ2AwckXZDmVy4D7k/VVgGVK8RmAetSwnoAmCZpWJrcn5ZiZmbWQANzNSzpTuAiYISkTooru64FBgFr0tXG69MVY+8FrpfUBRwCroiIygUCV1JckTaYYs6mMm+zCFguqYNi5NIGEBF7Jd0AbErlri+1ZWZmDZItwUTEpTXCi3ooew9wTw/72oFza8RfAmb3UGcxsLjuzpqZWZ/zJ/nNzCwLJxgzM8vCCcbMzLJwgjEzsyzqSjCS1tYTMzMzq+j1KjJJrwFeS3Gp8TAOf0p+CPDmzH0zM7OT2NEuU/4EcA1FMtnM4QTzC+DvM/bLzMxOcr0mmIj4MvBlSZ+MiK80qE9mZtYP1PVBy4j4iqQ/AMaV60TEskz9MjOzk1xdCUbScuBMYAvFrVyguEOxE4yZmdVU761iWoGJlbsfm5mZHU29n4PZBvybnB0xM7P+pd4RzAhgh6SNwMFKMCI+mKVXZmZ20qs3wVyXsxNmZtb/1HsV2cO5O2JmZv1LvVeRHeDw99q/GngV8KuIGJKrY2a5/OT63292F/q9sZ/f2uwu2Amg3hHM6eVtSRcDU7L0yMzM+oXjuptyRHwL+KPeykhaLOlFSdtKseGS1kh6Jv0cVtp3raQOSU9Lml6KT5a0Ne27Vem7liUNknR3im+QNK5UZ056jGckzTmeYzQzs1em3rspf7i0zJJ0I4dPmfVkCTCjKvYZYG1ETADWpm0kTQTagHNSndskDUh1bgfmARPSUmlzLrAvIs4CbgFuSm0NB+YD51OMsuaXE5mZmTVGvSOYPyst04EDwMzeKkTEI8DeqvBMYGlaXwpcXIrfFREHI+I5oAOYImkkMCQiHksf8lxWVafS1kpgahrdTAfWRMTeiNgHrOHIRGdmZpnVOwfzl330eG+KiN2pzd2S3pjio4D1pXKdKfbbtF4dr9R5PrXVJWk/cEY5XqNON5LmUYyOGDt27PEflZmZHaHeU2SjJd2X5lRekHSPpNF92A/ViEUv8eOt0z0YsTAiWiOitaWlpa6OmplZfeo9RfZ1YBXF98KMAr6dYsfqhXTai/TzxRTvBMaUyo0GdqX46BrxbnUkDQSGUpyS66ktMzNroHoTTEtEfD0iutKyBDief/lXAZWruuYA95fibenKsPEUk/kb0+m0A5IuSPMrl1XVqbQ1C1iX5mkeAKZJGpYm96elmJmZNVC9t4r5maS/AO5M25cCP++tgqQ7gYsovm65k+LKrhuBFZLmAj8BZgNExHZJK4AdQBdwVURUvhbgSoor0gYDq9MCsAhYLqmDYuTSltraK+kGYFMqd31EVF9sYGZmmdWbYD4OfJXicuAAHgV6nfiPiEt72DW1h/ILgAU14u3AuTXiL5ESVI19i4HFvfXPzMzyqjfB3ADMSZf9Vj5r8gWKxGNmZnaEeudg3lFJLlCchgLOy9MlMzPrD+pNMKdV3dZlOPWPfszM7BRUb5L4IvCopJUUczAfocZ8iZmZWUW9n+RfJqmd4gaXAj4cETuy9szMzE5qdZ/mSgnFScXMzOpyXLfrNzMzOxonGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8vCCcbMzLJoeIKR9DZJW0rLLyRdI+k6ST8txT9QqnOtpA5JT0uaXopPlrQ17btVklJ8kKS7U3yDpHGNPk4zs1NdwxNMRDwdEZMiYhIwGfg1cF/afUtlX0R8B0DSRKANOAeYAdwmaUAqfzswD5iQlhkpPhfYFxFnUXzN800NODQzMytp9imyqcCPIuLHvZSZCdwVEQcj4jmgA5giaSQwJCIei4gAlgEXl+osTesrgamV0Y2ZmTVGsxNMG3BnaftqSU9KWlz6Bs1RwPOlMp0pNiqtV8e71YmILmA/cEb1g0uaJ6ldUvuePXv64njMzCxpWoKR9Grgg8A3U+h24ExgErCb4ls0ofiCs2rRS7y3Ot0DEQsjojUiWltaWo6h92ZmdjTNHMH8CfB4RLwAEBEvRMShiHgZ+BowJZXrBMaU6o0GdqX46BrxbnUkDQSGAnszHYeZmdXQzARzKaXTY2lOpeJDwLa0vgpoS1eGjaeYzN8YEbuBA5IuSPMrlwH3l+rMSeuzgHVpnsbMzBqk7q9M7kuSXgv8MfCJUvjvJE2iOJW1s7IvIrZLWkHxdc1dwFURcSjVuRJYAgwGVqcFYBGwXFIHxcilLefxmJnZkZqSYCLi11RNukfEx3opvwBYUCPeDpxbI/4SMPuV99TMzI5Xs68iMzOzfsoJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyyaMplymZmx+vCr1zY7C70e9//5Pf7pB2PYMzMLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8uiKQlG0k5JWyVtkdSeYsMlrZH0TPo5rFT+Wkkdkp6WNL0Un5za6ZB0qySl+CBJd6f4BknjGn2MZmanumaOYN4XEZMiojVtfwZYGxETgLVpG0kTgTbgHGAGcJukAanO7cA8YEJaZqT4XGBfRJwF3ALc1IDjMTOzkhPpFNlMYGlaXwpcXIrfFREHI+I5oAOYImkkMCQiHouIAJZV1am0tRKYWhndmJlZYzQrwQTwoKTNkual2JsiYjdA+vnGFB8FPF+q25lio9J6dbxbnYjoAvYDZ1R3QtI8Se2S2vfs2dMnB2ZmZoVm3a7/wojYJemNwBpJP+ylbK2RR/QS761O90DEQmAhQGtr6xH7zczs+DVlBBMRu9LPF4H7gCnAC+m0F+nni6l4JzCmVH00sCvFR9eId6sjaSAwFNib41jMzKy2hicYSa+TdHplHZgGbANWAXNSsTnA/Wl9FdCWrgwbTzGZvzGdRjsg6YI0v3JZVZ1KW7OAdWmexszMGqQZp8jeBNyX5twHAt+IiH+WtAlYIWku8BNgNkBEbJe0AtgBdAFXRcSh1NaVwBJgMLA6LQCLgOWSOihGLm2NODAzMzus4QkmIp4F3lkj/nNgag91FgALasTbgXNrxF8iJSgzM2uOE+kyZTMz60ecYMzMLAsnGDMzy8IJxszMsnCCMTOzLJxgzMwsCycYMzPLwgnGzMyycIIxM7MsnGDMzCwLJxgzM8vCCcbMzLJwgjEzsyycYMzMLAsnGDMzy8IJxszMsnCCMTOzLBqeYCSNkfSQpKckbZf01yl+naSfStqSlg+U6lwrqUPS05Kml+KTJW1N+25V+h5mSYMk3Z3iGySNa/Rxmpmd6poxgukCPh0RvwdcAFwlaWLad0tETErLdwDSvjbgHGAGcJukAan87cA8YEJaZqT4XGBfRJwF3ALc1IDjMjOzkoYnmIjYHRGPp/UDwFPAqF6qzATuioiDEfEc0AFMkTQSGBIRj0VEAMuAi0t1lqb1lcDUyujGzMwao6lzMOnU1XnAhhS6WtKTkhZLGpZio4DnS9U6U2xUWq+Od6sTEV3AfuCMGo8/T1K7pPY9e/b0yTGZmVmhaQlG0uuBe4BrIuIXFKe7zgQmAbuBL1aK1qgevcR7q9M9ELEwIlojorWlpeUYj8DMzHrTlAQj6VUUyeWOiLgXICJeiIhDEfEy8DVgSireCYwpVR8N7Erx0TXi3epIGggMBfbmORozM6ulGVeRCVgEPBURXyrFR5aKfQjYltZXAW3pyrDxFJP5GyNiN3BA0gWpzcuA+0t15qT1WcC6NE9jZmYNMrAJj3kh8DFgq6QtKfZZ4FJJkyhOZe0EPgEQEdslrQB2UFyBdlVEHEr1rgSWAIOB1WmBIoEtl9RBMXJpy3xMZmZWpeEJJiK+R+05ku/0UmcBsKBGvB04t0b8JWD2K+immZm9Qv4kv5mZZeEEY2ZmWTjBmJlZFk4wZmaWhROMmZll4QRjZmZZOMGYmVkWTjBmZpaFE4yZmWXhBGNmZlk4wZiZWRZOMGZmloUTjJmZZeEEY2ZmWTjBmJlZFk4wZmaWhROMmZll0a8TjKQZkp6W1CHpM83uj5nZqaTfJhhJA4C/B/4EmAhcKmlic3tlZnbq6LcJBpgCdETEsxHx/4C7gJlN7pOZ2SlDEdHsPmQhaRYwIyL+Q9r+GHB+RFxdKjMPmJc23wY83fCONs4I4GfN7oQdN79+J6/+/tq9JSJaau0Y2OieNJBqxLpl04hYCCxsTHeaS1J7RLQ2ux92fPz6nbxO5deuP58i6wTGlLZHA7ua1Bczs1NOf04wm4AJksZLejXQBqxqcp/MzE4Z/fYUWUR0SboaeAAYACyOiO1N7lYznRKnAvsxv34nr1P2teu3k/xmZtZc/fkUmZmZNZETjJmZZeEEcwqQNEnSB5rdj/5C0qN1lLlG0msb0Z+qx32DpP9U2n6zpJWN7seJwM9F83kO5hQg6XKgtfwh09K+gRHR1fhe9W+SdlI853V/wE7SgIg41NN2nW2MA/53RJx7LPX6oxPpuTie1/IVPNYJ8572COYEI2mcpKckfU3SdkkPShqc9k2StF7Sk5LukzSsRv3ZkrZJ+oGkR9Il2tcDl0jaIukSSddJWijpQWCZpBZJ90jalJYLU1t/mOpskfSEpNMljUztbkmP856GPkEnAEm/TD8vkvRdSSsl/VDSHSr8FfBm4CFJD6Wy0yQ9JulxSd+U9PoU3ynp85K+B8yusf1dSa2p7IiUuJB0uaT7Jf1zuqHr/NS9G4Ez0+tzc/p92pbqvEbS1yVtTa/n+0pt3ZvaekbS3zXsycyrt+ficknfkvRtSc9JulrSp9Lzsl7S8FTuzPS8bJb0fyS9vfpB0vtpuaR16fn7jyl+kaSHJH0D2CppQOrHpvQe/kQqd8R7KpVdkra3SvrPqWxvvw/flPRt4EFJr5O0OD3WE5Kac5usiPByAi3AOKALmJS2VwB/kdafBP4wrV8P/M8a9bcCo9L6G9LPy4GvlspcB2wGBqftbwDvTutjgafS+reBC9P66ykua/808LkUGwCc3uznrAmv0S/Tz4uA/RQf4j0NeKz0PO4ERqT1EcAjwOvS9t8Cny+V+6+ltqu3v0sxEqq0s7P0mu4GzgAGA9uA1vT7s63q92lbWv808PW0/nbgJ8BrUlvPAkPT9o+BMc1+nvvgdertubgc6ABOB1rS63hF2ncLcE1aXwtMSOvnA+tqPM51wA/S6zACeJ7iH4yLgF8B41O5ecB/S+uDgHZgfK33FDAZWFN6jMp7ubffh05geNr+Hxz+u/EG4F8qv3+NXPrt52BOcs9FxJa0vhkYJ2koxS/Zwym+FPhmjbrfB5ZIWgHc28tjrIqI36T19wMTpd/dXWeIpNNTW1+SdAdwb0R0StoELJb0KuBbpX6eqjZGRCeApC0Uf8S+V1XmAoo7en8/PcevpkhGFXdXla/e7smaiPh5eux7gXcD3+ql/LuBrwBExA8l/Rg4O+1bGxH7U1s7gLdQ/KHszx6KiAPAAUn7Kf6hguKftHekUeYfAN8svTcG9dDW/en99Js0ap0C/CvF78dzqcy01O6stD0UmEDxofBu7ylJzwJvlfQV4J+AB+s4njURsbf0WB+U9F/S9mtI/zzW0U6fcYI5MR0srR+i+M+oLhFxhaTzgX8PbJE0qYeivyqtnwa8q5RwKm6U9E/AB4D1kt4fEY9Iem9qf7mkmyNiWb3964eqX6ta7ylRvPkv7aGNX/Wy3cXhU9mvqSpXPYF6tAnVWvfnq6jnOPqb8jG/XNp+meL4TwP+NSJ6eg+V9fRalF9LAZ+MiAeqK9d6T0l6JzAduAr4CPBxev99qH6sP4+Ipt7A13MwJ4n03+U+HZ7z+BjwcHU5SWdGxIaI+DzFHVzHAAcoht09eRAo32V6UqmtrRFxE8Vw/u2S3gK8GBFfAxYB//aVH12/VH7O1wMXSjoLQNJrJZ3dY83udlKcLgGYVbXvjyUNVzFHdzHFiLO31/oR4KOpD2dT/Efbn+8gfrTf+15FxC+A5yTNBlDhnT0Un5nmuM6gODW2qUaZB4Ar00gFSWenuZIj3lOSRgCnRcQ9wH/n8PtsJz3/PlQ/1ieVhl6SzqvroPuYE8zJZQ5ws6QngUkU8zDVbk6Tgtso/qD8AHiI4hTYFkmX1KjzV0BrmnjcAVyR4tekScYfAL8BVlO8ebZIegL4c+DLfXh8/clCYLWkhyJiD8U58jvTa7eeYg6kHl+g+KP0KMU597LvAcuBLcA9EdGeTpl9P71uN1eVvw0YIGkrxWm4yyPiIP3UUZ6Len0UmJveA9vp+TulNlKcyloP3BARtW6s+4/ADuDx9P78XxQjpYs48j01CvhuOu26BLg2tdHb70PZDcCrgCfTY91Q19H2MV+mbHYSUi+XnltjSbqO4sKPLzS7Lycaj2DMzCwLj2DMzCwLj2DMzCwLJxgzM8vCCcbMzLJwgjE7gaX7Uc1K6/8oaWJa/2xVuaPe4dms0TzJb3YCk7SE4o7AK6viv4yI1zenV2b18QjGrA9Juix9YPUH6Q67b5G0NsXWShqbyi2RdKukRyU9WxqlSNJXJe1It+l5Y6nt70pqlXQjMDh9cPaOtO+Xpfo36/BdeC9J8Zp3fm7082OnllPhfkNmDSHpHOBzFHeg/pmKW74vBZZFxFJJHwdupbitC8BIihtQvh1YBawEPgS8Dfh94E0Un/xeXH6ciPiMpKt7uEfWhynu8vBOik96b5L0SNp3HnAOsIvitjIXcuSNOc36jEcwZn3nj4CVkb5kLN3Z9l0UX4cAxW1d3l0q/62IeDkidlAkE4D3AndGxKF0u5F1x9iHd5fqv0Bxv7p/l/ZtjIjOiHiZ4vYy446xbbNj4gRj1nfE0e9oXN5fvg+YeihzPH3oyal4x2RrIicYs76zFvhIuqMu6RTZo0Bb2v9Rjn5K6hGgTcU3Go4E3tdDud9W7spbo/4lqX4LxYho4zEeh1mf8H8wZn0kIrZLWgA8LOkQ8ATFnaoXS/obYA/wl0dp5j6KU21bKb6F8IivZEgWUtwp9/GI+GhV/XdR3EU7KL4d8/+qxlf9muXmy5TNzCwLnyIzM7MsnGDMzCwLJxgzM8vCCcbMzLJwgjEzsyycYMzMLAsnGDMzy+L/A7cr3L8eao8AAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.countplot(df5[\"condition\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "conditions = [\n",
    "    (df5['condition'] == \"no stress\"),\n",
    "    (df5['condition'] == \"interruption condition\"),\n",
    "    (df5['condition'] == \"time pressure\" )\n",
    "    ]\n",
    "values = [1,0,-1]\n",
    "df5['condition2'] = np.select(conditions, values)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "df6 = pd.read_csv(\"frequency_domain_features_test.csv\")\n",
    "df7 = pd.read_csv(\"heart_rate_non_linear_features_test.csv\")\n",
    "df8 = pd.read_csv(\"time_domain_features_test.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((41033, 12), (41033, 7), (41033, 19))"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df6.shape, df7.shape, df8.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "df9 = df6.merge(df7,how='inner',on='uuid')\n",
    "df10 = df9.merge(df8,how='inner',on='uuid')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "conditions = [\n",
    "    (df10['condition'] == \"no stress\"),\n",
    "    (df10['condition'] == \"interruption condition\"),\n",
    "    (df10['condition'] == \"time pressure\" )\n",
    "    ]\n",
    "values = [1,0,-1]\n",
    "df10['condition2'] = np.select(conditions, values)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(41033, 37)"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df10.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "dftest = df10.drop(['uuid','condition','datasetId','SD1','SKEW_REL_RR','KURT_REL_RR',\n",
    "                    'SDSD_REL_RR', 'SDRR', 'HF_NU','SDSD'], axis = 1)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "dftrain = df5.drop(['uuid','condition','datasetId','SD1','SKEW_REL_RR','KURT_REL_RR',\n",
    "                    'SDSD_REL_RR', 'SDRR', 'HF_NU','SDSD'], axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>VLF</th>\n",
       "      <th>VLF_PCT</th>\n",
       "      <th>LF</th>\n",
       "      <th>LF_PCT</th>\n",
       "      <th>LF_NU</th>\n",
       "      <th>HF</th>\n",
       "      <th>HF_PCT</th>\n",
       "      <th>TP</th>\n",
       "      <th>LF_HF</th>\n",
       "      <th>HF_LF</th>\n",
       "      <th>SD2</th>\n",
       "      <th>sampen</th>\n",
       "      <th>higuci</th>\n",
       "      <th>MEAN_RR</th>\n",
       "      <th>MEDIAN_RR</th>\n",
       "      <th>RMSSD</th>\n",
       "      <th>SDRR_RMSSD</th>\n",
       "      <th>HR</th>\n",
       "      <th>pNN25</th>\n",
       "      <th>pNN50</th>\n",
       "      <th>KURT</th>\n",
       "      <th>SKEW</th>\n",
       "      <th>MEAN_REL_RR</th>\n",
       "      <th>MEDIAN_REL_RR</th>\n",
       "      <th>SDRR_REL_RR</th>\n",
       "      <th>RMSSD_REL_RR</th>\n",
       "      <th>SDRR_RMSSD_REL_RR</th>\n",
       "      <th>condition2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2661.894136</td>\n",
       "      <td>72.203287</td>\n",
       "      <td>1009.249419</td>\n",
       "      <td>27.375666</td>\n",
       "      <td>98.485263</td>\n",
       "      <td>15.522603</td>\n",
       "      <td>0.421047</td>\n",
       "      <td>3686.666157</td>\n",
       "      <td>65.018055</td>\n",
       "      <td>0.015380</td>\n",
       "      <td>199.061782</td>\n",
       "      <td>2.139754</td>\n",
       "      <td>1.163485</td>\n",
       "      <td>885.157845</td>\n",
       "      <td>853.763730</td>\n",
       "      <td>15.554505</td>\n",
       "      <td>9.063146</td>\n",
       "      <td>69.499952</td>\n",
       "      <td>11.133333</td>\n",
       "      <td>0.533333</td>\n",
       "      <td>-0.856554</td>\n",
       "      <td>0.335218</td>\n",
       "      <td>-0.000203</td>\n",
       "      <td>-0.000179</td>\n",
       "      <td>0.017080</td>\n",
       "      <td>0.007969</td>\n",
       "      <td>2.143342</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2314.265450</td>\n",
       "      <td>76.975728</td>\n",
       "      <td>690.113275</td>\n",
       "      <td>22.954139</td>\n",
       "      <td>99.695397</td>\n",
       "      <td>2.108525</td>\n",
       "      <td>0.070133</td>\n",
       "      <td>3006.487251</td>\n",
       "      <td>327.296635</td>\n",
       "      <td>0.003055</td>\n",
       "      <td>114.634458</td>\n",
       "      <td>2.174499</td>\n",
       "      <td>1.084711</td>\n",
       "      <td>939.425371</td>\n",
       "      <td>948.357865</td>\n",
       "      <td>12.964439</td>\n",
       "      <td>6.272369</td>\n",
       "      <td>64.363150</td>\n",
       "      <td>5.600000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-0.408190</td>\n",
       "      <td>-0.155286</td>\n",
       "      <td>-0.000059</td>\n",
       "      <td>0.000611</td>\n",
       "      <td>0.013978</td>\n",
       "      <td>0.004769</td>\n",
       "      <td>2.930855</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1373.887112</td>\n",
       "      <td>51.152225</td>\n",
       "      <td>1298.222619</td>\n",
       "      <td>48.335104</td>\n",
       "      <td>98.950472</td>\n",
       "      <td>13.769729</td>\n",
       "      <td>0.512671</td>\n",
       "      <td>2685.879461</td>\n",
       "      <td>94.280910</td>\n",
       "      <td>0.010607</td>\n",
       "      <td>118.939253</td>\n",
       "      <td>2.135350</td>\n",
       "      <td>1.176315</td>\n",
       "      <td>898.186047</td>\n",
       "      <td>907.006860</td>\n",
       "      <td>16.305279</td>\n",
       "      <td>5.182201</td>\n",
       "      <td>67.450066</td>\n",
       "      <td>13.066667</td>\n",
       "      <td>0.200000</td>\n",
       "      <td>0.351789</td>\n",
       "      <td>-0.656813</td>\n",
       "      <td>-0.000011</td>\n",
       "      <td>-0.000263</td>\n",
       "      <td>0.018539</td>\n",
       "      <td>0.008716</td>\n",
       "      <td>2.127053</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2410.357408</td>\n",
       "      <td>70.180308</td>\n",
       "      <td>1005.981659</td>\n",
       "      <td>29.290305</td>\n",
       "      <td>98.224706</td>\n",
       "      <td>18.181913</td>\n",
       "      <td>0.529387</td>\n",
       "      <td>3434.520980</td>\n",
       "      <td>55.328701</td>\n",
       "      <td>0.018074</td>\n",
       "      <td>127.318597</td>\n",
       "      <td>2.178341</td>\n",
       "      <td>1.179688</td>\n",
       "      <td>881.757865</td>\n",
       "      <td>893.460030</td>\n",
       "      <td>15.720468</td>\n",
       "      <td>5.748591</td>\n",
       "      <td>68.809562</td>\n",
       "      <td>11.800000</td>\n",
       "      <td>0.133333</td>\n",
       "      <td>-0.504947</td>\n",
       "      <td>-0.386138</td>\n",
       "      <td>0.000112</td>\n",
       "      <td>0.000494</td>\n",
       "      <td>0.017761</td>\n",
       "      <td>0.008660</td>\n",
       "      <td>2.050988</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1151.177330</td>\n",
       "      <td>43.918366</td>\n",
       "      <td>1421.782051</td>\n",
       "      <td>54.242160</td>\n",
       "      <td>96.720007</td>\n",
       "      <td>48.215822</td>\n",
       "      <td>1.839473</td>\n",
       "      <td>2621.175204</td>\n",
       "      <td>29.487873</td>\n",
       "      <td>0.033912</td>\n",
       "      <td>87.718281</td>\n",
       "      <td>2.221121</td>\n",
       "      <td>1.249612</td>\n",
       "      <td>809.625331</td>\n",
       "      <td>811.184865</td>\n",
       "      <td>19.213819</td>\n",
       "      <td>3.266724</td>\n",
       "      <td>74.565728</td>\n",
       "      <td>20.200000</td>\n",
       "      <td>0.200000</td>\n",
       "      <td>-0.548408</td>\n",
       "      <td>-0.154252</td>\n",
       "      <td>-0.000100</td>\n",
       "      <td>-0.002736</td>\n",
       "      <td>0.023715</td>\n",
       "      <td>0.013055</td>\n",
       "      <td>1.816544</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           VLF    VLF_PCT           LF     LF_PCT      LF_NU         HF  \\\n",
       "0  2661.894136  72.203287  1009.249419  27.375666  98.485263  15.522603   \n",
       "1  2314.265450  76.975728   690.113275  22.954139  99.695397   2.108525   \n",
       "2  1373.887112  51.152225  1298.222619  48.335104  98.950472  13.769729   \n",
       "3  2410.357408  70.180308  1005.981659  29.290305  98.224706  18.181913   \n",
       "4  1151.177330  43.918366  1421.782051  54.242160  96.720007  48.215822   \n",
       "\n",
       "     HF_PCT           TP       LF_HF     HF_LF         SD2    sampen  \\\n",
       "0  0.421047  3686.666157   65.018055  0.015380  199.061782  2.139754   \n",
       "1  0.070133  3006.487251  327.296635  0.003055  114.634458  2.174499   \n",
       "2  0.512671  2685.879461   94.280910  0.010607  118.939253  2.135350   \n",
       "3  0.529387  3434.520980   55.328701  0.018074  127.318597  2.178341   \n",
       "4  1.839473  2621.175204   29.487873  0.033912   87.718281  2.221121   \n",
       "\n",
       "     higuci     MEAN_RR   MEDIAN_RR      RMSSD  SDRR_RMSSD         HR  \\\n",
       "0  1.163485  885.157845  853.763730  15.554505    9.063146  69.499952   \n",
       "1  1.084711  939.425371  948.357865  12.964439    6.272369  64.363150   \n",
       "2  1.176315  898.186047  907.006860  16.305279    5.182201  67.450066   \n",
       "3  1.179688  881.757865  893.460030  15.720468    5.748591  68.809562   \n",
       "4  1.249612  809.625331  811.184865  19.213819    3.266724  74.565728   \n",
       "\n",
       "       pNN25     pNN50      KURT      SKEW  MEAN_REL_RR  MEDIAN_REL_RR  \\\n",
       "0  11.133333  0.533333 -0.856554  0.335218    -0.000203      -0.000179   \n",
       "1   5.600000  0.000000 -0.408190 -0.155286    -0.000059       0.000611   \n",
       "2  13.066667  0.200000  0.351789 -0.656813    -0.000011      -0.000263   \n",
       "3  11.800000  0.133333 -0.504947 -0.386138     0.000112       0.000494   \n",
       "4  20.200000  0.200000 -0.548408 -0.154252    -0.000100      -0.002736   \n",
       "\n",
       "   SDRR_REL_RR  RMSSD_REL_RR  SDRR_RMSSD_REL_RR  condition2  \n",
       "0     0.017080      0.007969           2.143342           1  \n",
       "1     0.013978      0.004769           2.930855           0  \n",
       "2     0.018539      0.008716           2.127053           0  \n",
       "3     0.017761      0.008660           2.050988           1  \n",
       "4     0.023715      0.013055           1.816544           1  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dftrain.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>VLF</th>\n",
       "      <th>VLF_PCT</th>\n",
       "      <th>LF</th>\n",
       "      <th>LF_PCT</th>\n",
       "      <th>LF_NU</th>\n",
       "      <th>HF</th>\n",
       "      <th>HF_PCT</th>\n",
       "      <th>TP</th>\n",
       "      <th>LF_HF</th>\n",
       "      <th>HF_LF</th>\n",
       "      <th>SD2</th>\n",
       "      <th>sampen</th>\n",
       "      <th>higuci</th>\n",
       "      <th>MEAN_RR</th>\n",
       "      <th>MEDIAN_RR</th>\n",
       "      <th>RMSSD</th>\n",
       "      <th>SDRR_RMSSD</th>\n",
       "      <th>pNN25</th>\n",
       "      <th>pNN50</th>\n",
       "      <th>KURT</th>\n",
       "      <th>SKEW</th>\n",
       "      <th>MEAN_REL_RR</th>\n",
       "      <th>MEDIAN_REL_RR</th>\n",
       "      <th>SDRR_REL_RR</th>\n",
       "      <th>RMSSD_REL_RR</th>\n",
       "      <th>SDRR_RMSSD_REL_RR</th>\n",
       "      <th>condition2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1868.532278</td>\n",
       "      <td>76.511189</td>\n",
       "      <td>570.643114</td>\n",
       "      <td>23.366245</td>\n",
       "      <td>99.478197</td>\n",
       "      <td>2.993254</td>\n",
       "      <td>0.122565</td>\n",
       "      <td>2442.168645</td>\n",
       "      <td>190.643094</td>\n",
       "      <td>0.005245</td>\n",
       "      <td>115.862444</td>\n",
       "      <td>2.209659</td>\n",
       "      <td>1.100715</td>\n",
       "      <td>934.665288</td>\n",
       "      <td>939.03173</td>\n",
       "      <td>11.801781</td>\n",
       "      <td>6.959924</td>\n",
       "      <td>3.933333</td>\n",
       "      <td>0.133333</td>\n",
       "      <td>-0.680262</td>\n",
       "      <td>-0.233075</td>\n",
       "      <td>0.000016</td>\n",
       "      <td>0.000288</td>\n",
       "      <td>0.012933</td>\n",
       "      <td>0.004578</td>\n",
       "      <td>2.825038</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>568.742845</td>\n",
       "      <td>26.301350</td>\n",
       "      <td>1553.971621</td>\n",
       "      <td>71.862973</td>\n",
       "      <td>97.509212</td>\n",
       "      <td>39.694850</td>\n",
       "      <td>1.835677</td>\n",
       "      <td>2162.409316</td>\n",
       "      <td>39.147940</td>\n",
       "      <td>0.025544</td>\n",
       "      <td>77.118903</td>\n",
       "      <td>2.186132</td>\n",
       "      <td>1.290615</td>\n",
       "      <td>817.062380</td>\n",
       "      <td>816.33879</td>\n",
       "      <td>20.558810</td>\n",
       "      <td>2.699200</td>\n",
       "      <td>24.600000</td>\n",
       "      <td>0.533333</td>\n",
       "      <td>-0.034454</td>\n",
       "      <td>-0.051689</td>\n",
       "      <td>0.000040</td>\n",
       "      <td>-0.002749</td>\n",
       "      <td>0.025148</td>\n",
       "      <td>0.013921</td>\n",
       "      <td>1.806517</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2101.871207</td>\n",
       "      <td>75.836461</td>\n",
       "      <td>655.175895</td>\n",
       "      <td>23.639042</td>\n",
       "      <td>97.829386</td>\n",
       "      <td>14.536877</td>\n",
       "      <td>0.524497</td>\n",
       "      <td>2771.583978</td>\n",
       "      <td>45.069921</td>\n",
       "      <td>0.022188</td>\n",
       "      <td>125.044377</td>\n",
       "      <td>2.051571</td>\n",
       "      <td>1.226663</td>\n",
       "      <td>876.762022</td>\n",
       "      <td>894.19889</td>\n",
       "      <td>13.853737</td>\n",
       "      <td>6.401942</td>\n",
       "      <td>7.066667</td>\n",
       "      <td>0.533333</td>\n",
       "      <td>-0.206953</td>\n",
       "      <td>-0.589940</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>-0.000172</td>\n",
       "      <td>0.015533</td>\n",
       "      <td>0.008149</td>\n",
       "      <td>1.906180</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>5757.544433</td>\n",
       "      <td>90.562305</td>\n",
       "      <td>592.913021</td>\n",
       "      <td>9.326123</td>\n",
       "      <td>98.817806</td>\n",
       "      <td>7.093235</td>\n",
       "      <td>0.111572</td>\n",
       "      <td>6357.550689</td>\n",
       "      <td>83.588517</td>\n",
       "      <td>0.011963</td>\n",
       "      <td>302.029812</td>\n",
       "      <td>2.080910</td>\n",
       "      <td>1.085143</td>\n",
       "      <td>1038.640693</td>\n",
       "      <td>998.91429</td>\n",
       "      <td>16.457194</td>\n",
       "      <td>12.986774</td>\n",
       "      <td>10.800000</td>\n",
       "      <td>1.866667</td>\n",
       "      <td>-0.820407</td>\n",
       "      <td>0.487198</td>\n",
       "      <td>-0.000238</td>\n",
       "      <td>-0.000464</td>\n",
       "      <td>0.016882</td>\n",
       "      <td>0.007587</td>\n",
       "      <td>2.225116</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>964.696325</td>\n",
       "      <td>70.256575</td>\n",
       "      <td>374.939530</td>\n",
       "      <td>27.305968</td>\n",
       "      <td>91.805057</td>\n",
       "      <td>33.468834</td>\n",
       "      <td>2.437457</td>\n",
       "      <td>1373.104689</td>\n",
       "      <td>11.202647</td>\n",
       "      <td>0.089265</td>\n",
       "      <td>72.579248</td>\n",
       "      <td>2.068728</td>\n",
       "      <td>1.252547</td>\n",
       "      <td>774.548508</td>\n",
       "      <td>778.90508</td>\n",
       "      <td>10.273114</td>\n",
       "      <td>5.020664</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>0.066667</td>\n",
       "      <td>1.738453</td>\n",
       "      <td>-0.005082</td>\n",
       "      <td>0.000054</td>\n",
       "      <td>-0.000006</td>\n",
       "      <td>0.013479</td>\n",
       "      <td>0.007811</td>\n",
       "      <td>1.725606</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           VLF    VLF_PCT           LF     LF_PCT      LF_NU         HF  \\\n",
       "0  1868.532278  76.511189   570.643114  23.366245  99.478197   2.993254   \n",
       "1   568.742845  26.301350  1553.971621  71.862973  97.509212  39.694850   \n",
       "2  2101.871207  75.836461   655.175895  23.639042  97.829386  14.536877   \n",
       "3  5757.544433  90.562305   592.913021   9.326123  98.817806   7.093235   \n",
       "4   964.696325  70.256575   374.939530  27.305968  91.805057  33.468834   \n",
       "\n",
       "     HF_PCT           TP       LF_HF     HF_LF         SD2    sampen  \\\n",
       "0  0.122565  2442.168645  190.643094  0.005245  115.862444  2.209659   \n",
       "1  1.835677  2162.409316   39.147940  0.025544   77.118903  2.186132   \n",
       "2  0.524497  2771.583978   45.069921  0.022188  125.044377  2.051571   \n",
       "3  0.111572  6357.550689   83.588517  0.011963  302.029812  2.080910   \n",
       "4  2.437457  1373.104689   11.202647  0.089265   72.579248  2.068728   \n",
       "\n",
       "     higuci      MEAN_RR  MEDIAN_RR      RMSSD  SDRR_RMSSD      pNN25  \\\n",
       "0  1.100715   934.665288  939.03173  11.801781    6.959924   3.933333   \n",
       "1  1.290615   817.062380  816.33879  20.558810    2.699200  24.600000   \n",
       "2  1.226663   876.762022  894.19889  13.853737    6.401942   7.066667   \n",
       "3  1.085143  1038.640693  998.91429  16.457194   12.986774  10.800000   \n",
       "4  1.252547   774.548508  778.90508  10.273114    5.020664   2.000000   \n",
       "\n",
       "      pNN50      KURT      SKEW  MEAN_REL_RR  MEDIAN_REL_RR  SDRR_REL_RR  \\\n",
       "0  0.133333 -0.680262 -0.233075     0.000016       0.000288     0.012933   \n",
       "1  0.533333 -0.034454 -0.051689     0.000040      -0.002749     0.025148   \n",
       "2  0.533333 -0.206953 -0.589940     0.000010      -0.000172     0.015533   \n",
       "3  1.866667 -0.820407  0.487198    -0.000238      -0.000464     0.016882   \n",
       "4  0.066667  1.738453 -0.005082     0.000054      -0.000006     0.013479   \n",
       "\n",
       "   RMSSD_REL_RR  SDRR_RMSSD_REL_RR  condition2  \n",
       "0      0.004578           2.825038          -1  \n",
       "1      0.013921           1.806517           1  \n",
       "2      0.008149           1.906180           1  \n",
       "3      0.007587           2.225116          -1  \n",
       "4      0.007811           1.725606           0  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dftest.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((369289, 28), (41033, 27))"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dftrain.shape, dftest.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df10.duplicated(subset='uuid').sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df5.duplicated(subset='uuid').sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = dftrain.drop(\"HR\",axis=1)\n",
    "y = dftrain[\"HR\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "X_train, X_val, y_train, y_val = train_test_split(X, y, test_size=0.1, random_state=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((332360, 27), (36929, 27))"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X_train.shape, X_val.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((332360,), (36929,))"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "y_train.shape, y_val.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "#cp = dftrain.corr()\n",
    "#mask = np.zeros_like(cp)\n",
    "#mask[np.triu_indices_from(mask)] = True\n",
    "#plt.figure(figsize=(30,30))\n",
    "#with sns.axes_style(\"white\"):\n",
    "#    sns.heatmap(cp,annot=True,linewidth=2,mask = mask,cmap=\"coolwarm\")\n",
    "#plt.title(\"Correlation Plot\")\n",
    "#plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "scrolled": true
   },
   "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>VLF</th>\n",
       "      <th>VLF_PCT</th>\n",
       "      <th>LF</th>\n",
       "      <th>LF_PCT</th>\n",
       "      <th>LF_NU</th>\n",
       "      <th>HF</th>\n",
       "      <th>HF_PCT</th>\n",
       "      <th>TP</th>\n",
       "      <th>LF_HF</th>\n",
       "      <th>HF_LF</th>\n",
       "      <th>SD2</th>\n",
       "      <th>sampen</th>\n",
       "      <th>higuci</th>\n",
       "      <th>MEAN_RR</th>\n",
       "      <th>MEDIAN_RR</th>\n",
       "      <th>RMSSD</th>\n",
       "      <th>SDRR_RMSSD</th>\n",
       "      <th>pNN25</th>\n",
       "      <th>pNN50</th>\n",
       "      <th>KURT</th>\n",
       "      <th>SKEW</th>\n",
       "      <th>MEAN_REL_RR</th>\n",
       "      <th>MEDIAN_REL_RR</th>\n",
       "      <th>SDRR_REL_RR</th>\n",
       "      <th>RMSSD_REL_RR</th>\n",
       "      <th>SDRR_RMSSD_REL_RR</th>\n",
       "      <th>condition2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-0.850461</td>\n",
       "      <td>-0.139308</td>\n",
       "      <td>-0.984299</td>\n",
       "      <td>0.123290</td>\n",
       "      <td>-0.183804</td>\n",
       "      <td>-0.404517</td>\n",
       "      <td>0.204009</td>\n",
       "      <td>-1.106199</td>\n",
       "      <td>-0.270674</td>\n",
       "      <td>0.126860</td>\n",
       "      <td>-0.739628</td>\n",
       "      <td>0.645106</td>\n",
       "      <td>0.332048</td>\n",
       "      <td>-0.158395</td>\n",
       "      <td>-0.087086</td>\n",
       "      <td>-1.207215</td>\n",
       "      <td>-0.424490</td>\n",
       "      <td>-1.054247</td>\n",
       "      <td>-0.874384</td>\n",
       "      <td>-0.548623</td>\n",
       "      <td>-0.319678</td>\n",
       "      <td>-0.556038</td>\n",
       "      <td>-0.446828</td>\n",
       "      <td>-1.121930</td>\n",
       "      <td>-0.724240</td>\n",
       "      <td>-0.523978</td>\n",
       "      <td>-1.797587</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.207041</td>\n",
       "      <td>-0.494574</td>\n",
       "      <td>1.859567</td>\n",
       "      <td>0.603678</td>\n",
       "      <td>0.948527</td>\n",
       "      <td>-0.631940</td>\n",
       "      <td>-0.786890</td>\n",
       "      <td>0.735448</td>\n",
       "      <td>0.206940</td>\n",
       "      <td>-0.878365</td>\n",
       "      <td>-0.043159</td>\n",
       "      <td>0.512342</td>\n",
       "      <td>0.642869</td>\n",
       "      <td>1.082557</td>\n",
       "      <td>1.114223</td>\n",
       "      <td>1.969545</td>\n",
       "      <td>-0.542877</td>\n",
       "      <td>2.435465</td>\n",
       "      <td>1.817293</td>\n",
       "      <td>-0.324967</td>\n",
       "      <td>-0.675614</td>\n",
       "      <td>-0.075092</td>\n",
       "      <td>-0.923057</td>\n",
       "      <td>0.908367</td>\n",
       "      <td>0.478176</td>\n",
       "      <td>0.072576</td>\n",
       "      <td>-0.484008</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3.465601</td>\n",
       "      <td>1.503677</td>\n",
       "      <td>-0.020234</td>\n",
       "      <td>-1.511359</td>\n",
       "      <td>-0.388739</td>\n",
       "      <td>0.458113</td>\n",
       "      <td>-0.557876</td>\n",
       "      <td>3.276715</td>\n",
       "      <td>-0.278136</td>\n",
       "      <td>0.319469</td>\n",
       "      <td>2.597732</td>\n",
       "      <td>-0.637816</td>\n",
       "      <td>-1.277965</td>\n",
       "      <td>0.611599</td>\n",
       "      <td>0.246970</td>\n",
       "      <td>0.633437</td>\n",
       "      <td>1.986375</td>\n",
       "      <td>0.336757</td>\n",
       "      <td>1.750001</td>\n",
       "      <td>-0.882233</td>\n",
       "      <td>0.591298</td>\n",
       "      <td>-0.544726</td>\n",
       "      <td>0.736833</td>\n",
       "      <td>0.827142</td>\n",
       "      <td>0.541158</td>\n",
       "      <td>-0.139616</td>\n",
       "      <td>0.829571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-0.070555</td>\n",
       "      <td>-0.218677</td>\n",
       "      <td>0.513966</td>\n",
       "      <td>0.139282</td>\n",
       "      <td>-0.803162</td>\n",
       "      <td>1.430990</td>\n",
       "      <td>0.814190</td>\n",
       "      <td>0.120531</td>\n",
       "      <td>-0.288248</td>\n",
       "      <td>0.719753</td>\n",
       "      <td>0.567114</td>\n",
       "      <td>0.025691</td>\n",
       "      <td>0.273998</td>\n",
       "      <td>-0.399367</td>\n",
       "      <td>-0.155129</td>\n",
       "      <td>0.638723</td>\n",
       "      <td>0.252095</td>\n",
       "      <td>0.255412</td>\n",
       "      <td>1.211666</td>\n",
       "      <td>-1.041454</td>\n",
       "      <td>-0.297102</td>\n",
       "      <td>0.103730</td>\n",
       "      <td>-0.029783</td>\n",
       "      <td>1.110613</td>\n",
       "      <td>1.349639</td>\n",
       "      <td>-0.959245</td>\n",
       "      <td>0.829571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.081083</td>\n",
       "      <td>0.875895</td>\n",
       "      <td>-0.574053</td>\n",
       "      <td>-0.831399</td>\n",
       "      <td>0.777681</td>\n",
       "      <td>-0.695688</td>\n",
       "      <td>-0.770841</td>\n",
       "      <td>-0.111165</td>\n",
       "      <td>-0.098269</td>\n",
       "      <td>-0.732786</td>\n",
       "      <td>-0.367277</td>\n",
       "      <td>-0.022202</td>\n",
       "      <td>-1.292749</td>\n",
       "      <td>0.370264</td>\n",
       "      <td>0.437421</td>\n",
       "      <td>-0.542377</td>\n",
       "      <td>-0.203448</td>\n",
       "      <td>-0.484830</td>\n",
       "      <td>-0.605217</td>\n",
       "      <td>0.209012</td>\n",
       "      <td>-0.781857</td>\n",
       "      <td>-0.685079</td>\n",
       "      <td>0.625590</td>\n",
       "      <td>-0.755545</td>\n",
       "      <td>-1.037874</td>\n",
       "      <td>1.456332</td>\n",
       "      <td>0.829571</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        VLF   VLF_PCT        LF    LF_PCT     LF_NU        HF    HF_PCT  \\\n",
       "0 -0.850461 -0.139308 -0.984299  0.123290 -0.183804 -0.404517  0.204009   \n",
       "1  0.207041 -0.494574  1.859567  0.603678  0.948527 -0.631940 -0.786890   \n",
       "2  3.465601  1.503677 -0.020234 -1.511359 -0.388739  0.458113 -0.557876   \n",
       "3 -0.070555 -0.218677  0.513966  0.139282 -0.803162  1.430990  0.814190   \n",
       "4  0.081083  0.875895 -0.574053 -0.831399  0.777681 -0.695688 -0.770841   \n",
       "\n",
       "         TP     LF_HF     HF_LF       SD2    sampen    higuci   MEAN_RR  \\\n",
       "0 -1.106199 -0.270674  0.126860 -0.739628  0.645106  0.332048 -0.158395   \n",
       "1  0.735448  0.206940 -0.878365 -0.043159  0.512342  0.642869  1.082557   \n",
       "2  3.276715 -0.278136  0.319469  2.597732 -0.637816 -1.277965  0.611599   \n",
       "3  0.120531 -0.288248  0.719753  0.567114  0.025691  0.273998 -0.399367   \n",
       "4 -0.111165 -0.098269 -0.732786 -0.367277 -0.022202 -1.292749  0.370264   \n",
       "\n",
       "   MEDIAN_RR     RMSSD  SDRR_RMSSD     pNN25     pNN50      KURT      SKEW  \\\n",
       "0  -0.087086 -1.207215   -0.424490 -1.054247 -0.874384 -0.548623 -0.319678   \n",
       "1   1.114223  1.969545   -0.542877  2.435465  1.817293 -0.324967 -0.675614   \n",
       "2   0.246970  0.633437    1.986375  0.336757  1.750001 -0.882233  0.591298   \n",
       "3  -0.155129  0.638723    0.252095  0.255412  1.211666 -1.041454 -0.297102   \n",
       "4   0.437421 -0.542377   -0.203448 -0.484830 -0.605217  0.209012 -0.781857   \n",
       "\n",
       "   MEAN_REL_RR  MEDIAN_REL_RR  SDRR_REL_RR  RMSSD_REL_RR  SDRR_RMSSD_REL_RR  \\\n",
       "0    -0.556038      -0.446828    -1.121930     -0.724240          -0.523978   \n",
       "1    -0.075092      -0.923057     0.908367      0.478176           0.072576   \n",
       "2    -0.544726       0.736833     0.827142      0.541158          -0.139616   \n",
       "3     0.103730      -0.029783     1.110613      1.349639          -0.959245   \n",
       "4    -0.685079       0.625590    -0.755545     -1.037874           1.456332   \n",
       "\n",
       "   condition2  \n",
       "0   -1.797587  \n",
       "1   -0.484008  \n",
       "2    0.829571  \n",
       "3    0.829571  \n",
       "4    0.829571  "
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sc = StandardScaler()\n",
    "\n",
    "X_train = pd.DataFrame(sc.fit_transform(X_train),columns=X.columns)\n",
    "X_val = pd.DataFrame(sc.transform(X_val),columns=X.columns)\n",
    "\n",
    "X_train.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "#param_grid = {'alpha': np.arange(0,2,0.05)}\n",
    "#\n",
    "#grid_tree = GridSearchCV(Lasso(), param_grid ,cv=10, scoring= 'neg_mean_absolute_error')\n",
    "#grid_tree.fit(X_train, y_train)\n",
    "#print(grid_tree.best_estimator_)\n",
    "#print(np.abs(grid_tree.best_score_))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "# linear regression\n",
    "lr = LinearRegression()\n",
    "\n",
    "\n",
    "# fit models on training data \n",
    "lr.fit(X_train, y_train)\n",
    "\n",
    "# make predictions on test data\n",
    "y_pred_lr = lr.predict(X_val)\n",
    "\n",
    "# calculate preformance measures linear regression\n",
    "rmse_lr = np.sqrt(mean_squared_error(y_val, y_pred_lr))\n",
    "mse_lr = mean_squared_error(y_val, y_pred_lr)\n",
    "mae_lr = mean_absolute_error(y_val, y_pred_lr)\n",
    "r2_lr = r2_score(y_val, y_pred_lr)\n",
    "adjr2_lr = 1-(1-r2_lr)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Linear Regression</td>\n",
       "      <td>1.418821</td>\n",
       "      <td>2.013052</td>\n",
       "      <td>0.949578</td>\n",
       "      <td>0.981004</td>\n",
       "      <td>0.98099</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              Method      RMSE       MSE       MAE        R2   Adj_R2\n",
       "0  Linear Regression  1.418821  2.013052  0.949578  0.981004  0.98099"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "resultsDf1 = pd.DataFrame({'Method':['Linear Regression'], 'RMSE': rmse_lr , 'MSE': mse_lr , 'MAE': mae_lr ,\"R2\": r2_lr, 'Adj_R2': adjr2_lr})\n",
    "resultsDf1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "#linear_coef = pd.DataFrame()\n",
    "#linear_coef['linear_Col'] = X_train.columns\n",
    "#linear_coef['linear_Coeff'] = abs(np.round((lr.coef_),3))\n",
    "#linear_coef.sort_values(by='linear_Coeff',ascending=False).head(33)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "# lasso regression\n",
    "lasso = Lasso(alpha = 0.0)\n",
    "\n",
    "# fit lasso on training data\n",
    "lasso.fit(X_train, y_train)\n",
    "\n",
    "# make prediciton on test data\n",
    "y_pred_lasso = lasso.predict(X_val)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "# calculate preformance measures lasso regression\n",
    "rmse_lasso = np.sqrt(mean_squared_error(y_val, y_pred_lasso))\n",
    "mse_lasso = mean_squared_error(y_val, y_pred_lasso)\n",
    "mae_lasso = mean_absolute_error(y_val, y_pred_lasso)\n",
    "r2_lasso = r2_score(y_val, y_pred_lasso)\n",
    "adjr2_lasso = 1-(1-r2_lasso)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Lasso Regression</td>\n",
       "      <td>1.418689</td>\n",
       "      <td>2.012677</td>\n",
       "      <td>0.949435</td>\n",
       "      <td>0.981008</td>\n",
       "      <td>0.980994</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             Method      RMSE       MSE       MAE        R2    Adj_R2\n",
       "0  Lasso Regression  1.418689  2.012677  0.949435  0.981008  0.980994"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "resultsDf2 = pd.DataFrame({'Method':['Lasso Regression'], 'RMSE': rmse_lasso , 'MSE': mse_lasso , 'MAE': mae_lasso ,\"R2\": r2_lasso, 'Adj_R2': adjr2_lasso})\n",
    "resultsDf2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "#lasso_coef = pd.DataFrame()\n",
    "#lasso_coef['lasso_Col'] = X_train.columns\n",
    "#lasso_coef['lasso_Coeff'] = abs(np.round((lasso.coef_),3))\n",
    "#lasso_coef.sort_values(by='lasso_Coeff',ascending=False).head(33)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "# ridge regression\n",
    "ridge = Ridge(alpha = 1.95)\n",
    "\n",
    "# fit ridge on training data\n",
    "ridge.fit(X_train, y_train)\n",
    "\n",
    "# make prediciton on test data\n",
    "y_pred_ridge = ridge.predict(X_val)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Ridge Regression</td>\n",
       "      <td>1.418801</td>\n",
       "      <td>2.012996</td>\n",
       "      <td>0.949477</td>\n",
       "      <td>0.981005</td>\n",
       "      <td>0.980991</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             Method      RMSE       MSE       MAE        R2    Adj_R2\n",
       "0  Ridge Regression  1.418801  2.012996  0.949477  0.981005  0.980991"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# calculate preformance measures ridge regression\n",
    "rmse_ridge = np.sqrt(mean_squared_error(y_val, y_pred_ridge))\n",
    "mse_ridge = mean_squared_error(y_val, y_pred_ridge)\n",
    "mae_ridge = mean_absolute_error(y_val, y_pred_ridge)\n",
    "r2_ridge = r2_score(y_val, y_pred_ridge)\n",
    "adjr2_ridge = 1-(1-r2_ridge)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n",
    "\n",
    "resultsDf3 = pd.DataFrame({'Method':['Ridge Regression'], 'RMSE': rmse_ridge , 'MSE': mse_ridge , 'MAE': mae_ridge ,\"R2\": r2_ridge, 'Adj_R2': adjr2_ridge})\n",
    "resultsDf3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "#ridge_coef = pd.DataFrame()\n",
    "#ridge_coef['ridge_Col'] = X_train.columns\n",
    "#ridge_coef['ridge_Coeff'] = abs(np.round((ridge.coef_),3))\n",
    "#ridge_coef.sort_values(by='ridge_Coeff',ascending=False).head(33)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Linear Regression</td>\n",
       "      <td>1.418821</td>\n",
       "      <td>2.013052</td>\n",
       "      <td>0.949578</td>\n",
       "      <td>0.981004</td>\n",
       "      <td>0.980990</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Lasso Regression</td>\n",
       "      <td>1.418689</td>\n",
       "      <td>2.012677</td>\n",
       "      <td>0.949435</td>\n",
       "      <td>0.981008</td>\n",
       "      <td>0.980994</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Ridge Regression</td>\n",
       "      <td>1.418801</td>\n",
       "      <td>2.012996</td>\n",
       "      <td>0.949477</td>\n",
       "      <td>0.981005</td>\n",
       "      <td>0.980991</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              Method      RMSE       MSE       MAE        R2    Adj_R2\n",
       "0  Linear Regression  1.418821  2.013052  0.949578  0.981004  0.980990\n",
       "0   Lasso Regression  1.418689  2.012677  0.949435  0.981008  0.980994\n",
       "0   Ridge Regression  1.418801  2.012996  0.949477  0.981005  0.980991"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "linear = pd.concat([resultsDf1, resultsDf2, resultsDf3],axis=0)\n",
    "linear"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "#pd.concat([linear_coef, lasso_coef, ridge_coef],axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Residual plots\n",
    "\n",
    "#for i in ['VLF', 'VLF_PCT', 'LF', 'LF_PCT', 'LF_NU', 'HF', 'HF_PCT', 'TP',\n",
    "#       'LF_HF', 'HF_LF', 'SD2', 'sampen', 'higuci', 'MEAN_RR', 'MEDIAN_RR',\n",
    "#       'RMSSD', 'SDRR_RMSSD', 'pNN25', 'pNN50', 'KURT', 'SKEW', 'MEAN_REL_RR',\n",
    "#       'MEDIAN_REL_RR', 'SDRR_REL_RR', 'RMSSD_REL_RR', 'SDRR_RMSSD_REL_RR',\n",
    "#       'condition2']:\n",
    "#    fig = plt.figure(figsize=(10,8))\n",
    "#    sns.residplot(x= X_val[i], y= y_val, color='green', lowess=True )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Train R2 score:  0.9999985183182631\n",
      "Validation R2 score:  0.9999912482134314\n"
     ]
    }
   ],
   "source": [
    "rf=RandomForestRegressor(random_state=1, n_estimators=20, max_depth=50)\n",
    "rf.fit(X_train,y_train)\n",
    "#rf.score(X_train,y_train)\n",
    "\n",
    "#this line doesn't make sense to me\n",
    "y_pred_rf_train=rf.predict(X_train)\n",
    "#\n",
    "\n",
    "y_pred_rf_val=rf.predict(X_val)\n",
    "print('Train R2 score: ', r2_score(y_train,y_pred_rf_train))\n",
    "print('Validation R2 score: ', r2_score(y_val,y_pred_rf_val))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Random Forest</td>\n",
       "      <td>0.030454</td>\n",
       "      <td>0.000927</td>\n",
       "      <td>0.007101</td>\n",
       "      <td>0.999991</td>\n",
       "      <td>0.999991</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Method      RMSE       MSE       MAE        R2    Adj_R2\n",
       "0  Random Forest  0.030454  0.000927  0.007101  0.999991  0.999991"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rmse_rf = np.sqrt(mean_squared_error(y_val, y_pred_rf_val))\n",
    "mse_rf = mean_squared_error(y_val, y_pred_rf_val)\n",
    "mae_rf = mean_absolute_error(y_val, y_pred_rf_val)\n",
    "r2_rf = r2_score(y_val, y_pred_rf_val)\n",
    "adjr2_rf = 1-(1-r2_rf)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n",
    "\n",
    "resultsDf1 = pd.DataFrame({'Method':['Random Forest'], 'RMSE': rmse_rf , 'MSE': mse_rf , 'MAE': mae_rf ,\"R2\": r2_rf, 'Adj_R2': adjr2_rf})\n",
    "resultsDf1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Random Forest</td>\n",
       "      <td>0.012589</td>\n",
       "      <td>0.000158</td>\n",
       "      <td>0.003191</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>0.999999</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Method      RMSE       MSE       MAE        R2    Adj_R2\n",
       "0  Random Forest  0.012589  0.000158  0.003191  0.999999  0.999999"
      ]
     },
     "execution_count": 46,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#This block doesn't make sense to me\n",
    "\n",
    "rmse_rf = np.sqrt(mean_squared_error(y_train, y_pred_rf_train))\n",
    "mse_rf = mean_squared_error(y_train, y_pred_rf_train)\n",
    "mae_rf = mean_absolute_error(y_train, y_pred_rf_train)\n",
    "r2_rf = r2_score(y_train, y_pred_rf_train)\n",
    "adjr2_rf = 1-(1-r2_rf)*(len(y_train)-1)/(len(y_train)-len(X_train.columns)-1)\n",
    "\n",
    "resultsDf1 = pd.DataFrame({'Method':['Random Forest'], 'RMSE': rmse_rf , 'MSE': mse_rf , 'MAE': mae_rf ,\"R2\": r2_rf, 'Adj_R2': adjr2_rf})\n",
    "resultsDf1\n",
    "#"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [],
   "source": [
    "#rf_coef = pd.DataFrame()\n",
    "#rf_coef['rf_Col'] = X_train.columns\n",
    "#rf_coef['rf_Coeff'] = abs(np.round((rf.feature_importances_),5))\n",
    "#rf_coef.sort_values(by='rf_Coeff',ascending=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Gradient Boost\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Train R2 score:  0.9999734386011124\n",
      "Test R2 score:  0.9999295694376167\n"
     ]
    }
   ],
   "source": [
    "\n",
    "gb=GradientBoostingRegressor(random_state=1,n_estimators=50 ,max_depth=50,learning_rate=.1)\n",
    "gb.fit(X_train,y_train)\n",
    "#gb.score(X_train,y_train)\n",
    "y_pred_gbtr=gb.predict(X_train)\n",
    "y_pred_gbval=gb.predict(X_val)\n",
    "print('Train R2 score: ',r2_score(y_train,y_pred_gbtr))\n",
    "print('Test R2 score: ',r2_score(y_val,y_pred_gbval))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Method</th>\n",
       "      <th>RMSE</th>\n",
       "      <th>MSE</th>\n",
       "      <th>MAE</th>\n",
       "      <th>R2</th>\n",
       "      <th>Adj_R2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Gradient Boost</td>\n",
       "      <td>0.086393</td>\n",
       "      <td>0.007464</td>\n",
       "      <td>0.044539</td>\n",
       "      <td>0.99993</td>\n",
       "      <td>0.99993</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           Method      RMSE       MSE       MAE       R2   Adj_R2\n",
       "0  Gradient Boost  0.086393  0.007464  0.044539  0.99993  0.99993"
      ]
     },
     "execution_count": 50,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rmse_gb = np.sqrt(mean_squared_error(y_val, y_pred_gbval))\n",
    "mse_gb = mean_squared_error(y_val, y_pred_gbval)\n",
    "mae_gb = mean_absolute_error(y_val, y_pred_gbval)\n",
    "r2_gb = r2_score(y_val, y_pred_gbval)\n",
    "adjr2_gb = 1-(1-r2_gb)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n",
    "\n",
    "\n",
    "resultsDf2 = pd.DataFrame({'Method':['Gradient Boost'], 'RMSE': rmse_gb , 'MSE': mse_gb , 'MAE': mae_gb ,\"R2\": r2_gb, 'Adj_R2': adjr2_gb})\n",
    "resultsDf2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "#AdaBoostRegressor"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "AdaBoost Regressor:\n",
      "\n",
      "Training Score: 0.9756740580889769\n",
      "Validate Score: 0.9755760352965355\n",
      "Root Mean Sq Error: 1.6088182614905226\n",
      "Mean Sq Error: 2.5882961985053874\n",
      "Mean Abs Error: 1.2815362314916658\n",
      "R2: 0.9755760352965355\n",
      "Adj R2: 0.975558164587151\n"
     ]
    }
   ],
   "source": [
    "ab = AdaBoostRegressor(n_estimators = 15, random_state = 1)\n",
    "ab.fit(X_train, y_train)\n",
    "ab_pred = ab.predict(X_val)\n",
    "\n",
    "print('AdaBoost Regressor:')\n",
    "print('')\n",
    "print('Training Score:', ab.score(X_train, y_train))\n",
    "print('Validate Score:', ab.score(X_val, y_val))\n",
    "\n",
    "\n",
    "\n",
    "rmse_ab = np.sqrt(mean_squared_error(y_val, ab_pred))\n",
    "mse_ab = mean_squared_error(y_val, ab_pred)\n",
    "mae_ab = mean_absolute_error(y_val, ab_pred)\n",
    "r2_ab = r2_score(y_val, ab_pred)\n",
    "adjr2_ab = 1-(1-r2_ab)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n",
    "\n",
    "print('Root Mean Sq Error:', rmse_ab)\n",
    "print('Mean Sq Error:', mse_ab)\n",
    "print('Mean Abs Error:', mae_ab)\n",
    "print('R2:', r2_ab)\n",
    "print('Adj R2:', adjr2_ab)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "#StackingRegressor\n",
    "estimators = [('bagging',BaggingRegressor()),('adaboost',AdaBoostRegressor()),('gradboost',GradientBoostingRegressor())]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Stacking Regressor:\n",
      "\n",
      "Training Score: 0.9999919056633009\n",
      "Validate Score: 0.9999834395774162\n",
      "Root Mean Sq Error: 0.04189233934102388\n",
      "Mean Sq Error: 0.0017549680954634973\n",
      "Mean Abs Error: 0.012078349992204121\n",
      "R2: 0.9999834395774162\n",
      "Adj R2: 0.9999834274603622\n"
     ]
    }
   ],
   "source": [
    "sr = StackingRegressor(estimators = estimators, final_estimator = RandomForestRegressor(n_estimators = 10, random_state = 1))\n",
    "sr.fit(X_train, y_train)\n",
    "sr_pred = sr.predict(X_val)\n",
    "\n",
    "print('Stacking Regressor:')\n",
    "print('')\n",
    "print('Training Score:', sr.score(X_train, y_train))\n",
    "print('Validate Score:', sr.score(X_val, y_val))\n",
    "\n",
    "\n",
    "\n",
    "rmse_sr = np.sqrt(mean_squared_error(y_val, sr_pred))\n",
    "mse_sr = mean_squared_error(y_val, sr_pred)\n",
    "mae_sr = mean_absolute_error(y_val, sr_pred)\n",
    "r2_sr = r2_score(y_val, sr_pred)\n",
    "adjr2_sr = 1-(1-r2_sr)*(len(y_val)-1)/(len(y_val)-len(X_val.columns)-1)\n",
    "\n",
    "print('Root Mean Sq Error:', rmse_sr)\n",
    "print('Mean Sq Error:', mse_sr)\n",
    "print('Mean Abs Error:', mae_sr)\n",
    "print('R2:', r2_sr)\n",
    "print('Adj R2:', adjr2_sr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [],
   "source": [
    "# creating a new dataframe of predictions, last column is true values\n",
    "columns = [y_pred_lr, y_pred_lasso, y_pred_ridge, y_pred_rf_val, y_pred_gbval, ab_pred, sr_pred, y_val]\n",
    "predictions = pd.DataFrame({'y_pred_lr': y_pred_lr, 'y_pred_lasso': y_pred_lasso,\n",
    "                            'y_pred_ridge': y_pred_ridge, 'y_pred_rf_val': y_pred_rf_val,\n",
    "                            'y_pred_gbval': y_pred_gbval, 'ab_pred': ab_pred, 'sr_pred': sr_pred,\n",
    "                            'y_val': y_val})                          "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>y_pred_lr</th>\n",
       "      <th>y_pred_lasso</th>\n",
       "      <th>y_pred_ridge</th>\n",
       "      <th>y_pred_rf_val</th>\n",
       "      <th>y_pred_gbval</th>\n",
       "      <th>ab_pred</th>\n",
       "      <th>sr_pred</th>\n",
       "      <th>y_val</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>47340</th>\n",
       "      <td>68.123206</td>\n",
       "      <td>68.116504</td>\n",
       "      <td>68.121723</td>\n",
       "      <td>67.917269</td>\n",
       "      <td>67.961004</td>\n",
       "      <td>69.025716</td>\n",
       "      <td>67.917153</td>\n",
       "      <td>67.923866</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>93858</th>\n",
       "      <td>78.008950</td>\n",
       "      <td>78.010202</td>\n",
       "      <td>78.013106</td>\n",
       "      <td>76.957500</td>\n",
       "      <td>76.943930</td>\n",
       "      <td>78.095026</td>\n",
       "      <td>76.961325</td>\n",
       "      <td>76.960355</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>201163</th>\n",
       "      <td>76.379754</td>\n",
       "      <td>76.365448</td>\n",
       "      <td>76.376156</td>\n",
       "      <td>76.087270</td>\n",
       "      <td>76.073863</td>\n",
       "      <td>75.922117</td>\n",
       "      <td>76.062999</td>\n",
       "      <td>76.086628</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>364634</th>\n",
       "      <td>74.927579</td>\n",
       "      <td>74.927992</td>\n",
       "      <td>74.928977</td>\n",
       "      <td>74.729004</td>\n",
       "      <td>74.723520</td>\n",
       "      <td>74.436961</td>\n",
       "      <td>74.747652</td>\n",
       "      <td>74.739285</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54177</th>\n",
       "      <td>96.590446</td>\n",
       "      <td>96.490199</td>\n",
       "      <td>96.585355</td>\n",
       "      <td>92.087715</td>\n",
       "      <td>92.004183</td>\n",
       "      <td>90.309333</td>\n",
       "      <td>92.114640</td>\n",
       "      <td>92.082130</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        y_pred_lr  y_pred_lasso  y_pred_ridge  y_pred_rf_val  y_pred_gbval  \\\n",
       "47340   68.123206     68.116504     68.121723      67.917269     67.961004   \n",
       "93858   78.008950     78.010202     78.013106      76.957500     76.943930   \n",
       "201163  76.379754     76.365448     76.376156      76.087270     76.073863   \n",
       "364634  74.927579     74.927992     74.928977      74.729004     74.723520   \n",
       "54177   96.590446     96.490199     96.585355      92.087715     92.004183   \n",
       "\n",
       "          ab_pred    sr_pred      y_val  \n",
       "47340   69.025716  67.917153  67.923866  \n",
       "93858   78.095026  76.961325  76.960355  \n",
       "201163  75.922117  76.062999  76.086628  \n",
       "364634  74.436961  74.747652  74.739285  \n",
       "54177   90.309333  92.114640  92.082130  "
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [],
   "source": [
    "# split the data\n",
    "X_predictions = predictions.drop(columns = 'y_val')\n",
    "y_predictions = predictions['y_val']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "X_train_stack, X_val_stack, y_train_stack, y_val_stack = train_test_split(X_predictions, y_predictions, test_size = .1, random_state = 1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Stacked Random Forest:\n",
      "\n",
      "Training Score: 0.9999966021614899\n",
      "Validate Score: 0.9999802659857451\n",
      "Root Mean Sq Error: 0.04600195318978735\n",
      "Mean Sq Error: 0.0021161796972753865\n",
      "Mean Abs Error: 0.00838958103206651\n",
      "R2: 0.9999802659857451\n",
      "Adj R2: 0.9999802284991508\n"
     ]
    }
   ],
   "source": [
    "\n",
    "stacked_rf=RandomForestRegressor(random_state=1, n_estimators=10)\n",
    "stacked_rf.fit(X_train_stack, y_train_stack)\n",
    "\n",
    "stacked_rf_pred = stacked_rf.predict(X_val_stack)\n",
    "\n",
    "print('Stacked Random Forest:')\n",
    "print('')\n",
    "print('Training Score:', stacked_rf.score(X_train_stack, y_train_stack))\n",
    "print('Validate Score:', stacked_rf.score(X_val_stack, y_val_stack))\n",
    "\n",
    "\n",
    "\n",
    "rmse_srf = np.sqrt(mean_squared_error(y_val_stack, stacked_rf_pred))\n",
    "mse_srf = mean_squared_error(y_val_stack, stacked_rf_pred)\n",
    "mae_srf = mean_absolute_error(y_val_stack, stacked_rf_pred)\n",
    "r2_srf = r2_score(y_val_stack, stacked_rf_pred)\n",
    "adjr2_srf = 1-(1-r2_srf)*(len(y_val_stack)-1)/(len(y_val_stack)-len(X_val_stack.columns)-1)\n",
    "\n",
    "print('Root Mean Sq Error:', rmse_srf)\n",
    "print('Mean Sq Error:', mse_srf)\n",
    "print('Mean Abs Error:', mae_srf)\n",
    "print('R2:', r2_srf)\n",
    "print('Adj R2:', adjr2_srf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 112,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Predicted HR</th>\n",
       "      <th>Actual HR</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>153685</th>\n",
       "      <td>80.850095</td>\n",
       "      <td>80.849998</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>219386</th>\n",
       "      <td>79.111344</td>\n",
       "      <td>79.112897</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>176059</th>\n",
       "      <td>69.159821</td>\n",
       "      <td>69.148603</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273714</th>\n",
       "      <td>82.375443</td>\n",
       "      <td>82.380744</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>141953</th>\n",
       "      <td>67.796638</td>\n",
       "      <td>67.815215</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>87362</th>\n",
       "      <td>68.936927</td>\n",
       "      <td>68.937307</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>253986</th>\n",
       "      <td>79.753597</td>\n",
       "      <td>79.744980</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>152942</th>\n",
       "      <td>74.346169</td>\n",
       "      <td>74.351216</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>216439</th>\n",
       "      <td>61.141498</td>\n",
       "      <td>61.154190</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>153321</th>\n",
       "      <td>85.044133</td>\n",
       "      <td>85.049572</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3693 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        Predicted HR  Actual HR\n",
       "153685     80.850095  80.849998\n",
       "219386     79.111344  79.112897\n",
       "176059     69.159821  69.148603\n",
       "273714     82.375443  82.380744\n",
       "141953     67.796638  67.815215\n",
       "...              ...        ...\n",
       "87362      68.936927  68.937307\n",
       "253986     79.753597  79.744980\n",
       "152942     74.346169  74.351216\n",
       "216439     61.141498  61.154190\n",
       "153321     85.044133  85.049572\n",
       "\n",
       "[3693 rows x 2 columns]"
      ]
     },
     "execution_count": 112,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# creating a new dataframe of predictions, last column is true values\n",
    "\n",
    "comp_predictions = pd.DataFrame({'Predicted HR': stacked_rf_pred, 'Actual HR': y_val_stack}) \n",
    "comp_predictions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Running models on Test Data "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [],
   "source": [
    "dftest = pd.DataFrame(sc.fit_transform(dftest), columns = dftest.columns)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>VLF</th>\n",
       "      <th>VLF_PCT</th>\n",
       "      <th>LF</th>\n",
       "      <th>LF_PCT</th>\n",
       "      <th>LF_NU</th>\n",
       "      <th>HF</th>\n",
       "      <th>HF_PCT</th>\n",
       "      <th>TP</th>\n",
       "      <th>LF_HF</th>\n",
       "      <th>HF_LF</th>\n",
       "      <th>SD2</th>\n",
       "      <th>sampen</th>\n",
       "      <th>higuci</th>\n",
       "      <th>MEAN_RR</th>\n",
       "      <th>MEDIAN_RR</th>\n",
       "      <th>RMSSD</th>\n",
       "      <th>SDRR_RMSSD</th>\n",
       "      <th>pNN25</th>\n",
       "      <th>pNN50</th>\n",
       "      <th>KURT</th>\n",
       "      <th>SKEW</th>\n",
       "      <th>MEAN_REL_RR</th>\n",
       "      <th>MEDIAN_REL_RR</th>\n",
       "      <th>SDRR_REL_RR</th>\n",
       "      <th>RMSSD_REL_RR</th>\n",
       "      <th>SDRR_RMSSD_REL_RR</th>\n",
       "      <th>condition2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-0.180651</td>\n",
       "      <td>0.729908</td>\n",
       "      <td>-0.655013</td>\n",
       "      <td>-0.670703</td>\n",
       "      <td>0.948685</td>\n",
       "      <td>-0.800531</td>\n",
       "      <td>-0.850120</td>\n",
       "      <td>-0.385830</td>\n",
       "      <td>0.199605</td>\n",
       "      <td>-0.880406</td>\n",
       "      <td>-0.346849</td>\n",
       "      <td>0.704120</td>\n",
       "      <td>-1.311071</td>\n",
       "      <td>0.705741</td>\n",
       "      <td>0.734375</td>\n",
       "      <td>-0.770174</td>\n",
       "      <td>-0.079551</td>\n",
       "      <td>-0.720737</td>\n",
       "      <td>-0.738994</td>\n",
       "      <td>-0.670233</td>\n",
       "      <td>-0.400379</td>\n",
       "      <td>0.099959</td>\n",
       "      <td>0.874270</td>\n",
       "      <td>-1.032491</td>\n",
       "      <td>-1.315621</td>\n",
       "      <td>2.165489</td>\n",
       "      <td>-1.798140</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-0.896694</td>\n",
       "      <td>-2.252038</td>\n",
       "      <td>1.047996</td>\n",
       "      <td>2.340531</td>\n",
       "      <td>0.468689</td>\n",
       "      <td>0.015174</td>\n",
       "      <td>0.128417</td>\n",
       "      <td>-0.531271</td>\n",
       "      <td>-0.212028</td>\n",
       "      <td>-0.464945</td>\n",
       "      <td>-0.704126</td>\n",
       "      <td>0.590695</td>\n",
       "      <td>1.727921</td>\n",
       "      <td>-0.239461</td>\n",
       "      <td>-0.195296</td>\n",
       "      <td>1.345525</td>\n",
       "      <td>-0.907455</td>\n",
       "      <td>1.780555</td>\n",
       "      <td>-0.333628</td>\n",
       "      <td>-0.309826</td>\n",
       "      <td>-0.139153</td>\n",
       "      <td>0.249628</td>\n",
       "      <td>-2.637482</td>\n",
       "      <td>1.208364</td>\n",
       "      <td>1.088565</td>\n",
       "      <td>-0.534311</td>\n",
       "      <td>0.832367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-0.052107</td>\n",
       "      <td>0.689836</td>\n",
       "      <td>-0.508612</td>\n",
       "      <td>-0.653765</td>\n",
       "      <td>0.546740</td>\n",
       "      <td>-0.543970</td>\n",
       "      <td>-0.620535</td>\n",
       "      <td>-0.214574</td>\n",
       "      <td>-0.195937</td>\n",
       "      <td>-0.533641</td>\n",
       "      <td>-0.262176</td>\n",
       "      <td>-0.058029</td>\n",
       "      <td>0.704495</td>\n",
       "      <td>0.240359</td>\n",
       "      <td>0.394666</td>\n",
       "      <td>-0.274422</td>\n",
       "      <td>-0.187973</td>\n",
       "      <td>-0.341509</td>\n",
       "      <td>-0.333628</td>\n",
       "      <td>-0.406093</td>\n",
       "      <td>-0.914328</td>\n",
       "      <td>0.064373</td>\n",
       "      <td>0.341957</td>\n",
       "      <td>-0.555479</td>\n",
       "      <td>-0.396722</td>\n",
       "      <td>-0.270132</td>\n",
       "      <td>0.832367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1.961771</td>\n",
       "      <td>1.564399</td>\n",
       "      <td>-0.616444</td>\n",
       "      <td>-1.542476</td>\n",
       "      <td>0.787696</td>\n",
       "      <td>-0.709407</td>\n",
       "      <td>-0.856399</td>\n",
       "      <td>1.649698</td>\n",
       "      <td>-0.091277</td>\n",
       "      <td>-0.742907</td>\n",
       "      <td>1.369913</td>\n",
       "      <td>0.083416</td>\n",
       "      <td>-1.560276</td>\n",
       "      <td>1.541415</td>\n",
       "      <td>1.188118</td>\n",
       "      <td>0.354574</td>\n",
       "      <td>1.091530</td>\n",
       "      <td>0.110338</td>\n",
       "      <td>1.017593</td>\n",
       "      <td>-0.748444</td>\n",
       "      <td>0.636937</td>\n",
       "      <td>-1.469479</td>\n",
       "      <td>0.004734</td>\n",
       "      <td>-0.307937</td>\n",
       "      <td>-0.541235</td>\n",
       "      <td>0.575273</td>\n",
       "      <td>-1.798140</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-0.678566</td>\n",
       "      <td>0.358449</td>\n",
       "      <td>-0.993948</td>\n",
       "      <td>-0.426080</td>\n",
       "      <td>-0.921862</td>\n",
       "      <td>-0.123201</td>\n",
       "      <td>0.472156</td>\n",
       "      <td>-0.941615</td>\n",
       "      <td>-0.287959</td>\n",
       "      <td>0.839242</td>\n",
       "      <td>-0.745989</td>\n",
       "      <td>0.024683</td>\n",
       "      <td>1.118718</td>\n",
       "      <td>-0.581154</td>\n",
       "      <td>-0.478939</td>\n",
       "      <td>-1.139501</td>\n",
       "      <td>-0.456370</td>\n",
       "      <td>-0.954729</td>\n",
       "      <td>-0.806555</td>\n",
       "      <td>0.679582</td>\n",
       "      <td>-0.072031</td>\n",
       "      <td>0.331337</td>\n",
       "      <td>0.533948</td>\n",
       "      <td>-0.932272</td>\n",
       "      <td>-0.483581</td>\n",
       "      <td>-0.748784</td>\n",
       "      <td>-0.482887</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        VLF   VLF_PCT        LF    LF_PCT     LF_NU        HF    HF_PCT  \\\n",
       "0 -0.180651  0.729908 -0.655013 -0.670703  0.948685 -0.800531 -0.850120   \n",
       "1 -0.896694 -2.252038  1.047996  2.340531  0.468689  0.015174  0.128417   \n",
       "2 -0.052107  0.689836 -0.508612 -0.653765  0.546740 -0.543970 -0.620535   \n",
       "3  1.961771  1.564399 -0.616444 -1.542476  0.787696 -0.709407 -0.856399   \n",
       "4 -0.678566  0.358449 -0.993948 -0.426080 -0.921862 -0.123201  0.472156   \n",
       "\n",
       "         TP     LF_HF     HF_LF       SD2    sampen    higuci   MEAN_RR  \\\n",
       "0 -0.385830  0.199605 -0.880406 -0.346849  0.704120 -1.311071  0.705741   \n",
       "1 -0.531271 -0.212028 -0.464945 -0.704126  0.590695  1.727921 -0.239461   \n",
       "2 -0.214574 -0.195937 -0.533641 -0.262176 -0.058029  0.704495  0.240359   \n",
       "3  1.649698 -0.091277 -0.742907  1.369913  0.083416 -1.560276  1.541415   \n",
       "4 -0.941615 -0.287959  0.839242 -0.745989  0.024683  1.118718 -0.581154   \n",
       "\n",
       "   MEDIAN_RR     RMSSD  SDRR_RMSSD     pNN25     pNN50      KURT      SKEW  \\\n",
       "0   0.734375 -0.770174   -0.079551 -0.720737 -0.738994 -0.670233 -0.400379   \n",
       "1  -0.195296  1.345525   -0.907455  1.780555 -0.333628 -0.309826 -0.139153   \n",
       "2   0.394666 -0.274422   -0.187973 -0.341509 -0.333628 -0.406093 -0.914328   \n",
       "3   1.188118  0.354574    1.091530  0.110338  1.017593 -0.748444  0.636937   \n",
       "4  -0.478939 -1.139501   -0.456370 -0.954729 -0.806555  0.679582 -0.072031   \n",
       "\n",
       "   MEAN_REL_RR  MEDIAN_REL_RR  SDRR_REL_RR  RMSSD_REL_RR  SDRR_RMSSD_REL_RR  \\\n",
       "0     0.099959       0.874270    -1.032491     -1.315621           2.165489   \n",
       "1     0.249628      -2.637482     1.208364      1.088565          -0.534311   \n",
       "2     0.064373       0.341957    -0.555479     -0.396722          -0.270132   \n",
       "3    -1.469479       0.004734    -0.307937     -0.541235           0.575273   \n",
       "4     0.331337       0.533948    -0.932272     -0.483581          -0.748784   \n",
       "\n",
       "   condition2  \n",
       "0   -1.798140  \n",
       "1    0.832367  \n",
       "2    0.832367  \n",
       "3   -1.798140  \n",
       "4   -0.482887  "
      ]
     },
     "execution_count": 62,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dftest.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([68.12320605, 78.00894965, 76.37975375, ..., 79.77525114,\n",
       "       72.20633914, 78.30965516])"
      ]
     },
     "execution_count": 63,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train linear predictions\n",
    "y_pred_lr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-2.97390795e+06, -2.55615037e+09,  1.38175752e+08, ...,\n",
       "       -4.10474304e+09,  2.74773720e+09,  1.05313137e+09])"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test linear predictions\n",
    "stacked_lr_pred = lr.predict(dftest)\n",
    "stacked_lr_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([68.11650427, 78.01020162, 76.36544841, ..., 79.77327981,\n",
       "       72.21687329, 78.32209409])"
      ]
     },
     "execution_count": 65,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train lasso predictions\n",
    "y_pred_lasso"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([65.53416161, 74.38782374, 69.21397052, ..., 60.68459432,\n",
       "       79.17132983, 79.95797628])"
      ]
     },
     "execution_count": 66,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test lasso predictions\n",
    "stacked_lasso_pred = lasso.predict(dftest)\n",
    "stacked_lasso_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([68.12172287, 78.01310619, 76.37615561, ..., 79.77349023,\n",
       "       72.20344123, 78.3108161 ])"
      ]
     },
     "execution_count": 67,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train ridge predictions\n",
    "y_pred_ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([65.52464556, 74.39682423, 69.22602176, ..., 60.69292488,\n",
       "       79.16308824, 79.95245058])"
      ]
     },
     "execution_count": 68,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test ridge predictions\n",
    "stacked_ridge_pred = ridge.predict(dftest)\n",
    "stacked_ridge_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([67.91726907, 76.95750044, 76.08726983, ..., 79.05699492,\n",
       "       71.65509243, 78.67101527])"
      ]
     },
     "execution_count": 69,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train random forest predictions\n",
    "y_pred_rf_val"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 113,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([64.73018515, 73.81165717, 69.2095931 , ..., 60.94626282,\n",
       "       77.42830137, 76.62098638])"
      ]
     },
     "execution_count": 113,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test random forest prediction\n",
    "stacked_rf_pred = rf.predict(dftest)\n",
    "stacked_rf_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([67.96100411, 76.94393044, 76.0738628 , ..., 79.03444023,\n",
       "       71.65660266, 78.64610007])"
      ]
     },
     "execution_count": 71,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train gradient boost predictions\n",
    "y_pred_gbval "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([64.7561422 , 73.81299698, 69.23239813, ..., 61.01824369,\n",
       "       77.41908973, 76.60754706])"
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test gradient boost predictions\n",
    "stacked_gb_pred = gb.predict(dftest)\n",
    "stacked_gb_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([69.02571595, 78.09502551, 75.92211663, ..., 78.09502551,\n",
       "       73.85115038, 78.24857146])"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train adaboost predictions\n",
    "ab_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([66.74921082, 73.85115038, 69.77614414, ..., 63.46033118,\n",
       "       78.24857146, 75.92211663])"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test adaboost predictions\n",
    "stacked_ab_pred = ab.predict(dftest)\n",
    "stacked_ab_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([67.91715341, 76.96132527, 76.06299855, ..., 79.05662618,\n",
       "       71.65729541, 78.67420658])"
      ]
     },
     "execution_count": 75,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# train stacking regressor predictions\n",
    "sr_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([64.73246265, 73.8089425 , 69.21373628, ..., 60.96377221,\n",
       "       77.4361271 , 76.59757333])"
      ]
     },
     "execution_count": 76,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test stacking regressor predictions\n",
    "stacked_sr_pred = sr.predict(dftest)\n",
    "stacked_sr_pred"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "# creating a new dataframe of stacked predictions\n",
    "\n",
    "stacked_predictions = pd.DataFrame({'stacked_lr_pred': stacked_lr_pred, 'stacked_lasso_pred': stacked_lasso_pred,\n",
    "                            'stacked_ridge_pred': stacked_ridge_pred, 'stacked_rf_pred': stacked_rf_pred,\n",
    "                            'stacked_gb_pred': stacked_gb_pred, 'stacked_ab_pred': stacked_ab_pred, \n",
    "                            'stacked_sr_pred': stacked_sr_pred})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>stacked_lr_pred</th>\n",
       "      <th>stacked_lasso_pred</th>\n",
       "      <th>stacked_ridge_pred</th>\n",
       "      <th>stacked_rf_pred</th>\n",
       "      <th>stacked_gb_pred</th>\n",
       "      <th>stacked_ab_pred</th>\n",
       "      <th>stacked_sr_pred</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-2.973908e+06</td>\n",
       "      <td>65.534162</td>\n",
       "      <td>65.524646</td>\n",
       "      <td>64.730185</td>\n",
       "      <td>64.756142</td>\n",
       "      <td>66.749211</td>\n",
       "      <td>64.732463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-2.556150e+09</td>\n",
       "      <td>74.387824</td>\n",
       "      <td>74.396824</td>\n",
       "      <td>73.811657</td>\n",
       "      <td>73.812997</td>\n",
       "      <td>73.851150</td>\n",
       "      <td>73.808942</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1.381758e+08</td>\n",
       "      <td>69.213971</td>\n",
       "      <td>69.226022</td>\n",
       "      <td>69.209593</td>\n",
       "      <td>69.232398</td>\n",
       "      <td>69.776144</td>\n",
       "      <td>69.213736</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9.992321e+08</td>\n",
       "      <td>60.597939</td>\n",
       "      <td>60.607264</td>\n",
       "      <td>60.167199</td>\n",
       "      <td>60.232024</td>\n",
       "      <td>61.429099</td>\n",
       "      <td>60.197427</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1.917144e+09</td>\n",
       "      <td>77.725435</td>\n",
       "      <td>77.714571</td>\n",
       "      <td>77.849928</td>\n",
       "      <td>77.832406</td>\n",
       "      <td>78.248571</td>\n",
       "      <td>77.847203</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   stacked_lr_pred  stacked_lasso_pred  stacked_ridge_pred  stacked_rf_pred  \\\n",
       "0    -2.973908e+06           65.534162           65.524646        64.730185   \n",
       "1    -2.556150e+09           74.387824           74.396824        73.811657   \n",
       "2     1.381758e+08           69.213971           69.226022        69.209593   \n",
       "3     9.992321e+08           60.597939           60.607264        60.167199   \n",
       "4     1.917144e+09           77.725435           77.714571        77.849928   \n",
       "\n",
       "   stacked_gb_pred  stacked_ab_pred  stacked_sr_pred  \n",
       "0        64.756142        66.749211        64.732463  \n",
       "1        73.812997        73.851150        73.808942  \n",
       "2        69.232398        69.776144        69.213736  \n",
       "3        60.232024        61.429099        60.197427  \n",
       "4        77.832406        78.248571        77.847203  "
      ]
     },
     "execution_count": 78,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stacked_predictions.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([64.72395885, 73.80934743, 69.20521636, ..., 60.95425468,\n",
       "       77.43361532, 76.62398021])"
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "final_prediction = stacked_rf.predict(stacked_predictions)\n",
    "final_prediction"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 115,
   "metadata": {},
   "outputs": [],
   "source": [
    "stacked_predictions_with_finals = pd.DataFrame({'stacked_lr_pred': stacked_lr_pred, 'stacked_lasso_pred': stacked_lasso_pred,\n",
    "                            'stacked_ridge_pred': stacked_ridge_pred, 'stacked_rf_pred': stacked_rf_pred,\n",
    "                            'stacked_gb_pred': stacked_gb_pred, 'stacked_ab_pred': stacked_ab_pred, \n",
    "                            'stacked_sr_pred': stacked_sr_pred, 'Prediction': final_prediction})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {
    "scrolled": true
   },
   "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>stacked_lr_pred</th>\n",
       "      <th>stacked_lasso_pred</th>\n",
       "      <th>stacked_ridge_pred</th>\n",
       "      <th>stacked_rf_pred</th>\n",
       "      <th>stacked_gb_pred</th>\n",
       "      <th>stacked_ab_pred</th>\n",
       "      <th>stacked_sr_pred</th>\n",
       "      <th>stacked_final</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-2.973908e+06</td>\n",
       "      <td>65.534162</td>\n",
       "      <td>65.524646</td>\n",
       "      <td>64.730185</td>\n",
       "      <td>64.756142</td>\n",
       "      <td>66.749211</td>\n",
       "      <td>64.732463</td>\n",
       "      <td>64.723959</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-2.556150e+09</td>\n",
       "      <td>74.387824</td>\n",
       "      <td>74.396824</td>\n",
       "      <td>73.811657</td>\n",
       "      <td>73.812997</td>\n",
       "      <td>73.851150</td>\n",
       "      <td>73.808942</td>\n",
       "      <td>73.809347</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1.381758e+08</td>\n",
       "      <td>69.213971</td>\n",
       "      <td>69.226022</td>\n",
       "      <td>69.209593</td>\n",
       "      <td>69.232398</td>\n",
       "      <td>69.776144</td>\n",
       "      <td>69.213736</td>\n",
       "      <td>69.205216</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9.992321e+08</td>\n",
       "      <td>60.597939</td>\n",
       "      <td>60.607264</td>\n",
       "      <td>60.167199</td>\n",
       "      <td>60.232024</td>\n",
       "      <td>61.429099</td>\n",
       "      <td>60.197427</td>\n",
       "      <td>60.177429</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1.917144e+09</td>\n",
       "      <td>77.725435</td>\n",
       "      <td>77.714571</td>\n",
       "      <td>77.849928</td>\n",
       "      <td>77.832406</td>\n",
       "      <td>78.248571</td>\n",
       "      <td>77.847203</td>\n",
       "      <td>77.848124</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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41028</th>\n",
       "      <td>8.125249e+08</td>\n",
       "      <td>80.327466</td>\n",
       "      <td>80.333911</td>\n",
       "      <td>81.035984</td>\n",
       "      <td>80.986996</td>\n",
       "      <td>83.057630</td>\n",
       "      <td>81.026176</td>\n",
       "      <td>81.027561</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41029</th>\n",
       "      <td>1.668889e+09</td>\n",
       "      <td>94.405159</td>\n",
       "      <td>94.405721</td>\n",
       "      <td>98.795314</td>\n",
       "      <td>98.730782</td>\n",
       "      <td>100.242714</td>\n",
       "      <td>98.806797</td>\n",
       "      <td>98.806229</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41030</th>\n",
       "      <td>-4.104743e+09</td>\n",
       "      <td>60.684594</td>\n",
       "      <td>60.692925</td>\n",
       "      <td>60.946263</td>\n",
       "      <td>61.018244</td>\n",
       "      <td>63.460331</td>\n",
       "      <td>60.963772</td>\n",
       "      <td>60.954255</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41031</th>\n",
       "      <td>2.747737e+09</td>\n",
       "      <td>79.171330</td>\n",
       "      <td>79.163088</td>\n",
       "      <td>77.428301</td>\n",
       "      <td>77.419090</td>\n",
       "      <td>78.248571</td>\n",
       "      <td>77.436127</td>\n",
       "      <td>77.433615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41032</th>\n",
       "      <td>1.053131e+09</td>\n",
       "      <td>79.957976</td>\n",
       "      <td>79.952451</td>\n",
       "      <td>76.620986</td>\n",
       "      <td>76.607547</td>\n",
       "      <td>75.922117</td>\n",
       "      <td>76.597573</td>\n",
       "      <td>76.623980</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>41033 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       stacked_lr_pred  stacked_lasso_pred  stacked_ridge_pred  \\\n",
       "0        -2.973908e+06           65.534162           65.524646   \n",
       "1        -2.556150e+09           74.387824           74.396824   \n",
       "2         1.381758e+08           69.213971           69.226022   \n",
       "3         9.992321e+08           60.597939           60.607264   \n",
       "4         1.917144e+09           77.725435           77.714571   \n",
       "...                ...                 ...                 ...   \n",
       "41028     8.125249e+08           80.327466           80.333911   \n",
       "41029     1.668889e+09           94.405159           94.405721   \n",
       "41030    -4.104743e+09           60.684594           60.692925   \n",
       "41031     2.747737e+09           79.171330           79.163088   \n",
       "41032     1.053131e+09           79.957976           79.952451   \n",
       "\n",
       "       stacked_rf_pred  stacked_gb_pred  stacked_ab_pred  stacked_sr_pred  \\\n",
       "0            64.730185        64.756142        66.749211        64.732463   \n",
       "1            73.811657        73.812997        73.851150        73.808942   \n",
       "2            69.209593        69.232398        69.776144        69.213736   \n",
       "3            60.167199        60.232024        61.429099        60.197427   \n",
       "4            77.849928        77.832406        78.248571        77.847203   \n",
       "...                ...              ...              ...              ...   \n",
       "41028        81.035984        80.986996        83.057630        81.026176   \n",
       "41029        98.795314        98.730782       100.242714        98.806797   \n",
       "41030        60.946263        61.018244        63.460331        60.963772   \n",
       "41031        77.428301        77.419090        78.248571        77.436127   \n",
       "41032        76.620986        76.607547        75.922117        76.597573   \n",
       "\n",
       "       stacked_final  \n",
       "0          64.723959  \n",
       "1          73.809347  \n",
       "2          69.205216  \n",
       "3          60.177429  \n",
       "4          77.848124  \n",
       "...              ...  \n",
       "41028      81.027561  \n",
       "41029      98.806229  \n",
       "41030      60.954255  \n",
       "41031      77.433615  \n",
       "41032      76.623980  \n",
       "\n",
       "[41033 rows x 8 columns]"
      ]
     },
     "execution_count": 82,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stacked_predictions_with_finals"
   ]
  },
  {
   "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
}