{ "nbformat": 4, "nbformat_minor": 0, "metadata": { "colab": { "name": "Quick regression example", "provenance": [], "include_colab_link": true }, "kernelspec": { "name": "python3", "display_name": "Python 3" } }, "cells": [ { "cell_type": "markdown", "metadata": { "id": "view-in-github", "colab_type": "text" }, "source": [ "\"Open" ] }, { "cell_type": "markdown", "metadata": { "id": "R5ApiQ53zuc1", "colab_type": "text" }, "source": [ "# Regression\n", "\n", "We'll use the famous [Boston house prices](https://www.cs.toronto.edu/~delve/data/boston/bostonDetail.html) dataset." ] }, { "cell_type": "code", "metadata": { "id": "IuDW61LQ7FHJ", "colab_type": "code", "outputId": "cfb1c567-46c9-42c3-f16f-78bef344db0b", "colab": { "base_uri": "https://localhost:8080/", "height": 204 } }, "source": [ "import pandas as pd\n", "from sklearn.datasets import load_boston\n", "\n", "boston = load_boston()\n", "df = pd.DataFrame(boston['data'], columns=boston['feature_names'])\n", "df.head()" ], "execution_count": 1, "outputs": [ { "output_type": "execute_result", "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
CRIMZNINDUSCHASNOXRMAGEDISRADTAXPTRATIOBLSTAT
00.0063218.02.310.00.5386.57565.24.09001.0296.015.3396.904.98
10.027310.07.070.00.4696.42178.94.96712.0242.017.8396.909.14
20.027290.07.070.00.4697.18561.14.96712.0242.017.8392.834.03
30.032370.02.180.00.4586.99845.86.06223.0222.018.7394.632.94
40.069050.02.180.00.4587.14754.26.06223.0222.018.7396.905.33
\n", "
" ], "text/plain": [ " CRIM ZN INDUS CHAS NOX ... RAD TAX PTRATIO B LSTAT\n", "0 0.00632 18.0 2.31 0.0 0.538 ... 1.0 296.0 15.3 396.90 4.98\n", "1 0.02731 0.0 7.07 0.0 0.469 ... 2.0 242.0 17.8 396.90 9.14\n", "2 0.02729 0.0 7.07 0.0 0.469 ... 2.0 242.0 17.8 392.83 4.03\n", "3 0.03237 0.0 2.18 0.0 0.458 ... 3.0 222.0 18.7 394.63 2.94\n", "4 0.06905 0.0 2.18 0.0 0.458 ... 3.0 222.0 18.7 396.90 5.33\n", "\n", "[5 rows x 13 columns]" ] }, "metadata": { "tags": [] }, "execution_count": 1 } ] }, { "cell_type": "code", "metadata": { "id": "i_9Iamfs8bOH", "colab_type": "code", "colab": {} }, "source": [ "dfX = df[['RM', 'AGE', 'DIS', 'LSTAT']]\n", "X = dfX.values\n", "y = boston['target']" ], "execution_count": 0, "outputs": [] }, { "cell_type": "markdown", "metadata": { "id": "mnNA450Hp829", "colab_type": "text" }, "source": [ "We could plot the data:" ] }, { "cell_type": "code", "metadata": { "id": "qITnC0Nk9kEA", "colab_type": "code", "colab": {} }, "source": [ "" ], "execution_count": 0, "outputs": [] }, { "cell_type": "markdown", "metadata": { "id": "FgqR9Hf4rp2c", "colab_type": "text" }, "source": [ "Let's do machine learning!\n", "\n", "First, we'll just use one feature:" ] }, { "cell_type": "code", "metadata": { "id": "RsUm4g5_sOEf", "colab_type": "code", "colab": {} }, "source": [ "x = X[:, 0].reshape(-1, 1)" ], "execution_count": 0, "outputs": [] }, { "cell_type": "code", "metadata": { "id": "UDZHP00WxTAt", "colab_type": "code", "colab": { "base_uri": "https://localhost:8080/", "height": 286 }, "outputId": "e11bb3e3-8647-4a45-e99c-cc5a2c4df9c9" }, "source": [ "plt.plot(x, y, 'o')" ], "execution_count": 29, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "[]" ] }, "metadata": { "tags": [] }, "execution_count": 29 }, { "output_type": "display_data", "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztnW2QXOWV3/9nWg30aHcZCU8oaBBi\nNykpy8rSmClgVxtXJGet2Lx4ildT7BZJuUI+uBKjdc16nHKMlCJBG8UBf9jyFmUnSwqMBZI9Flat\nIWWUpJYqcM0wmrCziJRfQLjBZryosUGN6Jk5+dB9R90997n3ubfve/9/VSrNdPe999zbc8/z3HP+\nzzmiqiCEEJJ/htI2gBBCSDTQoRNCSEGgQyeEkIJAh04IIQWBDp0QQgoCHTohhBQEOnRCCCkIdOiE\nEFIQ6NAJIaQgrEvyYB/60Id08+bNSR6SEEJyz+zs7C9VddTvc4k69M2bN2NmZibJQxJCSO4Rkdds\nPseQCyGEFAQ6dEIIKQh06IQQUhDo0AkhpCDQoRNCSEGwUrmIyKsAfg1gGcCSqo6LyEYAhwBsBvAq\ngNtV9XQ8ZpJ+mZ6r4eDTr+CNegOXjlQwuWcLJsaqaZvVF1GfU5LXiLaf20+t3kBJBMuqGKmUIQKc\nPtPEkAAr7f47I5Uy9t10FSbGqn0d321bAJFdu+m5GvY/tYDTZ5pr7E4CselY1Hbo46r6y47X/jOA\nt1X1gIhMAdigql/w2s/4+LhStpg803M1fPHbL6HRXF59rVIu4YGbt+XWqUd9TkleI9ruvh8/ykOC\nO665HEdma6GO73bMckkABZor5/xg2Gs3PVfD5OF5NJe7fWp5SHDwtu19fRciMquq436f6yfk8ikA\nj7R/fgTARB/7IjFy8OlX1tw4jeYyDj79SkoW9U/U55TkNaLt7vvxo7miePyF10Mf3+2YzWXtcuZB\n9ue2/15n7tid1L1m69AVwDMiMisi97Rfu1hV32z//HMAF7ttKCL3iMiMiMwsLi72aS4Jwxv1RqDX\n80DU55TkNaLt4W1bNkQUbPYX5Jhh7PPaJql7zdah/6GqfgTAJwB8VkQ+2vmmtuI2rldaVR9W1XFV\nHR8d9V25SmLg0pFKoNfzQNTnlOQ1ou3hbSuJhN5fkGOGsc9rm6TuNSuHrqq19v9vAfgOgGsA/EJE\nLgGA9v9vxWUk6Y/JPVtQKZe6XquUS6sJoTwS9TkleY1ou/t+/CgPCe689vLQx3c7ZrkkKA91DxJh\nr93kni2tmLyL3Unda74qFxFZD2BIVX/d/vnjAP4DgKMA7gZwoP3/d+M0lITHScYUSeUS9TkleY1o\ne/d+gqpcxq/YGOr4JtsBYN/RBdQbLWXKBeVwqUVn/5lWuYjIb6M1KwdaA8A3VfU/ishFAJ4AsAnA\na2jJFt/22hdVLoQQN9KU1eZBBWarcvGdoavqTwBsd3n97wF8LJx5hBDSoteh1uoNfPHbLwFAJA7V\nb7DwUu5kxaHbwpWihJBUiVN26QwWtXoDinODxfRcbfUzRVKB0aETQlIlTodqM1gUSQVGh04ISZU4\nHarNYFEkFRgdOiEkVeJ0qDaDxcRYFQ/cvA3VkQoEQHWkkqmEaBASbUFHCCG9xCm7nNyzxVXB0jtY\nTIxVc+nAe6FDJ4SkTlwOtYhrMLygQyeExErapZuLMvu2gQ6dEBIbcWvMSTdMihJCYqOIpZuzDGfo\nhJDI6exG5EYeF+3kATp0Qkik2HQjyuOinTxAh04IiRS/bkSdssG0E6ZFgw6dEBIpXuGUaofTZsI0\nepgUJYREiimcUh2p4Lmp3V3acCZMo4UOnRASKbZL+ZOscjg9V8POA8/iyqlj2Hng2a5qi0WCIRdC\nSKTYrs68dKTiqoKJOmE6SKEdOnRCSOTYrM60rbPSL0VqYOEHHTohJBWSqrNSpAYWftChE0KsyKvE\nMKnQThZgUpQQ4otNK7cs7NONIjWw8IMOnRDiiykOve/oQuT7jFq2WKQGFn4w5EII8cUUb643mpie\nq4VyjknGtgelhC5n6IQQX7zizfceOhFK212k5sxZgQ6dEOKLX7w5TPw7ydj2oCwsokMnhPgyMVbF\nhuGy52eCxr+Tim0nlXzNAqKqiR1sfHxcZ2ZmEjseISQ6bMriAi3HnCVp484Dz7rKFp3aMnlARGZV\nddzvc0yKEkKs6FwIZGpcIcDqe1lZYj9IC4sYciGEWDMxVsVzU7vx0B071sS/BUDv834WqicOUvKV\nDp0QEhi3+LcpeJv2TNiUfN21dbRwiVKGXAjJGVlZgt+r7TbFqtOeCbvVjNm1dRRHZmuFq8BIh05I\njshyKdikqieGwW3wKWIFRoZcCMkRWe7yk6cl9kVNlHKGTkiOyLojyssS+6JWYOQMnZAcMUiKjTgp\nagVGOnRCckRRHVHS5Ck8FATrkIuIlADMAKip6g0iciWAbwG4CMAsgD9R1Q/iMZMQAiTX5WcQyEt4\nKAhBYuifA/AygN9q//7nAB5U1W+JyF8C+AyAr0VsHyGkh7w7oqzILouIVchFRC4DcD2Ar7d/FwC7\nARxuf+QRABNxGEgIKQ6DVCgrDWxj6A8B+DMAK+3fLwJQV9Wl9u8/A+A6xIrIPSIyIyIzi4uLfRlL\nCMk3WZZdFgFfhy4iNwB4S1VnwxxAVR9W1XFVHR8dHQ2zC0JIQci67DLv2MTQdwK4SUQ+CeACtGLo\nXwUwIiLr2rP0ywDwmYkQ4klR9d9ZwXeGrqpfVNXLVHUzgE8DeFZV7wJwHMCt7Y/dDeC7sVlJCCkE\nUcouB6ULURD60aF/AcCfisiP0IqpfyMakwghRSUq/TeTq+6wYxEhJHcUoQtRENixiBCyStG030yu\nusOl/4QUnCKGJ1jTxh06dEIKTtG039NzNZz5YGnN66xpw5ALIbnFNoxiauicx/BEb4MPh5FKGftu\nuirXYaQooEMnJId8afolPPb8qdU+nqbORdNzNdfmzUA+wxNuTxsAsP78dQPvzAGGXAjJHdNztS5n\n7uAWRjn49CuuzlyAXIYnmAz1hg6dkJxhctLAWsdmcnSK9HuQhoHJUG/o0AnJGV6z0V7H5uXoNudw\nhSUbfHhDh04GlrwuHTc5abcwyuSeLSiXxLivWr2BvYdO4EvTL0VpYmwUtdNQVDApSgaSXrWEKamY\nRSb3bFmj9BAAd123yd12n8XgCuCx509h/IqNALLfDSnvDT7ihA6dDCRe2uysO4sgbegOPv0Kmiv+\n5T0UwP6nFvB+cyWXgxxpwZALGUjyrpaYGKtics8WXDpSwRv1Bg4+/YpryCjI+Zw+0yzUAqRBhDN0\nMpDkvS63bcjIdJ5B6HeQK1odmSzDGToZSLKulvBL2Nou53c7z6D0M8gVsY5MlqFDJwNJltUSNk7Q\nNGuu1RtdA0DneYah30GuaHVksg5DLmRgyapawiZhOzJcxukzTdfte8MvE2NVzLz2Nh59/pTvsaW9\n7/qZZiThkbznKvIGHTohGcPGCfr1pekdAB5/4XWrYyuA95srePCOHZEMdnnPVeQNhlwIyRg2y9vf\nabjPzjvpHACWA3QmizIkkvVcRdGgQyckY9g4QZsZbudnSmJeLepGVCGRLOcqighDLoRkDJuFQ26r\nRTvpHQDuvPZyqxi6Q5QhkazmKooIHTohGcTPCTrv7X9qYTU56tQ9r7oMAPdPbMN3XqzhvQ/WDgC9\n9dJ7BwPqyPMDHTohGcTWib77/rlWbAqgPCTGz55xcebOdg4C4Jarzw0m03M1TB6eR3O59alavYHJ\nw/MAWA4gi9ChE5IxbFeB7ju6sKZOS3NFse/oAibGqpieq2Hf0QXU2wnUIfFXxyiAQz98HeNXbMTE\nWBX7n1pYdearx1hW7H9qgQ49gzApSkjGsF2MUzcoXeqNZmtm/eR812csanQBaA0KzrFMWnfT6yRd\nOEMnJCbCxp6jWIzjVWWxJOIrY+TCn3xCh05IDAStt97p/IcMDrdXebLBsFp0w3DZsyDXiiqqPkW7\nnGONVMquTwIjlbJxW5IeDLkQEgNBapj01m4xzZ7fO7vUVc/lvhuvWtONqFwSXP/hS+ClOndK7ppw\nEqsAsO+mq1AekjXv77vpKo8jkLSgQyckBoKETdycvxv1RrOrSNfEWBUHb93etWjn4K3bcfzkorFJ\nUbkkq3XU3RABDt62ffUpYmKsioO39Ryj432SLRhyISQGgtQwCRKv7q3R4qZX33vohHH75nIr4blr\n6yiOzNa6BpJKucRVnDmHDp2QGNi1dRSPPX/Kc8GOQ9AmFH4DgN/+avUGjszWcMvVVRw/ubgmadsZ\nzx8ZLuPd95dWE6xsS5dt6NAJiZjpuRqOzNY8F+x04reMv5cLfRKSNvtrNJdx/OQinpvavcb2zm3d\nkq556b06iDCGTkjEuMXEFcDxk4uun3cKWG0YtlOOvPfBkmfHn96CWCbcZvr7ji5YDSyUNWYTOnRC\nIsa2m1AnE2NVzH3541ZO3YmDO7i1q5sYq+K5qd346YHrjd2KeuP503M142Ilv21JNqBDJyRivJyd\nWzu5TodsuwLTGTRs2tXZ1iS3rYHOeubZxdehi8gFIvJDEZkXkQUR2d9+/UoReUFEfiQih0TkvPjN\nJST7+DVm7tSj9zpkW5xBw0vv7gwUew+dwPnrhrBhuOxZk9wrjDJS8d6WZAObpOhZALtV9V0RKQP4\nGxH5awB/CuBBVf2WiPwlgM8A+FqMthKSCzrrmZvUJs7rthr0TjpnyF7hnc7kZr3RRKVc8mwtZ1LH\nbBguY+7LHw9kI0kH3xm6tni3/Wu5/U8B7AZwuP36IwAmYrGQkBwyMVbF5J4txk5BgtbsPIhcEQDW\nn3dOKz49V8OQYf8lEeuVqg6m0Mx9N3JVaF6wki2KSAnALIB/COAvAPwYQF1VnWLMPwPAZzCSK+Js\n3OCEUkzL+BWt2blNoaxOnJrmX5p+aY3O3aFSLhln/V5hFZtOSSTbWDl0VV0GsENERgB8B8BW2wOI\nyD0A7gGATZs2hbGRkMgJWjwrKPuf8pf/vREwbg60BoJ9RxfwTqPpum1JBA/cvM0Y7vFTp7BdXL4J\ntLBIVesichzA7wMYEZF17Vn6ZQBchbGq+jCAhwFgfHw86N8vIbHglUwM49B6V1faqFUc5+rmeL1m\n7l7SwmVV7D10AiPDZZSHpKuELtUpxcdG5TLanplDRCoA/gjAywCOA7i1/bG7AXw3LiMJiZooao47\n9CpVbKWHu7aOYtfWUdf3rvvtDZ6LgrxYtUGoThk0bGbolwB4pB1HHwLwhKp+T0T+DsC3ROR+AHMA\nvhGjnYRESpDiWSacWXnQxKbDo8+fMr736t83cNd1m1zrwVxQHrIaNJrLil+/v+SpbCHFwtehq+r/\nBTDm8vpPAFwTh1GExI1bvZMgIYneGHzUvFFv4P6JbQCAx194HcuqKInglqurGL9io/Wxl1VZTGuA\n4EpRMpD01juxDUk4i3XuPXQiNmcOtJ4UpudqOPTD11dj6cuqOPTD1wFgje1eJQP85IqkOIgGkEz1\ny/j4uM7MzCR2PJJt4pQNxkGQWXm5JFh/3jq802jFsoPcZuWS4OCt27Hv6IKx/duJ+7oX+vjZJgB+\neuB6eyNIphCRWVUd9/scy+eSVIhbNhgHtqs6qz2D05VTx4IdqO38TWoWt9edY33+iXmrfqSkmGTe\noedtFkfsiFo2mAR+ChhTx5+gDSyaK+obInEqKnbi/N5PboDkm0zH0G0qyZF8EqVsMCm8ZrleMfgw\nzrRWb2C4bL49TfdB2NwAKQaZnqHncRZH7IhCNhgFbk+AgPvyd5Myxs1h9u63Uh5Co7libVdJBOeX\nSzhj2MbrPuBqz8El0w49j7M4Yke/ssEocIvjTz45D0hLw+281hvb9wsBuu23XBLXlZummPyyKuo+\nWnPeB6SXTDv0rMziSPRkoRCU2xNgp8N16JwN28x+Xfe7rFh/XgkrzZUuTfnxk4vG+Lr4qGMUwOap\nYxiplLHvpqs4KyfZjqHbdlohJAxBZrg2n3U06iYH/d4Hy12a8iOzNezaOmpshuEytrhSbzQx+eQ8\nc0sk2w6dCZ7ikoWEd5AnPb/Pdp6PLY3mMr43/yYu8Eh+AlitqW6qrQ7YKWNI8cl0yAVggqeoZCHh\n7RbHLw9JVwwdsHsqDNN5CPCunOiwoopX24uCrpw6Ziy5y5g6ybxDJ8UkCwlvUxzf7TW/QcbL7upI\nBe+dXbJy3m50Ph14adqZWyJ06CQVspLwNj0BBn1KMJ1PdaSC56Z2hy7mVS5J19PB5J4tmDw83/UE\nAbSeLJhbInToJBWyIFsE1urFd20dxfGTi4GVN37n4/Y0cOaDJc8yuBuGy7jvxm71ivPz/qcWVrel\nyoU40KGTVEhLttjbWejd95dWpYq1eqOrRrmpvsz0XM3VoT5w87au189f553svP7Dl+DIbK1rEBC0\n5Ii99WA68corsVTGYMNqi6RQeDm0sGGPkghWVFdn8I//8HUsu2gKh8tDaC5rl5bdcdAbegYPoDWD\n79SiO5/tfP+Bm1s10W2ctNv5mVayknxhW22RDp0Ugt5Zs0OnQ/PSiKeFE2M32TZSKePs0oqVkzbt\nwzkGyS8sn0sGBq+Zd6O5jH1HFzDz2tuZc+bAuQbRJpWMmzLGJO/MgnKIpEumFxYRYoOfBrzeaHr2\n70wTZ7FQUHWPm5M27YNyxsGBM3SSW/pt0pwFllVx5dQxjAyXXYt3mRpCOy3qehU6vUlWlsoYLOjQ\nSaQkpbKIu0lzkiiwxmk7yhnAvWHFrq2jayo6Pvr8KVTKQ9gwXEb9TJMqlwGEDp1ERpJt5cIutc8L\nZ5daddBN8k7T+bdqrgsevGMHHfkAMjAOnfrc+LGtz2L6LoJ8R0VP9PmV7N176ITVtmSwGAiHnseG\nxHnERmVh+i5mXnu7K/7r9x0F7dOZR7wGLb/zL/qAR9wZCJWL18yRRIeNysL0XTz+wuuBvqNdW0f7\ntDb7eKlT3HoF2G5ListAzNCpz00Gm/osplnlsmGBm/MdTc/VsO/oQuiKhXnDT53iVtPFdltSXAZi\nhk59bjLYNCTxatLghiPPm3xyfmCceUnEarn+xFgVc1/+OB66YwebwBAAA7L0nzUussPmqWPG93qb\nJjvfUd615mGojlT6ShqTYsGl/x1koSFx2mTFGVQ96oY7crxavYGSCBrNZdeQQtERnAtNhU0ak2yQ\n9H03EDP0QSdLTyh+tmRhwVBJBMuqq/+H5aE7duDzT8wb91EuCaBwrc5osqkXFt7KLlHed7Yz9IGI\noQ86WVL5+MXZs7BgyHGc/ThzoHWuXvs4eOt2HLxte9e1MH3aL2lMskca991AhFwGnaypfLwaNBTF\nQW0YLgPwDjE516DzWphK4Jpm6EzsZ5c07jvO0AeAPKl8bGzaMFzGUDCxTKKUS4L7bmzVYXHTizu1\nWHYeeBZXTh3DzgPPYnqu5vn5O6+93PV1yhOzSxr3HWfoA0Dc/TujSPx0Vk40xZGBlt1nm8twaRiU\nCUoiOHjr9jWzb6+qiG4JTrfrOX7FxkwktokdafTNZVJ0QIgr2x5F4sdtH52t21SBdxrnqgfe61HH\nJE1sz5udhQaHqO67yGSLInI5gP8B4GK07rGHVfWrIrIRwCEAmwG8CuB2VT0d2FKSCF5x637wS/zY\n/DG77cNplOw4OOfG8CpKlRQiwKUXVlbllcuqq7JLoOWwvc45azkNEh9x3XcmbEIuSwA+r6ovishv\nApgVkf8J4F8A+IGqHhCRKQBTAL4Qn6kki5gW/DhhBBvdtMmR1eoNbJ46hiFBpkIsqnCdSdsWgTMV\n1spiToPkC9+kqKq+qaovtn/+NYCXAVQBfArAI+2PPQJgIi4jSTaZnqvBKzdpK9nyc2RZcuZA68mh\nl+m5Gj7/xLzVOZsSn0xwkn4JpHIRkc0AxgC8AOBiVX2z/dbP0QrJkAHi4NOvGJOXJtxm436VA7PG\n5ou6HbozM7fVitvUvCEkDNYqFxH5DQBHANyrqr+SjiJLqqoi4vrXLCL3ALgHADZt2tSftSRThIn5\nds7GOxNGF1bKxv6ZWeO5H7+NL02/hPsntgHwXwzl9gSSdGyVDAZWM3QRKaPlzB9T1W+3X/6FiFzS\nfv8SAG+5bauqD6vquKqOj44Wv4b1IBE05tsZVnBmtbV6Awqg3mji3bNLuVkY8fgLr6/+7FU4jKEU\nkiS+94+0puLfAPCyqv7XjreOAri7/fPdAL4bvXkkywQJldgs8W8uK1YitzIenPCKVx7BtgwuIVFh\nE3LZCeBPALwkIo5m7N8BOADgCRH5DIDXANwej4kkq3Q6Z7/ytr2qkLxL9IbErCcHWjr6r9y+nc6c\nJIqvQ1fVvwGMk5CPRWsOyRudseDf/fd/jTPNtXNsp65JJ7nvCareoZaMCXPIgJCXkCXJAf/p5g+3\nSsJ20FnXpJMs9gQteRSI6X3HJjT0xW+/tFqjhZAkYC2XjJKVhhRBCNJI5PjJxaTN82S4POT6dNEP\njgY9698bKQ506BnEdsVhEnYEHVRs5XhZi6H7OfOwIZSsnWdQ8jixGGQYcskgWWhI0SsrdAaVKEII\n03M1DAVsFp1X8rycP86/ARIPdOgZJAvFm+IaVPxWVWaRsEOPSYM+PVdzrYWeNbIwsSDBoEPPIFlo\nSBHXoBK2xVx5SFbVMv00twiz6V3XbQpcmsC0nD9Ps94sTCxIMOjQM0gWijfFNaiEcQYC4I5rLsd9\nN16F6kilr2JdCmCkslZGaaI6UsH9E9vwwM3bULIMEznflW2p4KzOerMwsSDBoEPPIFko3hTXoBLG\nGSiA782/icnD84lq1zvPd2Ksiq/cvt1qpt5oLmPvEyewY/8za8IqeZr1ZmFiQYJBlUtGSbN4k6Ns\naDSX1zRw6Ncmt7ZcNtQb0RXt8tpXdaRiVHQ4P9t0TFI9d5xOlZJpQdWQCKbnaplSkASRoZJsQIc+\nINjKz3olk8uqniGEoAQpFxCGztZ1QSs3urWAc7tuziAXBCesYhrQllVTkab6waqQ+YIhlwEgSCIu\niRjvxFgVz03txqsHro9kf05suzpSwV3XbUJ1pIK6jzO3CSWYrltYhc4b9cZqOM0tHp/VWDrJD3To\nA0AQJx1ljNdGnhckQenGSKWMHz/wSTx0xw68d3YJjz5/atUBe23zwM3buo59QXntrWC6brbJ0V6c\n/MHEWBUrls0wCAkCHfoAEMRJR6VssHkq+NL0S33FxstDgn03XbV6LJt9OdsAwNmlc6tDT59prrHP\ndN3CzNB7nwCoICFxQIc+AJichAJrZs5RKRv2HV3wfCqYnqvhsedPBdpnJ9WRCg7e1ipPa6NtF4tt\nep9aTNetOlJxrSDpZaujUnKeWmr1xhpNPBUkpF+YFB0AvJQlvXViolA2TM/VjLNlZ9Ybph+pg6C7\nvrpfmMIt2WlKyHa+PrlnCyafnEezR/j+Rvupw0nA2h67N+HcuY+SSNeAkkYiknVb8g8d+gDgpyzp\nrQrYr7LBK7HnzHr7iRVf2BN396qtXimXsGvrKHYeeNZKqbImPu4SLteO//2ceud5uj0VOPtwbEmz\nEFsWCsKR/mDIpeA4j/h7fbTTUSbjvPblhBT6iRXXG82uUJGpFd6G4TJuubqKI7O1rlj+3kMnjHHw\nZdXV/R58+hU0l73n4IrWLLxqERM3XZfeI6ShdsnTClZihjP0AuM26zLNKKNIxjmP7CYXuGG4vDrb\nC7vAyMFtBukWLth54FnXWbEXzn5tBznn6aD32rolQm2190mrXfK0gpWYoUMvMKZHfDf67SDUO3j0\nUimXujoXOU54/1MLaxYA2cSmge5QkVuYaHqu5ulATcdx9hu0TZ5fTNxtEPMaYJOMaZvOlaqbfMGQ\niyV5KXnq4OfMeum3g5CX0mTDcNlYi+ZXjaU1rzmO0QbTDNIZYLzwGjRq9YYxlOO3T7eYuLOs/5ar\nq6tx+pII/uB3NrqqinZtHU20KiPrthQDOnQL8lTyFLBzZr30+2jttf37Pd2Apudq2LH/GdzrEct2\nYtOAS6KyA7cZ5PRcDZ9/Yj50OKcTp0haEEwx8em5Go7M1lbPeVkVL556B7dcXV1TiO34ycVEY9pZ\nKAhH+kc0wUYD4+PjOjMzk9jxosLRDffiJofLAiZ7AfMjfr/n4nXMzv37hWZM9rht55yLUzgMaOnf\noyrk5SU57Dy+DYKWOsfNNsd+J7wy4lGHRgD8NKKSCSQ/iMisqo77fY4xdAvyljDysuuu6zbhyGyt\nyzEFebQ2xXV3bR3Fox4LhTr1537OvFwSvHd2CVdOHVsTO3akl53OtFZvYPLJeawAWO6nWLrBZsA9\n5n9BeQhLK+qrhAHg6aSdJz7nungVFWNMm3hBh25B3hJGJnudZg3jV2wMlWzz0ir7xeBt9efS9tRu\npWedxKfb00Dv4p8ocPt+O8NHDZ/G0g6VcgleD8JOAtVmP4xpEy8YQ7cgbwkjP3udaoc/PXA9npva\nbR0n9dIqeznqzmP7DYKCtc650VzG/qcWVn8P82QUtKCW2/cbpn2eE4t+xyMMZFsbhjFt4gcdugV5\nSxjFZa9X6MnkqEsiXcf2U46YJtqnzzRXk9BhnozuvPZya6cu0gqn7D10wqrbkHE/wOqAabJ5w3DZ\nKulaHalk9u+NZAeGXCzJW6H/OOz1Cj25aawr5VJXUarOpN/564bwTqOJoQDNIhzNeZhFScdPLuLO\nay/3jPMDrfg99Fwc26bbkIlOJ266Po4230/Dn9WnQZItOEMn1phm12c+aGnJTU8FvbLP02eaOLu0\nggfv2GGsC+6GM0N2nkCC8Ea9gfsntuGPr9vUpQPf+Tsbu2xef94615CP023IVpfe64S9npp639sw\nXMZIpZyLp0GSLShbtCRvlejisnd6ruYqDeycjffiJWkM0s6tV8roJ5V0296r9Z5XWzxHLjg9V/Pt\nKTpSKWPfTVdl+u+D5Atb2SJn6BbkdWFRHPZOjFWx/vy1kToncem2mtYr9uzmzMslQXmoO97tFnZw\nmzGXh6QVNnHBdB06r5eJzm5DfrH4zsYZhCQJHboFeatEF7e9Jgd9+kzTdRCxSWKWRM41obh1Ow7e\ntt03qesWxjh423YcvHW7MdFU8cpFAAAJPklEQVTodh381Cu9g4nfE0WW/zZIsWFS1IKiLCyKyl7b\n5GBn7Hny8LznApwV1TUrIG1CFqbk78RYFVdOHXNdydl7Hbyui1uYpmpx/ln92yDFhjN0C/LW/zFu\ne4MkB51O9+vP8547hOlZ6lcszfY6eLWac9Pp25z/yHA5V8XcSDGgQ7egaAuL+sUt1DFSce+x6ThL\nr4U1QW2zzRHYXofNF7k7dNPrnecPrK0MWS4J3n1/KTc5F1IcGHKxIIo+m2EJo1ZJw94btl/iWSPG\nFKbpXXjkh1NJsTeO3dtGD7C/Ds//5LTrsUyvO/t29tP7Hb13dmmNCsjNvn7Jm/KKxI+vbFFE/huA\nGwC8paq/135tI4BDADYDeBXA7apq/utvk2fZYhq4Vfjzkgd2bhfnjW6y65arqzh+ctH1uGHPxe+4\nnYStRLh56pjxvVdD7M8Uu4+yUmIU15Pkhyhli38F4J/3vDYF4Aeq+o8A/KD9O4mYMGqVJCSWJruO\nn1w01oiJohyBnxolbI7AS4YY5rolkXPJm/KKJINvyEVV/4+IbO55+VMA/mn750cA/C8AX4jQLoJw\nahWvG73TefYziw+roum3HIHX/gUInSPwKgkQJkxiWuYfZc4lb8orkgxhY+gXq+qb7Z9/DuDiiOwh\nHYQp22tzo3uVwbVxXmmVE/aSSyrsbHfj/oltRofeez3dBkJgbZz+gZu3xRr2yltJZ5IMfatctBWE\nNwbiReQeEZkRkZnFxf76Vg4aYdQqNo/7/T6up6X6mdyzxdhrNGibONvtO6+bWzhr8sl5TB6eXxPi\nAhCqRLEteVNekWQI69B/ISKXAED7/7dMH1TVh1V1XFXHR0f76yw/aISJO9vc6P0+rqdVTnhirIq7\nrtu0xqlH4chsrpvbQNh06ViURCw7byWdSTKEDbkcBXA3gAPt/78bmUUFIgq1SdC4s41UL4rH9bTK\nCft1XAp7zW2uW5D4dBKx7LyVdCbx4+vQReRxtBKgHxKRnwG4Dy1H/oSIfAbAawBuj9PIPNJvnLof\n/G70JJJ2cWI6v6DX3M35ezXKDlIPnbFskgY2Kpc7DW99LGJbCoWt2iQNolx4FGRGHGb27LeNX9lb\n0zUPM+C6DYTlIQEEXWGXKAZHLhoiYeBK0ZjIuqzMbaXj3kMnImsaHYUD9dvGb6GRg9s1DzPgmgZC\nt9f6cb5pPt2RfEOHHhN5kZX14zyCOMUwDtRvG9umzW7XPGotfZSONstPdyTbsDhXTORFVtaPhDGI\nUwzjQP22sXnaMV3zLFfQzPrTHckudOgxkVVZWW/ZWVPs2cZ5BHGKYRyo3zZ+zter8FdWB9zpuRqG\nDKUIsjDYkGzDkEuMZE1W5hZeEbivCjM5j85k3YWVMsolsUoIhlHW+G3j9n7n57wG0KCJYVOSMsrk\npfP9uHVEysJgQ7IPm0SnTJJqBtOMvNepm5yhWxKyPCT4jQvWoX6miUtHKti1ddSz4mKcKhen4bRX\nM+gweFWXdCsZHPZJzPT9lETwldu3Z2pyQJLFttoiHXqKJF0C1VTWFWiFhPwcrcnhOJ19ilrS1cvR\nus2mnesRlCTK7pJ8YuvQGXJJkaTVDCblja0D8kvWFVWdYTpvU7PosMnLvCijSHZhUjRFklYz9JsI\n9EtSJnk+Nj1Fo8J03qY66mEdcFYTtSQ/0KGnSNLSuX6VN34OJ6nzSaKJRyem877z2ssjdcBZVUaR\n/MCQS4qkUVOlH+WNnzIkqfNJOrTjdd5ehcLCHosOnISFSdGUKVrNjiTOxyZ5WLTrSgYbJkVzQtFm\nZEmcj1/ykLVQyKDCGDrJHX6xfDZQJoMKZ+gkFsKGPGy284vlsxYKGVTo0EnkhA15BNnOK7RDPTcZ\nVBhyIZETNuQRVaiEem4yqHCGTiInbMgjqlBJlB2ZCMkTdOgkcsKGPKIMlRRNPUSIDQy5kMgJG/Jg\nqISQ/uAMnURO2JAHQyWE9AdXihJCSMaxXSnKkAshhBQEOnRCCCkIdOiEEFIQ6NAJIaQg0KETQkhB\nSFTlIiKLAF5L7IDh+BCAX6ZtRALwPIvHoJzrIJ7nFao66rdBog49D4jIjI08KO/wPIvHoJwrz9MM\nQy6EEFIQ6NAJIaQg0KGv5eG0DUgInmfxGJRz5XkaYAydEEIKAmfohBBSEOjQOxCRkojMicj30rYl\nTkTkVRF5SUROiEhhq6WJyIiIHBaRkyLysoj8fto2RY2IbGl/j86/X4nIvWnbFQcisldEFkTkb0Xk\ncRG5IG2b4kJEPtc+z4Ug3yfL53bzOQAvA/ittA1JgF2qWnQt71cBfF9VbxWR8wAMp21Q1KjqKwB2\nAK0JCYAagO+kalQMiEgVwL8F8Luq2hCRJwB8GsBfpWpYDIjI7wH4VwCuAfABgO+LyPdU9Ud+23KG\n3kZELgNwPYCvp20L6R8RuRDARwF8AwBU9QNVradrVex8DMCPVTXri/fCsg5ARUTWoTU4v5GyPXHx\njwG8oKpnVHUJwP8GcLPNhnTo53gIwJ8BWEnbkARQAM+IyKyI3JO2MTFxJYBFAP+9HUb7uoisT9uo\nmPk0gMfTNiIOVLUG4L8AOAXgTQDvqOoz6VoVG38L4J+IyEUiMgzgkwAut9mQDh2AiNwA4C1VnU3b\nloT4Q1X9CIBPAPisiHw0bYNiYB2AjwD4mqqOAXgPwFS6JsVHO6R0E4An07YlDkRkA4BPoTVQXwpg\nvYj8cbpWxYOqvgzgzwE8A+D7AE4AWLbZlg69xU4AN4nIqwC+BWC3iDyarknx0Z7tQFXfQiveek26\nFsXCzwD8TFVfaP9+GC0HX1Q+AeBFVf1F2obExD8D8FNVXVTVJoBvA/iDlG2KDVX9hqperaofBXAa\nwP+z2Y4OHYCqflFVL1PVzWg9tj6rqoUc/UVkvYj8pvMzgI+j9YhXKFT15wBeFxGnw/THAPxdiibF\nzZ0oaLilzSkA14nIsIgIWt/nyynbFBsi8g/a/29CK37+TZvtqHIZPC4G8J3WPYF1AL6pqt9P16TY\n+DcAHmuHI34C4F+mbE8stAfmPwLwr9O2JS5U9QUROQzgRQBLAOZQ7BWjR0TkIgBNAJ+1TehzpSgh\nhBQEhlwIIaQg0KETQkhBoEMnhJCCQIdOCCEFgQ6dEEIKAh06IYQUBDp0QggpCHTohBBSEP4/epDC\njUXKb2QAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "tags": [] } } ] }, { "cell_type": "code", "metadata": { "id": "t_MvhboDp-3D", "colab_type": "code", "colab": {} }, "source": [ "import numpy as np\n", "from sklearn.linear_model import LinearRegression\n", "\n", "reg = LinearRegression()\n", "\n", "reg.fit(x, y)\n", "\n", "y_hat = reg.predict(x)" ], "execution_count": 0, "outputs": [] }, { "cell_type": "code", "metadata": { "id": "tMtjjoZ4xbov", "colab_type": "code", "colab": { "base_uri": "https://localhost:8080/", "height": 34 }, "outputId": "b5fd3351-c0b0-46e9-ba4e-c9e6eef2bf78" }, "source": [ "reg.coef_, reg.intercept_" ], "execution_count": 38, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "(array([9.10210898]), -34.67062077643857)" ] }, "metadata": { "tags": [] }, "execution_count": 38 } ] }, { "cell_type": "code", "metadata": { "id": "UCQH2zqqxbwW", "colab_type": "code", "colab": { "base_uri": "https://localhost:8080/", "height": 286 }, "outputId": "158d9cf0-cb31-405f-8e7f-8bf014253364" }, "source": [ "model_x = np.arange(3, 10, 0.1).reshape(-1, 1)\n", "model_y = reg.predict(model_x)\n", "\n", "plt.plot(x, y, 'o')\n", "plt.plot(model_x, model_y)" ], "execution_count": 42, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "[]" ] }, "metadata": { "tags": [] }, "execution_count": 42 }, { "output_type": "display_data", "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XmYU9X5wPHvm0wGBlCGTYQRxKoF\nV0BHRFErKOKKFHelaqXir9VWUVGwLmhRUNqqrV1EXOsGig6bCqKorYoKDqiIuIMEFAQGFcaZLOf3\nR5IxM5PlJrnJzSTv53l8ZiZzc+9JHN6c+55z3iPGGJRSShU+l9MNUEoplRsa8JVSqkhowFdKqSKh\nAV8ppYqEBnyllCoSGvCVUqpIaMBXSqkioQFfKaWKhAZ8pZQqEiVONyBa586dTa9evZxuhlJKtSjL\nli371hjTJdlxeRXwe/XqxdKlS51uhlJKtSgissbKcZrSUUqpIqEBXymlioQGfKWUKhIa8JVSqkho\nwFdKqSKhAV8ppYqEBnyllCoSGvCVUspJAR/896/gXZb1S+XVwiullCoqG1bA7Mvg6/egbixUHJzV\ny2nAV0qpXPP9CK/eDq/fDW06wZn/gX2HZ/2yGvCVUiqX1i4J9eo3fwL9RsGwSVDWISeX1oCvlFK5\nUPc9vHQLvH0ftO8Bo56BvY7JaRM04CulVLZ9ugjmXgHb1sGhl8CQG6BVu5w3QwO+Ukply44tsOCP\nsOJx6PxzuGgB9DzUseZowFdKqWxYWQXPXQ21W+HIq+GoceBp7WiTNOArpZSdvv86FOhXzYVufUO5\n+m4HOt0qQAO+UkrZwxhY/jgsmBCadnnsRDjs9+DOnzCbPy1RSqmWausamHs5fL4Yeh4Gw++Bzns5\n3apmNOArpVS6goHQNMuXbgEROPHPUDkaXPlZtUYDvlJKpWPTapjze/jqLdhrKJx8J5T3cLpVCWnA\nV0qpVAR88L+74LU7oLQt/PJeOPCsUA8/z2nAV0opq9ZXh8oifPMB7PdLOGEqtOvidKsssyXgi8iX\nwPdAAPAbYypFpCMwA+gFfAmcaYzZasf1VG5UVXuZumA162tq6V5exrhhvRnRv8LpZjXIpH3ZeG35\n1J50zxd5nremFgFM+PFSt+ALGowBtwgDf9aBLzfXWj5/0/YM7tOFxR9tSqt9N89dydYdPgDKyzxM\nHL5f9v8ufbXwymR44+/Qdhc46zHY5+TsXjMLxBiT/KhkJwkF/EpjzLdRj90BbDHGTBGR8UAHY8y1\nic5TWVlpli5dmnF7VOaqqr1MeOZ9an2BhsfKPG4mjzwgL4J+Ju3LxmvLp/ake75Yz7Mq0fmtnNdq\n+8Y9vQJfoHHM8riEqWf0zd7f5Zevh3L1Wz6Dg86HoX+CsvLsXCtNIrLMGFOZ7LhsDiWfCjwc/v5h\nYEQWr6VsNnXB6mb/QGt9AaYuWO1QixrLpH3ZeG351J50zxfreVYlOr+V81ptX9NgD+ALmuz8Xf74\nHcy7Eh46EYJ+OH82DP973gX7VNiVwzfAQhExwL3GmGlAV2PMhvDvvwa6xnqiiIwBxgD07NnTpuao\nTK2vqU3p8VzLpH3ZeG351J50z5fp/9tMX0cm7bP97/LjhTBvLHznhYGXwpA/hgZoWzi7evhHGGMO\nAk4ALhWRo6J/aUJ5o5i5I2PMNGNMpTGmskuXljP4Uei6l5el9HiuZdK+bLy2fGpPuufL9P9tpq+j\nfZknrfOnco2ktm+GZ8bA42eEqlmOfhGOv60ggj3YFPCNMd7w143As8AA4BsR6QYQ/rrRjmup3Bg3\nrDdlHnejx8o8bsYN6+1QixrLpH3ZeG351J50zxfreVYlOr/V826v91NV7U3YPo+7+dRHj0sy/7s0\nBj54Bv4xAD6YBb+4Fi55DXocktl580zGKR0RaQu4jDHfh78/DrgFmANcAEwJf52d6bVU7kQGwPJ1\nlk4m7cvGa8un9qR7vsjvJ85ZSU2tr9nvI7N2UpmlE5mdU+sL4BYhYAwV5WVs3V7HDl+w0bG+QCgX\nH6+dkcdtn6Xz3QaYfxWsng/d+8Pw2bDr/umfL49lPEtHRH5GqFcPoQ+Qx40xt4pIJ2Am0BNYQ2ha\n5pZE59JZOko5b9CUl/HGyIlXlJfx+vghls+TaLbQ2BnLY+Z4BfhiyklptDoNxkD1f2DB9RCog8F/\nhIG/y6tiZ1ZZnaWT8SszxnwO9I3x+GYgt/t3KaUyZtcgcqLZQt3Ly2J+qORsjGjLF6FiZ1+8Crsf\nAcP/Bp32zM21HZSfFX6UUo6xaxA50QeHY2NEwQC8+U/41+HgfTdU/+aCuUUR7EEDvlKqCbuCcaIP\njhH9K5g88gAqyssQQumirC/q27gK7j8uVK++15Fw6VtQeVHeVrbMhpaXrFJK2SJe+QW7BpHHDesd\nM4cf+eCIvlZW+evhf3fCa1Oh1U4wcjoccHqLKHZmNw34ShWhpgOq3ppaJjzzPvBTIM40GOfFTC/v\nslCxs40fwv6nwQl3QNvOubt+ntGAr1SRqar2ctXMFQSazNCLDKjaGZBz1otvqn4HLL4VlvwT2u0K\n5zwJvU/IfTvyjAZ8pYpIpGffNNhH5EvpjIx88d9QsbOtX8DBF8LQW6B1e6dblRc04CtVRJIVMute\nXpb3ZbHj+nEbvHgjLHsIOuwRmn2zx1FJn1ZMNOArVUQS9eDLPG4G9+mSMLeft1a/ECp29sPXcPjv\n4ejroLSN063KO8UzH0kpFXeqpFuEySMPYPFHm2wvHV1V7WXQlJfZY/x8Bk15OWG9nJRt/xZm/Qae\nOCtUtnj0Ijhukgb7ODTgK1VE4s2x/8uZoQ1E7C7VHBkz8NbUYvjpjiHjoG8MvP90qNjZyio4egKM\neRV2Oziz8xY4TekoVWAS5eCTTZW0u+RBovIKaaeItnlh/pXw8QtQcTAMvwe67pveuYqMBnylCkiy\n+fWRr/GC7bhhvZttI+hxp19+2NY7hmAQ3n0IFt4Y2oFq2G1w6P+BK72SzsVIUzpKFZB4PeqJc1Za\nP0nTGZsZFNS1bXOXzZ/BI8NDA7MV/eF3b8Jhl2qwT5EGfKUKSLyec02tj343L0yaO5+6YDW+YOMI\nn8mesRnX5Qn44fW/hYqdbVgBp/wNzp8DHfdIqz3FTlM6ShWQeDl4CAX9ZFMss7Hfb2uPq+GuI6UN\nS75ZGSqLsP5d6H0inPQX2Ll72u1Q2sNXqqAk6zknm2Jp5/66kfGEyO5UAHX+YIJnhPnrYPFkuPcX\nULMWTn8Qzn5cg70NNOArVUBG9K+gQ5vEm4F7a2rjzom3s059ohk6ca1bGgr0r06B/UfCZe+EvhZh\nZcts0ICvVIG56ZT9km4aHm9OvJ116uOlgWJ+4NRvhxeug+nHQt13cO5TMHIatOmY8nVVfBnvaWsn\n3dNWKXtUVXsbbfadSKp71VoVb2/caGUeN/cd+QNHfHgL1KyBQ34Dx9wErXe2vT2FzOqettrDV8oB\nWS03QKinXn3jcdx1Vr+G3no82aqQGSs9FG1ntnOj+TdHvDEaXCVw4XOhgVkN9lmjs3SUyjEri6Ps\nEr3IKl6PO1sbhzdd1RudSxjqWsokzwN04jv+7T+F//vtfeDJ0QbmRUx7+ErlWFqDmTZwYuPwEf0r\neH38EL6YchIV5WV0Yht/9/yN+0r/ymbTnhH1t/CfdhdpsM8R23r4IuIGlgJeY8zJIrIH8CTQCVgG\n/MoYU2/X9ZRqqbIx190KR7ccNIa79/2YPZdNog21TPWdyb2Bk/F4WjE5ix84qjE7UzqXA6uASALu\nduBOY8yTIvJvYDTwLxuvp1SLZHeBslQkqqOTtY1Ptq2DeWOp/GQhWzr25aLto3mjrnPL2lylQNgS\n8EVkN+Ak4FbgShERYAhwbviQh4GJaMBXinHDejfK4UP2UyvJZGVcIRiEZQ/AixPBBOD42+k44GIe\n0/o3jrGrh38XcA2wU/jnTkCNMcYf/nkdoB/jSuFwaiUO28sYf/spzP0DrHkdfnY0nHI3dOhlR1NV\nBjIO+CJyMrDRGLNMRI5O4/ljgDEAPXv2zLQ5SrUIiVIrTrBtXCHghzfvgVcmQ0mrUK36/qN0pWye\nsKOHPwgYLiInAq0J5fDvBspFpCTcy98NiDnR2BgzDZgGoYVXNrRHKZWidMcVovP+R+38NXeXTad8\n24fQ5+TQnPqdds1Wk1UaMp6WaYyZYIzZzRjTCzgbeNkYcx6wGDg9fNgFwOxMr6WUyo50pmxG8v6b\nar5jbMlMpteNw1fj5e1D7oKzH9Ngn4eyufDqWuBJEZkEVAP3Z/FaShWkrM2caSKdcYWpC1azj38V\nt5fex94uL08HjmKS7zzavr8Lr59kexOVDWwN+MaYV4BXwt9/Dgyw8/xKFZNcrsiNnNPyeet+YPQP\n93Jh6QLW04nz66/ltWBfALZleT2BSp+utFUqT02cs9KRFblJffYy2+8ewEUlL/BIYCjD6m5vCPaQ\nm/UEKj1aS0cphzStaBm9G1RVtZea2tiVLrO9Ijeu2q2w4HpY/ijfSgUX1N3IUtOn0SFC8k1YlHM0\n4CvlgKpqL+OeXoEv8NPEtJpaH+OeWgGQ1q5UWfXhHHjuatj+LRxxJcctOpA6SpsdZshOuknZQ1M6\nSjlg6oLVjYJ9RGTD8ES9+Jod9baXU47r+29gxq9g5q+g3S4wZjEcexOdy9vHPLxC0zl5TXv4qmjk\nasaLFYkCeqR98TYP2V4f4IoZy1m6ZguVu3fMzmsyBlY8AS9MAF8tDLkBBl0O7tD2iflYHkIlpwFf\nFYVcz3hJJlFA715exuA+XXh0ydqE53h0yVpmvP0VvmDoTsG211SzFuZeAZ+9BD0GwvC/Q5efNzok\nH8tDqOR0i0NVFOJt/pGt7f2Sqar2cuXM5QSb/PPzuISpZ/Rl6oLVSbcHjCft1xQMwjvTYdHE0M/H\nTgxtOejSzG++s7rFofbwVVFwqgY9xE4lAbhdQrBJHv+sAT0Y0b+CsTOWp309q68pul0Dd97CPe0e\noNOWd2HPY+CUu6Bca1sVGg34qig4VYM+XiqptccVc9D20SVrWfzRJtqXeeJOy0zGymuKtMvnq+O3\n7vlcXvcMtXWlLDtoMgcP/60WOytQeq+mioIT2/tB/LLDkbn3sXhratle78fjSj3oWn1NUxes5mf+\nz5hdegPXeGawKNifY+v+zB9W9dFgX8C0h6+KglODjOmmjHwBQ4c2HrbV+prl+eOpsPqafD9y7g8P\ncknpPLawM5fUj2VB8BAARMsiFDQN+KpoOFGDPl4qqdxCyqZmhw+rUyosD9SueRPmXMalJZ8yw380\nt/rP5TvaNWqvKlya0lEqi+KlkiYO34/yMk/C56YSfJPeSdR9D/OvhgePh0A9rx82nYny20bBXufR\nFz7t4SuVRclSSU0XL0VEgu/EOSstDd6Wt/FwfdX7PPHWVwSMwS3COYf2YNKIA+DTRaF59dvWwaG/\nhSHXM6hVOybvkj8L0VRuaMBXyiEj+lfw1NK1vP7Zlma/cwlcMWM5Vsdtt+7wNVqoFTCGuUs+ZOSa\nWzlo6/PQuTeMXgg9fqpYnm/bLKrs04CvVBqslmlItMIX4I0YwR5C5RMAywO2jRlOcL3NLZ4HKd+y\nHY66Gn5xTWiPWVXUNOArlSIrZRoiHwixBmyja9rbvc69C1v5k+chjne/w3vBPTjfN4Hnj/mdzVdR\nLZUGfKVSFG9u/dQFqxtq2cfLzUd4a2qxd7a74Qz3q1xf8iit8DHZdw7TAyeC6D9x9RP9a1AqRcnK\nNMT6QGhKCA20JlqAZdVuspHJJdM50v0BbwX7MN53MV+YbgCMOrRHxudXhUMDvlIpilf2IDKN0spi\nK0OoAnGZx530wyEeF0HOdy/kmpIZGITrfb/mscAxmPBsa5dA5e4d0zq3Kkw6D1+pFFRVe9le72/2\nuMclDXPYrc6fr6n1MXnkAVSUlyGEFk8lm5sfsad4ear0FiZ6HuGtYB+G1t3Bo4GhDcEeQgO+ju9/\nq/KKBnylUhBvpyp/VJnxccN6W8rPR44ZN6w33cvLGu4MPO74zy7Bz2XuZ3m+1QR+JusZW/9bfu27\nhvV0jnm8Y/vfqrykKR2lUhAvgBpDo+mWZR4XO3zBhOcywM1zV/KjL9iQ1qmp9eFxCW1L3Q1TMyP2\nl8+Z6pnGPq61LHIdwbU7zmMzsbcajNBSCSpaxgFfRFoDrwGtwud72hhzk4jsATwJdAKWAb8yxtRn\nej2lnJRop6paX4AJz7yHP2AadqFKJtagrS9o8EcF+1bUM7ZkFhe75/Et5SwZ8Hcufq1T0imdWipB\nNWVHSqcOGGKM6Qv0A44XkYHA7cCdxpi9gK3AaBuupZSjxg3rnTDlUusLWg72iUTOcKis4oXSa/m/\nkrk8FRzM0pOeZ+CJ5yftuZeXeZg88gBdSasaybiHb0J7JP4Q/tET/s8AQ4Bzw48/DEwE/pXp9ZRK\nl12bmAdsCOjJtGMH40ueYFTJS6wJ7sI59X9kSXA/vhiwDxB7E3EIBfqJw/fTQK9isiWHLyJuQmmb\nvYB/AJ8BNcaYyHSGdUDMv0ARGQOMAejZU7dUU9mR6SbmiVbO2u1oVzW3ee6nK1u5z38if/WfTi2t\nqYjq1esm4iodtgR8Y0wA6Cci5cCzQJ8UnjsNmAahTcztaI9STSVbHRtPVbXXcsXKTHXgO270/Idf\nul/n42AFv/NdwXKzFxA7H6/Fz1SqbJ2lY4ypEZHFwGFAuYiUhHv5uwFeO6+lVCrS2cTcSokEexhO\ndi1houdh2rOdu/wj+af/VOoJzcnv0MbDTadomkZlzo5ZOl0AXzjYlwFDCQ3YLgZOJzRT5wJgdqbX\nUipdqWxinsv0zS5s5VbPAwx1L2NF8Gec57uO1aZxarNNaYkGe2ULO3r43YCHw3l8FzDTGDNPRD4E\nnhSRSUA1cL8N11IFzK5B1VhiDXLGSpPksld/lvsV/ljyGKX4mOQ7jwcDxxPA3exIXTyl7GLHLJ33\ngP4xHv8cGND8GS1PNgORCsl0UDUZq4OcVgqfZaqHfMOUkukMcq/kzcC+TPD/hprWPQjEGSfQxVPK\nLlpaIYlIIPLW1GL4KRBVVeuQhJ0SDaraZUT/ikZlDKYuWN3s/6OV3nQbj4sObazVvInmIsho93Ms\nLL2WA12fM8E3mvN813HFmcez/KbjGDWwZ7OSDLp4StlJjMmfiTGVlZVm6dKlTjejkUFTXo6Zy60o\nL+P18UMcaFFh2mP8/JgrRwX4YspJtlwjVrrG4xLatS6hZoeP7uVl7Kj3xy1ZXFFexuA+XZi1zJvy\nXcDP5Svu8Eyjn+szFgX6c73vIr6mU8N5I39Lejep0iEiy4wxlcmO01o6SaQzu0OlLpVB1XTFuovw\nBU1DgPfW1OJxCR63NCqQVuZxN6xaHTTl5ZSCvQc/v3XP4bKSZ/meNvyh/jLmBA+DqL68t6aWfjcv\nbFgwpQFeZYumdJKIF3A0r2qvWCULPG6xNZ1h5UPaFzS0LS1pVLI4Euyrqr0pzdw5UD5jbukfudLz\nNM8HD2Vo3VTmBA+HGLU0a2p9jHtqhaYKVVZpDz8Jq7M7lA2a5nRszjYmKnwWbVutj+U3HdfosUg6\nyIrW1HFlydOMdj/HRjpwUf3VvBw8KOnzfEGTdCGYUpnQHn4SI/pXNNukQotS2W/qgtXNio5FAqBd\nxg3rTZmn+bTHpmLdvVmdvXOYayUvlI5nTMl8ngwM4bi6OywF+whNFaps0h6+BZpXzb5cjJWM6F/B\n0jVbeOKtrwjEmawQ7+4tWTt2YgcTSh7n3JKX+TLYlbPrr2dJcN9mx7lF2LV967h3GpoqVNmkAV/l\nhWwO2kavnBXiZ4oEOO3gioZ8/c1zVzYM6MYriNyhjYcBdW9xc8n9dKGGe/0ncaf/dH6kVczjzzm0\nB5W7d2Tc0yua7ZzlcQmD+3Rh0JSXdZaOygqdlqnyQqwpk9GzY6yeo+mURiBrK2c78h03eR7hVPcb\nrAr24FrfGN4ze8Y9fu9d2rKjPsj6mlral3nwBYINu1qVl3k4uW+3ZlM+U30PVHGyOi1TA77KG+nO\nQW/aG48QwCUQYwvaDBmGu95goudh2lHLPf5f8q/AcHxJbpib3l00Dea65kOlS+fhqxYnnbGSRLVv\nDPYH+13ZzCTPAxzrrqY6uBfX+MbwidnN0nObNqVpeWZd86GyTQO+yhm7V5FWVXu5auaKuAOwdhKC\nnO1ezISSx/EQ4E++UTwYOJ5ghhPdvDW19L9lYdzVvaADuco+GvBVTthdHC1yvlwE+93la6aUTOcw\n94e8HtiPCf7fsNZ0te38iYK9rvlQdtKAnyVaE6WxZDtONX2/BvfpwuKPNsV9/3JR1dJNgIvcz3NV\nyVPU4+Ea38XMDBxN/Dk7iTXN4SeaMQSh3H2x/90oe2nAz4Jsl/ptiRLlp6+vep/HlqxtCH7emloe\nXbK24ZhY71+2NyfpI2u53TONvq7PWRg4mOt9F7GRDhmd0xAK4pEPsUSvQUAHapXtNOBnQbr7pxay\neAGufZmnUbCPp9YX4Oa5K7O+v2wpPi4tmc3v3LPZRlsuq/8984IDSbdXH63pbJt4s3JA8/YqOzTg\nZ4HOtmguXk0iEeslcxLluu3QTz7lDs+9/Nzl5ZnAEdzi+xU17GTb+XfU+7m+6v2GVFX7Mg9ulxAI\nNl+ApXl7lQ0a8LMgF6V+c8HOcYh4O06NnbHczianpYwfuarkKS5yv8DXdODC+nG8Emy2iVvGtu7w\nNUpV1dT68LiE1qXuRguwImWSVWI6TpY6DfhZUAgVNrMxDhFrnn2uNguP53DXB0wpuY+erk084h/K\nHf6z+IE2Obu+L2jYpU0pK2/RfH0qdJwsPVotMwsKocJmLrYchNgVLAUYtGdH3JJ53jyendnOlJJp\nPF56G37cnFl3Azf6f53TYB9RzKm+dOXq77PQaA8/S1p6hU07xyES3XrHSvVEpmRma479UNdSJnke\noDPb+Jf/FO7yn0YdpbadPzLd0i1i6TW0tFRfPtBxsvRowFcx2TUOYeXWO/rDMVGphEx1ZhsTPQ9z\nsnsJHwZ3Z7Tvaj4wP0vrXIP27MjbX2xtVsMfQsG+zOPmtIMrku5/29JSffmiUMbJci3jlI6I9BCR\nxSLyoYisFJHLw493FJEXReST8NfMJjGrnIqVakknOMW79b5q5gr2GD+fQVNebrStX3YWVBl+6fov\nL7Yax1DXUu7wncnw+j+lHewBvtxcy1kDelARJ8DU+gIs/mhTs9TeqIE9W3SqL1/Y9fdZbOzo4fuB\nq4wx74rITsAyEXkRuBB4yRgzRUTGA+OBa224nsqBeLNqUg1O8QZkI6mOpj1+u2/Ju7GZWz33M8S9\nnGXBvbnGN4bPTOYB1ltTy6xlXiaPPICxM5bHnFq6vqa2xaf28pVdf5/FJuOAb4zZAGwIf/+9iKwC\nKoBTgaPDhz0MvIIG/BYl02BVVe1NWj4AGi9Ks7rvbDJCkPPcLzG+5AkEw02+C/hPYGjKxc7Kyzxx\nF3pF2q3pBWfoh2nqbJ2lIyK9gP7AW0DX8IcBwNeAfdWmVIswdcFqy4uq1tfUUlXtZXudP+Pr7iEb\neLJ0EpM8D/JucG+G1d/Bw4FhaVW2PLlvN8rLPHF/v76mVtMLqsWwbdBWRNoBs4ArjDHfSdSUOmOM\nEZGY//ZFZAwwBqBnz552NUflgVTSM+3LPBkP1roJcLF7PmNLZvEjHq72XcLTgaPIpCzC42+tJca4\nbIPu5WWaXlAthi0BX0Q8hIL9Y8aYZ8IPfyMi3YwxG0SkG7Ax1nONMdOAaRDa8cqO9qj8YDU9Eymx\nkEmw31e+5HbPNA5wfckLgUO4wXchmzIsdgYkDPYCDb14TS+olsCOWToC3A+sMsb8NepXc4ALwt9f\nAMzO9FqqZYmV6mgqMlOlJs06OaX4uKpkJrNLb2BX2cJv6y/n/3xjbQn2yRy+Z0cN8qpFsaOHPwj4\nFfC+iEQKo1wHTAFmishoYA1wpg3XUilyst5I5DpXxKmX08bjaqgemU4VzIPkY+7wTGMv13qeDhzF\nn3yj2Ea7zBqdgnfXbqOq2qtBX7UYdszS+R/xk6THZHp+lb5c1RuxspL2ypnLG6VHXAK3jTyw4edU\nqii04UfGlczgAvdC1tOJ8+uv5bVgX1teSyryqeS1FhJTVuhK2wKWi7r8VlfSRtoTLyBZLX18pOs9\nJnum053NPBIYylT/WWzHuemP+bCUXwuJKas04BewXNQbsfqhkmhQM3qlbTw78wM3lDzKGSWv8Vmw\nG2f6bmCp6ZNZ420QmWvvZA9bN9xRVmnAL2C5WBBkx4dKsgqHw1xv8yfPQ3TkO/7hH87f/CNtLXaW\nrshce6d72FpITFmlAb+AZbsuf1W1F1ecipCpfKjEC0xdqOFmz0Oc6H6blcHd+bXvGlaaXmm2NjkR\nSKVAZ6QXvaPeH7deEGQ/6OtKX2WVBvwClumCoERpikivNlawT/VDpXnAMpzm+i83eP5DGfXc7jub\n+wIn4s/in6vHJbRrXZLyNoqJ1hkEjMlJT78QNtxRuSEmSzXH01FZWWmWLl3qdDMUscsUl3ncDdUd\n423A7RbhL2f2TSnAVVV7Gff0CnwBQwWbmOyZzlHu93k72Jvxvov53HS35TU15Qr36NuXedhe78cX\nyM6/haabl2eDztIpbiKyzBhTmew47eGrmJINBMZLwwSNSRpoYgWnqacdQPWsqVxT8iQG4Xrfr3ks\ncAwmS5uyCXDuoT2ZNOIABk15OeU1AKnIRS5dV/oqKzTgq5jiBSlvuMhZunnjqmov455a0bBxiLem\nln899RzTOjzMCM/7vBo4kOt8o/HSJfMXkYABZrzzFZW7d7QlIJd5XNT6gjF/p7l0lS804KuYEtXB\nmfDM+zF3c7KSN544Z2VDsC/Bzxj3PC4veYbaH1pxpe//eCZ4JKkUO7NSfjkeX8Bw89yVGZdk9rgE\nf5yiOwIM7pOdDy9N46hU6SbmqplkZYrj7eaUaPemqmpvo9TJfvIls0tv4BrPTBYFD+LYuj/zTDD1\nypaZZt237vDRq1P6PfCK8jLatS6Jm/83wKxlXktrDVIRGWPx1tRi+GkqqN3XUYVFB21VI6nsKfvl\nlJNSPmcr6vlDyTNc4p7HFnbJPoqaAAAUM0lEQVTmBt+FLAgOyLTZjogMUMfb8arpsUFjbOuJxxs0\nz8UAsco/OmjbAuXDLfrEOSstBXt3CsVvIgPAlfIRt3vuY0/XBmb6f8Ek/3l8Fy521rbUzY/+IIFE\n9YjzTGTaZZnHxY44+fvoY8G+RVm62EqlQ1M6eSIfbtGrqr2WZ6vEmn8fz7aaLdxc8iBPt7qFUvyM\nqp/ANf5LGoK9xy3069G+RQX7iFpfIGmwj/WcZKuLk4k3EKwDxCoRDfh5ItE0yFy2waoKq4Hlk0Us\naj2eX7kX8aB/GMPqb+d/wQMafu0W4axDevDGZ1tSbW6LlmlPXLdVVOnQlE6eyIdbdKvXShZYqqq9\n3PvCO/xmx3ROc/+XQEkPTt9+E++anzc7NmgMiz/aZGnw1R2njINdOrTxpLzSNl2Z9sR1W0WVDg34\neSIf6qHEa0PbUjflbUotBZaqd9ex+NnpPOK6n3LXdv7mH8E/fvwldcTeCLx7eZnlD5qAMXjcQqnb\nxfb69LdDjKdNaQknHdiN+e9tyGrgj94aMRO62EqlSgN+nsiHeijx2nDrL+NPt2zk+6/pMO8i7na/\nxfvBXpzvm8Aqs3vCpwzu04Un3vrKcs/dFzD4AvYHewiNm8xa5qVVSXYznQatU6+coQE/T4zoX8HS\nNVsagp9bhNMOzm0PLu00gTGw/DFYcB2HBnZwm/8c7g+cSIDE+9mWeVzMWubNapomVbW+QEabqUOo\nRk+i8WfL4x9K2UwHbfNEVbW3UfALGJOVBTu22/ol/GcEzL4UdtmPC1vdybTAKRaCvZvWHnfawTXZ\n5ui5UF7miTlweu6hPeO2TwdWlZM04OeJfJilk9LU0GAAlvwb/nkYrFsGJ/0VLpzP2ScMSbpWNnL3\nUpNmnry8zMPkkQektBYg+rltSxN/WHRoE3u8IZrHJUwcvl/M1caTRhzQ8Dj8tGYh2WpkpbJNV9rm\niT3Gz485U0WALyyuaM2U5dWbGz+COb+HdW/D3sfByXdC+90aft1r/Pyk1yrzuBFMwjnsrnChnOgj\nPC5h6hmh8suprAqGn97L/rcsjDsoGykBfcWM5YnPFW6bzo5R+cDqSlvt4eeJfFhIk3RqaMAHr06F\ne4+EzZ/CyPvg3JmNgj1Yy1HX+gLU+hMvWDIG/npWv0Y96OhgH7kriu5BjxrYM27Pv3t5GVXV3oQz\ncCI98GSvwRi0ho1qcXTQNk9kc5aO1ZINCaeGrq+G2ZfBNx/A/qfB8bdDu9hVIMcN692woUkiyW4u\nu5eXxZx6eH3V+zy2ZG3DHVHAmIb3akT/Cip37xj3vUyUInOLNFwr1v+PeNLdMDwfSmmo4mJLD19E\nHhCRjSLyQdRjHUXkRRH5JPy1gx3XKlQj+lekVH3SqlTy8oP7dGmWfy/3BHhwt3lw3xDYsRnOfhxO\nfyBusG9gIVOYKAcf78OuqtrbKNhHRI93JHovE835j54tFDlHeVnyfD6kvkAuH0ppqOJjSw5fRI4C\nfgAeMcbsH37sDmCLMWaKiIwHOhhjrk10nmLO4WeL1bx8rHz4QFnFPTs9SOf6dXy5+xlc8vWpfLzN\nlbQ3Gu+a0co87pg19SE0aHrTKfvFPH+ycyer4Jno+QLceVa/ZteN7onH27Q91SqVWu1S2Smn1TKN\nMa+JSK8mD58KHB3+/mHgFSBhwFf2s1qyIXqWUDt2ML7kCUaVvITX15X/Hf4AF/+3bcPvk1V8TNTb\nFRoPdFbu3jGltEaic1uZtZMo3WSgITUTHeTbl3mInHrnshJ++NHfsIkLpJd6y4dSGqr4ZDOH39UY\nsyH8/ddA11gHicgYYAxAz549s9ic4mS1ZEMk0Ax2VXOr5366spX7/Cdyp/90OrzbgVpf43MkylvH\nu2as3muq5QES7U5lZQFX5FrxZuGsD2/hGH23E11BdOsOHx63UF7mYVutL+3cez6U0lDFJyezdEwo\nbxTzX6MxZpoxptIYU9mlS3b3MS0Ukd2j9hg/n0FTXk6Y97VaVXGf9j7u9PyDB0un8r1pw2n1N3Or\nfxQdyjuk3BvNZiXHccN6x53n36GNJ+H7Ennfxs5YnnAmT6w1EdF8AUPbViV8MeUkXh8/JK1xFq12\nqZyQzR7+NyLSzRizQUS6ARuzeK2i0bT3mSy9krRcgjHwwSyeNVcjru+4yz+Sf/hH4KOk0cyWVHqj\nqZRosDpTJfq41jE2DHdJqPcdmXLZ9H1p+r7FuhuIvN6xSebgQ2qpl0SvUWfpqFyybeFVOIc/L2rQ\ndiqwOWrQtqMx5ppE59BB2+RsHez7bj3MvwpWPwfdD+Ll3jdyw5umWQCKNaAb2Ty8Iln1zATBLtZ5\nIwufos8X6ziPS2jXuoStO3wJNzKPvC/x3rdYWw9aGXS2+n5bfY1KZSKng7Yi8gShAdrOIrIOuAmY\nAswUkdHAGuBMO65V7GwZ7DMG3n0YFt4QWkx13CQY+DuGuNwM+UXzw6N7o96a2kYBNtEdRrK7kUTl\nJKLPFes4X9BgTCh4Jkq/RN6XeO9P0JhmK5mTzcFPJfVi9TUqlQt2zdI5J86vjrHj/OonGQ/2bfkc\n5vwBvvwv9DoSTrkbOu2Z8CnRvfRYm5DEC2DJgp3VD694x1nZjjHyvqTyvjVNt0Rm6dTsSG2Qtqra\nG/dOQWfjKCdoaYUWJu3BvmAA3rgH/nk4rF8OJ98F58+xFOyjFwjFmwkTHcAig6PJgp3VchLpzlyJ\n3mjE6vsWPbALcN7AnrRtVZJWsI/czcSis3GUE7S0QguTaLAvbr78mw9hzmXgXQY/Pz5U2bK9tXRC\nshkrEZEAZqWgmSE0FjG4T5dmC69iBeF4ZSdae1xx6+IIoWAdeb+sDJLGSkE9umRtw++TDZBHS/S+\n6Wwc5RStllkgYgXanT2GGfu+wT6fTIPWO8MJd4Tq4KRQVjheFc9o0YOQVgY8o5932sEVLP5oU1qz\neYCYHy6JVuomYrXtVgZsE71vd8VYzatUJnI6aKuc17RH2Vc+5Xa5jz6rv4IDzggVO2vbKeXzxst9\nx5rdAqnlpmt9ARZ/tMnSbJemC7SaVssMGJN0xlAyVttu5bhEi8802CunaA6/QESCUGvquK7kMZ4p\nvYn2sp3R9VfDadPTCvYQP/f9lzP7xlx4lGpuOp3By+hxBWheLTNdVttu5ThdWKXykfbwC0T38jJ6\nfreUySXT6eX6hkf9xzDFfw7ty9ML9BGpLhBKNd+ezuBlst3B0l3MZKUkstWgrQurVD7SgF8IftzG\nY10fp9ePT/FFsCtn11/PkuC+tvUoU6l3Ey/QQfN8e7rti3dXEBlUtboK2UrbB/fpYmmMId75NMCr\nfKIBv6Vb/TzMG0uvH77hk70u4pKvjuOLbcGM89mZSBToIou33CLNatg3FW/WUaJxhUwXOWmQVoVM\nA35Ltf1beP4a+GAW7LIfnP0Ye1cczMtOtyuBSCC10gtPtEo3XtooXipGFzkpFaKDti2NMfDeU3DP\nIfDhHDj6OhjzClQc7HTLLEmWf7dyXLwdreLtQ6uLnJQK0R5+S7JtHcy7Ej5ZALsdAsP/Drvs43Sr\nUpJpOYXI4/FSL9naF1ipQqABvyUIBuHdh2DhjWACcPwUGDAGXO6kT803VmvapFMzSGfGKJWYBvx8\nt/mzULGzNf+DPX4RKnbWcQ+nW5W2ePl3q+UUkvXWddBVqfg04OergB+W/AMW3wbuVjD8Hug/KqWy\nCPnIai9ce+tK2U9r6eSjrz8IFTtbXw29T4KT/gI7d3O6VUqpPKW1dFoifx289mf431+hrAOc8RDs\nO6LF9+qVUvlBA36++OodmH0pfLsa+p4Dw26DNh2dbpVSqoBowHda/XZ46U/w1r9h5wo472nYe6jT\nrVJKFSAN+E76bDHM/QPUrIVDLoZjb4JWOzndKqVUgdKA74TaGlj4R6h+FDruCb9+HnY/3OlWKaUK\nnAb8XFs1D+ZfBds3wRFj4RfXgkeX/iulsk8Dfq78sBGeGwcfVkHXA+DcJ6F7f6dbpZQqIlkvniYi\nx4vIahH5VETGZ/t6eccYWP5EqNjZ6udgyA0wZrEGe6VUzmW1hy8ibuAfwFBgHfCOiMwxxnyYzevm\njZqvYN4V8Oki6HFoqNhZFy3kpZRyRrZTOgOAT40xnwOIyJPAqUBhB/xgEJbeD4smhnr4J9wRmoXj\n0mrUSinnZDvgVwBfRf28Djg0y9d01refwJzfw9o34WeDQ8XOOuzudKuUUsr5QVsRGQOMAejZs6fD\nrclAwA9v/h0WTw7Nujn1n9DvXC2LoJTKG9kO+F6gR9TPu4Ufa2CMmQZMg1DxtCy3Jzs2vBcqdrZh\nBewzHE78M+zU1elWKaVUI9kO+O8Ae4vIHoQC/dnAuVm+Zu74foTX7oD/3QVtOsGZj8C+pzrdKqWU\niimrAd8Y4xeRy4AFgBt4wBizMpvXzJm1S0K5+m8/hn7nwXGTtNiZUiqvZT2Hb4x5Dngu29fJmbof\n4KWb4e37oP1uMGoW7HWs061SSqmkHB+0bVE+fQnmXgHbvoIBF8MxN2qxM6VUi6EB34odW2Dh9bD8\nMei0N1z0AvQc6HSrlFIqJRrwk/lwNsy/GnZshiOvhqPGgae1061SSqmUacCP5/tv4LmrYNVc2PXA\nUK6+24FOt0oppdKmAb8pY2D547BgQmja5bET4bDfg1vfKqVUy6ZRLNrWNTD3cvh8MfQ8LFTsrPPe\nTrdKKaVsoQEfQsXO3rkPFt0cKoVw4p+hcrQWO1NKFRQN+Js+DpVF+Oqt0Hz6k++E8hZc00cppeIo\n3oAf8MHrd8Ort4OnDYz4N/Q9W4udKaUKVnEG/PXLYfZl8M37sO8IOHEqtNvF6VYppVRWFVfA99WG\nevSv/w3adoazHoV9TnG6VUoplRPFE/DXvBnK1W/+FPqPChU7K+vgdKuUUipnCj/g130fmn3zzn1Q\nvjv8qgr2HOx0q5RSKucKO+B/8iLMGwvb1sHA38GQ66G0rdOtUkopRxRmwN+xBV6YAO89CV36wOgX\nocchTrdKKaUcVVgB3xj4sAqeGwe1W+Goa+Coq6GkldMtU0opxxVOwP9uAzx3NXw0D7r1C+Xqd93f\n6VYppVTeKIyA//FCmPUbCNTB0Ftg4KVa7EwppZoojKjYac9Qjv6EO0LfK6WUaqZwAv6oWU63Qiml\n8pqWg1RKqSKhAV8ppYqEBnyllCoSGQV8ETlDRFaKSFBEKpv8boKIfCoiq0VkWGbNVEoplalMB20/\nAEYC90Y/KCL7AmcD+wHdgUUi8nNjTCDD6ymllEpTRj18Y8wqY8zqGL86FXjSGFNnjPkC+BQYkMm1\nlFJKZSZbOfwK4Kuon9eFH2tGRMaIyFIRWbpp06YsNUcppVTSlI6ILAJ2jfGrPxpjZmfaAGPMNGAa\nQGVlpcn0fEoppWJLGvCNMcemcV4v0CPq593CjyW0bNmyb0VkTRrXA+gMfJvmc53S0tqs7c0ubW92\nFXJ7d7dyULZW2s4BHheRvxIatN0beDvZk4wxXdK9oIgsNcZUJj8yf7S0Nmt7s0vbm13a3synZf5S\nRNYBhwHzRWQBgDFmJTAT+BB4AbhUZ+gopZSzMurhG2OeBZ6N87tbgVszOb9SSin7FNJK22lONyAN\nLa3N2t7s0vZmV9G3V4zRiTFKKVUMCqmHr5RSKoEWH/BFpLWIvC0iK8J1fW52uk1WiIhbRKpFZJ7T\nbUlGRL4UkfdFZLmILHW6PcmISLmIPC0iH4nIKhE5zOk2xSMivcPva+S/70TkCqfblYiIjA3/W/tA\nRJ4QkdZOtykZEbk83N6V+fj+isgDIrJRRD6IeqyjiLwoIp+Ev3bI9DotPuADdcAQY0xfoB9wvIgM\ndLhNVlwOrHK6ESkYbIzp10Kmtd0NvGCM6QP0JY/fZ2PM6vD72g84GNhBnIkQ+UBEKoA/AJXGmP0B\nN6G6WXlLRPYHLiZU3qUvcLKI7OVsq5p5CDi+yWPjgZeMMXsDL4V/zkiLD/gm5Ifwj57wf3k9MCEi\nuwEnAdOdbkuhEZH2wFHA/QDGmHpjTI2zrbLsGOAzY0y6iw9zpQQoE5ESoA2w3uH2JLMP8JYxZocx\nxg+8SqjoY94wxrwGbGny8KnAw+HvHwZGZHqdFh/woSE9shzYCLxojHnL6TYlcRdwDRB0uiEWGWCh\niCwTkTFONyaJPYBNwIPhlNl0EWnrdKMsOht4wulGJGKM8QJ/BtYCG4BtxpiFzrYqqQ+AI0Wkk4i0\nAU6kcSWAfNXVGLMh/P3XQNdMT1gQAd8YEwjfEu8GDAjfwuUlETkZ2GiMWeZ0W1JwhDHmIOAE4FIR\nOcrpBiVQAhwE/MsY0x/Yjg23wtkmIqXAcOApp9uSSDiPfCqhD9buQFsRGeVsqxIzxqwCbgcWEloI\nuhxoUQtBTWg6ZcaZi4II+BHhW/fFNM+F5ZNBwHAR+RJ4EhgiIo8626TEwr06jDEbCeWX87nU9Tpg\nXdRd3tOEPgDy3QnAu8aYb5xuSBLHAl8YYzYZY3zAM8DhDrcpKWPM/caYg40xRwFbgY+dbpMF34hI\nN4Dw142ZnrDFB3wR6SIi5eHvy4ChwEfOtio+Y8wEY8xuxphehG7hXzbG5G0PSUTaishOke+B4wjd\nIuclY8zXwFci0jv80DGESnzku3PI83RO2FpgoIi0EREh9P7m7aB4hIjsEv7ak1D+/nFnW2TJHOCC\n8PcXABlXJ85W8bRc6gY8LCJuQh9gM40xeT/VsQXpCjwb+rdNCfC4MeYFZ5uU1O+Bx8Jpks+BXzvc\nnoTCH6RDgUucbksyxpi3RORp4F3AD1TTMlawzhKRToCPUG2vvBrIF5EngKOBzuH6ZDcBU4CZIjIa\nWAOcmfF1dKWtUkoVhxaf0lFKKWWNBnyllCoSGvCVUqpIaMBXSqkioQFfKaWKhAZ8pZQqEhrwlVKq\nSGjAV0qpIvH/ABfz4KjuR28AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "tags": [] } } ] }, { "cell_type": "code", "metadata": { "id": "w4YF1l7etyVT", "colab_type": "code", "colab": { "base_uri": "https://localhost:8080/", "height": 34 }, "outputId": "d81d1bfc-67db-44c1-b088-dfd222391393" }, "source": [ "from sklearn.metrics import r2_score\n", "\n", "r2_score(y, y_hat)" ], "execution_count": 26, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "0.48352545599133423" ] }, "metadata": { "tags": [] }, "execution_count": 26 } ] }, { "cell_type": "code", "metadata": { "id": "boTJ6cwXx-PA", "colab_type": "code", "colab": {} }, "source": [ "" ], "execution_count": 0, "outputs": [] } ] }