07_SimpleStockLinearRegression.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "source": [
        "# Simple Stock Linear Regression"
      ],
      "metadata": {}
    },
    {
      "cell_type": "code",
      "source": [
        "library(quantmod)"
      ],
      "outputs": [
        {
          "output_type": "stream",
          "name": "stderr",
          "text": [
            "Warning message:\n",
            "\"package 'quantmod' was built under R version 3.6.3\"Loading required package: xts\n",
            "Warning message:\n",
            "\"package 'xts' was built under R version 3.6.3\"Loading required package: zoo\n",
            "Warning message:\n",
            "\"package 'zoo' was built under R version 3.6.3\"\n",
            "Attaching package: 'zoo'\n",
            "\n",
            "The following objects are masked from 'package:base':\n",
            "\n",
            "    as.Date, as.Date.numeric\n",
            "\n",
            "Loading required package: TTR\n",
            "Warning message:\n",
            "\"package 'TTR' was built under R version 3.6.3\"Registered S3 method overwritten by 'quantmod':\n",
            "  method            from\n",
            "  as.zoo.data.frame zoo \n"
          ]
        }
      ],
      "execution_count": 1,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:00.881Z",
          "iopub.execute_input": "2021-07-08T01:47:00.885Z",
          "iopub.status.idle": "2021-07-08T01:47:00.964Z",
          "shell.execute_reply": "2021-07-08T01:47:00.959Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "getSymbols('AAPL', src = 'yahoo', from = as.Date('2016-01-01'), to = as.Date('2019-01-01'))"
      ],
      "outputs": [
        {
          "output_type": "stream",
          "name": "stderr",
          "text": [
            "'getSymbols' currently uses auto.assign=TRUE by default, but will\n",
            "use auto.assign=FALSE in 0.5-0. You will still be able to use\n",
            "'loadSymbols' to automatically load data. getOption(\"getSymbols.env\")\n",
            "and getOption(\"getSymbols.auto.assign\") will still be checked for\n",
            "alternate defaults.\n",
            "\n",
            "This message is shown once per session and may be disabled by setting \n",
            "options(\"getSymbols.warning4.0\"=FALSE). See ?getSymbols for details.\n",
            "\n"
          ]
        },
        {
          "output_type": "display_data",
          "data": {
            "text/html": "'AAPL'",
            "text/markdown": "'AAPL'",
            "text/latex": "'AAPL'",
            "text/plain": "[1] \"AAPL\""
          },
          "metadata": {}
        }
      ],
      "execution_count": 2,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:00.969Z",
          "iopub.execute_input": "2021-07-08T01:47:00.980Z",
          "shell.execute_reply": "2021-07-08T01:47:02.152Z",
          "iopub.status.idle": "2021-07-08T01:47:02.158Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "head(AAPL)"
      ],
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": "           AAPL.Open AAPL.High AAPL.Low AAPL.Close AAPL.Volume AAPL.Adjusted\n2016-01-04   25.6525   26.3425  25.5000    26.3375   270597600      24.32322\n2016-01-05   26.4375   26.4625  25.6025    25.6775   223164000      23.71370\n2016-01-06   25.1400   25.5925  24.9675    25.1750   273829600      23.24963\n2016-01-07   24.6700   25.0325  24.1075    24.1125   324377600      22.26838\n2016-01-08   24.6375   24.7775  24.1900    24.2400   283192000      22.38613\n2016-01-11   24.7425   24.7650  24.3350    24.6325   198957600      22.74862"
          },
          "metadata": {}
        }
      ],
      "execution_count": 3,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.163Z",
          "iopub.execute_input": "2021-07-08T01:47:02.166Z",
          "iopub.status.idle": "2021-07-08T01:47:02.174Z",
          "shell.execute_reply": "2021-07-08T01:47:02.207Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "tail(AAPL)"
      ],
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": "           AAPL.Open AAPL.High AAPL.Low AAPL.Close AAPL.Volume AAPL.Adjusted\n2018-12-21   39.2150   39.5400  37.4075    37.6825   382978400      36.68959\n2018-12-24   37.0375   37.8875  36.6475    36.7075   148676800      35.74028\n2018-12-26   37.0750   39.3075  36.6800    39.2925   234330000      38.25717\n2018-12-27   38.9600   39.1925  37.5175    39.0375   212468400      38.00889\n2018-12-28   39.3750   39.6300  38.6375    39.0575   169165600      38.02837\n2018-12-31   39.6325   39.8400  39.1200    39.4350   140014000      38.39592"
          },
          "metadata": {}
        }
      ],
      "execution_count": 4,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.178Z",
          "iopub.execute_input": "2021-07-08T01:47:02.181Z",
          "iopub.status.idle": "2021-07-08T01:47:02.189Z",
          "shell.execute_reply": "2021-07-08T01:47:02.210Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "Stock <- AAPL\n",
        "Stock = data.frame(AAPL)\n",
        "head(Stock)"
      ],
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/html": "<table class=\"dataframe\">\n<caption>A data.frame: 6 × 6</caption>\n<thead>\n\t<tr><th></th><th scope=col>AAPL.Open</th><th scope=col>AAPL.High</th><th scope=col>AAPL.Low</th><th scope=col>AAPL.Close</th><th scope=col>AAPL.Volume</th><th scope=col>AAPL.Adjusted</th></tr>\n\t<tr><th></th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n</thead>\n<tbody>\n\t<tr><th scope=row>2016-01-04</th><td>25.6525</td><td>26.3425</td><td>25.5000</td><td>26.3375</td><td>270597600</td><td>24.32322</td></tr>\n\t<tr><th scope=row>2016-01-05</th><td>26.4375</td><td>26.4625</td><td>25.6025</td><td>25.6775</td><td>223164000</td><td>23.71370</td></tr>\n\t<tr><th scope=row>2016-01-06</th><td>25.1400</td><td>25.5925</td><td>24.9675</td><td>25.1750</td><td>273829600</td><td>23.24963</td></tr>\n\t<tr><th scope=row>2016-01-07</th><td>24.6700</td><td>25.0325</td><td>24.1075</td><td>24.1125</td><td>324377600</td><td>22.26838</td></tr>\n\t<tr><th scope=row>2016-01-08</th><td>24.6375</td><td>24.7775</td><td>24.1900</td><td>24.2400</td><td>283192000</td><td>22.38613</td></tr>\n\t<tr><th scope=row>2016-01-11</th><td>24.7425</td><td>24.7650</td><td>24.3350</td><td>24.6325</td><td>198957600</td><td>22.74862</td></tr>\n</tbody>\n</table>\n",
            "text/markdown": "\nA data.frame: 6 × 6\n\n| <!--/--> | AAPL.Open &lt;dbl&gt; | AAPL.High &lt;dbl&gt; | AAPL.Low &lt;dbl&gt; | AAPL.Close &lt;dbl&gt; | AAPL.Volume &lt;dbl&gt; | AAPL.Adjusted &lt;dbl&gt; |\n|---|---|---|---|---|---|---|\n| 2016-01-04 | 25.6525 | 26.3425 | 25.5000 | 26.3375 | 270597600 | 24.32322 |\n| 2016-01-05 | 26.4375 | 26.4625 | 25.6025 | 25.6775 | 223164000 | 23.71370 |\n| 2016-01-06 | 25.1400 | 25.5925 | 24.9675 | 25.1750 | 273829600 | 23.24963 |\n| 2016-01-07 | 24.6700 | 25.0325 | 24.1075 | 24.1125 | 324377600 | 22.26838 |\n| 2016-01-08 | 24.6375 | 24.7775 | 24.1900 | 24.2400 | 283192000 | 22.38613 |\n| 2016-01-11 | 24.7425 | 24.7650 | 24.3350 | 24.6325 | 198957600 | 22.74862 |\n\n",
            "text/latex": "A data.frame: 6 × 6\n\\begin{tabular}{r|llllll}\n  & AAPL.Open & AAPL.High & AAPL.Low & AAPL.Close & AAPL.Volume & AAPL.Adjusted\\\\\n  & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl>\\\\\n\\hline\n\t2016-01-04 & 25.6525 & 26.3425 & 25.5000 & 26.3375 & 270597600 & 24.32322\\\\\n\t2016-01-05 & 26.4375 & 26.4625 & 25.6025 & 25.6775 & 223164000 & 23.71370\\\\\n\t2016-01-06 & 25.1400 & 25.5925 & 24.9675 & 25.1750 & 273829600 & 23.24963\\\\\n\t2016-01-07 & 24.6700 & 25.0325 & 24.1075 & 24.1125 & 324377600 & 22.26838\\\\\n\t2016-01-08 & 24.6375 & 24.7775 & 24.1900 & 24.2400 & 283192000 & 22.38613\\\\\n\t2016-01-11 & 24.7425 & 24.7650 & 24.3350 & 24.6325 & 198957600 & 22.74862\\\\\n\\end{tabular}\n",
            "text/plain": "           AAPL.Open AAPL.High AAPL.Low AAPL.Close AAPL.Volume AAPL.Adjusted\n2016-01-04 25.6525   26.3425   25.5000  26.3375    270597600   24.32322     \n2016-01-05 26.4375   26.4625   25.6025  25.6775    223164000   23.71370     \n2016-01-06 25.1400   25.5925   24.9675  25.1750    273829600   23.24963     \n2016-01-07 24.6700   25.0325   24.1075  24.1125    324377600   22.26838     \n2016-01-08 24.6375   24.7775   24.1900  24.2400    283192000   22.38613     \n2016-01-11 24.7425   24.7650   24.3350  24.6325    198957600   22.74862     "
          },
          "metadata": {}
        }
      ],
      "execution_count": 5,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.194Z",
          "iopub.execute_input": "2021-07-08T01:47:02.196Z",
          "iopub.status.idle": "2021-07-08T01:47:02.250Z",
          "shell.execute_reply": "2021-07-08T01:47:02.246Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "x <- Stock$AAPL.Open\n",
        "y <- Stock$AAPL.Adjusted"
      ],
      "outputs": [],
      "execution_count": 6,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.256Z",
          "iopub.execute_input": "2021-07-08T01:47:02.259Z",
          "iopub.status.idle": "2021-07-08T01:47:02.263Z",
          "shell.execute_reply": "2021-07-08T01:47:02.296Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "relationship <- lm(y~x)\n",
        "print(relationship)"
      ],
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "\n",
            "Call:\n",
            "lm(formula = y ~ x)\n",
            "\n",
            "Coefficients:\n",
            "(Intercept)            x  \n",
            "     -1.717        1.001  \n",
            "\n"
          ]
        }
      ],
      "execution_count": 7,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.268Z",
          "iopub.execute_input": "2021-07-08T01:47:02.270Z",
          "iopub.status.idle": "2021-07-08T01:47:02.278Z",
          "shell.execute_reply": "2021-07-08T01:47:02.299Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "# Find Closing Price with Low Price of 180\n",
        "a <- data.frame(x = 180)\n",
        "result <- predict(relationship,a)\n",
        "print(result)"
      ],
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "       1 \n",
            "178.5278 \n"
          ]
        }
      ],
      "execution_count": 8,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.282Z",
          "iopub.execute_input": "2021-07-08T01:47:02.285Z",
          "iopub.status.idle": "2021-07-08T01:47:02.304Z",
          "shell.execute_reply": "2021-07-08T01:47:02.301Z"
        }
      }
    },
    {
      "cell_type": "code",
      "source": [
        "# Plot Chart\n",
        "plot(y,x, col = 'blue', main = \"Low & Closing Price Regression\", abline(lm(x~y)),cex=0.5, pch=16, xlab = 'Low Price', ylab = 'Adj Closing Price')\n",
        "\n"
      ],
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": "Plot with title \"Low & Closing Price Regression\"",
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAM1BMVEUAAAAAAP9NTU1oaGh8\nfHyMjIyampqnp6eysrK9vb3Hx8fQ0NDZ2dnh4eHp6enw8PD////UNI3wAAAACXBIWXMAABJ0\nAAASdAHeZh94AAAgAElEQVR4nO3diXayvBpA4YgiDlW4/6utzEkYHHgTEtjPWuf8tlVIK/tj\nVFUBYDG19gCALSAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIE\nEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIE\nEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCCk199A/fpXyLOjUsds+IO/\n8+sHSXr7eAbfjEE1kvT+20TaCaj0+vFMvxmd/ETDt8tf2vLzU/9MmiU6t36Qdsv634cz+CWk\nl8z+wbcTOH081y9GJz7NCOzyl7b8/NR3vZzN75+0JfXx2Qx+DEn9/TIRfQLi6yRC2q2fn/rX\nA18bbzf78a++kutrJfW8lltPS2YwOdvqv/dkkPB3E8iz14ap3Lh2jZBGlvP7uVzLVDsgx/of\n/fz1jXL77WkseUn9zVc3+qMfr6+fzc3j5anPQJvwa5qXcs3V7Ec1dyn/c399+/yo7/N8PeB0\ntUbYffWob73+/3kst/L6QF7DPrWzyRKVZM/xX7h/hHYnc6bd1I07GYMf+02sX9b+xTaHkIYh\ntVtm5crkUu+H3OuVT7n2ufR3fP2Dfno+T8b3qu9a20vtDPQJdztY9W5Kv9Rm/SZh8dfcYyKk\nogvpWE2m+UE75Uz/4m9mAuadrJl2U9fvZAx+9Dcxf1n7F9seQhqE1O35lMtAswqq/nEtqqVD\n+7c9fy1BSbestMz7aDMwJlxO8pVmfqqz65faRrXNlvRfj433Ua8L6zvcuh90j7prXyQjEyhX\nPWlh38maaTd1/U7G4Ed/E/OXtX+x7SEkO6Ry5fPaw8kv9ZJ4rDbf2sXK2qeolpZkYv1jf8ea\nsKp7y+tJ9kttcq+7Lar1X/nVLRkPqd1HKtcFef+D155Z8qiW62P1xetneb2saxPo/BXWneyZ\ntlM37mQMfuw3sX5Z6xfboI3+Wl8xn9xzu2WWVYvppVxw/ppl7m5sxb12HaqQxtc/g+9YEy4X\n1PPdelCz2OX1l2n9lX04Q+/g0T+onUjzqLzaP0vr3bjCWG1qj68eaNzJnunonYzBj/0m1i9r\n/WIbtNFf6yuDxbQ+LfSsvv8sF4TX0nAr9zkyfaut3LBLb1VJT/3c6FRI1oQvzZaO3kD32Pq/\nSTulyZDuxoi7+YzdNxn55iUf3smeqeoC0u5kDH7sN7F+WesX26CN/lpfGSymxq3ja8Ep/5e8\nNlwSfXHMqt35chsmSfU11XGwQ20v4PWtrF0utQN75vKm5kM6Zbn509mQzG++lvBTtzI17mTP\nVBuSNiVj8CO/if3LEtIODBbT7t/SMpvXUnJv1kZ341KC5p/uR70T3l/cMHXUzp5wkd/qI1un\nYiKk6TXS+PhHQkrGFtzmHqf2woZk9BHWkKwp6YMf+U3sX5aQdsB8clNr677eP7p3+0n6w6pl\npTr6qwX21+82/ennkewJV+5nY+k3l7fpfaTx8dc3TtY+knVBXmG0cSsK604j+0jF4E7m4Ed+\nE/uXJaQdMJ9c+3hTfdi3vptxELn8B/21qP7VB3q1gxDdlQ3lJPorG0YPB5Y/1P7FNpe39gCa\n+iok46jdrb7e72ZcVdc+4K/5hYw72TNt72zcyRj82G8yctRudOybsdFf6yva1n/zL3qjXm9k\nza2zueJpT1yady7l2pkYfQfInPC97rA54zsa0rvzSPbX9qOu+hd/Iw9o1xvGnYbnkeo763cy\nBj/6m1h/RULaASukbhlo0iiDKbeA7sq6RvTSPigt95P0E7XdUnR8NjMohhNud9GtKxsK7b/3\n5h7fhfSnX9nQTML4J6B7QHlILbfvZM20u7NxJ2PwY7+J9csS0g7YIb029xP9zEjSHwA2H/go\nr2lLqjueb8ZPqikc9QvNRiZc7VWcrvpd7OXtca4umvsupOqiONUekK8uvDNfudRPIGvOLxl3\nMmfa39m4kz74sd/E+mUJCevLrctitzvTeBFSwJptysfJvpxvczONHyEFTHuB4PAY9qZmGj9C\nClj3+gQ18rYQm5pp/AgpZPmlPEmVnL2uGlaZafQICRBASIAAQgIEEBIggJAAAYQECCAkQAAh\nAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAh\nAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAADyEpIDI/LOXy4awwC0BIFREh\nAYso7f9/eKBThIRYKOM/PzzSJUJCJJT13x8e6hAhIQ5qcOOHx7pDSIiCGrn1w4OdISTEQI3e\n/OHRrhASIqAmbv/wcEcICeFTk1/88Hg3CAmhs65lICTgB/YiSkjA9wZLKCEBXxsuoIQEfGtk\n+SQk4EtjiychAd8ZXToJCfjK+MJJSMAXpl4KS0jA5yaXTEICPja9YBISMOFwOJjfmFkuCQkY\ndzhYJc0tloQEjLNDml0qCQkYZ4WkRrb1jJ9+i5CwDyMdTZZESMB71ekjQgIWqRdIQgK+M36Y\ngX0k4BvmuueTxZGQgIEmpLqmj5ZGQgIG6pDq//9sYSQkoKZvzfUhfbgsEhJQMfaLDp0PH01I\nQGUY0txROhshAZWRkNT8ZUE6QgJqRjLtYYZPN+8ICRhXLYWEBCxSB0RIwMdGYmkLYh8J+NDI\nakd9c+y7ICSgO2ukfUd9vE3XPeBrhIT4zFYxOPmq3j1igJCwB2/WL2ZIU28COYeQsAdvN9S0\nO/y09BES9mA8JOsy1eE3P0dI2IWJjqYPe3+JkLBbY818e9i7e5yXhwQ4C2Ds9BFrJOBbg2S+\nfO2EhpCAlnZZ0Lc1ERJQU+3+0XevjW0e/MP8vn9IgLMATP17EhMS8CvtVXyEBAx91ITxKj72\nkQDbXBX2ZUG/HbHTJuD4IQHOAnthhGSW0v3o8zc5meI/pOtRqfTudBZATw/JWju1X7757KNP\neAypvjj9pCqZk1kAQ+YKaSSkb94taIrvkDKV5UXxzNTVxSyAWaMhffVuQVN8h5SovLydq6OL\nWQBz7FqqrxcfZ6j4Dql98eH8ixAJCQ4cjBNEv19XN8Z3SOc2pMTFLIAZB70kbbNOgteQ0sv1\nrm6vm3k2f7SBkODAMKTv3+RkiteQatXNJHcxC2COFpJ1Y3FLPs8jPR7Xa5pWhxyy2Y4ICU6M\ndDR8K67fcGUDdkQ/0mCukbYTktK5mQXQ6POJNiQuEYJz78OwQ1o6Ry4Rwva8T8NaIy2fJZcI\nYXvexmEfvls+Sy4RwvbMxVGdhjUvcIgxJC4RggdmG/pXdTfxh8QlQvDNOJzQXd1gv9ZvYU5c\nIoStM8oZhqTdacFMuEQIW6eXczhMlBRRSFwihHVsLaSgZoFdGelo5G2/FyAkbJxWiDI7Ejlc\n107ay0MCnAV2ojve3X1mi7GZJzUbQsK29asgVYwceJCaDSFh2yZPHBES8IVDeyXD2MfFGv9d\ngpCwB3NrH5E1EyFhDwgJWEypuc03QgI+8W5xYh8JeKM87C154nUKIWHLJq4HkkdI2DJCAgQQ\nEiBAmR25K4qQsGFq9pMvhefk4SEBzgLbp9rP4mu/QUjA14ZLESEB3xpbiNhHAr7jeRkiJGyS\n70WIkLBFcp9p+cUM3T8kwFlgy7pXlXudo/uHBDgLbFh72JuQ3M8Cm9V8QAMhERIWaBceQiIk\n/K5bdgiJkPCzftEhJELCL8pulPW1N4SEbfDzsqNJhIS4NQF5ev3eJEJCtLSAfL0QdhIhIVaH\nMWsNhpAQq9GQ1iqJkBArY+eou7nSYAgJ0TJfRO77zJGJkLABzQLDPpL/WWBD1jxc1yAkRG/d\nA981QkLc1j7u3SAkRG3tw94tQkLUAumIkBC1MCoqCAlRU6F0REiImArigF2FkBAtZW3YcfX3\nCrNA9OztOq7+XmMWiNormcF2HSGtMQvEbPRqb0JaYxaIWdOMFQ77SCvMAjEL5mBdh5AQHxXC\n9d4mQkIMjHBCXDwICREwNuWCXDoICRHQQwpz4SAkREALKdBlg5AQg9A7IiREJdglg5AQDxXu\ngkFIiEbIiwUhIRZBLxVeQ/q7pKqUZn+uZoHNWv9NIOd4DCk/qt7JySywXX1HQZbkMaRMJbdH\ndet5T1TmYhbYrHaRIKQiUY/u9kMlLmaBreqWCEIyj13OH8gkJBi0BSLMjlgjIXwBnz7q+N1H\nuj+rW+wj4QvTS0NAayefh79P2lG7Y+5kFtie2Y6CKcnveaSsOo+UpBfOI+FDM8vCbkMKaRaI\nw9yiQEijk9W5mQWiM78khNPRCiFdj0qld6ezwFbEsyB4P4/UHHGYPWgX0d8PUsbWLmFfXmfw\nHVKmsrwonpm6upgFojWyv6PU5I/C4zukRFXHvXN1dDELRGtYS/BXBRl8h9T+K8MlQjAMalHT\nPwqR75DObUhcIgTDZEfsI9mPU+nlele318084xIhzIluCfAaUneOSKmES4QwLb4FwOd5pMfj\nek3T6pBDNttRhH9HSIrw+Q/nygbPs0C4Ynz6CQmhifLZJySEJdILLQkJQYn1qSckhCTaZ56Q\nEJB4n3hCQjgift4JCcGI+WknJIQi6medkODf6FWocT/phATv7NdFlF9FevqoQ0jwzgopitcb\nvUNI8E4v59BYdUACCAl+6LEMOyIkJwhpcyZq2UpHhAQ/5kPyPx5phAQvpnrZyrYdIcGPNx3F\nXhIhYU2E5BQhbZp90G4LZ5IICb712XSroug7IiR4pm3IqQ2siVqEBK+0naLDdg5+ExI8O4xZ\ne1ACCAle2QkRkkuEtF2EtOAhAc4CqzE36rbRESHBt7Gdo/hrIiR4NnKMYQPbd4QE95pM2k05\nQvr1IQHOAv70e0Pd0QVCIiR8axjSyCppxfFJICQ4pxektnPBt4GQ4J7e0ejRhvgRErzQ0tli\nR4QEP/R09KTWHJMkQoJ7E+ugLa2WCAnOmYcX+ngIyTVCipwRyGFI+8E6AxRHSJA3crp1LCT2\nkVwjpLjpW3EjGW0nnx4hQV4Ty2hD2+yIkOBCt0LaRUQlQoIzej7bTahGSHDGXA1tOyVCgivm\n9tyWt+sKQoIzdULtm0BuegepICQ40+8cEZLUQwKcBZzSNup2sWVHSHBhcLh76x0RElwYnDgi\nJJGHBDgLuDQ4BUtIIg8JcBZwaqyktcfkFCHBCbX9q4IMhAQXqmdwN9cHFYQEJ+onUOto8yUR\nEiTVxajui12cQyoREgTVybTbddp7fhOSyEMCnAVcqC+v625p315vTH4QEhbSj2/v6AyshZCw\njH5RqrbuISQXDwlwFhBihqT07684Ku/8h3Q9KpXenc4CHhkh7feJ8xiSqh54UpXMySywAm0f\nacfPm++QMpXlRfHM1NXFLLCqPT9tvkNKVF7eztXRxSywpl0/a75DUkr7QnwW8Mo8nrDvJ813\nSOc2pMTFLOCTeYR758+Z15DSy/Wubq+beTZ/tGHnT0ok9NcaKbW3490WryHVqptJ7mIW8Kh7\nrdF+LvGe4fM80uNxvaZpdcghm+2IkGJgvYyckDw8JMBZYClC0oUTktK5mQVEtemobgtv7RGt\naJWQ3oZCSBFo95F4skqEhN+0ryPnuaqscNTug603npzgtW9swlNV8xjSX0JIcTNf83rY+eEF\n0+KQ7mnZRPr84IF5qk7V/di0i5IezmGmpF3mtTSkU71yUcknJRU3VV3YQEhRasPpDnZ3x72H\nd1tpiOtZGNJVnfIyi6s6f/TY50mlOSHF6XAw89G/tO+21hhXszCk8lURxlXdb11UciekOI2F\nNOyGkH54SLVZ91VIxeP4/oQrIQVGv4bB6mfYzR47WhrSsVkjPeZfqGc6E1Jk+p2jblOu72eX\n3QzI7CPdk/mXji+ZBdY3DGnqVRO7rWrpUbu0OSt0khrQcBZY3SAkVYw3s8/9o5LIeSSV3oSG\nMzoLrK7fH6pfxlf0lwjZ9yMkhw8JcBb40cFaPQ1/uNLA1kVI+M7g+Lf903WGtbalIeVZ+S4m\n717xumgWCMrBWiWtPZ5ALAzpmbTvwfDRJUK/zALhqHaPRo6AY2lIJ3Uu10V5plKpEdmzwGq6\nTg4j1h1acASubDBviCCkEJgXLxDSLIFr7Uo5IW2M3gsdvbcwpEyd/l7/+Tu9+XiJBbPAGvpg\nRldHdGQTeT0SVzZszpuOCMm2+DzSrbyy4SR6pR0hradtZDokDtaN4oQsNH0krJG+Q0jQdJH0\nwYyuktYeZ3gWhFS/qM/Jm6MS0krskIqxkDCCkKCzt+wKo6SClKawaQfT1IpI//7aYwzQwpBS\n0dNHo7OAV5O7Rvr6ae1BBkjqEiFZhLSaiX0j8+frDS9YAm9+4gAhreZdSJQ0bmFIeVpdIiSN\nkNYzvo9k/pySBhZv2nHUbpPedERINkKCbvxAg3kHQhrB4W9ozH5U863hPdYZXcgICVoYwxXR\n2CrJ+wgjsCSkZ5Yo6bc9sWYB9/RShtt0bMp9ZkFIz/oT+GTf9sScBTxoQxnfPSKkzywI6axO\neZGfPvxkpJ9mAQ/MVQ8h/WZBSPX7NTxVIjkecxbw4E1InIH9yMKrv/v/iCIkjwhJAiHt3vQ+\n0uxHLsNASDj01wSZnxM78U75GEFIqBkhmd/Ee4tCMqw8KvymC4WQFiGkfdNKqW4oKx227D7E\nJUL7Zq9/3twDUwhp3/RMxg8sENJHCGnn7B2iQTWE9BFCQsk67G39aJUhxYWQ9m7kPCy+R0g7\nNzwPS0i/IKSdM0PicPevCGnnmnD6gCjpJ4S0G3of1skjpZ2HJaSfyL2L0Enw3YsJSd7gjFH/\nI1UMLhXCl+RCUoKv8CMkecNTr8Xgq3Yrb5UBxm3ppt05ub/+/56ovyKV+0BmQpI3GVJ/sI6V\n0e8WhpSpR/XfhzoVuTrKjImQnLD2kdqv1ExIdPUxqU+jqD91TGJE9izgQt/RzBqJNdTnFoaU\ndGukhJAi0iWiisE+0si98NbiTbt2Hykrbq/Nu/VGhW/0q6AP7oUPLD3YcGoPfpcrpOuKo8IX\n9MsY5mqio48tPiF7T18ZpeVqSV1khjSYBcQctLc06Q8zkMtyXNmwJ9qKqPysCUKSQ0h7YoRU\nf82OkAxC2oX+qNzoOwVZN/C9pSFdjryLUPjadPoV0chBBlZNSywM6cLbccXgoKm+ocYvYyCk\nny0+ISt2yHtqFljODqk7DVtw2bcQqUuEPvJ3SauVV5r9iY8Kc4yQ6idtsJZiH2mJhSGl6vNP\nvsyP2ksu5i+CICRxWkfWd+hHwsKQnsnpzcqll6nkVl+Z96wuKZIdFeYNOiq0lNYb1mbIvbDv\n7ePaC1xLj/lXARKSuDYZNfweHQnwGJL6fOeKkOQNO2KvSJDHE7KskVZgpEI37ngMqXzJxbO6\nxT6SH9bhBDbjHFr4iX1ffT7SSbv3cfZoHyEtYJ4Y0koiJId8hlT8ZdV5pCS9cB7JmT6XiZCo\nyQUuWt2Yw0RIXVCsl5wIJyRnn6O5L9alCoOOCMmNpSFdj0XxPKrjx6dlqym8mysh/axb7fRf\nEZIHC0O6l00k5Trkm5IIyaFu861/UTmHHNxbGNJJ3YqHOn70DkJffAo6IS1j7A/pFRWE5IjA\n1d+P8pzQB7s1fwkhefI2JFqSJhBSqu6fvZ4iT9Xp2T5KelTQjYZESQ4t3rR73MurfT58c8ib\nem0KEpIP3Q7R2HsFEZK45QcbVPl+dqpcK33geVJpTkiO9REV7WeIEZJjiw9/11fNHW+fPvqi\nkjshOaVt1hUTb15HR9L8n5B9HN+fcCWkJbS9Iq2jtUe1cWtc2XAmJGdGDjEQkg+LQ7qV13Sn\nH2/Z/TILfOLQX7ZAR/4JfhqFIEL62nRCROTDwpCu3ecjib6/HSF9bTKktQe2EwtDOnaf2Cf2\n+bH2LPCRJho6WonUG0TylsUra6IZZkRMXoitkWbfzGTJLPCJ8RVSwblXX9hH2oaxXSTtB3CN\no3bxGStj8igDIfmx/DxSynkkv8bTGKyQ9B+M3RmiwnnPBs+ziNdoSF8dqmMtJY+QokNIIVr0\nvnbO3viHkOa830d6PwFCEkZI8Zs6zjD7CB8D2xM27SJlHeD+LiSII6Q49cWMdERI/i0K6Xmu\nTsPmR+lPZCakd2ZDWntwe7QkpGei0vK/d6WSp9yQCkJ6by6ktce2S0tCOqpz/eEsfyfZi78J\n6b3pfaSVB7ZTC0K6l28f1EiV6LUNhPQVMlrfgpDOqv+wsKfsxXaE9B06Wt3CDxrrv8l5pNW8\n+lFUtLIFISWEFARWRyFYtGnXv7vqvT5+J4WQ3jsYL5Sgo5UtCOnRH/R+Jhxs8KxLhxVSEJYc\n/s5Ucilfaf64JLywzzdCCsuiKxsu3RWrZ8EhFYT0AfN8LJehrm3ZtXbPrHqf1YvsdQ2E9Iku\nIFZGIeCi1WhxFjYkhBQdLq4LESHFho6CREixoaIgEVJsCClIhBSdvqFD8y7Fa48IhBQnbWXE\naikMhBShAyEFZ+HLKIw35DqLnZclpFnfh0RsrkmGJHfBHSHNMQ80fNgRJbkluGn3lHtNEiHN\n+OGIHSE5J7mPlIm9JomQNNaqRxFSiHjL4tBZZ41Uf7XqdxNxNkCUCCl0ZkiqoIogLd20S6uP\nvvxLZF+QREg9IyRFRYFaGFLWfRhzJjOe4Sx2btgRJYVoYUhK2TdEEFLH2q4jpEAtDCnp1kiJ\nzHiGs9g54zgDIQVr8aZd8vf6zz3R3r5YACH1jOMMHGkI1dKDDafmmJ3o29oRkk5bJ1FRsBaf\nkL2lZUb3ifv+iJA0B9Paw8Eorv4O2GHsA5DWHhRGCYX0yDjYIK7thpAiIBHS83JUHLWT1nfD\nPlIEFoeU347lCyhkd5IIyTrEwLoodAtDutVH7XinVXF9SHQUgyUh3c+vhpLsIXtVgzGLvToM\nrT0kzFr0QWOvisrTsYQkbaQjSgrbopdRZO0NseFYs9ir0ZAoKWSskUI0SIiQQiewj/RHSIKG\np46KgmtVw8dRu7DY547afugocELnkYQvttt5SN22HCuiaHBlQyjMNVH3vt6EFAf/19pdj+9X\nYDsMyT68wKuP4uLx6u/6mETzAqb5t3ggpB3+BeLmO6RMZXn5Ic7q6mIWEaOjuPkOKVF5eTtX\nRxeziBknX6PmO6T2nNP8uacdhjQ4B8v+UVR8h3RuQ5o9OEFIHPqOjNeQ0sv1rm6vm3k2f7SB\nkAgpMl5D6t4kXKkkdzGLmI1s2RFSPHy++cnjcb2maXXIIZvtaHchmSdjD/o3EQfeRSgAbT0c\nsotXOCE5+4yY8Okbc3QUpxVCuibqOHs6dp9rJM4hRc3rPlKqkmtx+eBzm3cd0tqjwS88hvSo\nL7JT57x4plwipKGj+HkM6VyeO8rqM7FcIqShow3wfvV387kVXCLUI6IN8B7Srd6m4xKhHiFt\ngNdNu3N7GjY/c4mQhpDi5zGkPNE+cXb+BbW7CKkrp/yzEFLkvJ5H6l6Pnrz5DPQ9hNSthFT3\n1cojwgLhXNngeRarazfnVPfluuPBIoS0liakPfyqe0BIq+EAw5YQ0nooaUMIaTUc894SQlqD\n/qqJtccCEYS0gjogOtoSQloB16huDyGtgI62h5DWQEWbQ0irIKStIaQ1KK4I2hpCWsHGf71d\nIiT/tv3b7RQhecd23RYRkmdKcaRhiwjJo+ZVE4S0QYTkTxsQIW0QIfnRXVzXfoFtISQvuCho\n6wjJi/pV5XS0XYTkBSujrSMkPyhp4wjJB3aQNo+QPOBIw/YRkgeEtH2E5ByHvveAkFwjol0g\nJNfIaBcIyTVC2gVCcoxXH+0DIbm1nd8EswjJJVZGu0FIDrF7tB+E5I4ipP0gJEeMV/Jh8wjJ\nhfYkLB3tBiE5wCVB+0NIDhDS/hCSA3S0P4TkgCKj3SEkebyX6g4RkjjFmdgdIiRpvCnxLhGS\nsHrodLQ3hCQr3pFjEUISFe3AsRAhSYp13FiMkOSoOIcNCYQkJspBQwghSYlxzBBDSEIiHDIE\nEZKM+EYMUYQkIroBQxghSYhtvBBHSAIiGy4cIKTl4hotnCCkpTgNi4KQFotoqHCIkJaJZ6Rw\nipB+Vb3kKIaBwgevIf1dUlVKsz9Xs/CmehFsBOOEHx5Dyo+qd3IyC494NTl0HkPKVHJ7VLee\n90RlLmbhESFB5zGkRD262w+VuJiFP7wFJAweQzJOuMyffQk+JN5LFSbWSD+hI5j87iPdn9Wt\nqPeRqoAICSafh79P2lG7Y+5kFu41h70pCQa/55Gy6jxSkl7iPY/UBkRI0HFlw5cICWPCCUnp\n3MxCQhcQHUGzQkjXRB2vbmfhDvtGGOUzpEeqkmtxifgSIT7TEhM8hvSoCsrUOS+eqZpdJwUa\n0oGQMMFjSOfy3FFWn4nN1dHFLByjI0zxfomQSrUvpGfhRpcOIWGK95Bu9TZdRJcIdfkoNu0w\nxeum3bm9nCE/R3KJUBlNnw8hYYrPF/Yl3facml8hBRPSwUBImOT1PFLW5pPMro+CD2ntYSE8\n4VzZ4HkWHzE74moGTCKkWXZJwDhCmkdJ+AghzToQEj5CSLMICZ8hpFl0hM8Q0jwqwkcIaVyb\nDhnhI4Q0dBhadTyIACENjHRESHiDkAboCN8jJBMV4SeEZCAj/IaQDHSE3xCSgTUSfkNInUP/\n7viEhC8RUmsQECHhc4TUGq6J6AgfI6QWm3RYgJA67B3hdzsPyeqFkvCjfYdk99IlREj4DiH1\nL5cgJPyMkA7dDeuInacxYBN2HZK5QmLnCL/bc0hTHRESvkZIBR1hOULiBBIE7DmkYrwjQsL3\ndhtSFwwdQcBeQ+qboSMI2HdIY5zPGltESHQEAbsM6aB/jCUZQcAeQxrmQ0hYiJBYJUEAIbWr\nJK74xgL7C6ld8wzWQ4SE3+0upMNkSLx0Ar/ba0jsGUEUIdERBOwrJI7UwZFdhTTeEa8sx3KE\n1JdQBMgAAAnASURBVB14cDNP7MPOQ9JfkuRmntiHHYU0s3tESFho2yHpecweZ6AjLLPpkOxW\nJjsCFtp1SDJzAfYbkswMgMamQyomOyIkyNp2SB026+DWHkKaPV4HSNhBSKMdERJE7TUkyRkA\n+wqpICI4soOQCtZCcG4PIXEtHZzbR0hs08GxXYSklPFhSLITB4qNh9QeYCisN+GSmTrQ8x/S\n9ahUenc6i4Z+sJuQ4JTHkFT1wJOqZE5mYRqEVBASHPEdUqayvCiembq6mMXUS/m0H9ARXPAd\nUqLy8nauji5moa9wuIwBHvkOSSntC/FZaCFplzL8OjHgY75DOrchJS5mYYek2JKDH15DSi/X\nu7q9bubZ/NGGRftIbNPBP68h1aqbSS47i74aSsIKfJ5Hejyu1zStDjlksx19P4vxgwyEBF82\ncmUDIWFd4YSkdN8+eCIk6TECE3yGlJ+VOjUXB0kf/h7dR/p6KsCPPIaUJ9XKJq0n4uY8Eisk\nrMNjSNVlQfk1OVUTkQzJOHdESFiBx5CS+oHP5PiUDamNxu6IkOCN96u/Xyul00k0pLaaQUeE\nBG88hnRU7cmj44mQsC0eQ7qqc3PrqU6EhE3xefg76+q5vzlV9PEsDv2L9ugIK/J6QvaRtree\nZ5GQml76nOgI6wjnyoZfZtE3M1wfERI82kZIo5yOEDDEHdJUSbyeD57FHBLrIgQj4pCoCOGI\nNKTD6PFuOsJa4gyJQwwITMQhjebkfmzAiIhDYoWEcMQZUr+PZH2wpfuhAWMiDanEygjh2FBI\nlIT1EBIgIN6QONCAgEQbEh0hJLGGREUISqQhsTpCWOIMiY4QmLhCapKhIoQmqpDabAgJodlA\nSO6HA7wTZUiKjhCYqEJq9pGU9iUhIQhxhWT/kI4QiPhC8jB74FvRhURHCFFsIdERghRZSHSE\nMEUV0vefdg74EVNIisN0CFVEISlOHCFY8YRUfo+QEKhoQqq+RUgIVCwhNd+hI4QpkpA4XIew\nxRESHSFwMYTE6SMEL4KQyAjhCz8kOkIEgg+JjhCD0EOiI0Qh8JDoCHEIOyQ6QiSCDomOEIuQ\nQ6IjRCPckDgNi4gEGxIZISahhkRHiEqgIdER4hJoSO5nAUgiJEAAIQECCAkQQEiAgEBDAiLz\nw1IuH47vGTn8FZg0kw5oHo5nFOczwKQ3MWmv83A8ozifASa9iUl7nYfjGcX5DDDpTUza6zwc\nzyjOZ4BJb2LSXufheEZxPgNMehOT9joPxzOK8xlg0puYtNd5OJ5RnM8Ak97EpL3Ow/GM4nwG\nmPQmJu11Ho5nFOczwKQ3MWmv83A8ozifASa9iUl7nYfjGcX5DDDpTUza6zyAzSMkQAAhAQII\nCRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAc5Duh5VkuXVzSzp\nbspO+ed3Pp+Qn5U6P+rbwoPWJi096tJfMznhUWuTlh61Pj3hUWuTdvG3tmbmcNqlrPoFkvKv\nc6puHuWn/JD+MyXV5KrFXXjQ2qTFR/2SJ/XkpEfdT1p61Pr0hEetTdrF39riOKSHOr+W9Ks6\nl/+iJY/ikag/8Sk/VCozzUZWTjSrJio8aH3S0qMupfWyIj1qbdLSo9amJz1qbdIu/tYWxyGl\n7fZAufjcX7du6iI+5avUNBuJypspSw9an7T0qItynPXSLj1qbdLSo9amJz1qbdIO/tY2Pwcb\nyichVc9C/t+GOqSr6DSbKSeFs0GXk5Yf9VOd6qVdftTdpKVHrU1PetTapN0sIQYvIeXq1Pwz\n3P1Hcsqpup9fe6mS0y3/fSz/9k4GXU9aftQn9Wx3rQvtP6KTlh61Nj3pUWuTdrKEmLyEdC1X\n2k6WyWrKab0neRKc7mtLxsmTq01afNQXdSschdRPWnrU2vTkQ+om7WAJsfkI6ZmUa2sXIbVT\nvr3WTZnk6vuaJtVmtYNB95OWHXW1TeQmJGPSsqPWpic9amPS4kvIYG7uJt3Kk+pfAgfLZDPl\n5gvZw73F2cWTq026JjfqY3kmwE1I2qRr0n/ranpu/tbaUKVHbfAQ0qkefiL/dzoZfxjxZyBx\nMuh20g2pSZ+rY1711IRHrU+6IX1Kppyem7+1Pj2XJ5Kch/Q8np7VjfqgzFPuUFI35YaLZ0B8\n0P2khzcXTrL/aHvhUeuT7r4lM2l9evH8rUdn427SlXu3h3ep/l27K6ljJ/2U65Mzcs9AO72j\n/KC1SQuPWl/ahUetT9rV3zp197dO5Uc9wnFIz/5IifCJa23KWfm3z+vzeRKqyw/ytNyRcXJl\nQz1p6VFX3F3ZULTneqX/1t305P/W3aSd/K1NjkM6axsFR9FDkNqU8/oCNrnTBEk/UtlB65MW\nH3Wp2XiRHnU/aelR69MTHrU2aSd/a5PjkPSt67y6uNfVlI+Shzb76ckOejBp0VEXXUjiozYn\nLTnq3N3fOnf7tzbxeiRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABC\nAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABC\nAgQQEiCAkAABhAQIICRAACGF5ouP3q4/sjA5P/VvORgS3uPvHpqvQ3ql9NS+5WJMeIu/e2i+\nCqn8//zk9FOG8RFCCs3XIRW5SlwNBp8ipNDoIV2P1Wdx/6lz+dVd3cv/nNWfddfyv0rlR5U2\n38oSdXq2U0icfpo3aoQUGi2kU7UHdCqKpPreud6E69c/+hpJqVS9fl59q3pckr9upe0U4Bgh\nhaYP6aaSR/FI1K24vP5X/iSpvnsx7/qs9pFeueTNt27lzaq6e3nrtQ919/5b7A4hhaYPKa0C\neMVQPMuVyt9rnfMoVzfP/q7NUbu8vP3XPjotb1arqVSVceWvTT44Rkih6UNqbpX/Ob2KyNTj\ntTJ6ahtq+nkk7c76FBq+Br9f/IlDMxrS/ZVQciyOx2Yrz74rIa2NP3FoRkMq1PHvtdOTlYfm\n8pG7TobkerRo8JcOzXAfqdzDydT59dX99f/nsbta24H9PhKHGTwhpNCMHbUrjzQoVR1S0NOY\nCOlaHqvLyqN21RReX3OwwTlCCo22W9OdR3o5qmP1ncS46+C2dR7pZF+LB0cIKTT68YFrUl3Z\nULpUZ2MvxmV1UyG9Vkcq7a5sUGc6co+QAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGE\nBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGE\nBAggJEAAIQECCAkQQEiAAEICBPwDNY5i9XXCGoIAAAAASUVORK5CYII="
          },
          "metadata": {}
        }
      ],
      "execution_count": 9,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false,
        "execution": {
          "iopub.status.busy": "2021-07-08T01:47:02.308Z",
          "iopub.execute_input": "2021-07-08T01:47:02.311Z",
          "iopub.status.idle": "2021-07-08T01:47:02.353Z",
          "shell.execute_reply": "2021-07-08T01:47:02.368Z"
        }
      }
    }
  ],
  "metadata": {
    "kernel_info": {
      "name": "ir"
    },
    "language_info": {
      "name": "R",
      "codemirror_mode": "r",
      "pygments_lexer": "r",
      "mimetype": "text/x-r-source",
      "file_extension": ".r",
      "version": "3.6.1"
    },
    "kernelspec": {
      "name": "ir",
      "language": "R",
      "display_name": "R"
    },
    "nteract": {
      "version": "0.28.0"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}