{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 6. Linear Model Selection and Regularisation – Conceptual\n", "\n", "Excercises from **Chapter 6* of [An Introduction to Statistical Learning](http://www-bcf.usc.edu/~gareth/ISL/) by Gareth James, Daniela Witten, Trevor Hastie and Robert Tibshirani." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1. We perform best subset, forward stepwise, and backward stepwise selection on a single data set. For each approach, we obtain p + 1 models, containing 0, 1, 2, . . . , p predictors. Explain your answers:\n", "\n", "\n", "### (a) Which of the three models with k predictors has the smallest training RSS?\n", "\n", "Best subset selection will exhibit the smallest training RSS because all possible combinations of predictors are considered for a given k. It is possible that the most effective model is also found by forward stepwise or backward stepwise, but not possible that they will find a model more effective than best subset, as the former consider a subset of the models considered by best subset selection.\n", "\n", "### (b) Which of the three models with k predictors has the smallest test RSS?\n", "\n", "Best subset selection will certainly yield the smallest test rss, this model might also be discovered by forward stepwise or backward stepwise but this is not certain.\n", "\n", "### (c) True or False:\n", "\n", "- i. The predictors in the k-variable model identified by forward stepwise are a subset of the predictors in the (k+1)-variable model identified by forward stepwise selection.\n", "\n", " - True\n", "\n", "- ii. The predictors in the k-variable model identified by backward stepwise are a subset of the predictors in the (k + 1)- variable model identified by backward stepwise selection.\n", "\n", " - True\n", "\n", "- iii. The predictors in the k-variable model identified by backward stepwise are a subset of the predictors in the (k + 1)- variable model identified by forward stepwise selection.\n", "\n", " - False\n", "\n", "- iv. The predictors in the k-variable model identified by forward stepwise are a subset of the predictors in the (k+1)-variable model identified by backward stepwise selection.\n", "\n", " - False\n", "\n", "- v. The predictors in the k-variable model identified by best subset are a subset of the predictors in the (k + 1)-variable model identified by best subset selection.\n", "\n", " - False" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2. For parts (a) through (c), indicate which of i. through iv. is correct. Justify your answer.\n", "\n", "\n", "### (a) The lasso, relative to least squares, is:\n", "\n", "- i. More flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Incorrect. Least squares uses all features, whereas lasso wil either use all features or set coefficients to zero for some features – so lasso is either equivalent to least squares or less flexible.\n", "\n", "- ii. More flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", "\n", " - Incorrect\n", "\n", "- iii. Less flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Correct. As λ is increased the variance of the model will decrease quickly for a small increase in bias resulting in improved test MSE. At some point the bias will start to increase dramatically outweighing any benefits from further reduction in variance, at the expense of test MSE.\n", "\n", "- iv. Less flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", "\n", " - Incorrect\n", "\n", "### (b) Repeat (a) for ridge regression relative to least squares.\n", "\n", "- i. More flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Incorrect. Same as for lasso\n", "\n", "- ii. More flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", "\n", " - Incorrect. Same as for lasso\n", "\n", "- iii. Less flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Correct. Same as for lasso\n", "\n", "- iv. Less flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", " - Incorrect. Same as for lasso\n", "\n", "### (c) Repeat (a) for non-linear methods relative to least squares.\n", "\n", "- i. More flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Incorrect\n", "\n", "- ii. More flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", "\n", " - Correct. If the data contains non-linear relationships then the reduction in bias will outweigh the increase in variance as the model better fitst the data. At some point the model will become too flexible and will start to overfit the data at which point increased variance will begin to outwiegh any further reduction in bias.\n", "\n", "- iii. Less flexible and hence will give improved prediction accuracy when its increase in bias is less than its decrease in variance.\n", "\n", " - Incorrect\n", "\n", "- iv. Less flexible and hence will give improved prediction accuracy when its increase in variance is less than its decrease in bias.\n", "\n", "\n", " - Incorrect\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3. Suppose we estimate the regression coefficients in a linear regression model by minimizing for a particular value of s. For parts (a) through (e), indicate which of i. through v. is correct. Justify your answer.\n", "\n", "\n", "### (a) As we increase s from 0, the training RSS will:\n", "\n", "- iv. Steadily decrease.\n", "\n", "As we increase s we increase the models flexibility, in the training setting this will lead to overfitting and a monotonic decrease in the RSS error.\n", "\n", "### (b) Repeat (a) for test RSS.\n", "\n", "- ii. Decrease initially, and then eventually start increasing in a U shape.\n", "\n", "The increased flexibility of our model improves our models performance in the test setting up until the point where our model starts to overfit the data.\n", "\n", "\n", "### (c) Repeat (a) for variance.\n", "\n", "- iii. Steadily increase.\n", "\n", "Reduced constraint on sum of coefficients means monotonic increase in flexibility – increase variance\n", "\n", "### (d) Repeat (a) for (squared) bias.\n", "\n", "- iv. Steadily decrease.\n", "\n", "Reduced constraint on sum of coefficients means monotonic decrease in bias.\n", "\n", "### (e) Repeat (a) for the irreducible error.\n", "\n", "- v. Remain constant.\n", "\n", "The irreducible error represents the inherent noise in our data, because this noise is random it contains no useful information and remains constant regardless of the models flexibility." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4. Suppose we estimate the regression coefficients in a linear regression model by minimizing for a particular value of λ. For parts (a) through (e), indicate which of i. through v. is correct. Justify your answer.\n", "\n", "### (a) As we increase λ from 0, the training RSS will:\n", "\n", "- iii. Steadily increase.\n", "\n", "As lambda increases a heavier penalty is placed on a higher total value of coefficients so the model becomes less flexible and so is less able to fite variability in the training data.\n", "\n", "### (b) Repeat (a) for test RSS.\n", "\n", "- ii. Decrease initially, and then eventually start increasing in a U shape.\n", "\n", "As lambda increases the models bias is increased reducing any overfitting of the training data and so improving the test RSS accuracy. At some point the bias will increase to the point that our model is less able to represent true relationships in the data.\n", "\n", "\n", "\n", "### (c) Repeat (a) for variance.\n", "\n", "- iv. Steadily decrease.\n", "\n", "As lambas increases a higher penalty is placed on model flexibility and so variance decreases.\n", "\n", "\n", "### (d) Repeat (a) for (squared) bias.\n", "\n", "- iii. Steadily increase.\n", "\n", "As lambda increases coefficient estimates are reduced and so bias increases.\n", "\n", "### (e) Repeat (a) for the irreducible error.\n", "\n", "- v. Remain constant.\n", "\n", "Irreducible error is not effected by the model." ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## 5. It is well-known that ridge regression tends to give similar coefficient values to correlated variables, whereas the lasso may give quite different coefficient values to correlated variables. We will now explore this property in a very simple setting.\n", "\n", "## Suppose that n = 2, p = 2, x11 = x12, x21 = x22. Furthermore, suppose that y1+y2 =0 and x11+x21 =0 and x12+x22 =0, so that the estimate for the intercept in a least squares, ridge regression, or lasso model is zero: βˆ0 = 0.\n", "\n", "### (a) Write out the ridge regression optimization problem in this setting.\n", "\n", "![6_5.jpg](./images/6_5.jpg)\n", "\n", "### (b) Argue that in this setting, the ridge coefficient estimates satisfy βˆ 1 = βˆ 2 .\n", "\n", "![6_5bi.jpg](./images/6_5bi.jpg)\n", "![6_5bii.jpg](./images/6_5bii.jpg)\n", "\n", "### (c) Write out the lasso optimization problem in this setting.\n", "\n", "### (d) Argue that in this setting, the lasso coefficients βˆ1 and βˆ2 are not unique—in other words, there are many possible solutions to the optimization problem in (c). Describe these solutions.\n", "\n", "![6_5cd.jpg](./images/6_5cd.jpg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 6. We will now explore (6.12) and (6.13) further." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### (a) Consider (6.12) with p = 1. For some choice of y1 and λ > 0, plot (6.12) as a function of β1. Your plot should confirm that (6.12) is solved by (6.14)." ] }, { "cell_type": "code", "execution_count": 52, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEKCAYAAADn+anLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAFsBJREFUeJzt3X+MHPd53/H3c6SOYZlrfKbOV4YnhYLDEJHSmlavilu2gWzFsWy0oBREjlw0FBqVNCAJsAG3hZT+ERetYbeprTpBI1SMBJOBI0VITIhwBDuyIlcIEP84uaosWWDEynJ0EkMy1EVmGUaH4z39Y+eoE3Vzt7u82dnbeb+Axe5+d/b4zO7MfjjfmflOZCaSJC1lqO4CJEn9y5CQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklRqfd0FXIxLL700t23bVncZUvuOHGnd79hRbx1qtCeffPKvMnOsnWnXdEhs27aNqampusuQ2nftta37r3+9zirUcBHxg3antbtJklTKkJAklTIkJEmlDAlJUilDQpJUak0f3dSN+fnk1JlZZufOMbx+HZs3DTM0FHWXJUl9qVEhMT+fHDl+mr0Hp5ieOcvE6Eb275lkx/iIQSFJS2hUd9OpM7PnAwJgeuYsew9OcerMbM2VSVJ/alRIzM6dOx8QC6ZnzjI7d66miiSpvzUqJIbXr2NidOOb2iZGNzK8fl1NFUlSf2tUSGzeNMz+PZPng2Jhn8TmTcM1VyZJ/alRO66HhoId4yMcum2XRzdJUhsaFRLQCoqxkQ11lyFJa0KjupskSZ0xJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSqcpCIiJ+JCK+FRH/JyKejYj/WLRfERHfjIjnI+L3I2K4aN9QPD9avL6tqtokSe2pckvideB9mfkuYCdwfUS8B/gvwN2ZuR2YAW4tpr8VmMnMnwTuLqaTJNWospDIlv9XPL2kuCXwPuAPivYDwA3F493Fc4rXr4sIx8uQpBpVuk8iItZFxFPACeBR4P8Cf52Zc8Uk08DW4vFW4CWA4vXXgM1V1idJWl6lIZGZ5zJzJzABXAP89FKTFfdLbTXkhQ0RsS8ipiJi6uTJk6tXrCTpLXpydFNm/jXwdeA9wNsiYmFgwQngleLxNHAZQPH6jwGvLvG37s3MycycHBsbq7p0SWq0Ko9uGouItxWPNwI/DzwHPA78UjHZLcDDxePDxXOK1/8kM9+yJSFJ6p0qhwrfAhyIiHW0wuihzPxyRHwPeDAi/jPwv4H7iunvA343Io7S2oK4ucLaJEltqCwkMvNp4N1LtL9Aa//Ehe1/C9xUVT2SpM55xrUkqZQhIUkqZUhIkkoZEpKkUoaEJKmUISFJKmVISJJKVXky3cCYn09OnZlldu4cw+vXsXnTMENDDlArqfd6/XtkSKxgfj45cvw0ew9OMT1zlonRjezfM8mO8RGDQlJP1fF7ZHfTCk6dmT3/hQBMz5xl78EpTp2ZrbkySU1Tx++RIbGC2blz57+QBdMzZ5mdO1dTRZKaqo7fI0NiBcPr1zExuvFNbROjGxlev66miiQ1VR2/R4bECjZvGmb/nsnzX8xCH+DmTcM1Vyapaer4PXLH9QqGhoId4yMcum2XRzdJqlUdv0eGRBuGhoKxkQ11lyFJPf89srtJklTKkJAklTIkJEmlDAlJUilDQpJUypCQJJUyJCRJpQwJSVKpykIiIi6LiMcj4rmIeDYiPla0fzIiXo6Ip4rbhxa9566IOBoRRyLiA1XVJklqT5VnXM8Bn8jM70TECPBkRDxavHZ3Zv63xRNHxJXAzcBVwI8DX4uIn8pMh1uVpJpUtiWRmccy8zvF49PAc8DWZd6yG3gwM1/PzO8DR4FrqqpPkrSynuyTiIhtwLuBbxZNd0TE0xFxf0SMFm1bgZcWvW2a5UNFklSxykMiIn4U+EPg45n5Q+Ae4J3ATuAY8NmFSZd4ey7x9/ZFxFRETJ08ebKiqiVJUHFIRMQltALii5n5JYDMPJ6Z5zJzHtjPG11K08Bli94+Abxy4d/MzHszczIzJ8fGxqosX5Iar8qjmwK4D3guMz+3qH3LosluBJ4pHh8Gbo6IDRFxBbAd+FZV9UmSVlbl0U27gF8BvhsRTxVtvwZ8JCJ20upKehH4KEBmPhsRDwHfo3Vk1O0e2SRJ9aosJDLzT1l6P8Mjy7znU8CnqqpJktQZz7iWJJUyJCRJpQwJSVKpKndcN9r8fHLqzCyzc+cYXr+OzZuGGRpaaheNpKZaC78ThkQF5ueTI8dPs/fgFNMzZ5kY3cj+PZPsGB/puwVAUj3Wyu+E3U0VOHVm9vwXDzA9c5a9B6c4dWa25sok9Yu18jthSFRgdu7c+S9+wfTMWWbnPO1DUsta+Z0wJCowvH4dE6Mb39Q2MbqR4fXraqpIUr9ZK78ThkQFNm8aZv+eyfMLwEJf4+ZNwzVXJqlfrJXfCXdcV2BoKNgxPsKh23b19VELkuqzVn4nDImKDA0FYyMb6i5DUh9bC78TdjdJkkoZEpKkUoaEJKmUISFJKmVISJJKGRKSpFKGhCSplCEhSSplSEiSShkSkqRShoQkqZQhIUkqVVlIRMRlEfF4RDwXEc9GxMeK9rdHxKMR8XxxP1q0R0T8ZkQcjYinI+LqqmqTJLWnyi2JOeATmfnTwHuA2yPiSuBO4LHM3A48VjwH+CCwvbjtA+6psDZJUhsqC4nMPJaZ3ykenwaeA7YCu4EDxWQHgBuKx7uBg9nyDeBtEbGlqvokSSvryT6JiNgGvBv4JjCemcegFSTAO4rJtgIvLXrbdNF24d/aFxFTETF18uTJKsuWpMarPCQi4keBPwQ+npk/XG7SJdryLQ2Z92bmZGZOjo2NrVaZkqQlVHpluoi4hFZAfDEzv1Q0H4+ILZl5rOhOOlG0TwOXLXr7BPBKlfX1k/n55NSZ2b6+jKGkcoO6DlcWEhERwH3Ac5n5uUUvHQZuAT5T3D+8qP2OiHgQ+FngtYVuqUE3P58cOX6avQenmJ45e/6C6DvGRwZiIZMG3SCvw1V2N+0CfgV4X0Q8Vdw+RCsc3h8RzwPvL54DPAK8ABwF9gO3VVhbXzl1Zvb8wgUwPXOWvQenOHVmtubKJLVjkNfhyrYkMvNPWXo/A8B1S0yfwO1V1dPPZufOnV+4FkzPnGV27lxNFUnqxCCvw55x3QeG169jYnTjm9omRjcyvH5dTRVJ6sQgr8OGRB/YvGmY/Xsmzy9kC/2ZmzcN11yZpHYM8jpc6dFNas/QULBjfIRDt+0auCMjpCYY5HXYkOgTQ0PB2MiGusuQ1KVBXYftbpIklWorJCLid9tpkyQNlna3JK5a/CQi1gH/cPXLkST1k2VDIiLuiojTwD+IiB8Wt9O0htJ4eLn3SpLWvmVDIjM/nZkjwG9k5t8tbiOZuTkz7+pRjZKkmrTb3fTliNgEEBH/KiI+FxE/UWFdkqQ+0G5I3AP8TUS8C/j3wA+Ag5VVJUnqC+2GxFwxttJu4POZ+XlgpLqyJEn9oN2T6U5HxF20RnX9Z8XRTZdUV5YkqR+0uyXxy8DrwK9m5l/Suqzob1RWlSSpL7QVEkUwfBH4sYj458DfZqb7JCRpwLV7xvWHgW8BNwEfBr4ZEb9UZWGSpPq1u0/iPwD/KDNPAETEGPA14A+qKkySVL9290kMLQRE4VQH75UkrVHtbkl8JSK+CjxQPP9lWtekliQNsGVDIiJ+EhjPzH8XEb8I/FNa163+M1o7slWj+fnk1JnZgbvIidQPXL9aVtqS+O/ArwFk5peALwFExGTx2r+otDqVmp9Pjhw/zd6DU0zPnD1/ucQd4yONXJCl1eT69YaV9itsy8ynL2zMzClgWyUVqS2nzsyeX4ABpmfOsvfgFKfOzNZcmbT2uX69YaWQ+JFlXtu43Bsj4v6IOBERzyxq+2REvBwRTxW3Dy167a6IOBoRRyLiA+2V31yzc+fOL8ALpmfOMjt3rqaKpMHh+vWGlULi2xGx98LGiLgVeHKF934BuH6J9rszc2dxe6T4e1cCN9O6uNH1wG8XQ3+oxPD6dUyMvjmnJ0Y3Mrzej026WK5fb1gpJD4O/OuI+HpEfLa4/S/g3wAfW+6NmfkE8GqbdewGHszM1zPz+8BR4Jo239tImzcNs3/P5PkFeaHPdPOm4Zork9Y+1683LLvjOjOPA/8kIt4L/EzR/EeZ+ScX8W/eERF7gCngE5k5Q2ssqG8smma6aFOJoaFgx/gIh27b1fijL6TV5vr1hrbOk8jMx4HHV+Hfuwf4T0AW958FfpXWYbVv+WeX+gMRsQ/YB3D55ZevQklr19BQMDayoe4ypIHk+tXS07OmM/N4Zp7LzHlgP290KU0Dly2adAJ4peRv3JuZk5k5OTY2Vm3BktRwPQ2JiNiy6OmNwMKRT4eBmyNiQ0RcAWynNaCgJKlG7Q7L0bGIeAC4Frg0IqaBXweujYidtLqSXgQ+CpCZz0bEQ8D3gDng9sxs3rFmktRnKguJzPzIEs33LTP9p4BPVVWPJKlzjuQqSSplSEiSShkSkqRShoQkqZQhIUkqZUhIkkoZEpKkUoaEJKlUZSfTqT953V41kct99wyJBvG6vWoil/uLY3dTg3jdXjWRy/3FMSQaxOv2qolc7i+OIdEgXrdXTeRyf3EMiQbxur1qIpf7i+OO6wbxur1qIpf7i2NINIzX7VUTudx3z+4mSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklaosJCLi/og4ERHPLGp7e0Q8GhHPF/ejRXtExG9GxNGIeDoirq6qLklS+6rckvgCcP0FbXcCj2XmduCx4jnAB4HtxW0fcE+FdUmS2lRZSGTmE8CrFzTvBg4Ujw8ANyxqP5gt3wDeFhFbqqpNktSeXg/LMZ6ZxwAy81hEvKNo3wq8tGi66aLtWI/r0xK8qpf6ictjb/XL2E1LfcO55IQR+2h1SXH55ZdXWZPwql7qLy6Pvdfro5uOL3QjFfcnivZp4LJF000Aryz1BzLz3syczMzJsbGxSouVV/VSf3F57L1eh8Rh4Jbi8S3Aw4va9xRHOb0HeG2hW0r18qpe6icuj71X5SGwDwB/BuyIiOmIuBX4DPD+iHgeeH/xHOAR4AXgKLAfuK2qutQZr+qlfuLy2HtVHt30kczckpmXZOZEZt6Xmacy87rM3F7cv1pMm5l5e2a+MzP/fmZOVVWXOuNVvdRPXB57r192XKtPeVUv9ROXx94zJLQir+qlfuLy2FuO3SRJKmVISJJKGRKSpFKGhCSplCEhSSplSEiSShkSkqRSniehVedQzuqEy0t/MyS0qhzKWZ1weel/djdpVTmUszrh8tL/DAmtKodyVidcXvqfIaFV5VDO6oTLS/8zJLSqHMpZnXB56X/uuNaqcihndcLlpf8ZElp1DuWsTri89De7myRJpQwJSVIpQ0KSVMqQkCSVMiQkSaU8ukl9wUHeBoPf4+CpJSQi4kXgNHAOmMvMyYh4O/D7wDbgReDDmTlTR33qLQd5Gwx+j4Opzu6m92bmzsycLJ7fCTyWmduBx4rnagAHeRsMfo+DqZ/2SewGDhSPDwA31FiLeshB3gaD3+NgqiskEvjjiHgyIvYVbeOZeQyguH/HUm+MiH0RMRURUydPnuxRuaqSg7wNBr/HwVRXSOzKzKuBDwK3R8TPtfvGzLw3Myczc3JsbKy6CtUzDvI2GPweB1MtO64z85Xi/kREHAKuAY5HxJbMPBYRW4ATddSm3nOQt8Hg9ziYer4lERGbImJk4THwC8AzwGHglmKyW4CHe12b6rMwyNvW0b/D2MgGf1jWKL/HwVPHlsQ4cCgiFv7938vMr0TEt4GHIuJW4C+Am2qoTZK0SM9DIjNfAN61RPsp4Lpe1yNJKucZ11qzPLu3Wn6+AkNCa5Rn91bLz1cL+ulkOqltnt1bLT9fLTAktCZ5dm+1/Hy1wJDQmuTZvdXy89UCQ0Jrkmf3VsvPVwvcca01ybN7q+XnqwWGhNashbN729X0Qzo7nf9OP18NJkNCjdD0QzqbPv/qnvsk1AhNP6Sz6fOv7hkSaoSmH9LZ9PlX9wwJNULTD+ls+vyre4aEGqHbQzrn55OTp1/n5Zm/4eTp15mfz16Uu6JO6/KQVnXLHddqhG4O6ezXnb3d1OUhreqWWxJqjE4viNOvO3u7rcsLAqkbbklIJbrd2dvp+QidTu9OaPWSISGVWNjZu/gHeaWdvZ12BXXTddRNXVK37G6SSnSzs7fTrqBuuo7cCa1ecktCKtHNzt5Ou4K66TpyJ7R6yZCQltHp+EWddgV123XkuErqFbubpFXUaVeQXUfqd25JSKuo064gu47U7/ouJCLieuDzwDrgdzLzMzWXJHWk064gu47Uz/qquyki1gH/A/ggcCXwkYi4st6qJKm5+iokgGuAo5n5QmbOAg8Cu2uuSZIaq99CYivw0qLn00WbJKkG/RYSS+2te9PwlhGxLyKmImLq5MmTPSpLkpqp30JiGrhs0fMJ4JXFE2TmvZk5mZmTY2NjPS1OkpomMvtjfHyAiFgP/DlwHfAy8G3gX2bmsyXTnwR+0OU/dynwV12+dxA0ef6bPO/Q7Pl33lt+IjPb+l92Xx0Cm5lzEXEH8FVah8DeXxYQxfRdb0pExFRmTnb7/rWuyfPf5HmHZs+/8975vPdVSABk5iPAI3XXIUnqv30SkqQ+0uSQuLfuAmrW5Plv8rxDs+ffee9QX+24liT1lyZvSUiSVtC4kIiImyLi2YiYj4jJC167KyKORsSRiPhAXTX2SkR8MiJejoinituH6q6pahFxffH9Ho2IO+uup5ci4sWI+G7xXU/VXU/VIuL+iDgREc8sant7RDwaEc8X96N11liVknnvan1vXEgAzwC/CDyxuLEYSPBm4CrgeuC3iwEHB93dmbmzuA30UWUOIAnAe4vvugmHgX6B1rq82J3AY5m5HXiseD6IvsBb5x26WN8bFxKZ+VxmHlnipd3Ag5n5emZ+HzhKa8BBDQ4HkGyQzHwCePWC5t3AgeLxAeCGnhbVIyXz3pXGhcQymjq44B0R8XSxeTqQm96LNPU7XpDAH0fEkxGxr+5iajKemccAivt31FxPr3W8vg9kSETE1yLimSVuy/2vccXBBdeiFT6Le4B3AjuBY8Bnay22egP5HXdgV2ZeTau77faI+Lm6C1JPdbW+990Z16shM3++i7etOLjgWtTuZxER+4EvV1xO3QbyO25XZr5S3J+IiEO0ut+eWP5dA+d4RGzJzGMRsQU4UXdBvZKZxxced7K+D+SWRJcOAzdHxIaIuALYDnyr5poqVawkC26ktVN/kH0b2B4RV0TEMK0DFQ7XXFNPRMSmiBhZeAz8AoP/fS/lMHBL8fgW4OEaa+mpbtf3gdySWE5E3Aj8FjAG/FFEPJWZH8jMZyPiIeB7wBxwe2aeq7PWHvivEbGTVpfLi8BH6y2nWp0OIDlgxoFDEQGt9f73MvMr9ZZUrYh4ALgWuDQipoFfBz4DPBQRtwJ/AdxUX4XVKZn3a7tZ3z3jWpJUyu4mSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCWgURcUdE/HkxwvCni7bSEYeltcJDYKWLFBEbaJ25vR14jdb4UFcBfw+YB/4n8G8zc+CH59bgadzJdFIFxoG/zMxXI+Kqou21zJwBKE5gk9Yku5ukizcEZETcBTwNPJyZ8zXXJK0KQ0JaJZn5aVrDvWyLiL111yOtBkNCWkWZ+SrwAHB13bVIq8GQkFbHZRHxU8XjnwWWuvqhtOYYEtLqmAMejojvAj8O/E5E3FiMwPmPaY04/NVaK5S64CGw0kWKiG3AlzPzZ2ouRVp1bklIkkq5JSFJKuWWhCSplCEhSSplSEiSShkSkqRShoQkqZQhIUkq9f8BCppdlTAPTQ0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import seaborn as sns\n", "\n", "\n", "def rr(y1, β1, λ):\n", " return np.power(y1 - β1, 2) + λ*(β1**2)\n", "\n", "y1 = 5\n", "λ = 1\n", "β = list(range(-10, 15))\n", "results = [rr(y1, β1, λ) for β1 in β]\n", "βR = y1/(1+λ)\n", "\n", "ax = sns.scatterplot(x=β, y=results)\n", "ax.axvline(x=βR, color='r')\n", "plt.xlabel('β1')\n", "plt.ylabel('Cost');\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### (b) Consider (6.13) with p = 1. For some choice of y1 and λ > 0, plot (6.13) as a function of β1. Your plot should confirm that (6.13) is solved by (6.15)." ] }, { "cell_type": "code", "execution_count": 62, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "np.absolute(y1) <= λ/2\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEKCAYAAADn+anLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAFXBJREFUeJzt3X9s3Hd9x/HXywnOQjDDTV0v5AdBLERrkQid17Flm4LYaKkmpUEra7eRwKoEqa0EG2xr2R+gSRVljJ/TUi1ZqyYTtKoGUaNSAaVrVRENqIOi/sqyRlBaN8EJadRaaYjl3Ht/3PfM1fHXvjvf977fu+/zIVl3/vjOfl/d8yvfz09HhAAAmE1f3gUAAIqLkAAApCIkAACpCAkAQCpCAgCQipAAAKQiJAAAqQgJAEAqQgIAkGpx3gUsxMUXXxxr167Nuwx0kyNHqrfr1+dbB5CjgwcP/iIihhp5bFeHxNq1azU6Opp3GegmmzZVbx99NM8qgFzZ/lmjj6W7CQCQipAAAKQiJAAAqQgJAEAqQgIAkKqrZzcBQNlUKqFTZyY1OXVe/YsXafmyfvX1ObOfR0gAQJeoVEJHxie0fe+oxk6f1arBpdq9dUTrhwcyCwq6mwCgS5w6MzkdEJI0dvqstu8d1akzk5n9TEICALrE5NT56YCoGTt9VpNT5zP7mYQEAHSJ/sWLtGpw6WvaVg0uVf/iRZn9TEICALrE8mX92r11ZDooamMSy5f1Z/YzGbgGgC7R12etHx7Qvhs3MrsJAHChvj5raGBJ535ex34SAKDrEBIAgFSEBAAgFSEBAEiVWUjYXm37EduHbT9t+2NJ+2dsv2j7UPJxdd1zbrV91PYR21dmVRsAoDFZzm6akvSJiPix7QFJB20/lHztSxHxL/UPtn2ppOskXSbpzZK+Z/vtEZHdUkIAwJwyu5KIiOMR8ePk/oSkw5JWzvGUzZLujYhzEfFTSUclXZFVfQCA+XVkTML2WknvkvTDpOlm20/Yvsv2YNK2UtILdU8b09yhAgDIWOYhYfsNkr4h6eMR8YqkOyS9TdIGScclfaH20FmeHrN8vx22R22Pnjx5MqOqAQBSxiFh+3WqBsTXIuKbkhQR4xFxPiIqknbrV11KY5JW1z19laRjM79nROyKiJGIGBkaGsqyfAAovSxnN1nSnZIOR8QX69pX1D1si6Snkvv7JV1ne4ntt0paJ+lHWdUHAJhflrObNkr6kKQnbR9K2j4l6XrbG1TtSnpO0kclKSKetn2fpGdUnRl1EzObACBfmYVERHxfs48zPDjHc26TdFtWNQEAmsOKawBAKkICAJCKkAAApOLQIQDISaUSOnVmsmOnzLWCkACAHFQqoSPjE9q+d1Rjp89On1e9fnigUEFBdxMA5ODUmcnpgJCksdNntX3vqE6dmcy5stciJAAgB5NT56cDombs9FlNThVreRghAQA56F+8SKsGl76mbdXgUvUvXpRTRbMjJAAgB8uX9Wv31pHpoKiNSSxf1p9zZa/FwDUA5KCvz1o/PKB9N25kdhMA4EJ9fdbQwJK8y5gT3U0AgFSEBAAgFSEBAEhFSAAAUhESAIBUhAQAIBUhAQBIRUgAAFIREgCAVIQEACAVIQEASEVIAABSERIAgFSEBAAgFVuFA0AbVCqhU2cmC302RCsICQBYoEoldGR8Qtv3jmrs9NnpU+bWDw90fVDQ3QQAC3TqzOR0QEjS2Omz2r53VKfOTOZc2cJlFhK2V9t+xPZh20/b/ljSfpHth2w/m9wOJu22/VXbR20/YfvyrGoDgHaanDo/HRA1Y6fPanLqfE4VtU+WVxJTkj4REb8l6d2SbrJ9qaRbJD0cEeskPZx8Lknvl7Qu+dgh6Y4MawOAtulfvEirBpe+pm3V4FL1L16UU0Xtk1lIRMTxiPhxcn9C0mFJKyVtlrQnedgeSdck9zdL2htVP5D0JtsrsqoPANpl+bJ+7d46Mh0UtTGJ5cv6c65s4ToycG17raR3SfqhpOGIOC5Vg8T2JcnDVkp6oe5pY0nb8U7UCACt6uuz1g8PaN+NG5nd1Czbb5D0DUkfj4hX7NT/aLN9IWb5fjtU7Y7SmjVr2lUmACxIX581NLAk7zLaLtPZTbZfp2pAfC0ivpk0j9e6kZLbE0n7mKTVdU9fJenYzO8ZEbsiYiQiRoaGhrIrHgCQ6ewmS7pT0uGI+GLdl/ZL2pbc3ybp/rr2rcksp3dLernWLQUAyEeW3U0bJX1I0pO2DyVtn5J0u6T7bN8g6XlJ1yZfe1DS1ZKOSnpV0kcyrA0A0IDMQiIivq/Zxxkk6b2zPD4k3ZRVPQCA5rHiGgCQipAAAKQiJAAAqQgJAEAqQgIAkIqQAACkIiQAAKkICQBAKkICAJCKM64BYBaVSujUmcme2/q7WYQEAMxQqYSOjE9Mn1tdO0Ro/fBA6YKC7iYAmOHUmcnpgJCq51Vv3zuqU2cmc66s8wgJAJhhcur8dEDUjJ0+q8mp8zlVlB9CAgBm6F+8aPq86ppVg0vVv3hRThXlh5AAgBmWL+vX7q0j00FRG5NYvqw/58o6j4FrAJihr89aPzygfTduZHZT3gUAQBH19VlDA0vyLiN3dDcBAFIREgCAVIQEACAVIQEASEVIAABSERIAgFSEBAAgFSEBAEhFSAAAUhESAIBUhAQAIBV7NwHoeRxF2rrMriRs32X7hO2n6to+Y/tF24eSj6vrvnar7aO2j9i+Mqu6AJRL7SjSLTsPaOPnHtGWnQd0ZHxClUrkXVpXyLK76W5JV83S/qWI2JB8PChJti+VdJ2ky5Ln7LRdvtM9ALQdR5EuTGYhERGPSXqpwYdvlnRvRJyLiJ9KOirpiqxqA1AeHEW6MHkMXN9s+4mkO2owaVsp6YW6x4wlbRewvcP2qO3RkydPZl0rgC7HUaQL0+mQuEPS2yRtkHRc0heS9tlGkGbtMIyIXRExEhEjQ0ND2VQJoGdwFOnCdHR2U0SM1+7b3i3pgeTTMUmr6x66StKxDpYGoEdxFOnCNHQlYfs/G2lr4PusqPt0i6TazKf9kq6zvcT2WyWtk/SjZr8/AMymdhTpysHXa2hgCQHRhEavJC6r/ySZefTbcz3B9j2SNkm62PaYpE9L2mR7g6pdSc9J+qgkRcTTtu+T9IykKUk3RQSjSgCQszlDwvatkj4laantV2rNkiYl7ZrruRFx/SzNd87x+Nsk3TZntQCAjpqzuykiPhsRA5I+HxFvTD4GImJ5RNzaoRoBADlpdHbTA7aXSZLtv7L9RdtvybAuAEABNBoSd0h61fY7Jf29pJ9J2ptZVQCAQmg0JKYiIlRdGf2ViPiKpIHsygIAFEGjs5smkkHsD0n6w2R20+uyKwsAUASNXkn8uaRzkv46In6u6pYZn8+sKgBAITQUEkkwfE3Sr9v+U0m/jAjGJACgxzW64vqDqq6AvlbSByX90PafZVkYACB/jY5J/KOk34mIE5Jke0jS9yT9V1aFAcBsOGWusxoNib5aQCROifOxAXRY7ZS52iFCtR1d1w8PEBQZafQP/bdtf8f2h21/WNK3JD2YXVkAcCFOmeu8+fZu+k1JwxHxd7Y/IOkPVN276X9UHcgGgI7hlLnOm+9K4suSJiQpIr4ZEX8bEX+j6lXEl7MuDgDqccpc580XEmsj4omZjRExKmltJhUBQApOmeu8+Qauf22Ory2d42sA0HacMtd5811JPG57+8xG2zdIOphNSQCQjlPmOmu+K4mPS9pn+y/1q1AYkdSv6vGjAIAeNmdIRMS4pN+3/R5J70iavxUR/515ZQCA3DW0mC4iHpH0SMa1AAAKhlXTAIBUjW7L0TPY9wUAGleqkGDfFwBoTqm6m9j3BQCaU6qQYN8XAGhOqUKCfV8AoDmlCgn2fQGKp1IJnZw4pxdPv6qTE+dUqUTeJaFOqQau2fcFKBYmkxRfqa4kJPZ9AYqEySTFl1lI2L7L9gnbT9W1XWT7IdvPJreDSbttf9X2UdtP2L48q7oAFAeTSYovyyuJuyVdNaPtFkkPR8Q6SQ8nn0vS+yWtSz52SLojw7oAFASTSYovs5CIiMckvTSjebOkPcn9PZKuqWvfG1U/kPQm2yuyqg1AMTCZpPg6PXA9HBHHJSkijtu+JGlfKemFuseNJW3HZ34D2ztUvdrQmjVrsq0WQKaYTFJ8RZndNNv/EbPOg4uIXZJ2SdLIyAhz5YAuV5tMgmLq9Oym8Vo3UnJ7Imkfk7S67nGrJB3rcG0AgBk6HRL7JW1L7m+TdH9d+9ZkltO7Jb1c65YCAOQns+4m2/dI2iTpYttjkj4t6XZJ9yVnZD8v6drk4Q9KulrSUUmvSvpIVnW1gu3FAZRVZiEREdenfOm9szw2JN2UVS0LwYpQAGVWuhXXzWJFKIAyIyTmwYpQAGVGSMyDFaFA49jRtfcQEvNgRSjQmNr43ZadB7Txc49oy84DOjI+QVB0uaIspissVoQCjUkbv9t340YWy3UxQqIBrAgF5sf4XW+iuwlAWzB+15sICQBtwfhdb6K7CUBbMH7XmwgJAG3D+F3vobsJAJCKK4mMsCkggF5ASGSATQEB9Aq6mzLApoAAegUhkQEWFQHoFYREBlhUBKBXEBIZYFERegE7ukJi4DoTLCpCt2PyBWq4kshIbVHRysHXa2hgCW8sdBUmX6CGkABwASZfoIaQAHABJl+ghpAAcAEmX6CGgWsAF2DyBWoIiYJgrycUDTu6QiIkCoHphgCKijGJAmC6IYCiIiQKgOmGAIqKkCgAphsCKKpcQsL2c7aftH3I9mjSdpHth2w/m9wO5lFbHphuiE5gLya0Is+B6/dExC/qPr9F0sMRcbvtW5LP/yGf0jqL6YbIGpMj0KoidTdtlrQnub9H0jU51tJx7PWELDE5Aq3KKyRC0ndtH7S9I2kbjojjkpTcXpJTbUDPYXIEWpVXd9PGiDhm+xJJD9n+30afmITKDklas2ZNVvUBPaU2OaI+KJgcgUbkciUREceS2xOS9km6QtK47RWSlNyeSHnurogYiYiRoaGhTpUMdDUmR6BVHb+SsL1MUl9ETCT33yfpnyTtl7RN0u3J7f2drq3bsJUHGsXkCLQqj+6mYUn7bNd+/tcj4tu2H5d0n+0bJD0v6docausazFZBs9iLCa3oeEhExE8kvXOW9lOS3tvperpV2myVfTdu5A8BgLYp0hRYNIHZKgA6gZDoUmzlUW6snkanEBJditkq5VUbj9qy84A2fu4Rbdl5QEfGJwgKZILzJLoUs1XKi/EodBIh0cWYrVJOjEehk+huAroM41HoJEKiZBjw7H6MR6GT6G4qERbg9QbGo9BJXEmUCNtF9w62lkenEBIlwoAngGYREiXCgGcxMU6EIiMkSoQBz+JhYRyKjoHrEmHAs3hYGIeiIyRKppUFeJxbkR3GiVB0hATmxLTZbHGsKIqOMQnMiWmz2WKcCEXHlQTmRHdIthgnQtEREpgT3SHZY6NGFBndTZgT3SHNY90DeglXEphTK90hZZ4NxUA/eg1XEphXM/sElX1xGAP96DWEBNqq7H8kGehHryEk0FZl/yPJ/ljoNYQE2qrVP5JFHextti4G+tFrGLhGW9X+SM4cuJ3rj2RRB3tbqYt1D+g1jijGv9haMTIyEqOjo3mXgRmand10cuKctuw8cMFajEw2udu0qXr76KPzPrSjdQEdZPtgRIw08liuJNB2zS4Oa3UcI+uptmUfXwEkQgIF0Mqq7k50UbHaHCjgwLXtq2wfsX3U9i1514PstTLY28pU20olNHm+onNTFQahgQYV6krC9iJJ/ybpTySNSXrc9v6IeCbfypClVgZ7m+0Kql15/PLFV3Ru6rw+ufMAg9BAA4p2JXGFpKMR8ZOImJR0r6TNOdeEDmhmVbfU/FTb2pXHuSREGl3k12xdQK8pWkislPRC3edjSRvwGs12BTEIDbSmUN1Nkmb7Z9prOo5t75C0Q5LWrFnTiZpQQM12BbESGmhN0a4kxiStrvt8laRj9Q+IiF0RMRIRI0NDQx0tDsXSTFdQ7cpjSRIKDEIDjSnalcTjktbZfqukFyVdJ+kv8i0JvaB25TG18o2KkPbduJFBaKABhQqJiJiyfbOk70haJOmuiHg657LQI/r6rP5F1YtnVkwDjSlUSEhSRDwo6cG86wAAFG9MAgBQIIQEACAVIQEASEVIAABSdfV5ErZPSvpZi0+/WNIv2lhOtynz6y/za5fK/fp57VVviYiGFpp1dUgshO3RRg/d6EVlfv1lfu1SuV8/r7351053EwAgFSEBAEhV5pDYlXcBOSvz6y/za5fK/fp57U0q7ZgEAGB+Zb6SAADMo3QhYfta20/brtgemfG1W5OztY/YvjKvGjvF9mdsv2j7UPJxdd41Za3MZ6jbfs72k8nvejTverJm+y7bJ2w/Vdd2ke2HbD+b3A7mWWNWUl57S+/30oWEpKckfUDSY/WNti9VdWvyyyRdJWlncuZ2r/tSRGxIPnp6Y8W6M9TfL+lSSdcnv/cyeU/yuy7DNNC7VX0v17tF0sMRsU7Sw8nnvehuXfjapRbe76ULiYg4HBFHZvnSZkn3RsS5iPippKOqnrmN3sEZ6iUSEY9JemlG82ZJe5L7eyRd09GiOiTltbekdCExh7Ker32z7SeSy9OevPSuU9bfcU1I+q7tg8kxwGU0HBHHJSm5vSTnejqt6fd7T4aE7e/ZfmqWj7n+1Tjv+drdaJ7/FndIepukDZKOS/pCrsVmryd/x03YGBGXq9rddpPtP8q7IHRUS+/3wh061A4R8cctPG3e87W7UaP/LWzvlvRAxuXkrSd/x42KiGPJ7Qnb+1Ttfnts7mf1nHHbKyLiuO0Vkk7kXVCnRMR47X4z7/eevJJo0X5J19lekpyxvU7Sj3KuKVPJm6Rmi6qD+r1s+gx12/2qTlTYn3NNHWF7me2B2n1J71Pv/75ns1/StuT+Nkn351hLR7X6fu/JK4m52N4i6V8lDUn6lu1DEXFlRDxt+z5Jz0iaknRTRJzPs9YO+GfbG1TtcnlO0kfzLSdbJT9DfVjSPttS9X3/9Yj4dr4lZcv2PZI2SbrY9pikT0u6XdJ9tm+Q9Lyka/OrMDspr31TK+93VlwDAFLR3QQASEVIAABSERIAgFSEBAAgFSEBAEhFSABtYPtm2/+X7DD82aQtdcdhoFswBRZYINtLVF25vU7Sy6ruD3WZpN+QVJH075I+GRE9vz03ek/pFtMBGRiW9POIeMn2ZUnbyxFxWpKSBWxAV6K7CVi4Pklh+1ZJT0i6PyIqOdcEtAUhAbRJRHxW1e1e1trennc9QDsQEkAbRcRLku6RdHnetQDtQEgA7bHa9tuT+78rabbTD4GuQ0gA7TEl6X7bT0p6s6T/sL0l2YHz91Tdcfg7uVYItIApsMAC2V4r6YGIeEfOpQBtx5UEACAVVxIAgFRcSQAAUhESAIBUhAQAIBUhAQBIRUgAAFIREgCAVP8Pt/fr56r8fwoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def lasso(y1, β1, λ):\n", " return np.power(y1 - β1, 2) + λ*(np.absolute(β1))\n", "\n", "y1 = -1\n", "λ = 3\n", "β = list(range(-10, 15))\n", "results = [lasso(y1, β1, λ) for β1 in β]\n", "\n", "if y1 > λ/2:\n", " print('y1 > λ/2')\n", " βL = y1 - λ/2\n", "if y1 < -λ/2:\n", " print('y1 < -λ/2')\n", " βL = y1 + λ/2\n", "if np.absolute(y1) <= λ/2:\n", " print('np.absolute(y1) <= λ/2')\n", " βL = 0\n", "\n", "ax = sns.scatterplot(x=β, y=results)\n", "ax.axvline(x=βL, color='r')\n", "plt.xlabel('β1')\n", "plt.ylabel('Cost');" ] }, { "cell_type": "code", "execution_count": 75, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHr1JREFUeJzt3X9wnNV97/HPdyWtvJaMpcoyv2wwdalTx+OMI4WfnY4TmkBTXzxgaJsUOw0/bF9uLk0nEKelvgl12omBJJPchJEhkPAzJcEwUGhiA41vpoYANgZfI8cQh+baBLBwLFuWZa2kPfcPeR92pd3Vo/31PLv7fs1kxtI+u3uy2M93z/d7zveYc04AAEhSJOgBAADCg6AAAPAQFAAAHoICAMBDUAAAeAgKAAAPQQEA4CEoAAA8BAUAgKc+6AFM1owZM9ycOXOCHgYAVJTt27e/55xrn+i6igsKc+bM0bZt24IeBgBUFDP7jZ/rSB8BADwEBQCAh6AAAPAQFAAAHoICAMBTcauPAKCWJBJOB/vjig+PKFpfp7amqCIRK9n7ERQAIKQSCac97/bpuvu2af+hAc1qjemuFZ2ad/K0kgUG0kcAEFIH++NeQJCk/YcGdN1923SwP16y9yQoAEBIxYdHvICQtP/QgOLDIyV7T4ICAIRUtL5Os1pjab+b1RpTtL6uZO9JUACAkGpriuquFZ1eYEjWFNqaoiV7TwrNABBSkYhp3snT9Nj1F7L6CAAwGhjapzWW7/3K9k4AgNAjKAAAPAQFAICHoAAA8BAUAAAeggIAwENQAAB4CAoAAA+b1wAgJMp9dkImBAUACIEgzk7IhPQRAIRAEGcnZEJQAIAQCOLshEwICgAQAkGcnZAJQQEAQiCIsxMyodAMACEQxNkJmRAUACAkyn12QsYxBPruAIBQCTQomNlsM/uZme02s9fM7G+DHA8A1Lqg00fDkr7gnHvZzKZJ2m5mTzvnugMeFwDUpEBnCs65t51zL5/4c5+k3ZJOD3JMAFDLQlNTMLM5khZJeiHDYyvNbJuZbevp6Sn30ACgZoQiKJhZs6SNkj7vnDsy9nHn3J3OuU7nXGd7e3v5BwgANSLomoLMrEGjAeFB59yjQY8HADIJQwfTcgg0KJiZSbpb0m7n3DeCHAsAZBOWDqblEHT66EJJyyV9zMxeOfG/TwY8JgCQNBoMevoGtb/3mN45fFztzaMby4LqYFoOgc4UnHP/Kam6wiyAipdIOPUOxPV273GtemC7NztYv2yhbt+0Rzv29QbSwbQcgp4pAECoJFNFr+477AUEaXR2sGbjTq1ePFdSMB1My4GgAAApkofdTI3WZTzfoCXWEFgH03IgKABAiuRhN70DQxnPN5jVGtNj119YlUVmiaAAAGmSh910bdmr9csWjjvf4NTpMbVPa6zKgCCFYJ8CAIRJ8rCb6+7bpts37dG6pQt01owmTW2s04ym6g0GSQQFAEgRlsNugkJQAIAxwnDYTVCoKQAAPAQFAICH9BGAmlErTe0KQVAAUBMmampHwBhF+ghATUjuVE5tW5FsapcMGJfdsVUXrv+ZLrtjq/a826dEwgU86vIjKACoCcmdyqmSTe1yBYxaQ1AAUBOSO5VTJZva5QoYtYagACB0kucYvHXomHr6BouSxknuVB7btqKtKZozYNQac66ycmadnZ1u27ZtQQ8DQImU8pSzbMXkWjhZzcy2O+c6J7yOoAAgTHr6BnXZHVvT0jnJzqSl3GVc7auP/AYFlqQCCJWg8vu13NoiFTUFAKFCfj9YBAUAgUstLNdFlLUgjNIjfQQgUJmKvPddfY4evf4CDQ0nqjK/H2bMFAAEKtPGsRX3vCiT6fTWqVV9ylkYERQABIqNY+FCUAAQKArL4UJQABCoXDuNC1WKndHVjkIzgMA11ke0bukCTY3W6Vh8RI31hX9frYVdyqVAUABQsEJ2Ax/sj2vFPS8WfQdzts6npd4ZXekICgAKUug38vjwiNqbG7V2yXy1xBrUOzCkri17Cy40U8DOD0EBQEEK/UYei9bpi5fM002P7PSCym1XLFQsWlihOVnAHjsDoYCdG4VmAAUp9Bv5cMJ5ASH53Jse2anhAovCpSxgVzNmCgAKUug38qHhRMagMjSc8H7Op2YRiZjmnTxNj11/YdV2Pi2FwGcKZnaPmR0ws11BjwXA5BX6jXyifQqFnJ+c7HzKzmj/Aj9Pwcz+RNJRSfc55xZMdD3nKQDhU8jqo0yF6g3LOzRv5jTV10cCO1+h2lTMeQrOuZ+b2ZygxwEgf4WcRRCJmM5ub9ZD156rA32DOtgf17eeeV1/9/F5mnfyNFYRlVngQQEADg0M6dPfeyHt5t/9dp8eu/5CVhGVWeA1BT/MbKWZbTOzbT09PUEPB0CR5ZoNsIqovCpipuCcu1PSndJoTSHg4QAVId88f6H1gXyem2s2wCqi8qqIoABgcvLdZVzI7uRCnpucDYx9bnI2wPnJ5ROG1Uc/lLRY0gxJ70r6snPu7mzXs/oImFi+K3YKWelT6CqhQmYomFglrT76VNBjAKpNvit2/D4v9QYei9ZpOOE0EB/R2iXz1bVlr3bs6/X9nknMBsIh8KAAoPjyXbHj53mpaaL25sZxfYvWL1uo2zft0Y59vawSqkAVsfoIwOTku2LHz/N6B+J65/Bxff3KD+nWKxbq+1vfTOtbtGbjTq1ePJdVQhUq8JrCZFFTAPyZKEef7fHU35uZ6kyKRCLezX3320e06oHtGWcGST+/abFi0XrqAiFSMTUFAKWRK0c/0UqhtqZoxsfbmqNeQJDenxmsXTJfq+7fLml0dhGL1lMfqFCkj4Aql+mc4mxnIBzsj0tKPyNh0ewWrV0yX/2Dw4oPJ9TenH6z339owJtFkDKqfMwUgAAVYxlmrtfINiNobqzLucoouQpp0ewW3XjxPK3ZmH4Azq0/fT9dNKs1ptNaYtq65qMsJa0CzBSAgBTSEtrva2SbEQwOu5ztqpOrkFYvnusFhOTzb3pkp2646GzvORuWd6jOpBHn5FRZNUqMR1AAAjJRCqcYr5Ft38FAfFjrly3MusoouQqprSma8flzZzZr65qP6ofXnafhkYSWdT2vP7l1iy6/47lJBzaEC+kjICDFaAk90Wtk23fw28PH1bVlr9Yuma+WWINmtcZ06vSYl/ZJ9ht658jxjM+PNYzOKHa9dVhrH9+V9/nMCB9mCkBAJjpxrJDXaDhxOE18eEQPXXuuPjF/pvfYdz+9SE3ROn3pzz4gSbr7P3/tNZ5LFYmYTjlpStZ9C/HhEU2N5q5NoPIwUwACMlETuMm8xjef3qNlHbPV1hTVzGmNGhwa0afuesF73a6rOrRu6QJFG0xvHRrUlx79v2mPtcYaMr5+rg6l0fo6HYuPcNZBlWHzGhCgYqw+Gh5OaM+BPq26f3vOFULrli7QzJMa9a1nXtfm7gPe8/M92jKRcPqvg/1698jxtDYXqUdpIjz8bl7jvxoQoGIcLH9oYMgLCNL7K4RWL57rXbP/0ICmRuu06v7tWtYxO+35+aZ7IhHTnLYm/dGp0/TgtefqkdXna+2S+frWM6/rjZ6jFJsrFOkjoMJlKza3pKSEZrXG1DswlLbRLPWxfNM9kYhpaET66yxHaVJsrjzMFIAKl63YfCw+4v15/bKF6tqyV7NaY5o5rbGoR1sWYxUVwoOZAlChkvWIRCKhDcs70moKd63o1MyTotpy42K9+V6/bt+0Rz1HB3XXik6dNj1W1KMt823TjXDyVWg2s0bn3OBEvysHCs3A+PYVn5g/U//45/NVd2JVUKaOp6VqQVHIMZwoH7+FZr9B4WXn3Icn+l05EBRQ7fzcyAs9+rLYOEoz/IrSOtvMTpF0uqSYmS2SlPyvfJKkqQWPEkAav9+6w5bH5yjN6jFRTeFiSX8jaZakb6T8vk/SP5RoTEDNytbL6LHrL1RbUzTt8Bvy+CiFnEHBOXevpHvNbJlzbmOZxgTUrFwzgLE1hK6rOrT6ge1574YGMvG1+sg5t9HM/lzSByVNSfn9P5VqYEAtSebkR5zT9//mI/r2s2+k7UaWlDaDSO5I/tGq8+WcI4+PovG1T8HMuiT9paT/qdG6wpWSzizhuICqkOnUs0zX/NfBfu1667De7j0uSfrKpfO1aHaL15soufEs1ebuA3LOFbQbGhjL7z6FC5xzC81sp3PuFjP7uqRHSzkwoNL5LRr3DsT17pHjXgvqZO+irqs+rBEnfeWJXVrWMTtjDaGB/kIoMr9/o5J/E4+Z2WmShiSdVZohAeHm59u/5P8QnYH4iNdQLnndTY/s1HDCyTmnzd0H1LVl77hDcb75Fx9SPbMDFJnfmcKTZtYi6TZJL0tykr5XslEBITWZjVrZisaJRMI76yBaX6cR5zJf55ymNNRrVmtMO/b16vZNe7R2yXy1NUXVMjWqnr7jGoiPSE0l/7+NGuJrpuCcW+ec6z2xAulMSR9wzq0t7dCA8JnMEZrJ9g+LZrdow/IOPbzyPP141fk62B9PO1NZThl7FyWc1Bpr0IarOrzAsO7Jbg0OJ3TTj19VQ12EJagoOt+9j8zsAklzks8xMznn7ivRuIBQmsymsbamqH543bkaGnH6XX9cB/vjOnX6FK0c0+b6q091j1teun7ZQv3zU93658sW6tSWKVq3dIGmRuvUOzDk9TGaOa2RJagoOl9BwczulzRX0iuSkn/7nSSCAmrKZJu/HTk+nNaoruuqDrU3N6Y9f3P3AX350g965yUnb/w79vXqy/9tRKdOj+mU6VPSUlYblnfotJQzlYFi8TtT6JQ031XaMW1AFvn26pnMEZoH++PjDr9Z/cB2rVu6QJ/9wUvedbNaY5KT1j3ZnTHY5DoSEyg2v0Fhl6RTJL1dwrEAZVFIV8/J3KCzpZrObJvqzTaSy0+nNkYytr9OBht6C6FcJmqI928aTRNNk9RtZi9K8tplO+cuLe3wgOLL1V/Iz43X7w06W6qpoS7i1QiOxUd08klTNH1Ko6af0shsAIGbaKZwe6kHYGaXSPqWpDpJ33POfa3U74nalUg4xYdH9PUrP6TegSF1bdmrHft6S9JhNFOqacPyDp3cHFVdZJqGRhJqqItoZvP7u5HbpzV6qa23Dw8QHFB2EzXE+z9+XsTMnnfOnT/ZNzezOknflfRxSfslvWRmTzjnuif7WsBEMqWN1i9b6K3mKfbyzkypptZYg97oOZo1dcWBNQhasfbIT5n4kozOkfQr59yvnXNxSf8qaWmRxgSkyZQ2WrNxp2678kN64Jpz5eSy7k7OVzLVlOxPdGhgKOc+h8nsgwBKoVhBId9/SadL2pfy8/4Tv0tjZivNbJuZbevp6cnzrVDrshV+e4/Ftfj2Lbr8jue0592+ogcGP2NIpq7CdngOak/Q3bQyzYfH/Yt0zt3pnOt0znW2t7eXYVioRsnCb6pZrTHvW3g5vpVnG0MydTXR40Cp+W2d/Tkza811SZ7vv1/S7JSfZ0n6bZ6vBeSULPymNpVbv2yhurbs9a5J/Vbut+21n+Z4ucaQuvR0oseBUvO7T+EUjRaBX5Z0j6RNYzayLc/z/V+SdLaZnSXpLUl/JenTeb4WkNPYwq+Z6StP7PIOs5He/1aeLPh+8+k9WtYxW21NUQ3Eh3Xa9JjqT7SrzlYUPvmkRg3EMy8rnWifAxvVEDS/DfH+UdLZku7W6JnNb5jZv5jZ3BOP78rnzZ1zw5I+J2mTpN2SfuScey2f1wL8SBZ+T50e0/GhEX32wrMyfis/2B/XN5/eo89ccJbWPdmtK7qe16e/94L2HHi/5pCtKPzqvsNes7tMNYqxxeexN/yJHgdKyXdDPOecM7N3JL0jaVhSq6RHzOxp59wX8x2Ac+7fJf17vs8H8nGwP64V97yo9uZGr+fQ6Eay0ZtwfHhEyzpma83G9HMOVt2/XY9ef4FmTpuStSg8NVrn/Xkym+KAMPDbEO8GSZ+R9J5Gz1G4yTk3ZGYRSW9IyjsoAEFIJBJpDei+9pNfase+Xm1d81GpSV7aJtNN/9jgiBJNLuuO5d6BobTrWTmESuJ3pjBD0uXOud+k/tI5lzCzJcUfFlA6iYTTe/1xrwFdsuB873NvqqE+cqJgnFD7tMaMN/033+tXU2N9xh3Lt12xULf+dE/a9ckaRT4N+IBy8xUUnHP/K8dju4s3HKD0MnUvXbNxpx669lwdPT6sFfe8qP2HBvSVJR8Yd87BHX/9YT3w/G/0hyc3jysKN9RHdPT4sHqOjrYHS9YoWmMN7FJGxfBdUwCqRbZagJlpxT0veI+d2tqkbz/7elqa6Tv/8YY+dc6Z3r6Bsc3xZjS5cSuHCm3AB5QTQQFVw2+KJlstIGJK+11LrEGbuw9oc/eBtOf/wyfnZ903kKmDKruUUUmC3tEMFEVyz0Dq2cfZWlZk2yAWi6bvJu4dGMq4u3hqY92k0j7sUkYlsUo7TK2zs9Nt27Yt6GEgINlmAz19g7rsjq3jvv1nS9Gkvk5DfUT1EdNAfEQjCaevPtWtzd0H9In5M3XDRX+YVlPIpxZA51OEgZltd851TnQd6SNUjFw318mmaJJpnkyvuWF5h9YtXaBIJKLWWEPBu4vZpYxKQvoIFSGRcHrnyHH1Dw5r7ZL5WjS7Ja2BXb4pmkxF4FX3b1ckElH7tEbV10eKsruYXcqoFMwUEFrJFE8ikdB7KctIUw/H2bGvV4lEQg31EW24qkOrHsh8xnE2FIGBdAQFhFJqWmftkvla92T3uJYUX7xknr6/9U0vYLQ3N2rd0gU6a0aTpjbWaUbTxN/Is61EogiMWkX6CKGUmtZpiTWovblRN148T+ue7NZf3vkLrX18l5ob63XLpR/0ZhA79vXqsz94SVfd/YJM5itFk20lUmusYVItsYFqwUwBoZRM6yya3aLfa4rqhovOHtec7r8/+LIeXnmer/RPtlVLmdppR+tMew70paWrWC2EWsFMAaEUra/TJ+bP1I0Xz9Ntm36pM9qmZrz5jzhNWGCeaA9DJGJqa4rqyPFh/cWG57X9//WOa4PBOcmoFQQFhEbqKWZOTrdc+kHd+9ybWtYxW/URy3jzn9IQmfCksmxtJt7rH8x4TUusgeIzahbpI4RCpv0CP/jsR3TNH/++vvDjV9Xe3Kjbrliomx7ZmZbSmdHUqBlNjTn3AGRbYZRsgZ08PyF5TXInM8Vn1CKCAkIh07f5fb8b0NrHd2n/oQHtPzSgW3+6R+uWLtDc9iZNiaavLsrVWC7bCqNkC+z2aY1p13Rt2av1yxZ6NQzOSUYtIX2EUMj0bX5qtC7td8nVRQf6BnXwqP/8fltTVBuu6khLMa1ftlDffvYNLyWUugppx75e3fvcm3ro2nO1dc1H9dj1F1JkRs1gpoBQyPRt/lh8JOM3/IP9cX3+4Vd8t56OREyntkzRuqULNDVap96BId2+aY96jg6mtcCmFQXATAEhkWm/wJltU8f9bv2yherasnfShd+WWFSnTJ+iL/z4Va26f7t6jg6OSwnRigKgSypCJNNeAkl658hx/bZ3QAf74+raslc79vXm7IA6mdfnxo9aQZdUhFKuG3PqATWp18WidWpqrNfnH36loMJvpgNwAKQjKKBs/J4rkOm6+64+R49ef4GGhhN8ywdKiJoCyibbJrKxO4UzXbfinhdlMvL9QIkRFFA2fttU084aCA5BAWXj9yCcUp9pnNpOgw6oQDqCAkoi0403W5vqtqbouL5H9119Ts5+RoWMK1dzPKDWsSQVRZeroCwp47LTTNeffFKjBuLFXT7a0zeoy+7YOm5D3GSXtwKVxu+SVGYKKLqxheL25ka9c/i49vce08H+uNqaomkF42wF6JGEil5Ypl4B5MaSVOSUz4av1BvvotktuvHieeOay6UuQ811o+7pGyzqZjOO3wRyY6aArLLl33/Xn7tIm1ooXr147rgT08YuQ81WWB5JuKLn/nPVNQAEGBTM7Eoze83MEmY2YZ4L5ZctrfPqvsM5b9SpN14/B9a0xhrUNaaLaddVHfrqU91FP/0stfEdHVCB8YJMH+2SdLmkDQGOATlkS+u0TG3w/nzdfdvGFWnHdhydKF1zaGBI3372da1dMl8tsQb1DgzJTNrcfWDcexcj90+7CyC7wIKCc263JJnxDS2ssuXfp8catGh2i3bs6816o07eeBMJp7tWdI5bWZSarokPj2hz94G0ILBheQe5fyAA1BSQVVtT1Ls5S++3rv7aT3Zr9eK53u9y3aj9pGsy1RQ2bt837r3J/QOlV9KZgpk9I+mUDA/d7Jx7fBKvs1LSSkk644wzijQ6TCQSMc1oiqaldW7ftEc79vXqmj/+fd836onSNckaROps4u8+Pk9ntzdz6A1QZoFvXjOzLZJudM752pHG5rXyyrbZ6+GV5ykWrVNLrDg3as46AEqLzWsoikxLONcvW6hb/u01vXtksGjvw6lnQDgEVmg2s8sk/W9J7ZKeMrNXnHMXBzUeZJasCfxo1fne6WfJFFL32320hwCqTJCrjx6T9FhQ7w//IhGTc05XdD2f9nvaQwDVh/QRfCm0nTXtqoHKQFCAL4W0h6BdNVA5Al99NFmsPgpOviuEaFcNBM/v6iO6pMK3fNtD0K4aqBykj1BypT5eE0DxEBRQcrSrBioH6SP4lm9NYWzXVHYsA+FFUIAvuc5d9hsYKCoD4Uf6CL5kO3Cn0ENvAIQLQQG+sIIIqA0EBfjCCiKgNhAU4AsriIDaQKEZvrCCCKgNBAX4xgoioPqRPgIAeJgpVBGOtARQKIJClSh0cxkASKSPqgabywAUA0GhSrC5DEAxEBSqBJvLABQDQaFKsLkMQDFQaK4SbC4DUAwEhSrC5jIAhSIoVDH2LQCYLIJClWLfAoB8UGiuUuxbAJAPgkKVYt8CgHwQFKoU+xYA5IOgUKXYtwAgHxSaqxT7FgDkg6BQxdi3AGCyAksfmdltZvZLM9tpZo+ZWUtQYwEAjAqypvC0pAXOuYWSXpf09wGOZUKJhFNP36DeOnRMPX2DSiRc0EMCgKILLH3knNuc8uMvJF0R1FgmwkYwALUiLKuPrpb0k6AHkQ0bwQDUipLOFMzsGUmnZHjoZufc4yeuuVnSsKQHc7zOSkkrJemMM84owUhzS90Itmh2i1YvnquWWIPiwyNKJByzBQBVo6RBwTn3p7keN7PPSFoi6SLnXNYkvXPuTkl3SlJnZ2fZk/nJjWDtzY268eJ5WrNxJ2kkAFUpyNVHl0haI+lS59yxoMbhR3Ij2A0Xne0FBIk0EoDqE+Q+he9IapT0tJlJ0i+cc6sDHE9WyY1gTY119BMCUNWCXH30B0G9dz4iEVOsoV6zWmNpgYF+QgCqSVhWH1UE+gkBqHa0uZgE+gkBqHYEhUminxCAalYTQYGzigHAn6oPCrSoAAD/qr7QTIsKAPCv6oMCZxUDgH9VHxQ4qxgA/Kv6oMDeAgDwr+oLzewtAAD/qj4oSOwtAAC/qj59BADwj6AAAPAQFAAAHoICAMBDUAAAeCzH0cihZGY9kn4T9DhSzJD0XtCDCAk+i3R8Hun4PNKV+/M40znXPtFFFRcUwsbMtjnnOoMeRxjwWaTj80jH55EurJ8H6SMAgIegAADwEBQKd2fQAwgRPot0fB7p+DzShfLzoKYAAPAwUwAAeAgKBTKz28zsl2a208weM7OWoMcUJDO70sxeM7OEmYVuZUW5mNklZrbHzH5lZl8KejxBMrN7zOyAme0KeixBM7PZZvYzM9t94t/J3wY9prEICoV7WtIC59xCSa9L+vuAxxO0XZIul/TzoAcSFDOrk/RdSX8mab6kT5nZ/GBHFagfSLok6EGExLCkLzjn/kjSeZL+R9j+bhAUCuSc2+ycGz7x4y8kzQpyPEFzzu12zu0JehwBO0fSr5xzv3bOxSX9q6SlAY8pMM65n0v6XdDjCAPn3NvOuZdP/LlP0m5Jpwc7qnQEheK6WtJPgh4EAne6pH0pP+9XyP7hI3hmNkfSIkkvBDuSdDVxyE6hzOwZSadkeOhm59zjJ665WaNTwwfLObYg+Pk8alymY/1Y5gePmTVL2ijp8865I0GPJxVBwQfn3J/metzMPiNpiaSLXA2s8Z3o84D2S5qd8vMsSb8NaCwIGTNr0GhAeNA592jQ4xmL9FGBzOwSSWskXeqcOxb0eBAKL0k628zOMrOopL+S9ETAY0IImJlJulvSbufcN4IeTyYEhcJ9R9I0SU+b2Stm1hX0gIJkZpeZ2X5J50t6ysw2BT2mcjux8OBzkjZptJD4I+fca8GOKjhm9kNJz0uaZ2b7zeyaoMcUoAslLZf0sRP3i1fM7JNBDyoVO5oBAB5mCgAAD0EBAOAhKAAAPAQFAICHoAAA8BAUgCIxs5+aWa+ZPRn0WIB8ERSA4rlNo2vQgYpFUAAmycw+cuL8jClm1nSiL/4C59yzkvqCHh9QCHofAZPknHvJzJ6Q9FVJMUkPOOdq/gAZVAeCApCff9Joj6Pjkm4IeCxA0ZA+AvLze5KaNdr3akrAYwGKhqAA5OdOSWs1en7G+oDHAhQN6SNgksxshaRh59xDJ85jfs7MPibpFkkfkNR8olPsNc65musSi8pGl1QAgIf0EQDAQ1AAAHgICgAAD0EBAOAhKAAAPAQFAICHoAAA8BAUAACe/w/gBeOThHa9qAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "np.random.seed(1)\n", "x1 = np.random.normal(0, 1, 100)\n", "y_hat = x1 + np.random.normal(0, 1, 100)*0.2\n", "\n", "ax = sns.scatterplot(x=x1, y=y_hat)\n", "plt.xlabel('x1')\n", "plt.ylabel('y_hat');" ] }, { "cell_type": "code", "execution_count": 76, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHjFJREFUeJzt3XuQXOV55/HvMyPNMBoJayyJiyXZUrkU1lqWDasJjqNNhRiMyYZChbELnHDxJVzCsuxW4lqRwnLVQrLFLZvKVkyBTFiDceIlkSm0mARjB5IqYbBGtkslRGSwHcJwHQsJD7oNo372j5ke9Yy6p/v0nHPe9/T5faoozeXQ83RPz3nO+7zP+x5zd0RERFrVFToAEREpFiUOERFJRIlDREQSUeIQEZFElDhERCQRJQ4REUlEiUNERBJR4hARkUSUOEREJJF5oQPIwtKlS33VqlWhwxARKZQdO3b83N2XNTuuIxPHqlWrGBoaCh2GiEihmNlLrRynUpWIiCSixCEiIokocYiISCJKHCIikogSh4iIJKLEISIiiXRkO65Ip6lUnL0HxhgbP0rPvG6W9PfQ1WWhw5KSUuIQiVyl4ux5Y5SrHhhieN8hVgz08ZUrBjnt5EVKHhKESlUikdt7YGwqaQAM7zvEVQ8MsffAWODIpKyUOEQiNzZ+dCppVA3vO8TY+NFAEUnZKXGIRK5nXjcrBvqmfW3FQB8987oDRSRlp8QhErkl/T185YrBqeRRneNY0t8TODIpK02Oi0Suq8s47eRFPHzdenVVSRSUOEQKoKvLWLaoN3QYIoBKVSIikpASh4iIJKLEISIiiShxiIhIIkocIiKSiLqqRBrQxoJSJHm+X5U4ROrQxoJSJHm/X1WqEqlDGwtKkeT9flXiEKlDGwtKkeT9flXiEKlDGwtKkeT9flXiEKlDGwtKkeT9fjV3z+SBQxocHPShoaHQYUjBqatKiiSN96uZ7XD3wWbHBR9xmNn5ZrbHzF40sxtnOe6TZuZm1vRJiaShurHg8oEFLFvUq6QhUcvz/Ro0cZhZN/Bl4LeAtcCnzWxtneMWATcAz+YboYiIzBR6xHEW8KK7/9Tdx4BvABvqHHcLcDtwOM/gRMqgUnFGRo/wyr6DjIweoVLpvPK1pCt04lgOvFzz+fDk16aY2ZnASnd/NM/ARMqgunDsoru2sf62J7norm3seWNUyUNmFTpx1CvCTb1jzawL+DPgD5s+kNnVZjZkZkMjIyMphihlU6YrcC10lHaE3nJkGFhZ8/kK4NWazxcBpwNPmRnAKcBWM7vQ3ae1Tbn7ZmAzTHRVZRm0dK6ybTWihY7SjtAjju3AGjNbbWY9wKXA1uo33f1td1/q7qvcfRXwDHBc0hBJS9muwLXQUdoRNHG4+zhwPfA48DzwkLs/Z2Y3m9mFIWOTcirbFbgWOko7QpeqcPfHgMdmfO1LDY49O4+YpLyqV+C1yaOTr8C7uozTTl7Ew9et10JHaVnoUpVIVMp4Ba6FjpJU8BGHSEx0BS7SnBKHyAzVK3ARqU+lKhERSUSJQ0REElHiEBGRRDTHIZIB3ctDOpkSh0jKyrZtiZSPSlUiKSvbtiVSPkocIikr27YlUj5KHCIp08aB0umUOERSVsZtS6RcNDkukjJtWyKdTolDJAPatkQ6mUpVIiKSiBKHiIgkosQhIiKJKHGIiEgimhwXaYH2nhI5RolDpAntPSUynUpVIk10wt5TlYozMnqEV/YdZGT0CJWKlzoOmRuNOESaKPreU7GMmGKJQ+ZOIw4ptHpXsGlf1RZ976lYRkyxxCFzpxGHFFajK9jeeV1ccd/3U7uqre49NfPnFGXvqVhGTLHEIXOnEYcUVqMr2Jf2Hkz1qrZ276ltG3+Th69bX6jySiwjpljikLlT4pDCanQFu6Cn+7ivzfWqtrr31PKBBSxb1FuYpAHx7NabRRyabA9DpSoprOoVbG3yWDHQx8Gx6UmiCFe1Wa4TiWW33rTj0GR7OObeeRl6cHDQh4aGQochGZt54jhv7Ul88bfXUgF+NnKA//3dFxh550j0JxOdANszMnqEi+7adtyFw8PXrdfOxG0ysx3uPtjsOI04pLBqr2ArlQo/PzDG79z77NTJ957L1nHq4hNY3Bf3Ku9GczU6Ac5Ok+3haI5DCq0699DV1cU1X9sx7eR7zYM7OFoh6qQBxTwBxjC3oMn2cJQ4ZM5iOIkU8eRbVbQTYLW0dtFd21h/25NcdNc29rwxmvvvPZZJ/zIKnjjM7Hwz22NmL5rZjXW+/wdmttvMdprZd83sAyHilPpiOYkU7eRbq2gnwFgW8hW9TbrIgs5xmFk38GXgY8AwsN3Mtrr77prDfggMuvtBM/t94HbgkvyjlXpiqc8XeZFeLF1PrQo1umvUeaZ5oPyFnhw/C3jR3X8KYGbfADYAU4nD3Z+sOf4Z4LJcI5RZxVIiKtrJd6YinQAbtUFnObpT51lcQpeqlgMv13w+PPm1Rj4P/F29b5jZ1WY2ZGZDIyMjKYYos4mpRFTkRXpFEqK0Fkt5TCaEHnHU+8uuWxw3s8uAQeA36n3f3TcDm2FiHUdaAcrsilYi0g2Z5i7E6C6Wka1MCJ04hoGVNZ+vAF6deZCZnQvcBPyGux/JKTZpQd4nkbmc+FXuSE/epbUQ5TFpLHSpajuwxsxWm1kPcCmwtfYAMzsTuAe40N3fDBCjNJFXiWiuHVztljtiaDfOQ8zPs2idZ50u6IjD3cfN7HrgcaAbuM/dnzOzm4Ehd98K3AEsBP7GzAD+1d0vDBa0BDPXDq52yh1lGaXE/jyL3vzQaUKXqnD3x4DHZnztSzUfn5t7UBKl2U78I6NHpk4oA33z2Xfo3eNOMO2UO2JpN85aEZ5nbJ1nZZ4vC544RFrV6MR/tOJTm92dt/Ykbjjnl7j2wR3HXTm3M5FflknZsjzPtMQ+Qsta6DkOkZbVq3Pfc/k6/vhbu6dOehevWzmVNGD6PEY7K41jajfOUlmeZ1pqR2hnrlzMpgvWcuDIOK//4nBUc0NZ0YhDCqNenbtSqfDt3cd6Jhb3zZ/1yjlpuaNo7cbtKsvzTEt1hHbmysV84eOnsXHLzmmv28kn9nJorHNLWEocUigzT/wjo0emla/2H3o31bbNskzKluV5pqU6Qrv27A9OJQ04NsK9ZcPpfPar2zu2hKVSlRTazPLVlh0vc/dl61Jt24xxRXoWrbONnmfMbbqhVN93S/p7Zr19caeucNeIQwqt3pXyQN/8jr5ybnditp0uoLJPAjdSfd+9/ovDdUe4+w+9O/V5JzYZaMQhhTfzSnnevK7oRghpamchY7uLJ7VHVGNdXcYpJ55wXMPGHZ88g7uf+snUcWk3GcQwAtSIQzLRCT3usT6Hdlpn212noTbd2c0c8c6f18U7h8cZeWdiZ6S0mwxiGQEqcUjqYnlzz0XMz6GdhYztJgDtEdXczIaNpf2eWak0loWaKlVJ6uZa3kg6FM9i6B5ziaadfZvaXaehPaKSy7KZIpYRoEYckppqaefg2DibLljL3U/9hB++vB9o/c2d9Eo/q5FB3n+gScpi7bTOtrtOQ226cYllBKjEIamodwK/7eIzuPPxPfzw5f0tv7mTDsWzGrrn+QfaavKby5xLowQATNvnq95jxrZHVJnFslBTiUNSUe8EvnHLTjZdsJZbHt3d8ps76ZV+ViODPP9AW0l+aYysZiaAmOdxpL5YRoBKHJKKRifwD50y8SZv9c2d9Eo/q5FBkj/QuXZftZL8shhZxTLRKsnEMAJU4pBUNDqB9/XMy3RvqCxHBq38gaZx1d4zr5vz1p7ExetWsrhvPvsPvcuWHS9PS35ZjKximWiV4lHikFSkdQJPOhTPc2RQTxpX7QN984/bCv7uy9Yx0Dd/6pgsRlazPWasa1gkDubeefvODA4O+tDQUOgwSifmk01W9fxX9h1k/W1PHvf1bRt/k+UDC1p6jJHRI1P3E6laMdB33BzHv+w9wEt7D7Kgp5uDY0f5wJIFrFrS33b8jV6TNcsW8sLIO5r7KCEz2+Hug82O04hDUhND7bWRmLuvWi0ZHRmvsOmRXdNO5nPRaLQW89xHzBcnZaLEkSK9qZPJ8/XKsvvqgc+dddxIIEmJrpXkk9XJvF6yj3XuQ11g8dDK8ZS0u4lcWeX9emV5h7vqSOCSzc+w6ZFdHBmvJPr/W1mdnefJPNa7Aca8mr9sWkocZnbcJU29r5WZ3tTJ5P16ZbV1RhrPo9EtbYGprVTMLLeTeazbjMQ6EiqjVktV3wP+QwtfKy29qZPJ+/XKauFUWs+j2eK889aexN2XrZvWedXKybyVcmC9Y2JYZDZTLNttSJPEYWanAMuBPjM7E6i+c04EWmsZKQm9qZMJ8XplMXmf1fOYOZKp3lf9oWs+gru3dDJvZU5gtmNCT4TPlOdqfs1Xzm7WdlwzuxL4DDAI1Pa3jgJfdfdvZhpdm0K042riLplOeb2K3ubbyjExyeOE3invzXak0o7r7vcD95vZxe6+JbXoIpXFJnKd/kZrV6e8Xlk9j+qcxswTulm6W5kUrcSaR8t3zO3IsWhpjsPdt5jZbwP/Fjih5us3ZxVY3rLYRE5ml/brFaq8kMXvvdvgtovPYOOWndN2G+5O8HRaKaOpxHq8oiXTEFrtqrobuAT4L0zMc3wK+ECGceVOXVHF1mnt0F1dXdz/9M/YdMFa/u/Vv8qmC9Zy/9M/o6ur9Q76VrqjYu2gCinWduSYtLTliJntdPczav5dCHzT3c/LPsTk2pnjSKOmLOEUrVbfTFp19na7qopWMkyT5jjS23Kk+td40MzeB+wFVrcbXIw0ZC+2TisvzOXGSzMfp1niVIl1uk6Zf8tSq4njUTNbDNwB/ABw4N7MogogljtrSXs6MfG3c+OlEKOHThyxKJnOLvHuuJMrxk9w97dTCcDsfODPgW7gXne/tc7PewBYx8RI5xJ3/5fZHrPddtxO/AMoi7zLCyHeK7OV45b097D/0Biv7T/MNTMWCWZZYilzWacTtVqqajlxmNmvAauoGaW4+wPtBjj5mN3Aj4GPAcPAduDT7r675pjrgDPc/VozuxS4yN0vme1xta16OeV1Mg91smw0D/fsH32Utw6+y+tvH57aPbcq63meTptbKrtWE0erXVVfA+4E/iPwK5P/zW1P5wlnAS+6+0/dfQz4BrBhxjEbgPsnP/5b4BxL0swuuapUfGp/pZHRI7l2NVXLC8sHFrBsUW9mJ/FGHXj7D41l+twbdfscdbjqgSEW9HTnPs/TaXNL0ppW5zgGgbWe/l2flgMv13w+DHy40THuPm5mbwNLgJ+nHIvMUVnKFvVOlssW9mZeJmo0D+fuDO87xP5D7+Y+z9OJc0vSXKtN4buAUzL4+fX+omYmp1aOwcyuNrMhMxsaGRlJJThJpixrYepd+d9wzpqppAHZPPdGu+hW47n7qZ9w28Vn5LomI891ICFHszJds00O/x8TJ+lFwG4z+z5wpPp9d79wjj9/GFhZ8/kK4NUGxwyb2TzgPcBbMx/I3TcDm2FijmOOcUkb5lq2KEpzQr0r/9VL+4OVbGrjufPxPdyy4XRWL+1nQW83S/vbK9m1+rvIq3W1LKPZomhWqroz45+/HVhjZquBV4BLgd+ZccxW4EomtnH/JPAPGZTMJAVzKVsU6cRQ72TpeOYlm9leozRP3kl/F7O1rjZKQEkvErR/VFxmLVW5+z/O9l/1ODP7Xjs/3N3HgeuBx4HngYfc/Tkzu9nMqqOZvwSWmNmLwB8AN7bzsyR7cylbFK3MNXMifml/b+Ylm9leozQbA9L6XTTaBmZ8vJJ4exhNwsclrXuOn9D8kPrc/THgsRlf+1LNx4eZ2BtLIjeXskW7J4aYyltLFvbwV1d9mG4z+nq6Wdx3LJY04szr5JnWz2mUgB665iOJRw+ahI9LWolDpSMB2l9x286JIZbyVqM4Fvf1pBpnXifPtH5OowT07tFK4sSknR3i0vpWm5KIOkCSaafMFUt5q1kcacWZVwdTWj+n0bqT+d1diXefbdRRFtv8V1m0NOIws+uBr7v7vkaHpBdS8cVyJVwk7ZS5Yql7N4sjzfuS59HBlNbPqR0lLFvYyw3nrGH10n7mdRkPfO4srrjv+4lGD9o/Kh6tlqpOAbab2Q+A+4DHZ3Q2XZ56ZAWmDpD2JD0xxFL3bhZHmnE2e43SmvNJ6yS9ZGEPf33Vh6k4/Mm3dvPt3W9OJYqt16/n0Fj4uSlJrqVSlbt/EVjDRIfTZ4AXzOx/mtkHJ7+/K7MICyiWK+FOF8vis2Zx5BXn+HiF4X0HeWnvAXa9+gtuenhnsJtZVUfdn7jraX799qf43Xuf5cpfW82ZKxdPXUgdrZD59jCSjZYnx93dzex14HVgHBgA/tbMnnD3/55VgEUUy5Vwp4tl8VmzOPKIs1Jx9rw5yjVfO7blyW0Xn8GfPbGHP7nojMxHujNHOt1dHDfq3rhlJ5suWDsVoy6kiqvVTQ5vMLMdwO3ANuDfufvvM7HV+cUZxldIuh1nfvLY2LCVye1mcWQd594DY1Mn5GqMG7fs5OJ1KzM/Qddbr/Ha/sMsWzg9WQ3vO8TivvmALqSKrtURx1LgE+7+Uu0X3b1iZhekH1ax6Q5inaUIpcdGMS7p78n8BF0vsV7z4A5u2XA6n/3q9qnjVgz0TW3EqAupYmspcdQuyKvzvefTC6dzqAOkcxSh9NgoxpMW9WZ+gm6UtFYv7Z+KacVAH/dcvo6l/T1TN57ShVRxpbUAUKRjFWHxWb0Y77l8He97T1/mJ+hGSWtBb7dG3R0q8a1ji0B3AJS0xbS1SSOhYizbuqUivBfa1eodADXiEGlBEUqPoWIs05xe2ZJkI9pyRApFW7nEKa/b9oYWyzY3oWnEIYWhqz0JrQgddnnQiEMKQ1d7ElqjjRtj6rDLgxKHFIau9iQ0Le6doFKVFEYe6yli6piJKRaZUKZGgNkocUhhZL2eIs85lGZJQfM58SpCh13WtI6jhIp8JZtl7COjR7jorm3HjWjS3g6/laSQVywitVpdx6E5jpKptyFdqK2325Fl22decyitTPJrPkdipsRRMupMaiyvjplWkoK6dyRmShwloyvZxvLqmGklKah7R2KmyfGSKcJOr6Hk1THTyiS/unckZpocLxl168ShyA0KIel1y5Y2OZS6dCUbB7V0JqeLnnhojqOEyrIhnXQWNXbEQyMOkRlUDomTGjviocQhUkPlkHipsSMeKlWJ1FA5JF5qUY6HRhwiNVQOiZcaO+KhxCGFkcfcg8ohcVM3WhyClarM7L1m9oSZvTD570CdY37ZzL5nZs+Z2U4zuyRErBJeXntsqRwi0lywBYBmdjvwlrvfamY3AgPuvnHGMb8EuLu/YGbvA3YAH3L3/bM9thYAdp48d4tVV1Vx6Xc3N0VYALgBOHvy4/uBp4BpicPdf1zz8atm9iawDJg1cUjnyXPuQeWQYlJHXH5CdlWd7O6vAUz+e9JsB5vZWUAP8JMcYpPIaLdYaUYdcfnJNHGY2XfMbFed/zYkfJxTga8Bn3X3SoNjrjazITMbGhkZSSN8iYjmHoqrUnFGRo/wyr6DjIweyezeL+qIy0+mpSp3P7fR98zsDTM71d1fm0wMbzY47kTgW8AX3f2ZWX7WZmAzTMxxzC1yiY1aMYuptny0bGEvN5yzhtVL+1nQ283S/nS3u1FHXH5Clqq2AldOfnwl8MjMA8ysB3gYeMDd/ybH2CRC2mNruryu5OeiWj5atrCXL3z8NDY9souz73yKT9z1dOpdcRqV5idkV9US4CHg/cC/Ap9y97fMbBC41t1/z8wuA/4P8FzN//oZd//RbI+trirpdEWZCH5l30HW3/Yk91y+jlse3Z3L/dzVVdW+6Luq3H0vcE6drw8Bvzf58YPAgzmHJhK9RhPBWbQnz0W1fLS4b34u8w/qiMuH9qoSSSiGElFRJoKr5aODY0fVFddBlDik0PI+iee1gr2ZorQnV5sa/v3K93DPZes0/9AhtFeVFFaIOv/MEtGyhb28/vZh+nu76Zs/L7eaeiv3LY9FV5fx3v5eFvf1qCuuQyhxSGGFqPPXlojOXLmYL3z8NDZu2Zn7BHUR25M1/9A5VKqSVOVZOgpR568tEV179genkkb1Z+e5UlntyRKKEoekJu/6f4g6f+1agbw6hURio8Qhqcl7r6AQC75qS0QrBvoKMUEtzcXQKVckmuOQ1ORdOgpV56+WiCoVL8wEtTRWlMWUMVHikNSE2Cso5IRrESeo5XhFWUwZE5WqJDVl3CtIE9TFV5TFlDHRiENSoytwKSLtqpucRhySKl2BS9GUcaQ8VxpxiEipaaScnBKHiJSeVrUno1KViIgkosQhIiKJqFQl0dDd20SKQYlDoqDVuxOUPKUIVKqSKOS9z1WMYrlJlEgzShwSBa3eVfKU4lDikCgU5VaoWVLylKJQ4pAoaPWukqcUh7l3Xv10cHDQh4aGQochCZV9YlgNAhKame1w98Fmx6mrSqJR9tW72vpCikKJQyQiZU+eUgya4xARkUQ04hCR3JR9HqtTKHGUnP6QJS+a/O8cKlWVmFYqS560wLFzKHGUmP6QJU9a4Ng5lDhKTH/IkictcOwcwRKHmb3XzJ4wsxcm/x2Y5dgTzewVM/uLPGPsdPpDljxpd4DOEWzluJndDrzl7rea2Y3AgLtvbHDsnwPLJo+/vtlja+V4azRZKXlTM0bcirByfANw9uTH9wNPAcclDjNbB5wM/D3Q9AlJ67RSWfKmBY6dIWTiONndXwNw99fM7KSZB5hZF/CnwOXAObM9mJldDVwN8P73vz/9aDuU/pBFJKlME4eZfQc4pc63bmrxIa4DHnP3l81mvwp2983AZpgoVSWJU0REWpdp4nD3cxt9z8zeMLNTJ0cbpwJv1jnsI8Cvm9l1wEKgx8zecfcbMwpZRESaCFmq2gpcCdw6+e8jMw9w99+tfmxmnwEGlTRERMIKuY7jVuBjZvYC8LHJzzGzQTO7N2BcIiIyC93ISUREgNbbcbVyXEREElHiEBGRRJQ4REQkESUOERFJRIlDREQSUeIQEZFElDhERCQRJQ4REUlEiUNERBIJuVeVSNR00yGR+pQ4ROrQ3RFFGlOpSqSOvQfGppIGwPC+Q1z1wBB7D4wFjkwkPCUOkTrGxo9OJY2q4X2HGBs/GigikXgocYjU0TOvmxUDfdO+tmKgj5553YEiEomHEodIHUv6e/jKFYNTyaM6x7GkvydwZCLhaXJcpI6uLuO0kxfx8HXr1VUlMoMSh0gDXV3GskW9ocMQiY5KVSIikogSh4iIJKLEISIiiShxiIhIIkocIiKSiLl76BhSZ2YjwEuh46ixFPh56CAiotdjOr0ex+i1mC7v1+MD7r6s2UEdmThiY2ZD7j4YOo5Y6PWYTq/HMXotpov19VCpSkREElHiEBGRRJQ48rE5dACR0esxnV6PY/RaTBfl66E5DhERSUQjDhERSUSJIydmdoeZ/bOZ7TSzh81sceiYQjKzT5nZc2ZWMbPoukbyYGbnm9keM3vRzG4MHU9IZnafmb1pZrtCxxIDM1tpZk+a2fOTfyf/NXRMtZQ48vMEcLq7nwH8GPijwPGEtgv4BPBPoQMJwcy6gS8DvwWsBT5tZmvDRhXUV4HzQwcRkXHgD939Q8CvAv85pveHEkdO3P3b7j4++ekzwIqQ8YTm7s+7+57QcQR0FvCiu//U3ceAbwAbAscUjLv/E/BW6Dhi4e6vufsPJj8eBZ4HloeN6hgljjA+B/xd6CAkqOXAyzWfDxPRiUHiYWargDOBZ8NGcoxu5JQiM/sOcEqdb93k7o9MHnMTE8PQr+cZWwitvB4lVu9WgmpxlGnMbCGwBfhv7v6L0PFUKXGkyN3Pne37ZnYlcAFwjpegD7rZ61Fyw8DKms9XAK8GikUiZGbzmUgaX3f3b4aOp5ZKVTkxs/OBjcCF7n4wdDwS3HZgjZmtNrMe4FJga+CYJBJmZsBfAs+7+/8KHc9MShz5+QtgEfCEmf3IzO4OHVBIZnaRmQ0DHwG+ZWaPh44pT5ONEtcDjzMx8fmQuz8XNqpwzOyvge8Bp5nZsJl9PnRMga0HLgc+Onm++JGZ/afQQVVp5biIiCSiEYeIiCSixCEiIokocYiISCJKHCIikogSh4iIJKLEIZIzM/t7M9tvZo+GjkWkHUocIvm7g4kefZFCUuIQyYiZ/crk/VdOMLP+yfsqnO7u3wVGQ8cn0i7tVSWSEXffbmZbgT8G+oAH3V03KpLCU+IQydbNTOxLdRi4IXAsIqlQqUokW+8FFjKxT9kJgWMRSYUSh0i2NgObmLj/ym2BYxFJhUpVIhkxsyuAcXf/q8l7jD9tZh8F/gfwb4CFkzsEf97dS7U7sBSbdscVEZFEVKoSEZFElDhERCQRJQ4REUlEiUNERBJR4hARkUSUOEREJBElDhERSUSJQ0REEvn/e/FB4NC5YrcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "np.random.seed(1)\n", "x1 = np.random.normal(0, 1, 100)\n", "y_hat = np.random.normal(0, 1, 100)*0.2\n", "\n", "sns.scatterplot(x=x1, y=y_hat)\n", "plt.xlabel('x1')\n", "plt.ylabel('y_hat');" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.5" } }, "nbformat": 4, "nbformat_minor": 2 }