{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Regression 102" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is our first basic regression example.\n", "\n", "In this example, we won't be doing anything complicated or using complex databases yet." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "from sklearn import linear_model\n", "from sklearn.preprocessing import PolynomialFeatures\n", "lr = linear_model.LinearRegression() #In sklearn we define an object as our ML operator\n", "import pylab as plt\n", "%matplotlib inline \n", "#again, this is only in jupyter notebooks" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(40L,)\n", "[ 0.00000000e+00 2.50000000e-01 1.00000000e+00 2.25000000e+00\n", " 4.00000000e+00 6.25000000e+00 9.00000000e+00 1.22500000e+01\n", " 1.60000000e+01 2.02500000e+01 2.50000000e+01 3.02500000e+01\n", " 3.60000000e+01 4.22500000e+01 4.90000000e+01 5.62500000e+01\n", " 6.40000000e+01 7.22500000e+01 8.10000000e+01 9.02500000e+01\n", " 1.00000000e+02 1.10250000e+02 1.21000000e+02 1.32250000e+02\n", " 1.44000000e+02 1.56250000e+02 1.69000000e+02 1.82250000e+02\n", " 1.96000000e+02 2.10250000e+02 2.25000000e+02 2.40250000e+02\n", " 2.56000000e+02 2.72250000e+02 2.89000000e+02 3.06250000e+02\n", " 3.24000000e+02 3.42250000e+02 3.61000000e+02 3.80250000e+02]\n" ] } ], "source": [ "#Let's create some data\n", "x = np.arange(0,20,0.5) #Data from 0 to 10 in steps of 1\n", "y = np.power(x,2)\n", "print x.shape\n", "print y" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEACAYAAACwB81wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF0dJREFUeJzt3X+s3Xd93/HnK82CLmWNAixOldAQFkTMtihhwu2UShxG\n7RuYlhCQsoxqM6OdkVIoGtYUG4R86WiDKzkd6hSpEFoMIspSNhNnKr5xGg6TpzVhJRkBm9RV5yzJ\nsGEDKiKsKVne++N8bQ6Xc5x77j3nnl/Ph3Tk7/mc7/mezzcn930/9/19fz7fVBWSpNl33rg7IEna\nGAZ8SZoTBnxJmhMGfEmaEwZ8SZoTBnxJmhMDB/wk5yX5WpKDzfOLkjyQ5Ikky0ku7Np3d5LjSY4l\n2TbMjkuSBrOWEf4HgKNdz3cBD1bV64CHgN0ASV4P3AxsBt4K3Jkk6+uuJGmtBgr4SS4D3gbc1dV8\nI7C/2d4PvL3ZvgG4p6qer6oTwHFgy7p6K0las0FH+L8H/Guge3rupqo6BVBVJ4GLm/ZLgae69num\naZMkjcGqA36SfwScqqrHgHOlZlyrQZIm0PkD7HsdcEOStwELwN9M8jngZJJNVXUqySXAd5r9nwFe\n1fX+y5q2n5DEXxCStAZVNdB10VWP8KvqQ1X1C1X1GuAW4KGq+mfA/cC7m922A/c12weBW5JckOQK\n4ErgkT7HntnHnj17xt4Hz8/zm8fzm+Vzq1rbOHmQEX4/HwfuTfIe4Ek6lTlU1dEk99Kp6HkOuLXW\n2ktJ0rqtKeBX1VeArzTb3wN+pc9+twO3r7l3kqShcabtiLVarXF3YaQ8v+k2y+c3y+e2Vhl3liWJ\nmR5JGlASalQXbSVp0iwvL7Nt2zvZtu2dLC8vj7s7E88RvqSptLy8zE03bef06b0ALCzcxoED+1lc\nXBxzzzbGWkb4BnxJU2nbtndy+PANdKrBAfazdetBHnjgP4yzWxvGlI4kqa9h1OFL0obbuXMHR45s\n5/TpzvOFhdvYuXP/ud8050zpSJpay8vL7Nv3SaDzC2Be8vdgDl+S5oY5fElSXwZ8SZoTBnxJmhMG\nfEkTzxm1w+FFW0kTbd5n1PZjlY6kmTPvM2r7sUpHktTXIDcxf0mSh5M8muSbSX6nad+T5OkkX2se\n13e9Z3eS40mOJdk2ihOQNNt27tzBwsJtwH5gfzOjdse4uzWVBkrpJHlpVf0oyc8A/wXYSeduVz+s\nqjtW7LsZuBt4I50bmD8IvHZl/saUjqQXM88zavtZS0pnoLV0qupHzeZL6Px18P0zn91j9xuBe6rq\neeBEkuPAFuDhQT5TkhYXFw3yQzBQDj/JeUkeBU4C7ao62rz0viSPJbkryYVN26XAU11vf6ZpkySN\nwaAj/BeAa5P8HPBAkjcBdwK/VVWV5GPAPuDXBznu0tLS2e1Wq+W9KCVphXa7TbvdXtcx1lyWmeQj\nwI+qal9X2+XA/VV1dZJdQFXV3ua1Q8Ceqnp4xXHM4Us6y3z96oy0LDPJK8+ka5IsAFuBx5Jc0rXb\nO4BvNNsHgVuSXJDkCuBK4JFBOidpvpyZZHX48A0cPnwDN9203Zm1QzRISufngf1JQucXxeeq6k+T\nfDbJNcALwAngvQBVdTTJvcBR4DngVofyks5l375PNjNqO5OsTp/utDnKH45VB/yqehx4Q4/2f36O\n99wO3L62rkmShslbHEqaGN62cLRcS0fSRPGi7eq4eJokzQkXT5Mk9WXAl6Q5YcCXNBbexWrjmcOX\ntOG8i9X6edFW0lTwLlbr50VbSVJfTryStOGcYDUepnQkjYUTrNbHHL4kzQlz+JKkvgz4kjQnDPiS\nRsoJVpPDHL6kkXGC1eiM+haHL0nycJJHk3wzye807RcleSDJE0mWz9wGsXltd5LjSY4l2TZIxyRN\nv5+8g1Un8J+pzNHGW3XAr6r/C7y5qq4Frgb+YZLrgF3Ag1X1OuAhYDdAktcDNwObgbcCdza3R5Qk\njcFAOfyq+lGz+ZLmvd8HbgTOzJjYD7y92b4BuKeqnq+qE8BxYMt6OyxpeuzcuYOFhdvohIb9zQSr\nHePu1twaKOAnOS/Jo8BJoF1VR4FNVXUKoKpOAhc3u18KPNX19meaNklzYnFxkQMHOuvkbN160Pz9\nmA20tEJVvQBcm+TngOUkLWDlFdeBr8AuLS2d3W61WrRarUEPIWlCLS4uGuSHoN1u026313WMNVfp\nJPkIcBr4NaBVVaeSXAJ8uao2J9kFVFXtbfY/BOypqodXHMcqHUka0KirdF55pgInyQKwFXgUOAi8\nu9ltO3Bfs30QuCXJBUmuAK4EHhmkc5Kk4Rkkh//zwJebHP6fAQer6k+BvcDWJE8AbwE+DtDk9+8F\njgJ/AtzqUF6abU6ymmxOvJI0FE6y2liulilpbLyL1cZytUxJUl/e8UrSUHgXq8lnSkfS0HgXq41j\nDl+S5oQ5fElSXwZ8SQOz3n46mdKRNBDr7SeDOXxJI2e9/WQwhy9J6ss6fEkDsd5+epnSkTQw6+3H\nzxy+JM0Jc/iSpL4M+JLOyZr72WFKR1Jf1txPrlHf4vCyJA8l+WaSx5O8v2nfk+TpJF9rHtd3vWd3\nkuNJjiXZNkjHJI3fvn2fbIL9dqAT+M9crNX0GaQs83ngg1X1WJKXAX+e5HDz2h1VdUf3zkk2AzcD\nm4HLgAeTvNbhvCSNx6oDflWdBE42288mOQZc2rzc68+KG4F7qup54ESS48AW4OH1dVnSRrHmfras\n6aJtklcD1/Dj4P2+JI8luSvJhU3bpcBTXW97hh//gpA0BRYXFzlwoLN0wtatB83fT7mBZ9o26Zwv\nAB9oRvp3Ar9VVZXkY8A+4NcHOebS0tLZ7VarRavVGrRbkkZkcXHRID8B2u027XZ7XccYqEonyfnA\nfwK+VFWf6PH65cD9VXV1kl1AVdXe5rVDwJ6qenjFe0zrS9KANmLi1R8CR7uDfZJLul5/B/CNZvsg\ncEuSC5JcAVwJPDLg50naINbbz75Vp3SSXAf8KvB4kkeBAj4EvCvJNcALwAngvQBVdTTJvcBR4Dng\nVofy0mRaWW9/5Mh28/UzyIlXklzjfgq5lo4kqS/Xw5dkvf2cMKUjCXCN+2njeviSNCfM4UtaFUsw\n55MjfGnOuOTxbDClI+lFWYI5G0zpSJL6sixTmjOWYM4vUzrSHLIEc/qZw5ekOWEOX5LUlwFfmmHW\n26ubKR1pRllvP9vM4Us6y3r72TbSHH6Sy5I8lOSbSR5P8ptN+0VJHkjyRJLlrpuYk2R3kuNJjiXZ\nNkjHJEnDNUgO/3ngg1X1d4B/APxGkquAXcCDVfU64CFgN0CS1wM3A5uBtwJ3Jhnot5Gktdu5cwcL\nC7cB+4H9Tb39jnF3S2O05pROki8C/655vKmqTjX3t21X1VU9bmL+JWDJm5hLG8d6+9m1YTn8JK8G\n2sDfBZ6qqou6XvteVb08ye8D/7Wq7m7a7wL+pKr+44pjGfAlaUAbUoef5GXAF4APVNWzdG5m3s3o\nLW0wyy+1GgOtpZPkfDrB/nNVdV/TfCrJpq6Uznea9meAV3W9/bKm7acsLS2d3W61WrRarUG6Jc21\nleWXR45st/xyBrXbbdrt9rqOMVBKJ8lngf9dVR/satsLfK+q9ia5DbioqnY1F20/D/wicClwGHjt\nyvyNKR1pfSy/nE9rSemseoSf5DrgV4HHkzxKJ3XzIWAvcG+S9wBP0qnMoaqOJrkXOAo8B9xqZJek\n8XHilTTlnFE7n5xpK80pyy/njwFfkuaEyyNLc8ASTK2VI3xpipiv1xmmdKQZZwmmzjClI0nqa6CZ\ntpLGa+fOHRw5sp3TpzvPOytg7h9vpzQ1TOlIU8YSTIE5fGmmGNh1LgZ8aUZYjaMXY8CXZoTVOHox\nVulIkvqySkeaQFbjaBRM6UgTyou2Ohdz+JI0J8zhS1PKBdG0EVYd8JN8OsmpJF/vatuT5OkkX2se\n13e9tjvJ8STHkmwbdselWXGmBPPw4Rs4fPgGbrppu0FfI7HqlE6SXwaeBT5bVVc3bXuAH1bVHSv2\n3QzcDbyRzs3LH6TH/WybfU3paK5Zgqm1GGlKp6qOAN/v9bk92m4E7qmq56vqBHAc2DJIxyRJwzWM\nHP77kjyW5K4kFzZtlwJPde3zTNMmaYWdO3ewsHAbsB/Y35Rg7hh3tzSD1hvw7wReU1XXACeBfevv\nkjRfFhcXOXCgk8bZuvWgSyhoZNY18aqqvtv19FPA/c32M8Crul67rGnraWlp6ex2q9Wi1Wqtp1vS\nxOpXW7+4uGiQ1zm1223a7fa6jjFQHX6SVwP3V9Xfa55fUlUnm+1/Bbyxqt6V5PXA54FfpJPKOYwX\nbTXnXBBNw7SWi7arHuEnuRtoAa9I8j+BPcCbk1wDvACcAN4LUFVHk9wLHAWeA241qmve7dv3ySbY\nd6pxTp/utBnwtVFWHfCr6l09mv/oHPvfDty+lk5JkobPxdOkDeKCaBo319KRRqDfxVkXRNOwuHia\nNAG8OKuNYMCXJoBLJWgjuFqmJKkvL9pKQ+bFWU0qUzrSCHhxVqNmDl8aA4O7xsGAL20wK3I0LgZ8\naYNZkaNxsUpHktSXVTrSOliRo2liSkdaJZdL0CQxhy+NiBdnNWkM+NKIeHFWk8aLtpKkvlYd8JN8\nOsmpJF/varsoyQNJnkiynOTCrtd2Jzme5FiSbcPuuLSRdu7cwcLCbcB+YH9zcXbHuLslDWSQEf4f\nASsTlruAB6vqdcBDwG6A5p62NwObgbcCdyYZ6E8PaZyWl5fZtu2dbNv2TpaXl1lcXOTAgU4aZ+vW\ng+bvNZUGvYn55XRuYn518/xbwJuq6lSSS4B2VV2VZBdQVbW32e9LwFJVPdzjmObwNVG8QKtpMI4c\n/sVVdQqgqk4CFzftlwJPde33TNMmTbyfvNl4J/CfKbuUptmwJ16taai+tLR0drvVatFqtYbUHUma\nDe12m3a7va5jrDelcwxodaV0vlxVm3ukdA4Be0zpaBqY0tE02IiUTprHGQeBdzfb24H7utpvSXJB\nkiuAK4FHBvwsaeRWXpwFvECrmbXqEX6Su4EW8ArgFLAH+CLwx8CrgCeBm6vqB83+u4FfA54DPlBV\nD/Q5riN8jYUjeU0zZ9pKA3D2rKaZM20lSX25PLLmlksba96Y0tHc6LWMsUsba1qZw5f68AKtZo0B\nX+rDC7SaNV60lST1ZcDXzOk1mcrljSVTOpox58rVe4FWs8QcvuaeuXrNC3P4kqS+DPiaWubqpcGY\n0tFUMleveWcOX3PDXL3mnTl8SVJfjvA1FVamaQCXStBcM6WjmdQvXw+Yq9fcGlvAT3IC+GvgBeC5\nqtqS5CLg3wOXAyfo3A3rr3u814CvczJfL/20cebwX6BzM/Nrq2pL07YLeLCqXgc8BOwe0mdJktZg\nWAE/PY51I51iaJp/3z6kz9IMs7ZeGp1hpXT+CvgB8P+AP6iqu5J8v6ou6trne1X18h7vNaUjwNp6\naRBrSekM6xaH11XVt5P8LeCBJE8AK6N436i+tLR0drvVatFqtYbULU2Tffs+2QT7Tq7+9OlO2+Li\n4tmHNK/a7TbtdntdxxhKwK+qbzf/fjfJF4EtwKkkm6rqVJJLgO/0e393wNf86FVqKam3lYPhj370\nowMfY90BP8lLgfOq6tkkPwtsAz4KHATeDZwZst233s/S7FiZvjlyZDsf/vD7OXLkNm8qLo3IMEb4\nm4ADSao53uer6oEk/w24N8l7gCeBm4fwWZoRvdI3X/nKQQ4c2N816ncilTRM6w74VfU/gGt6tH8P\n+JX1Hl/Tb5ALrubqpdFxpq1G6lyzZF0aQVo7l1bQxDnXLFlLLaW1G2dZpgQMVnlj+kbaWAZ8DY2V\nN9JkM+BraKy8kSabAV9rYuWNNH28aKuBWXkjjZ9VOhq6XiN5K2+k8bNKR0PV6yLsmZF8P6ZvpMll\nwNdZK0fn/Vav3LlzB0eObLfyRpoyBnwBvUfzV111Vc99FxcXrbyRppA5/Dm02rz8tdd+im996y+9\nCCtNIHP4elGD5OVf+cpNHDjwEUfy0oxwhD/jeuXle1XY7Ny5w5JKaYo4wp9jvdI05uUldTPgT5nV\nBvYzwXtllQ18ioWF3mvbWFIpzbaRB/wk1wP/FjgP+HRV7R31Z86KXitPrjawn3nfSublpTlWVSN7\n0AnyfwlcDvwN4DHgqhX71Cw6dOhQbd36jtq69R116NChc7b3a1tY2FTwmYLP1MLCprr22jc1z6t5\nfObs+3q19zpGd18kTa8mdg4Wkwd9w0AHh18CvtT1fBdw24p9RvYfZFz6Bdpe7R/72Md67tsriL/8\n5X974MDe7xePpOm2loA/6pTOpcBTXc+fBraM+DPH7lwplpXtd9zxb1adjrn88ks4ffqn8+/nuuBq\nXl7SGRNx0XZpaensdqvVotVqja0vk6LX8gW33965uGpgl+ZPu92m3W6v6xgjrcNP8kvAUlVd3zzf\nRefPkL1d+9Qo+zAOgywf/OEPv5/f/u3f71n/7sqTkvqZuOWRk/wM8ATwFuDbwCPAP62qY137zFzA\nh/43COlXVmlglzSIiQv4cLYs8xP8uCzz4yten8mAL0mjNJEB/0U7YMCXpIGtJeCfN6rOSJImiwFf\nkuaEAV+S5oQBX5LmhAFfkuaEAV+S5oQBX5LmhAFfkuaEAV+S5oQBX5LmhAFfkuaEAV+S5oQBX5Lm\nhAFfkuaEAV+S5sS6An6SPUmeTvK15nF912u7kxxPcizJtvV3VZK0HsMY4d9RVW9oHocAkmwGbgY2\nA28F7kwy0EL9s2K9Nx2edJ7fdJvl85vlc1urYQT8XoH8RuCeqnq+qk4Ax4EtQ/isqTPr/9N5ftNt\nls9vls9trYYR8N+X5LEkdyW5sGm7FHiqa59nmjZJ0pi8aMBPcjjJ17sejzf//mPgTuA1VXUNcBLY\nN+oOS5LWZmg3MU9yOXB/VV2dZBdQVbW3ee0QsKeqHu7xPu9gLklrMOhNzM9fz4cluaSqTjZP3wF8\no9k+CHw+ye/RSeVcCTzS6xiDdliStDbrCvjA7ya5BngBOAG8F6Cqjia5FzgKPAfcWsP6U0KStCZD\nS+lIkibbRMy0PdcErmmW5Pok30ryF0luG3d/hi3JiST/PcmjSXqm7KZFkk8nOZXk611tFyV5IMkT\nSZa7qtCmTp/zm5mfuySXJXkoyTebwpLfbNpn4jvscX7vb9oH+g4nYoSfZA/ww6q6Y9x9GZYk5wF/\nAbwF+F/AV4FbqupbY+3YECX5K+DvV9X3x92X9Uryy8CzwGer6uqmbS/wf6rqd5tf2BdV1a5x9nOt\n+pzfzPzcJbkEuKSqHkvyMuDP6cwH+hfMwHd4jvP7JwzwHU7ECL8xaxdvtwDHq+rJqnoOuIfOFzRL\nwmT9P7RmVXUEWPmL60Zgf7O9H3j7hnZqiPqcH8zIz11Vnayqx5rtZ4FjwGXMyHfY5/zOzG1a9Xc4\nST+svSZwTbOVk8+eZvYmnxVwOMlXk/zLcXdmBC6uqlPQ+YEDLh5zf0Zh1n7uSPJq4Brgz4BNs/Yd\ndp3fmTL3VX+HGxbwB5zANfV/Ys6J66rqDcDbgN9o0gazbPz5z+GauZ+7Jt3xBeADzUh45Xc21d9h\nj/Mb6Dtcb1nmqlXV1lXu+ing/lH2ZYM8A/xC1/PLmraZUVXfbv79bpIDdNJYR8bbq6E6lWRTVZ1q\ncqjfGXeHhqmqvtv1dOp/7pKcTycYfq6q7muaZ+Y77HV+g36HE5HSab6IM7oncE2zrwJXJrk8yQXA\nLXQmpM2EJC9tRhsk+VlgG9P/vYWfzIceBN7dbG8H7lv5hinzE+c3gz93fwgcrapPdLXN0nf4U+c3\n6Hc4KVU6n6WTkzo7getM3m2aNSVSn6Dzi/XTVfXxMXdpaJJcARyg8yfy+cDnp/n8ktwNtIBXAKeA\nPcAXgT8GXgU8CdxcVT8YVx/Xo8/5vZkZ+blLch3wn4HH6fw/WcCH6Mzwv5cp/w7PcX7vYoDvcCIC\nviRp9CYipSNJGj0DviTNCQO+JM0JA74kzQkDviTNCQO+JM0JA74kzQkDviTNif8Pnqm0l2mF7XsA\nAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.scatter(x,y)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "LinearRegression(copy_X=True, fit_intercept=True, n_jobs=1, normalize=False)" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#lr.fit(x,y)#<-- This doesn't work!!\n", "#We need to reshape from 1D to 2D vectors\n", "\n", "x = x[:,None]\n", "y = y[:,None]\n", "lr.fit(x,y)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHEtJREFUeJzt3X+Q3PV93/HnC4OSw8ZELtFRS0Qm2BjJU2phR+4MbbUt\n2T2TTgVIHZXJ2JUNtTqhIGZ8TSVZ4+E0MZXl9Ih/tGRGxLHOLikh0VwQU6y9o3jxqBMbYsAQS2B5\nYqkgm8NgG6PmEkT87h/7vcty2j20u9/d73d3X4+ZHX33c9/93ufLcvvez/vzSxGBmZkNtrOyroCZ\nmWXPwcDMzBwMzMzMwcDMzHAwMDMzHAzMzIyUgoGk8yX9iaQjkr4j6QOSlkqakvSMpLKk82vO3yHp\naHJ+KY06mJlZ69JqGXwOeCAiVgH/GHga2A48GBHvBh4CdgBIWg1sAlYBVwN3SlJK9TAzsxa0HQwk\nvRX4ZxHxJYCIeC0iXgauASaS0yaAa5Pj9cA9yXnHgKPA2nbrYWZmrUujZXAx8KKkL0l6TNJeSecC\nwxExAxARzwPLkvOXA8/WvP5EUmZmZhlJIxicDVwB/PeIuAL4f1RTRAvXufC6F2ZmOXV2Ctd4Dng2\nIv4ieb6fajCYkTQcETOSLgReSH5+Ario5vUrkrLTSHIAMTNrQUQ01RfbdssgSQU9K+nSpOgq4DvA\nAeAjSdlm4L7k+ABwvaQlki4G3gk8ssj1+/Zx2223ZV4H35vvz/fXf49WpNEyANgK3C3pHOCvgI8C\nbwLulXQDcJzqCCIi4rCke4HDwCngpmi19mZmlopUgkFEfBv4tTo/+vUG5+8Gdqfxu83MrH2egZyh\nQqGQdRU6pp/vDXx/va7f768VynOGRpIzSGZmTZJEdLsD2czMep+DgZmZORiYmZmDgZmZ4WBgZmY4\nGJiZGQ4GZmaGg4GZmeFgYGZmOBiYmRkOBmZmhoOBmZnhYGBmZjgYmJkZDgZmZoaDgZmZ4WBgZmak\nFAwkHZP0bUmPS3okKVsqaUrSM5LKks6vOX+HpKOSjkgqpVEHMzNrXVotg58DhYhYExFrk7LtwIMR\n8W7gIWAHgKTVwCZgFXA1cKekprZnMzOzdKUVDFTnWtcAE8nxBHBtcrweuCciXouIY8BRYC1mZpaZ\ntIJBANOSHpX075Oy4YiYAYiI54FlSfly4Nma155IyszMLCNnp3SdKyPih5J+GZiS9AzVAFFr4XMz\nM8uJVIJBRPww+fdHkv6MatpnRtJwRMxIuhB4ITn9BHBRzctXJGV1jY2NzR8XCgUKhUIaVTYz6xuV\nSoVKpdLWNRTR3hd2SecCZ0XESUlvBqaAXcBVwI8jYo+kbcDSiNiedCDfDXyAanpoGnhX1KmIpHrF\nZma2CElERFMDc9JoGQwDk5Iiud7dETEl6S+AeyXdABynOoKIiDgs6V7gMHAKuMmf+GZm2Wq7ZdBJ\nbhmYWSvK5TJ7x8cB2DI6ysjISMY16q5WWgYOBmbWV8rlMpuvu449s7MAbBsaYmJycqACgoOBmQ28\njaUS66en2Zw8nwAOFIvsn5rKslpd1Uow8NpEZmaW2jwDM7Nc2DI6yuZDh6A2TTQ6mnGt8s9pIjPr\nO+5Adp+BmdnAc5+BmZm1xMHAzMwcDMzMzMHAzMxwMDAzMxwMzMwMBwMz6xHlcpmNpRIbSyXK5XLW\n1ek7nmdgZrnnxeea40lnZtaXvPhcczzpzMzMWuKF6sws97z4XOc5TWRmPWHQF59rhvsMzMws2z4D\nSWdJekzSgeT5UklTkp6RVJZ0fs25OyQdlXREUimtOpiZWWvS7EC+FThc83w78GBEvBt4CNgBIGk1\nsAlYBVwN3CmpqQhmZmbpSiUYSFoB/AbwBzXF11AdAUby77XJ8Xrgnoh4LSKOAUeBtWnUw8zMWpNW\ny+D3gN8GahP8wxExAxARzwPLkvLlwLM1551IyszMLCNtDy2V9K+AmYh4QlJhkVNb6gkeGxubPy4U\nChQKi/0KM7PBU6lUqFQqbV2j7dFEkv4L8CHgNWAIOA+YBN4PFCJiRtKFwNciYpWk7UBExJ7k9QeB\n2yLim3Wu7dFEZmZNymQ0UUR8IiJ+JSJ+FbgeeCgiPgzcD3wkOW0zcF9yfAC4XtISSRcD7wQeabce\nZmbWuk4uR/FpoCjpGeCq5DkRcRi4l+rIoweAm/z138zSVC6XKZU2Uipt9AqnZ8iTzsys47o5e7hc\nLnPddZuZnd0DwNDQNiYnJwZqxrJnIJtZ7nR7+elSaSPT0+uhZo3TYvEAU1P7O/L78qiVYOCF6sys\no/aOj7Nndnb+o5nZWfaOjw/UN/Ve4GBgZn1ldHQLhw5tnlvglKGhbYyOTiz+InMwMLPO6vTy0+Vy\nmfHxvUA1EIyMjDA5OVFTNlj9Ba1yn4GZdVynOpDdWVyfO5DNbKC4s7g+b3tpZmYtcZ+BmfUsdxan\nx2kiM+tp9TqQB53TRGbWtxZbYuI8XuE8XsmoZv3BLQMzy71Go4aArs5u7hUeTWRmfanRqKHzeIX1\n09M1pXCgWGT/1FQ2Fc0Jp4nMzKwlHk1kZrm32KihTs5uHiROE5lZZurNTG40Omix8m4tj90r3Gdg\nZj2j3tLWt+zcye23f8HLS7TJwcDMesbGUum0zt/b3vbLHP/x7+LlJdrjDmQzM2uJO5DNLBP1lra+\n5eO3cvvt27y8RAbaThNJ+gXg68CS5HFfRHxC0lLgj4GVwDFgU0S8nLxmB3AD8Bpwa0TUHRTsNJFZ\nf7v99tu5647PAfCxj9/Kzp07vbxECjLrM5B0bkT8taQ3Af8HGAXWAy9FxGckbQOWRsR2SauBu4Ff\nA1YADwLvqvep72Bg1r+8F0HnZNZnEBF/nRz+QnLNnwDXUO0TIvn32uR4PXBPRLwWEceAo8DaNOph\nZu0pl8tsLJXYWCqdtv5P2sbH9yaBYDNQDQpzLQLrvlSCgaSzJD0OPA9UIuIwMBwRMwAR8TywLDl9\nOfBszctPJGVmlqG5oZ7rp6ero3yuu67jAcHyI5UO5Ij4ObBG0luBsqQCsDC/01K+Z2xsbP64UChQ\nKBRaq6SZLWrv+Dh7ZmfnB3UyO8ve8fFU0jb1+gG8F0F6KpUKlUqlrWukOpooIn4m6QHg/cCMpOGI\nmJF0IfBCctoJ4KKal61IyuqqDQZm1nsW9g0cOrR5vm/AG9enY+EX5V27djV9jTRGE10AnIqIlyUN\nAWVgF1ACfhwRexp0IH+AanpoGncgm2Wu3ozgNJaD9j7F3ddKB3IaLYN/CExIEtU+iK9ExP9O+hDu\nlXQDcBzYBBARhyXdCxwGTgE3+RPfLHsjIyNMTE7Or/Mz4XV+BoqXozCzjvIQ0u7z2kRmlrl6ncWe\nSNZdDgZmlim3AvLBwcDMMuXO4nzwqqVmZtYSr1pqZqnxRLLe5TSRmbWk2e0prXvcZ2BmXeGO4nxz\nn4GZdcViK452c+VTS4+DgdkA6tQHdhornzqYZCQicvuoVs/M0nTw4MEYHhqKfRD7IIaHhuLgwYOL\nnl8sbohiccP8eQcPHoyhoeGAfQH7YmhoOA4ePBgbisXYBxHJYx/EhmKxY3Wz+pLPzuY+b5t9QTcf\nDgZm6WvmA7vRh/7czxYGiXaDQbuvt6pWgoGHlpr1iXK5PL/I3JaUFpl7fd9Ade/68fG9jIyMzD9q\n1dvkfmJ0tO16WOc5GJj1gYXLT28+dKjh8tOd/MBud+VTB5PseGipWR/YWCpVO2yT5xPAgWKR/VNT\ndc9v1IpYOEcA6PoQ0k60cAaN5xmYDahmg0E9jeYOAJ5E1mMcDMwGVBq7lHmRuf6R1U5nZpYx71Jm\n7XLLwGwANdqAxktM9AeniczsNAs7ZKFxp7AXmesPmQQDSSuALwPDwM+BuyLi85KWAn8MrASOAZsi\n4uXkNTuAG4DXgFsjom4vl4OB9Yq8joCp15fw9suu4PHHP4b7BvpXVn0GrwEfj4gnJL0F+JakKeCj\nwIMR8RlJ24AdwHZJq4FNwCpgBfCgpHf5U996VTNj/Ltt7/g4e2Zn5z/2mZ3ltuPfzbJKllNtL1QX\nEc9HxBPJ8UngCNUP+WuojnAj+ffa5Hg9cE9EvBYRx4CjwNp262GWldoP3M3AntnZ+VZCHr1j5QqG\nhrZR/bOcSDag2dL0dbygXH9JdTSRpHcA7wW+AQxHxAxUA4akZclpy4E/r3nZiaTMzFJ2xbp1bJ2e\nnn++FfjPGzey4/3vr+kbaL6TOM+tIWtNasEgSRH9KdU+gJOSFqZ9WkoDjY2NzR8XCgUKhUKrVTTr\niDwvofDww4/xM0a5mScBOMnlPPzwY+zcubOtD+566ae94+MOBhmpVCpUKpW2rpFKMJB0NtVA8JWI\nuC8pnpE0HBEzki4EXkjKTwAX1bx8RVJWV20wMMuj/I/x/0ec5L8mxxPA97OsjHXAwi/Ku3btav4i\nzS5zWu9BdTTRHQvK9gDbkuNtwKeT49XA48AS4GLgeySjmupct+mlW80GVTP7DqTxu7zvQH7RwhLW\naQwtvRL4OvAU1VRQAJ8AHgHupdoKOE51aOlPk9fsAG4ETuGhpWZtW2zCWKfmDuR1OK150plZS/rh\nQ83rClktr01k1iSPijGrcsvABloaSz9nIQ/7Dlh+uWVgNgAW9g8cOrSZyckJJicn2po7YIPNLQMb\naGnsA9BJ9foz3D9gb8QtA7MmLTZHIOuO5XK5zIfXr+d3X30VgA8//DBfOXCgq3WwAdLsWNRuPvA8\nA8tIHsbRr1uzJvZBRPLYB7FuzZqOzR1o5ODBg7GhWIwNxaLnEvQIWphn4JaBWR15WG7h2PHn6paN\njIx0rX/Ao60Gh4OBWRM6lTqqd923rbyUrT/+KdW5mbCVc7hk5aVANb3VjQ/kPARF6w4HA7M66i0+\nd8u6dR35ltyob2D37k+yfv313PxqdVHfV5ecYPfuT7b1u8waajav1M0H7jOwDC3MlW8oFk/L4W8o\nFtv+PY36BubqsHC9oW7KQ9+JNQ/3GZilZ2EqplMb1jTqG6hXh27L/4qslhYHA7MzlNa+BQtnDy/W\nN5AHWQck6462t700y1I3t16c+5Z8oFjkQLH4hv0F9eo2N3t4eno909Prue66zWzceDV/s+TN3Mx7\nuJn38DdL3uy+Aeu+ZvNK3XzgPgNbRJ7z2Y3qVixuSOYHzHUR7JvvD8iyb8D6C+4zsEGS52GPjeoG\n59U936kYy5qDgVkHvPTii3XLduwe5dChzXPdDgwNbWN0dKLLtTM7nYOB9aw8b0T/M85lK+fwuk5h\nzu3q7GGzZnjVUutpWS8mV1uP2hFC4+N7mZ6+mLfwJAAnuZxi8fteWdS6wttemmWg3v7DO3fewu23\nf8GbzVgmWgkGqQwtlfRFSTOSnqwpWyppStIzksqSzq/52Q5JRyUdkVRKow5mWRkf35t86G8GqkHh\n4YcfY3Kyus9AsXjAgcByL615Bl8CFv6fvh14MCLeDTwE7ACQtBrYBKwCrgbulNRUBDPLSrlcplTa\nSKm08Q3nNYyMjDA1tZ+pqf0OBJZ7qaWJJK0E7o+Iy5PnTwPrImJG0oVAJSIuk7Sd6hjYPcl5XwXG\nIuKbda7pNJHNy7p/oF46aHKyOhLI+w9bnuRtp7NlETEDEBHPS1qWlC8H/rzmvBNJmVlDeVhX//Xp\noOogpvHxvUxN7fcIIet53Rxa2tJX/LGxsfnjQqFAoVBIqTrWS7KYYLZwhNBiPGnMslSpVKhUKm1d\no5PBYEbScE2a6IWk/ARwUc15K5KyumqDgVm3LEwJHTq0mZ07b+HQoW25mDCWdcrM8mXhF+Vdu3Y1\nf5Fm169o9ADeATxV83wPsC053gZ8OjleDTwOLAEuBr5H0ndR55oprNJh/aCT6xDVWxcoz2sI5XlN\nJssHWlibKK1A8EfAD4C/Bf4v8FFgKfAg8AwwBfxSzfk7kiBwBCgtct1O/vcaaL24yXkn6txoc/lG\nwSAPOrXJjvWPVoJBKmmiiPjNBj/69Qbn7wZ2p/G7rXl56IxtRSfy8o06hUdHt3gNIRso3s9gANV2\nxm4G9syvqJkPndqjoNk5AnmdNLZldLS6DhMwQXVNpi05WZPJelizTYluPnCaqCPynGboVD68UTqo\nUXne9WKaz7qHrPoMOvUYpGDQzT/uPHdAphWoFnb0LtYH0O1OYX+QW6e1Egy8hHUOdDuH3++bnNcb\nFnrZZZc1PL+bcwR6tb/GBkCz0aObDwakZZDntE23NdtqOdNhoWvWXJmLdJDfa+sG3DKwXtdMq6Ve\nC2BuraCFLrhgmMnJT3rJCLNGmo0e3XwwIC2DPOfw8+RM+wHy3Cns99q6AbcMelO/5/DT0Ew/QJ63\nlvR7bXnlnc4sdxYuEDcyMkKptJHp6fUwv1TdBGvW3MXTT3/PS0ebLZC3JazNmuZ+ALNsuGVgmTnT\nFkCxeIDR0S3eQMbsDLllYJk706WVm2kBQL77Acz6QrM9zt18MCCjifrFYiNl+mEkkFmvwKOJLEuN\ndiMDuj4SyJu/mDXHwcBatjDn/9KLL552zksvvlh3mWi4i6Gh+ruGtbs8hJd8MGueg4G1pF7O/+1v\nX8ZWzgFOAbCVc7iEc7mgzus7ORIoi/2SzXqdg4G9oXqjfup92//JT36Hn7GVm3kSgJNczgUXfH/R\njWLO45Wu34+Znc7BYEDV+4CvVw6n5/sbjfpZufJCZmf/Bydrhn/OfeNf2A8AdCyVs2V0lM2HDs3l\no6obwXjzF7PFNdvj3M0HHk3Utnqrejaz0cuaNeuaGvVzpnsDdHr1Tu8ZYIOMXhpNJOmDwGepbr35\nxYjYk1Vdek273+ob7ftbPX59+fHjv1O3DouN+slDbr6bexSY9YNMgoGks4D/BlwF/AB4VNJ9EfF0\nFvXptnof5mmkbeoN32z0oX+mqqmf9Ef9OJVjljPNNiXSeAD/BPhqzfPtwLY656XdespcvfTKpz71\nqbbTNvUmcb3tbZc0neJpJ/XTyn8Lp3LM0kcPpYmWA8/WPH8OWJtRXbqqmqL5EG/hbgBOzn6IO+74\nUttpm3oafatfLMXTzdSPUzlm+ZH70URjY2Pzx4VCgUKhkFld0vDiizO8lfv5/PxY/AqvvlpvJH59\ni6VtFg7f3L27Wt7ow73eB3G98n6Zzdsv92G2UKVSoVKptHeRZpsSaTyopokO1jwfmDTRujVrThtF\n895LLkklbdOJdE6/7MzVL/dhdiZoIU2UVTB4E/A9YCWwBHgCWFXnvE78d8pUoyGV3fyAT6O+zcq6\nf8Ab0dsgaSUYZJImioi/k3QzMMXfDy09kkVduq3RKJpm0ja9xmsFmfWAZqNHNx/0YcsgIvtvyc1I\nI72Sh2/lThPZIKFXWgaDLo1VObvVEdovG7j3y32YdYq3vewxC1Mu24aGUk25dCLQdLrOZvZ6rWx7\n6WDQYzaWSqyfnq7ZIRgOFIvsn5pq+9qd/ND2sE6z7vEeyNaWTu4D0A8d4Wb9zMGgx3hNHzPrBKeJ\nelCnUi7O7Zv1B/cZWNuc2zfrfQ4GZmbWUjA4q1OVse4rl8tsLJXYWCpRLpezro6Z9RC3DPqE8/1m\nNsctgwFWOyx0M7AnGRbaSW6JmPUPDy21lnjxObP+4mDQJ7o9/6CTE9TMrPscDPqEF2Izs3a4A7nP\neYKa2eDxPAN7nV5c4dTM2udgYK/TyRVOzSy/PLTUzMxa0lYwkPRvJP2lpL+TdMWCn+2QdFTSEUml\nmvIrJD0p6buSPtvO77fFbRkdraaGqLYKtg0NscUrnJpZHe22DJ4CrgMeri2UtArYBKwCrgbulDTX\nZPl94MaIuBS4VJITzYtoZ2LX3AijA8UiB4pFd/CaWUOp9BlI+howGhGPJc+3U92QeU/y/KvAGHAc\neCgiVifl1wPrIuK3Glx3oPsMPGLHzFqRpz6D5cCzNc9PJGXLgedqyp9LyqyOLJaYMLPB9IaTziRN\nA8O1RUAAOyPi/k5VzMzMuucNg0FEFFu47gngoprnK5KyRuUNjY2NzR8XCgUKhUIL1elN3uLSzM5E\npVKhUqm0dY00+wz+U0R8K3m+Grgb+ADVNNA08K6ICEnfALYCjwL/C/h8RBxscN2B7jMAT+wys+Z1\nfdKZpGuBLwAXAD8FnoiIq5Of7QBuBE4Bt0bEVFL+PmAf8IvAAxFx6yLXH/hgYGbWLM9ANjOzXI0m\nMjOzHuJgYGZmDgZmZuZgYGZmOBiYmRkOBmZmhoOBmZnhYGBmZjgYmJkZDgZmZoaDgZmZ4WBgZmY4\nGJiZGQ4GZmaGg4GZmeFgYGZmOBiYmRkOBmZmhoOBmZnRZjCQ9BlJRyQ9IWm/pLfW/GyHpKPJz0s1\n5VdIelLSdyV9tp3fb2Zm6Wi3ZTAFvCci3gscBXYASFoNbAJWAVcDd0qa25z594EbI+JS4FJJI23W\noWdVKpWsq9Ax/Xxv4Pvrdf1+f61oKxhExIMR8fPk6TeAFcnxeuCeiHgtIo5RDRRrJV0InBcRjybn\nfRm4tp069LJ+/h+yn+8NfH+9rt/vrxVp9hncADyQHC8Hnq352YmkbDnwXE35c0mZmZll6Ow3OkHS\nNDBcWwQEsDMi7k/O2Qmcioj/2ZFamplZRyki2ruA9BHgY8C/jIi/Tcq2AxERe5LnB4HbgOPA1yJi\nVVJ+PbAuIn6rwbXbq5yZ2YCKCL3xWX/vDVsGi5H0QeC3gX8+FwgSB4C7Jf0e1TTQO4FHIiIkvSxp\nLfAo8O+Azze6frM3Y2ZmrWmrZSDpKLAEeCkp+kZE3JT8bAdwI3AKuDUippLy9wH7gF8EHoiIW1uu\ngJmZpaLtNJGZmfW+XM9AlnSbpOckPZY8Pph1ndIg6YOSnk4m3m3Luj5pk3RM0rclPS7pkazr0y5J\nX5Q0I+nJmrKlkqYkPSOpLOn8LOvYjgb31xd/e5JWSHpI0nckPSVpa1LeF+9fnfu7JSlv+v3LdctA\n0m3AKxFxR9Z1SYuks4DvAlcBP6Dad3J9RDydacVSJOmvgPdFxE+yrksaJP1T4CTw5Yi4PCnbA7wU\nEZ9JAvrSiNieZT1b1eD++uJvL5nbdGFEPCHpLcC3gGuAj9IH798i9/dvafL9y3XLINFvnchrgaMR\ncTwiTgH3UH3z+onojf+3zkhEHAIWBrZrgInkeIIenjzZ4P6gD/72IuL5iHgiOT4JHKE6ObYv3r8G\n9zc3d6up968X/mBvTtY++oNebcotsHBCXj9OvAtgWtKjkj6WdWU6ZFlEzED1DxJYlnF9OqGv/vYk\nvQN4L9XVEob77f2rub9vJkVNvX+ZBwNJ08nCdXOPp5J//zVwJ/CrydpHzwM93WQdIFdGxBXAbwD/\nMUlD9Lv85ltb01d/e0kK5U+pjmw8yenvV0+/f3Xur+n3r615BmmIiOIZnnoXcH8n69IlJ4BfqXm+\nIinrGxHxw+TfH0mapJoaO5RtrVI3I2k4ImaSvO0LWVcoTRHxo5qnPf23J+lsqh+UX4mI+5Livnn/\n6t1fK+9f5i2DxSRv0pwNwF9mVZcUPQq8U9JKSUuA66lO0usLks5NvqUg6c1Aif5438Trc7AHgI8k\nx5uB+xa+oMe87v767G/vD4HDEfG5mrJ+ev9Ou79W3r+8jyb6MtUc2M+BY8B/mMvz9bJkmNfnqAbj\nL0bEpzOuUmokXQxMUm12nw3c3ev3J+mPgALwD4AZqkur/BnwJ8BFVJdZ2RQRP82qju1ocH//gj74\n25N0JfB14Cmq/08G8AngEeBeevz9W+T+fpMm379cBwMzM+uOXKeJzMysOxwMzMzMwcDMzBwMzMwM\nBwMzM8PBwMzMcDAwMzMcDMzMDPj/xbl41yijcOQAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "y_noise = 200*np.random.rand(40,1) - 100\n", "y_n = y + y_noise\n", "poly = PolynomialFeatures(5, include_bias = False)\n", "x_predict = np.arange(0,30,0.5)[:, None] \n", "x_p = poly.fit_transform(x_predict)\n", "x_5 = poly.fit_transform(x)\n", "lr.fit(x_5,y_n)\n", "predict_y_n = lr.predict(x_p)\n", "plt.scatter(x, y, c ='blue')\n", "plt.scatter(x, y_n, c ='red')\n" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPX59/H3HSASFRBEg0oFJCDgRtispZUoTZD6FAUU\npVpw6eqGGhEQf8JPQQWLC33EFQFb69IqFR8gAbTRi98PV0CQRXEBBQGRigoFk5D7+WNOYMgCycwk\nM5N8Xtc1FzPfOXPmzkk49/mux9wdERGp31LiHYCIiMSfkoGIiCgZiIiIkoGIiKBkICIiKBmIiAjV\nSAZmNt3MtprZigreyzWzEjNrEVY2xszWmdkaM8sJK+9mZivM7CMzezD6H0FERKJVnZrBDKBf2UIz\naw1kAxvCyjoDQ4DOQH9gmplZ8PYjwNXu3hHoaGbl9ikiIrWrysnA3RcD31Tw1gPAyDJlFwDPuXux\nu68H1gG9zKwV0MTd3wm2exq4sNpRi4hITEXVZ2BmA4Av3H1lmbdOAL4Ie70pKDsB2BhWvjEoExGR\nOGoY6QfNLA24jVATkYiIJLGIkwHQHmgLvB/0B7QGlppZL0I1gRPDtm0dlG0CflRBeYXMTAsniYhE\nwN3t0FvtV91mIgseuPsH7t7K3U9y93aEmnwy3f0rYA5wiZmlmlk7IAN42923AN+aWa8ggQwDXj7E\nD5Twj3HjxsU9hroQo+JUnIn+SJY4I1GdoaV/A/6X0Aigz83syrLnbfYnitXAC8BqYB5wje+P8Fpg\nOvARsM7d8yKKXEREYqbKzUTu/qtDvH9Smdf3APdUsN17wGlV/V4REal5moEcA1lZWfEO4ZCSIUZQ\nnLGmOGMrWeKMhEXavlQbzMwTOT4RkURkZngNdyCLiEgdpGQgIiJKBiIiomQgIiIoGYiICEoGIiKC\nkoGIiKBkICIiKBmIiAhKBiIigpKBiIigZCAiIigZiIgISgYiIkJ090AWEUkIhYWFfPjhh3zwwQes\nWrWKtWvX0rVrV4YOHUr79u3jHV5SUM1ARJJWUVERU6ZMIT09nSFDhvDSSy/RqFEjBg0axJYtW/jJ\nT35Cr169uP/++9mxY0e8w01ourmNiCSlJUuW8Ic//IH09HQefvhhOnToUG6b4uJi/vWvfzF9+nQ+\n+OADFi5cyHHHHReHaGtXjd7cxsymm9lWM1sRVjbZzNaY2XIze9HMmoa9N8bM1gXv54SVdzOzFWb2\nkZk9WJ1gRUT+85//8Pvf/56LLrqI2267jfz8/AoTAUDDhg3Jzs7m2WefZejQoZx99tls2LChliNO\nDtVpJpoB9CtTtgA4xd27AuuAMQBm1gUYAnQG+gPTzKw0Sz0CXO3uHYGOZlZ2nyIiFdq1axfnn38+\n33//PatXr+aSSy5h/6mlcmbG2LFjue666zj77LP56KOPaiHa5FLlZODui4FvypQtcveS4OWbQOvg\n+QDgOXcvdvf1hBJFLzNrBTRx93eC7Z4GLowifhGpJ3bu3En//v1p164df/nLX2jWrFm19zFixAjG\njRtHVlYWK1eurIEok1csRxNdBTwbPD8BWBL23qagrBjYGFa+MSgXEanU999/T//+/encuTOPPfYY\nKSmRj3256qqraNCgARdffDErVqwgNTU1hpEmr5iMJjKzsUCRuz97yI1FRKrhu+++47zzzuPUU0+N\nOhGUGjZsGB06dGDKlCkxiLBuiLpmYGZXAL8Azg0r3gT8KOx166CssvJKjR8/ft/zrKwssrKyoglX\nRJKIuzNs2DA6d+7MtGnTYpIIINSHMHXqVHr27MnQoUNp27ZtTPYbLwUFBRQUFES1j2oNLTWztsAr\n7n5a8Po8YApwtrtvD9uuC/AMcCahZqCFQAd3dzN7E7gBeAeYC0x197xKvk9DS0XqscmTJ/Piiy/y\nxhtvcNhhh8V8/xMnTuStt95izpw5Md93PEUytLTKycDM/gZkAUcDW4FxwG1AKlCaCN5092uC7ccA\nVwNFwAh3XxCUdwdmAo2Bee4+4iDfqWQgUk8VFBRw6aWX8vbbb3PiiSfWyHf88MMPnH766dx3330M\nGDCgRr4jHmo0GcSDkoFI/bR582a6d+/OzJkzycnJOfQHovDqq69y9dVXs2rVKo444oga/a7aUqOT\nzkREakNRURFDhgzhD3/4Q40nAoC+ffvyk5/8hIkTJ9b4dyUy1QxEJKGMHTuWpUuXMnfu3Jh1GB/K\nxo0bOe2009iwYQNNmzY99AcSnGoGIpLUli1bxpNPPsmMGTNqLREAtG7dmp///Of89a9/rbXvTDRK\nBiKSEIqKirjqqquYPHkyrVq1qvXv/+Mf/8i0adOor60RSgYikhCmTJnCsccey7Bhw+Ly/eeccw7F\nxcUsXrw4Lt8fb+ozEJG4+/DDD+nduzfvvvtuXCeATZ06lSVLlvDss8m9mIKGlopI0ikpKaFPnz4M\nGTKE66+/Pq6x7Nixg3bt2rF27VrS09PjGks01IEsIknn0UcfpaSkhGuuuSbeoXDUUUdx0UUXMX36\n9HiHUutUMxCRuNm2bRtdunTh9ddfp0uXLvEOB4ClS5cycOBAPv30Uxo0aBDvcCKimoGIJJU77riD\nyy67LGESAUC3bt047rjjmDdvXrxDqVWqGYhIXKxYsYLs7GzWrl1L8+bN4x3OAZ5++mmeffZZ5s+f\nH+9QIqKagYgkBXfnxhtvZNy4cQmXCACGDBnCkUceSVFRUbxDqTWqGYhIrXvppZcYN24cy5Yto2HD\nWN5wUUBDS0UkCezZs4cuXbrwxBNP0Ldv33iHUyepmUhEEt4DDzzAGWecoUSQYFQzEJFa89VXX9Gl\nSxfeeust2rdvH+9w6iw1E4lIQrvxxhspKSlh6tSp8Q6lTlMyEJGEtX79erp3787q1auTeqmHZKBk\nICIJ64orruDEE0/kzjvvjHcodV6NdiCb2XQz22pmK8LKmpvZAjP70MzyzaxZ2HtjzGydma0xs5yw\n8m5mtsLMPjKzB6sTrIgkpw8++ID58+eTm5sb71CkEtUZTTQD6FembDSwyN1PBl4DxgCYWRdgCNAZ\n6A9MM7PSLPUIcLW7dwQ6mlnZfYpIHXP77bczatQomjVrduiNJS6qnAzcfTHwTZniC4BZwfNZwIXB\n8wHAc+5e7O7rgXVALzNrBTRx93eC7Z4O+4yI1EFLlixh6dKlCbEqqVQu2nkGx7r7VgB33wIcG5Sf\nAHwRtt2moOwEYGNY+cagTETqIHdn9OjRjB8/nsaNG8c7HDmIWM8Dj3lv7/jx4/c9z8rKIisrK9Zf\nISI1ZNGiRWzZsiVut7KsLwoKCigoKIhqH9UaTWRmbYBX3P304PUaIMvdtwZNQP9y985mNhpwd58U\nbJcHjAM2lG4TlF8K9HH3P1byfRpNJJKk3J3evXtz3XXX8atf/Sre4dQrtbEchQWPUnOAK4Lnw4GX\nw8ovNbNUM2sHZABvB01J35pZr6BDeVjYZ0SkDlmwYAE7duzgkksuiXcoUgVVbiYys78BWcDRZvY5\noSv9e4G/m9lVhK76hwC4+2ozewFYDRQB14Rd4l8LzAQaA/PcPS82P4qIJAp354477mDcuHFJe7ew\n+kaTzkQk5ubNm8ett97KihUrSEnRepi1TauWikjcuTvjxo1j/PjxSgRJRL8pEYmpuXPnUlhYyKBB\ng+IdilSDkoGIxExprWDcuHGqFSQZ/bZEJGZeeeUV9u7dy4UXamGBZKNkICIx4e7cdddd3HHHHaoV\nJCH9xkQkJvLy8tizZ49qBUlKyUBEolZaK7j99ttVK0hS+q2JSNRee+01/v3vf3PRRRfFOxSJkJKB\niETtzjvvZOzYsZptnMSUDEQkKq+//jqbNm1i6NCh8Q5FoqBkICJRueuuu7jtttto2DDWK+JLbVIy\nEJGILVmyhE8++YRf//rX8Q5FoqRkICIRmzBhAqNGjaJRo0bxDkWipFVLRSQiS5cuZcCAAXzyyScc\ndthh8Q5HwmjVUhGpNRMnTmTkyJFKBHWEagYiUm2rVq2ib9++fPrppxx++OHxDkfKUM1ARGrF3Xff\nzY033qhEUIeoZiAi1fLxxx9z1lln8cknn9C0adN4hyMVUM1ARGrcvffey7XXXqtEUMfEpGZgZmOA\ny4G9wErgSuAI4HmgDbAeGOLu34ZtfxVQDIxw9wWV7Fc1A5EE8vnnn5OZmcm6deto0aJFvMORSsSl\nZmBmbYDfApnufjrQEBgKjAYWufvJwGvAmGD7LsAQoDPQH5hmZtUKWkTiY/LkyfzmN79RIqiDYjF/\n/DugEDjCzEqANGAToZN/n2CbWUABoQQxAHjO3YuB9Wa2DugFvBWDWESkhmzevJm//e1vrF27Nt6h\nSA2Iumbg7t8AU4DPCSWBb919EZDu7luDbbYAxwYfOQH4ImwXm4IyEUlgf/rTnxg2bBjHHnvsoTeW\npBN1zcDMTgJuItQ38C3wdzO7DCjb2B9R4//48eP3Pc/KyiIrKyuiOEUkctu2bWPGjBmsWLEi3qFI\nBQoKCigoKIhqH1F3IJvZECDb3X8bvP418GPgXCDL3beaWSvgX+7e2cxGA+7uk4Lt84Bx7l6umUgd\nyCKJYezYsWzfvp1HH3003qFIFUTSgRyLZHAG8FegJ/ADMAN4BzgR+Le7TzKzUUBzdx8ddCA/A5xJ\nqHloIdChorO+koFI/H3zzTdkZGTw3nvv0bZt23iHI1UQSTKIupnI3d83s6eB9wgNLV0GPA40AV4w\ns6uADYRGEOHuq83sBWA1UARcozO+SOL685//zIABA5QI6jjNQBaRSn3//fecdNJJ/M///A8dO3aM\ndzhSRZqBLCIx9cgjj/Dzn/9ciaAeUM1ARCq0a9cu2rdvz6JFizj11FPjHY5Ug2oGIhIzjz32GD/9\n6U+VCOoJ1QxEpJzdu3dz0kknMX/+fLp27RrvcKSaVDMQkZh4/PHHOfPMM5UI6hHVDETkAHv27KF9\n+/a88sordOvWLd7hSARUMxCRqE2fPp1u3bopEdQzqhmIyD4//PADGRkZvPTSS/Ts2TPe4UiEVDMQ\nkajMmDGD0047TYmgHlLNQEQAKCwspEOHDjz//PP8+Mc/jnc4EgXVDEQkYk899RSdOnVSIqinVDMQ\nEfbs2UOHDh34xz/+wZlnnhnvcCRKqhmISESefPJJzjjjDCWCekw1A5F6bvfu3WRkZDBnzhy6d+8e\n73AkBlQzEJFqe+yxx+jZs6cSQT2nmoFIPbZr1y4yMjK0BlEdo5qBiFTLI488Qu/evZUIRDUDkfpq\n586dtG/fnldffVXLVNcxqhmISJU9+OCDnHvuuUoEAsSoZmBmzYAngVOBEuAq4CPgeaANsB4Y4u7f\nBtuPCbYpBka4+4JK9quagUgN2L59OyeffDJLliyhQ4cO8Q5HYiyeNYOHgHnu3hk4A1gLjAYWufvJ\nwGvAmCDILsAQoDPQH5hmZtUKWkSiM2nSJAYPHqxEIPtEXTMws6bAMndvX6Z8LdDH3beaWSugwN07\nmdlowN19UrDdfGC8u79Vwb5VMxCJsU2bNnH66aezYsUKTjjhhHiHIzUgXjWDdsDXZjbDzJaa2eNm\ndjiQ7u5bAdx9C3BssP0JwBdhn98UlIlILbjzzju5+uqrlQjkAA1jtI9uwLXu/q6ZPUCoiajsJX1E\nl/jjx4/f9zwrK4usrKzIohQR1q1bx4svvsiHH34Y71AkhgoKCigoKIhqH7FoJkoHlrj7ScHrnxJK\nBu2BrLBmon+5e+cKmonygHFqJhKpeUOHDuXUU09l7Nix8Q5FalBcmomCpqAvzKxjUNQXWAXMAa4I\nyoYDLwfP5wCXmlmqmbUDMoC3o41DRA5u+fLlFBQUMGLEiHiHIgkoFs1EADcAz5hZI+BT4EqgAfCC\nmV0FbCA0ggh3X21mLwCrgSLgGl3+i9Qsd2fkyJHcfvvtHHnkkfEORxKQZiCL1APz58/npptuYuXK\nlTRq1Cje4UgN0wxkESmnuLiYW265hcmTJysRSKWUDETquBkzZnDMMcfwy1/+Mt6hSAJTM5FIHbZz\n5046duzInDlz6NGjR7zDkVqiZiIROcB9993Hueeeq0Qgh6SagUgdVbrsxNKlS2nTpk28w5FaFEnN\nQMlApI664ooraNWqFffee2+8Q5FaFkkyiNU8AxFJIEuWLGHhwoWsXbs23qFIklCfgUgds3fvXq6/\n/nomTZpEkyZN4h2OJAklA5E65qmnnqJx48Zcdtll8Q5Fkoj6DETqkG+++YbOnTszf/58MjMz4x2O\nxIk6kEXqueuvv57i4mIeeeSReIcicaQOZJF6bMWKFTz//POsWbMm3qFIElKfgUgdMH/+fM7t04cT\njj6ad999N97hSBJSMhBJcvn5+VxywQU027GDG9auZfjAgeTn58c7LEkySgYiSW7q3XdDURGzCd1I\nZNLu3Tw+ZUq8w5Iko2QgkuRWfPghfYDT4x2IJDV1IIsksblz57I3JYW3Gzdm1p49AIxKS2NWbm6c\nI5Nko5qBSILLz89ncE4Og3Ny9vUF5Ofnc8G55zLkoou44YYbePqf/2ROdjZzsrOZNXs2/fr1i3PU\nkmw0z0AkgeXn5zN84EAm7d4NhK76rx87lj9PnEjX3bvZCXyclqYEIAeI66QzM0sB3gU2uvsAM2sO\nPA+0AdYDQ9z922DbMcBVQDEwwt0XVLJPJQOpd/Lz8/d1AG//+muuXLaM4cF7s4C7WrTgkn//m6eA\nlcBcYE52Ni8uqPC/kdRD8Z50NgJYDTQNXo8GFrn7ZDMbBYwBRptZF2AI0BloDSwysw4664uUrwnk\nppRvyS1x5wngCaBl7YYndVhM+gzMrDXwC+DJsOILCF3IEPx7YfB8APCcuxe7+3pgHdArFnGIJLvH\np0xh0u7dDAeGA1eUlHBTSgqzCP0nGpWWRsuMDL5v0IAdYWW/U4exRClWNYMHgJFAs7CydHffCuDu\nW8zs2KD8BGBJ2HabgjIRKeM0oMsZZzCnZagOcN3ZZ/Pkk0/yzAsv8MyjjwIwKzdX/QUStaiTgZmd\nD2x19+VmlnWQTSNqBho/fvy+51lZWWRlHewrRJLb73JzGb54MYR1GM+65x769evH9u3bOf300/nr\nX//KOeecw6BBg/Z9Lj8/nylTHgcgN/d3Sg71TEFBAQUFBVHtI+oOZDO7G7icUGdwGtAEmA30ALLc\nfauZtQL+5e6dzWw04O4+Kfh8HjDO3d+qYN/qSpB6J7wD+XdhV/2XXnoprVq14sEHHyy3/cCBw9m9\nexIAaWmjmD17lhJCPRb3JazNrA+QG4wmmgxsd/dJQQdyc3cv7UB+BjiTUPPQQqDCDmQlA5GQmTNn\nMnnyZN577z3S0tIOeC8nZzALFw6AsDFH2dlzWLDgxVqPUxJDvEcTlXUv8IKZXQVsIDSCCHdfbWYv\nEBp5VARcozO+SOVWrVrFLbfcQkFBwb5EEN4s9PXX2+MZntQRmnQmksB27txJr169GDlyJFdeeSVQ\nvlkoNfVGoBGFhfcBaiaSBGgmijUlA6nP3J3hw4ezefNmzELTd3Jzf8eUKY+XaxbKzHyCli3T922j\nRFC/JVozkYhEYcaMGbz++ut89dUe9uyZDMDixcPp1KlTuW1btkxXH4FERclAJAGtWLGCUaNG0aFD\nVz7//HJKawGhEadPkJY2qnT0KWlpo8jNnVXZrkSqRKuWiiSY559/njPPPItWrdqxZ09Rufdbtkxn\n9uzQiKHs7DnqH5CYUJ+BSAJ55ZVXuPDCwZSUnAcMVuewRER9BiJJqHSYqLuzdOn/UlJyBvBPIIXC\nQoLO4TkA5OYqEUjNUDIQqYbKZgdHs7/9w0TnAV8DdxLegqvOYakNSgYiVVR2eenhixdHfVOZKVMe\nDxJBM2AxcDUpKaMpKTkMUOew1B4lA5EqCl9eGoDdu3l8ypRqJ4Pys4fXAE8Rqhms4owzVqlZSGqd\nkoFILSo7e7hhw+uAN2hMGxpyBYWpm7jnnueUAKTWKRmIVFGFy0tX4aYyZWsCoUQwHFhNcXEJhwPT\nWA/ASFJrJniRQ9DQUpFqqG4HctmaQEpKLiUlU4CzgbM5DOcxNh1wj2Pdz1iipaGlIjWsX79+B00A\nZW8ys7+DOHS6LylZidkNuKcC59E45QUoqYXARQ5ByUAkRsrWAipeR6g5DRvuISPjJFq3/g99+tzB\nqIkTq930JBJrSgYiUai8P6CidYQ+A+7i5ptHcu+99+7bR48ePfY1Pel+xhIvSgYiEaqoP6Cs0DpC\n/8Xtt9/N+++/xW23/RdnnXUWg3NygP39DkoAEm9KBiIRqqg/ICXlJkqCPoDSCWPbtm1j/frVzJ8/\nl+Li4phPXBOJBSUDkWo4+O0mT6Ndu+Mo/mYkAFffeD2vvvoqL774IgUFBZxyyikMzsmJycQ1kVhT\nMhCpoopuN5maOpLCQva9/u6L/3BfYSG7gRHjx9Pp1FN5++23Ofroo+MXuEgVRJ0MzKw18DSQTmiQ\n3BPuPtXMmgPPA22A9cAQd/82+MwY4CqgGBjh7hpULQlp4sSJPHH/QwA0bN7qgGahsiuKFn7djiuX\nLaMbcAnw05ISmh5zzAGJINKJayI1LRY1g2LgZndfbmZHAu+Z2QLgSmCRu082s1HAGGC0mXUBhgCd\ngdbAIjProNllkmgmTpzI5NtvZ2rw+oZ/bwNWHrBN+Iqig7KzyQduAe4DDJiTcuD9o/r168es2bM1\nekgSTtTJwN23AFuC5zvNbA2hk/wFQJ9gs1lAATAaGAA85+7FwHozWwf0At6KNhaRaIX3Cax68zWm\nsv+28wDX8xDfcxpw4IqiX375Jeu3b+eDlBQmlpRg7L/qr2jWshKAJJqY3vbSzNoCXYE3gXR33wr7\nEsaxwWYnAF+EfWxTUCb1WH5+PoNzchick0N+fn6t7SM/P5+cnMHk5Axm4sSJDBw4nIULB7Bw4QC+\n/35nue2bNDnigNtNZmdnM2PGDDIzM7ngggv45yuv8GZ2NnOys5k1ezYAwwcOZMDChQxYuJDhAwdG\n/POJ1KSYdSAHTUT/INQHsNPMyjb7RNQMNH78+H3Ps7KyyMrKijRESVCxuE9AVfZR9godOKBD+NVX\nS9cNCtUFvucVbmD/TWVuAG4dNZKxY8cC8N5779G7d2/cnblz59KjRw8AfvGLX+z7jEYPSW0oKCig\noKAgup24e9QPQkklj1AiKC1bQ6h2ANAKWBM8Hw2MCtsuDzizkv26JKe8vDwflJ3tg7KzPS8v76Db\nDsrO9pngHjxmgg/Kzq7W91W0j96ZmftimDBhgh+Tmuozg/eOSU31zMzeDjN9/8d+XOb1TG/fvou3\naXGMt2lxjE+YMMHd3b/66iv//e9/7+np6T59+nTfu3dvjf5sItUVnDurdR6PVc3gKWC1uz8UVjYH\nuAIoHX7xclj5M2b2AKHmoQzg7RjFIQmgJu4IVl0rgdXvv89vgxlgNy9axP3u+6/QCwu57eNVZT7V\nu9yksYcf3n9zma+++opbb72VJ598kssvv5w1a9bQvHnzg8ah0UOSNKqbPco+gN7AXmA5sAxYCpwH\ntAAWAR8CC4Cjwj4zBviYUO0h5yD7rsnkKTWkKlfD4TWHCRMmeHpa2r6r9vS0tEPWJsrKy8s7YB9H\np6QcEMOPg/LwmI5v0szT0tKD2sBMT0tL9wkTJnh29iDPzh60L4Yvv/zSc3NzvXnz5n7NNdf4559/\nXu3YqlpLEokFIqgZxKSZqKYeSgaJ6VAnt0Mlg7In7vS0NJ8wYULUJ8zwuLq2b39ADLngTbF939mU\nRp6Z2dvz8vLKnfzd3ffu3euLFi3yiy66yI866ii/7rrr/IsvvogoLpHapmQgNa6iE3nZk/ehtqmp\ndvTwE3v79qd5UxqFnfwbOjT2IznFj+QUT009qsKk89RTT/kpGRl+RFqat23b1h9++GHfsWPHQb9T\nV/2SaCJJBlqOQso52N28qnJT+NqaWBU+J6BPn25MnPjnA+8oxg1cxwoAdnI6mZlv0rJlOgC5uaGY\n3Z2VK1cye/ZsZs2axYbPPuNsYCDw+JYttG/fnmbNmlX6/fHuGxGJmepmj9p8oJpBrauNq/qq1C6q\nso/w9v6UlKPLjATK9ZSU5gf0B+Tl5fnevXt99erVPm3aNL/44ov9mGOO8ZNOOslvuukmP7tHD3+q\nGj+bRgpJokI1A4nWoa78YzE6JtKaw8FuJFNS8miZrU/jlFM6kpY2k127vicj48dMmDCBiy++mKOP\nPpqsrCzOP/987rvvPtq0aQOE5gTEdBamSBJRMqhjqnvD9uqKVRNQdZdkKLtiqNnNwHbgXUKT2A8H\n/gA8BmwDPuHjjxvTqVMnTj+9E926dSMzM5OuXbtWuoJodROdho1KXWKhGkViMjNP5PgSTdk27FFp\naVHP5I1kH9X1ww8/sG3bNrZt28bXX3/N9u3bWbx4MXl5r1FUVEhGRhvef38V27cfQWj5t23ATkKT\n2k8EjqJBg48455zebN78DY0bH87NN/+BoUOHYmbViqW6ybSmk69IJMwMd6/WH7+SQR0yOCcntP5N\n8HoWMCc7mxcXVG+F8IqWbYjkhOfubN++nfXr17N+/Xo2bNjAxo0b2bx5M5s3b+bLL79k69at7Nmz\nhyZNmlBYuJfU1FSOPz6dNWs+Ye/evsCRNGo0jx/96EQ+/fRC4NeElrn6J5mZT4Z1CP9OJ2KRQCTJ\nQM1EUk54E05VRsy4O59++inLly9n7dq1BzwaNWpE27Ztadu2LW3atKF169b07NmTjRs38vLLi2jd\nuhN9+57J3Xf/X3bvvh+AHTtyKSl5hNL+gKKiWTRr9gRpaU+we3cGAGlpo7nnnllRJwBd2YuEKBnU\nITXRhl1Rh/JDEyfy3XffsWTJEpYuXcry5ctp2rQpXbt2pUuXLvTt25drr72Wk08+ucLlGvLz8/nN\nb27e1/7/xhsHLhBXvjN4/43lSzuQc3Njkwg0NFQkRM1EdUysr3QH5+Tws4ULOQp4jdCqgt82asR5\n/ftz1lln0b17dzIzM2nZsuVB91P2jmGffJLL/jsFnEWo87f09S2kpDxFSckDQGiNoNmzoz/5V/Sz\nxaJZTSTR011SAAALzElEQVTRqJlIyo3SiSQ5FBUV8cYbbzBv3jze++gjZgM9gFOAksMO45+zZ9O/\nf/+D7iN8GOjxxzdh9qxZB7ljWNkF4v7K2LG5vP566HaSsagFiMghVHdiQm0+0KSzqFRncteuXbt8\n9uzZPmzYMG/RooX37NnT//u//9vffPNNnzt37iGXXJgwYcK+pZ6HDx9+wISwJjQsNzmrCQ0PuUBc\nTYvF5DeRRITWJpJwh5ohW1RU5PPnz/fLL7/cmzVr5uecc45PnTr1kKtyll3cbcKECd402P9M8Kbg\nkLtvNvCRNKtwxdDqnvxrYh0grS0kdVEkyUDNRPXQ8uXLmTFjBs899xxt27blsssu409/+hPp6ekV\nbn+wNYAWLx5O84Y/lLtX8HXksZM/AbCTnx/0jmGVfWfZ4a010dmr+xGLBKqbPWrzQR2vGVTlqrTs\nNtW5kg1vBnkYvEmjRp6RkeE/+tGP/I477vB169ZV+rnSq/YJEyYcYg2gmd6iYdohm4GGDx9e7o5h\nVYl7ZtB80yczU+sAiVQRaiZKHpEsBX1UauoBt26sShv3Qw895G2OP94bNWzoP/vZzzwvL8+Li4vL\nfU/VT/4V3xqybDPR8OHDI+4DqKh5q32LFge9raWaeET2UzJIIlVZ8bLsNhXdrauiq+Ndu3b59OnT\nvXv37t6qVSvPyDjF+/Q5f98JM7qTf8WrgYZ3IB/qyj+SY9M7MzPqxChSX0SSDNRnUId89tln3Hrr\nrbz88hyaNWtB//7ZrFr1BVu2jOTjj+Htt4czduz1B7T5v/rqoSZ8VW3YZ79+/Q7aB1AdFU6eu+ce\nYP+yGGd8/TVXLlt20PsqiEg1VDd7xOpB6D7Ja4GPgFGVbBP7lFmLDta+H2kz0dENG+57fUxqqt91\n113erdtP/ZhjjvPDDz/cGzY83GFy5e37LdpXcOWf60eS7UeS7TC43JV/vIZ9RnNrTZH6jGRpJgJS\ngI+BNkAjYDnQqYLtauI41YqqnOzLNq1UdD/esuP3U1OP8iM5xY+gk6ekpDo0cGjtcIWbtThk+375\nZDA4pu39tUVzBEQqF0kyiMtyFGb2Y2Ccu/cPXo8Ogp9UZjuPR3xVFT7kMjf3dwD7Xhd+/dkBzRiz\ngBmZmaS2bAeUH6KZmnoj0IjCwvuA0BIMZZt0UlJyKSm5EdgC/I3Qcs6/ByYGzw+9rEPZfTZL+R0P\nlRQm5ZIMWmROpGLJtBzFCcAXYa83Ar3iFEtEyt5s5fXXLyX8ZN7UfsOVZT6zbPkHfOcjgPJt9YWF\njxJ+It+9G+6//65g/5cAf6ek5AdgCqGR+iuAi4GTCSUCqGr7fo8ePcKS1imwbFkMj0zt0RwBkdhR\nB3KEpkx5/IDbLpY9mX/n93IDnwBFANxACt/5DRxsZc4DOUVFe4CZQC6h1YHOxewN3E8CXiU1dS0w\nksLC0CcqO/mX7detaIlq3a1LpH6LVzLYROgWVaVaB2XljB8/ft/zrKwssrKyajKualrJkeQAsJMd\nZV7Dd9zAdawIXn8BnBb22QOv4lNT1+J+C0VFG4ClmP0/Gjduwe7d71JcPB5oGTTz3BJ2sn+Od999\nlyfuHwnAb28ewdixY8ud/A8mVrexFJH4KSgooKCgILqdVLeTIRYPoAH7O5BTCXUgd65gu1j2qcRU\n6Xo8vwC/Fjwt6HydGdYR27Bhs32jclJTj/LU1GPKjdLJyvqld+36Ez///PO9ZcuWnpZ2hJ94YoZP\nnTrVS0pKKuxULhVpJ6rW4xGp20iWDmQAMzsPeIjQyKLp7n5vBdt4vOI7lME5OfyfhQvZCcwFFgId\ngauATsA7QH6XLhyRnoG7M2zYQL777jtmznye3bt3cfzxR/Pll1+yceNGunfvzvnnn88vf/lLOnXq\nVOX79kayHn887nEsIrUrmTqQcfc8Qr2fSSsFuD549AJ6AxuA14HVwNZPP8U/+4wGDRqwYcMHtGzZ\nkuOOa0mbNmfQo0cPevbsSZcuXWjYsPZ+DRXduUyTtUREHcgRKjtLdl1qKuuB+woL6UnoivulGr7i\nronbXIpI/aTbXkahomWWa3vce3XH2quZSKTui6SZSMmgHtJkLZG6TckgyekkLSKxoGSQxNR8IyKx\nomSQxCIZJioiUpFIkkFKTQUj5eXn5zM4J4fBOTnk5+fHOxwRkX00tLSWlG0GKntDdw0TFZF4UjNR\nLalKM5A6kEUkFpJqBrKUpyWZRSRelAxqiZqBRCSRqZmoFqkZSERqg4aWioiIhpaKiEhklAxERETJ\nQERElAxERAQlAxERQclARESIMhmY2WQzW2Nmy83sRTNrGvbeGDNbF7yfE1bezcxWmNlHZvZgNN8v\nIiKxEW3NYAFwirt3BdYBYwDMrAswBOgM9AemmVnpmNdHgKvdvSPQ0cySfuZVQUFBvEM4pGSIERRn\nrCnO2EqWOCMRVTJw90XuXhK8fBNoHTwfADzn7sXuvp5QouhlZq2AJu7+TrDd08CF0cSQCJLhDyQZ\nYgTFGWuKM7aSJc5IxLLP4CpgXvD8BOCLsPc2BWUnABvDyjcGZSIiEkeHXKjOzBYC6eFFgANj3f2V\nYJuxQJG7P1sjUYqISI2Kem0iM7sC+C1wrrv/EJSNBtzdJwWv84BxwAbgX+7eOSi/FOjj7n+sZN9a\nmEhEJAK1ej8DMzsPGAmcXZoIAnOAZ8zsAULNQBnA2+7uZvatmfUC3gGGAVMr2391fxgREYlMVDUD\nM1sHpALbg6I33f2a4L0xwNVAETDC3RcE5d2BmUBjYJ67j4g4ABERiYmEXsJaRERqR0LPQDazcWa2\n0cyWBo/z4h1TODM7z8zWBhPoRsU7nsqY2Xoze9/MlpnZ2/GOp5SZTTezrWa2IqysuZktMLMPzSzf\nzJrFM8YgporiTKi/TTNrbWavmdkqM1tpZjcE5Ql1PCuI8/qgPNGO52Fm9lbwf2aVmd0dlCfa8aws\nzmofz4SuGZjZOOB7d78/3rGUZWYpwEdAX+BLQn0gl7r72rgGVgEz+xTo7u7fxDuWcGb2U2An8LS7\nnx6UTQK2u/vkIME2d/fRCRhnQv1tBnN4Wrn7cjM7EngPuAC4kgQ6ngeJ8xIS6HgCmNnh7v4fM2sA\n/A+QS2gOVcIcz4PE+XOqeTwTumYQSNRO5F7AOnff4O5FwHOE/qgTkZGAv2t3XwyUTVAXALOC57NI\ngEmJlcQJCfS36e5b3H158HwnsIbQJNCEOp6VxFk61yhhjieAu/8neHoYof8/35BgxxMqjROqeTwT\n7gRRgeuCtY+ejHeVrIyyE+sSeQKdAwvN7B0z+228gzmEY919K4ROHMCxcY7nYBLyb9PM2gJdCa0K\nkJ6oxzMszreCooQ6nmaWYmbLgC1AgbuvJgGPZyVxQjWPZ9yTgZkttNDCdaWPlcG/vwSmAScFax9t\nARKmCplkert7N+AXwLVBs0eySNR2zIT82wyaXv5BaATfTsofv4Q4nhXEmXDH091L3D2TUA3rZ2aW\nRQIezzJxnm1mfYjgeEY1zyAW3D27ips+AbxSk7FU0ybgxLDXrYOyhOPum4N/t5nZbEJNXIvjG1Wl\ntppZurtvDdqXv4p3QBVx921hLxPib9PMGhI6wf7F3V8OihPueFYUZyIez1Lu/p2ZzQN6kIDHs1QQ\n51ygh7u/HvZWlY5n3GsGBxMc7FKDgA/iFUsF3gEyzKyNmaUClxKabJdQzOzw4CoMMzsCyCGxjqNx\nYNvmHOCK4Plw4OWyH4iTA+JM0L/Np4DV7v5QWFkiHs9ycSba8TSzlqVNK2aWBmQDy0iw41lJnMsj\nOZ6JPproaUJtiiXAeuD3pe11iSAYrvUQoaQ63d3vjXNI5ZhZO2A2oepsQ+CZRInTzP4GZAFHA1sJ\nLVnyT+DvwI8ILV8yxN13xCtGqDTOc0igv00z6w28Aawk9Lt24DbgbeAFEuR4HiTOX5FYx/M0Qh3E\npYMv/uLufzKzFiTW8awszmqfOxM6GYiISO1I6GYiERGpHUoGIiKiZCAiIkoGIiKCkoGIiKBkICIi\nKBmIiAhKBiIiAvx//0vQoE7U6RoAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.scatter(x, y, c ='blue')\n", "plt.scatter(x, y_n, c ='red')\n", "plt.plot(x_predict, predict_y_n, c ='black', label = '3rd Degree')" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([[ 9.20733482e+01, -2.59681982e+01, 3.02899345e+00,\n", " -1.50657823e-01, 2.80193672e-03]])" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lr.coef_" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[ 9.13742192e+01 -2.57459858e+01 3.00133841e+00 -1.49176871e-01\n", " 2.77335009e-03]]\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEACAYAAABGYoqtAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VPWd//HXZxJCglyKXOUiFeRqdQURi3VLbAW0bPFS\nf4q7YrReKwp96M8allbRxVLosgr+qq6ulahUFi9UbDVBxcHFqqBgQcEIWhDYElQUFcIt+fz+mJM4\nCRkgySRnJnk/H495cOab7znzmZNwPvO9nTF3R0REpCaRsAMQEZHUpSQhIiIJKUmIiEhCShIiIpKQ\nkoSIiCSkJCEiIgnVO0mYWT8zW2VmK4N/d5rZRDNrb2aLzazYzIrMrF3cPpPNbL2ZrTOzUXHlQ8xs\ntZl9YGb31Dc2ERGpH0vmOgkziwBbgNOAG4DP3H2mmd0KtHf3fDMbBMwDTgV6AC8Bfd3dzexN4AZ3\nX2FmzwOz3b0oaQGKiEitJLu76SzgQ3ffDJwLFATlBcB5wfZYYL67H3D3jcB6YJiZdQXauPuKoN6j\ncfuIiEgIkp0kLgb+EGx3cfcSAHffBnQOyrsDm+P22RqUdSfWCqmwJSgTEZGQJC1JmFkLYq2EJ4Oi\n6v1Yuv+HiEiayUzisc4B3nb3T4PnJWbWxd1Lgq6k7UH5VqBn3H49grJE5QcxMyUcEZE6cHerTf1k\ndjddAjwR93wRcHmwnQc8G1c+zsyyzOw44HhgedAltdPMhpmZAZfF7XMQd0/5x+233x56DIpTMSpO\nxVnxqIuktCTMrBWxQetr4opnAAvM7KfAJuAiAHdfa2YLgLXAfuB6/yb6CcBcIBt43t0LkxGfiIjU\nTVKShLvvBjpVK9tBLHHUVH86ML2G8reBE5MRk4iI1J9WXDeg3NzcsEM4IoozedIhRlCcyZYucdZF\nUhfTNRYz83SMW0QkTGaGhzhwLSIiTYyShIiIJKQkISIiCSlJiIhIQkoSIiKSkJKEiIgkpCQhIiIJ\nKUmIiEhCShIiIpKQkoSIiCSkJCEiIgkpSYiISEJKEiIikpCShIhILdx11128+uqrYYfRaJQkRERq\nYcGCBWRnZ4cdRqPR90mIiByh0tJSjj76aD7//PO0TBT6PgkRkQa0evVq+vfvn5YJoq6SkiTMrJ2Z\nPWlm68zsPTM7zczam9liMys2syIzaxdXf7KZrQ/qj4orH2Jmq83sAzO7JxmxiYgky9tvv80pp5wS\ndhiNKlktidnA8+4+EPgH4H0gH3jJ3fsDS4DJAGY2CLgIGAicA9xnZhXNn/uBK929H9DPzEYnKT4R\nkXpTkqgDM2sL/KO7PwLg7gfcfSdwLlAQVCsAzgu2xwLzg3obgfXAMDPrCrRx9xVBvUfj9hERCZ2S\nRN0cB3xqZo+Y2Uoze9DMWgFd3L0EwN23AZ2D+t2BzXH7bw3KugNb4sq3BGUiIqHbs2cPxcXFnHTS\nSWGH0qiSkSQygSHA79x9CLCLWFdT9elHmo4kImlr9erV9OvXj5ycnLBDaVSZSTjGFmCzu78VPH+a\nWJIoMbMu7l4SdCVtD36+FegZt3+PoCxReY2mTp1auZ2bm0tubm793oWIyCGkY1dTNBolGo3W6xhJ\nWSdhZkuBq939AzO7HWgV/GiHu88ws1uB9u6eHwxczwNOI9ad9CLQ193dzN4AJgIrgD8Dc9y9sIbX\n0zoJEWlUV111FYMHD2bChAlhh1JndVknkYyWBMQu7PPMrAXwEXAFkAEsMLOfApuIzWjC3dea2QJg\nLbAfuD7uij8BmAtkE5stdVCCEBEJw9tvv81VV10VdhiNTiuuRUQOY8+ePbRv354dO3ak9ZiEVlyL\niDSANWvW0Ldv37ROEHWlJCEichjpOGidLEoSIiKHoSQhIiIJNeckoYFrEZFD2Lt3L+3bt+fTTz+l\nVatWh98hhWngWkQkydasWUOfPn3SPkHUlZKEiMghNOeuJlCSEBE5pJUrVypJiIhIzZp7S0ID1yIi\nCZSWltKhQwc+++yzJrGQTgPXIiJJ9NZbb3HCCSc0iQRRV0oSIiIJvP766wwfPjzsMEKlJCEikoCS\nhJKEiEiN3F1JAiUJEZEabdy4ETOjV69eYYcSKiUJEZEaVLQizGo1GajJUZIQEamBuppilCRERGqg\nJBGjxXQiItXs3r2bjh07NplFdBVCW0xnZhvN7K9mtsrMlgdl7c1ssZkVm1mRmbWLqz/ZzNab2Toz\nGxVXPsTMVpvZB2Z2TzJiExGprbfeeovvfOc7TSpB1FWyupvKgVx3H+zuw4KyfOAld+8PLAEmA5jZ\nIOAiYCBwDnCffTMydD9wpbv3A/qZ2egkxScicsTU1fSNZCUJq+FY5wIFwXYBcF6wPRaY7+4H3H0j\nsB4YZmZdgTbuviKo92jcPiIijUZJ4hvJShIOvGhmK8zsqqCsi7uXALj7NqBzUN4d2By379agrDuw\nJa58S1AmItJo3J033nhDSSKQmaTjfM/d/25mnYDFZlZMLHHES+pI89SpUyu3c3Nzyc3NTebhRaSZ\n2rhxI5FIhGOPPTbsUOotGo0SjUbrdYykJAl3/3vw7ydm9kdgGFBiZl3cvSToStoeVN8K9IzbvUdQ\nlqi8RvFJQkQkWZrSIrrqH6DvuOOOWh+j3t1NZtbKzFoH20cBo4A1wCLg8qBaHvBssL0IGGdmWWZ2\nHHA8sDzoktppZsOCgezL4vYREWkUGo+oKhktiS7AQjPz4Hjz3H2xmb0FLDCznwKbiM1owt3XmtkC\nYC2wH7g+btHDBGAukA087+6FSYhPROSIvf7664wbNy7sMFKGFtOJiAR27dpF586d+eyzz8jOzg47\nnKTTN9OJiNTDX/7yF4YMGdIkE0RdKUmIiASi0ahmSlajJCEiElCSOJjGJEREiI1HdOnShe3bt9Oq\nVauww2kQGpMQEamj119/nZNPPrnJJoi6UpIQEUFdTYkoSYiIoCSRiMYkRKTZqxiPKCkp4aijjgo7\nnAajMQkRkTqoGI9oygmirpQkRKTZU1dTYkoSItLsKUkkpjEJEWnWmst4BGhMQkSk1jQecWhKEiLS\nrKmr6dCUJESkWYtGo4wYMSLsMFKWxiREpNmq+P6I7du3N4vuJo1JiIjUwmuvvabxiMNQkhCRZquo\nqIizzz477DBSmpKEiDRbhYWFShKHkbQkYWYRM1tpZouC5+3NbLGZFZtZkZm1i6s72czWm9k6MxsV\nVz7EzFab2Qdmdk+yYhMRqe7jjz9m+/btnHLKKWGHktKS2ZKYBKyNe54PvOTu/YElwGQAMxsEXAQM\nBM4B7jOzioGU+4Er3b0f0M/MRicxPhGRSkVFRYwaNYpIRB0qh5KUs2NmPYAfAf8VV3wuUBBsFwDn\nBdtjgfnufsDdNwLrgWFm1hVo4+4rgnqPxu0jIpJU6mo6MslKoXcDtwDx81K7uHsJgLtvAzoH5d2B\nzXH1tgZl3YEtceVbgjIRkaTav38/L7/8MqNGjTp85WYus74HMLMxQIm7v2NmuYeomtSFDVOnTq3c\nzs3N1YpJETlib7zxBn369KFLly5hh9KgotEo0Wi0Xseo92I6M/s1cClwAMgB2gALgaFArruXBF1J\nr7j7QDPLB9zdZwT7FwK3A5sq6gTl44AR7v6zGl5Ti+lEpM6mTJmCu/PrX/867FAaVSiL6dz9X939\nWHfvDYwDlrj7eOA54PKgWh7wbLC9CBhnZllmdhxwPLA86JLaaWbDgoHsy+L2ERFJGq2POHL17m46\nhN8AC8zsp8RaCRcBuPtaM1tAbCbUfuD6uGbBBGAukA087+6FDRifiDRD27dvZ8OGDQwfPjzsUNKC\n7t0kIs3K448/zjPPPMMzzzwTdiiNTvduEhE5DE19rR21JESk2SgvL6dr166sWLGCXr16hR1Oo1NL\nQkTkEFauXEnHjh2bZYKoKyUJEWk2/vSnP3HOOeeEHUZaUZIQkWZj4cKFnHee7vZTG0oSItIsfPTR\nR2zbto3TTz897FDSipKEiDQLCxcuZOzYsWRkZIQdSlpRkhCRZmHhwoWcf/75YYeRdjQFVkSavJKS\nEgYMGMC2bdto2bJl2OGERlNgRURq8Oyzz3L22Wc36wRRV0oSItLkqaup7tTdJCJN2s6dO+nZsydb\nt26lTZs2YYcTKnU3iYhU8+c//5nvf//7zT5B1JWShIg0aepqqh91N4lIk1VaWkrXrl3ZsGEDnTp1\nCjuc0Km7SUQkzksvvcTJJ5+sBFEPShIi0mQ99dRTXHDBBWGHkdbU3SQiTdLu3bvp1q0b69at45hj\njgk7nJSg7iYRkcCiRYs47bTTlCDqqd5JwsxamtmbZrbKzN4zs18H5e3NbLGZFZtZkZm1i9tnspmt\nN7N1ZjYqrnyIma02sw/M7J76xiYizde8efO49NJLww4j7SWlu8nMWrn7bjPLAF4DbgbGAp+5+0wz\nuxVo7+75ZjYImAecCvQAXgL6urub2ZvADe6+wsyeB2a7e1ENr6fuJhFJ6JNPPqFv375s2bKF1q1b\nhx1Oygitu8nddwebLYNjfg6cCxQE5QVAxTd9jAXmu/sBd98IrAeGmVlXoI27rwjqPRq3j4jIEVuw\nYAFjxoxRgkiCpCQJM4uY2SpgGxB197VAF3cvAXD3bUDnoHp3YHPc7luDsu7AlrjyLUGZiEitPP74\n4+pqSpLMZBzE3cuBwWbWFigys1ygen9QUvuHpk6dWrmdm5tLbm5uMg8vImlqw4YNfPTRR4wcOTLs\nUEIXjUaJRqP1OkbSp8Ca2a+AUuBKINfdS4KupFfcfaCZ5QPu7jOC+oXA7cCmijpB+ThghLv/rIbX\n0JiEiNTojjvuYMeOHcyePTvsUFJOKGMSZtaxYuaSmeUAI4FVwCLg8qBaHvBssL0IGGdmWWZ2HHA8\nsDzoktppZsPMzIDL4vYRETksd1dXU5Ilo7vpGKAguLBHgMfc/eVgjGKBmf2UWCvhIgB3X2tmC4C1\nwH7g+rhmwQRgLpANPO/uhUmIT0SaieXLlxOJRBg6dGjYoTQZWnEtIk3GjTfeSOfOnfnVr34Vdigp\nqS7dTUoSItIklJaW0rNnT1asWMFxxx0XdjgpSbflEJFm68knn+TUU09VgkgyJQkRaRIeeOABrrvu\nurDDaHKUJEQk7f31r39l8+bNjBkzJuxQmhwlCRFJew888ABXX301mZlJWR8scTRwLSJp7auvvqJX\nr168++67dOvWLexwUpoGrkWk2Zk3bx5nnnmmEkQDUZIQkbTl7hqwbmBKEiKStt544w127drFD3/4\nw7BDabKUJEQkbT3wwANce+21RCK6lDUUDVyLSFoqKSlhwIABrF+/no4dO4YdTlrQwLWINBtz5szh\nkksuUYJoYGpJiEja+fLLL+nduzfLly+nd+/eYYeTNtSSEJFm4cEHH+Sss85SgmgEakmISFrZu3cv\nffr04bnnnmPw4MFhh5NW1JIQkSZv3rx5nHDCCUoQjUQtCRFJG+Xl5QwaNIj77ruPH/zgB2GHk3bU\nkhCRJm3RokW0adOGM888M+xQmg0lCRFJC+7OjBkzuPXWWzGr1YdhqYd6Jwkz62FmS8zsPTNbY2YT\ng/L2ZrbYzIrNrMjM2sXtM9nM1pvZOjMbFVc+xMxWm9kHZnZPfWMTkabj5ZdfZseOHZx//vlhh9Ks\nJKMlcQC4yd1PAIYDE8xsAJAPvOTu/YElwGQAMxsEXAQMBM4B7rNvPhbcD1zp7v2AfmY2OgnxiUia\nc3fy8/P5t3/7NzIyMsIOp1mpd5Jw923u/k6w/TWwDugBnAsUBNUKgPOC7bHAfHc/4O4bgfXAMDPr\nCrRx9xVBvUfj9hGRZuypp57C3bnwwgvDDqXZSerXOJnZt4GTgTeALu5eArFEYmadg2rdgdfjdtsa\nlB0AtsSVbwnKRaQZ279/P1OmTOF3v/udbuQXgqQlCTNrDTwFTHL3r82s+hzVpM5ZnTp1auV2bm4u\nubm5yTy8iKSIRx55hJ49e3LWWWeFHUraiUajRKPReh0jKeskzCwT+BPwgrvPDsrWAbnuXhJ0Jb3i\n7gPNLB9wd58R1CsEbgc2VdQJyscBI9z9ZzW8ntZJiDQDu3fvpl+/fixcuJBTTz017HDSXpjrJH4P\nrK1IEIFFwOXBdh7wbFz5ODPLMrPjgOOB5e6+DdhpZsOCgezL4vYRkWbo3nvv5bvf/a4SRIjq3ZIw\ns+8BrwJriHUpOfCvwHJgAdCTWCvhInf/IthnMnAlsJ9Y99TioPwUYC6QDTzv7pMSvKZaEiJN3Oef\nf06/fv1YtmwZ/fv3DzucJqEuLQndlkNEUtLEiRPZu3cv//mf/xl2KE1GXZJEUmc3iYgkw8qVK1mw\nYAHvvfde2KE0e5pPJiIppaysjOuuu47p06fToUOHsMNp9pQkRCSlPPjgg7Rs2ZK8vLywQxE0JiEi\nKaSkpIQTTzyRJUuW8J3vfCfscJocDVyLSFq77LLL6Nq1KzNnzgw7lCZJA9cikrai0ShLly7VYHWK\n0ZiEiIRu586dXHHFFdx33320bt067HAkjrqbRCR048ePp3Xr1tx///1hh9KkqbtJRNLO/PnzWbFi\nBStXrgw7FKmBWhIiEpqPP/6YoUOH8sILL3DKKaeEHU6TF+YN/kREaqWsrIzx48dz0003KUGkMCUJ\nEQnFjBkzALjllltCjkQORWMSItLoXnjhBe69916WL1+u76xOcUoSItKoiouLycvL45lnnqFnz55h\nhyOHoe4mEWk0X3zxBWPHjmX69OmcccYZYYcjR0Czm0SkUZSVlTFmzBj69+/P7NmzD7+DJJ1mN4lI\nyvrFL35BWVkZs2bNCjsUqQWNSYhIg5s+fTovvPACy5YtIzNTl510ot+WiDSoOXPm8PDDD/Pqq69y\n9NFHhx2O1FJSupvM7GEzKzGz1XFl7c1ssZkVm1mRmbWL+9lkM1tvZuvMbFRc+RAzW21mH5jZPcmI\nTUTC8/DDDzNr1ixefvllunXrFnY4UgfJGpN4BBhdrSwfeMnd+wNLgMkAZjYIuAgYCJwD3GdmFQMp\n9wNXuns/oJ+ZVT+miKSJJ554gttuu40XX3yRXr16hR2O1FFSkoS7LwM+r1Z8LlAQbBcA5wXbY4H5\n7n7A3TcC64FhZtYVaOPuK4J6j8btIyJpZO7cudx0000UFRXRr1+/sMORemjIMYnO7l4C4O7bzKxz\nUN4deD2u3tag7ACwJa58S1AuImnC3Zk+fToPPfQQr7zyCgMGDAg7JKmnxhy4TurChqlTp1Zu5+bm\nkpubm8zDi0gtlZWVMXHiRJYtW8Zrr72mMYgUEI1GiUaj9TpG0hbTmVkv4Dl3Pyl4vg7IdfeSoCvp\nFXcfaGb5gLv7jKBeIXA7sKmiTlA+Dhjh7j+r4bW0mE4khZSWlnLppZeyY8cO/vjHP9KuXbvD7ySN\nLuzFdBY8KiwCLg+284Bn48rHmVmWmR0HHA8sd/dtwE4zGxYMZF8Wt4+IpKgNGzYwfPhwsrOzKSws\nVIJoYpI1BfYPwF+IzUj62MyuAH4DjDSzYuCHwXPcfS2wAFgLPA9cH9csmAA8DHwArHf3wmTEJyIN\nY+HChZx++ulcffXVPP7447Rs2TLskCTJdO8mEam1/fv3k5+fz9NPP82CBQsYNmxY2CHJEdB3XItI\ng3vnnXe48sorOeaYY3j77bfp0KFD2CFJA9IN/kTkiOzZs4cpU6YwatQobrzxRp577jkliGZALQkR\nOaz/+Z//4ZprrmHQoEGsXr2arl27hh2SNBIlCRFJqLi4mPz8fFauXMmsWbO48MILww5JGpm6m0Tk\nINu3b2fChAmcccYZDB8+nOLiYiWIZkpJQkQqbdq0iUmTJjFgwACysrJ4//33+cUvfkF2dnbYoUlI\nlCREhDVr1jB+/HiGDBlCy5Yteffdd7n77rs1MC1KEiJNSVFRET8ZNYqfjBpFUVHRQc/j7dq1i7lz\n53LGGWcwevRoTjjhBD788ENmzpyp+y5JJS2mE2kiioqKyDv/fGaUlgLw86wsWgC/3bcPgFtzcnjk\n6adp06YN8+bN47//+785/fTTueqqqxgzZgwtWrQIMXppDHVZTKckIZJGioqKeHDWLACuuflmRo/+\n5nu5fjJqFGNffJG84Plw4DrgUmL3zLkTWNayJcf37cvFF19MXl4ePXv2bNw3IKHSimuRJqx6SyFv\n2TJunDKFpUtXArDv008r624DPgEeAG4GuhG7k+aJvXuzfM2aRo5c0pmShEgjOVQroKafA8ya9SAA\nN998DQ/OmsWM0tLKlgKlpUy67U52lj8IfE4k8mf+AvwS2AHsBbYS4S7K6QBMpAV9so9u4HcpTY2S\nhEgjqKkVcPZFFxF97nkAcn/8I55/4onK8YN/eeUVvoq0Zt++ewBYtiyP/t1aA/AVsBJ4AdhVng3c\nAeykvLwV5ZzGZ+whQhvK2EoZY7id1QB8zUl07Pi3RnzX0hRoTEKkBof71H+kx6hoCez79G9csWpV\nZSvg/wIPAXOC5z8H7oHKnxcA13IsGXSnnF3s4SgyMlZSXraXFkAPYDOwn6G0ogXGUexiC5Afd5T/\nSyTye8rL7wYgJ+dWFi4sqNN7kaZBYxIiSXC4vv+bb74GqNoV9NZbb/HQf8wG4OqbJjF06FDOPz+P\n0tIZALS1q7gi7jVeI5YgxhP7kvdjgCixL1n5AFgK7OVjuvAxxwBriWCZbdlbdgkt2Mx2IuxnN215\nrTLRTAR2Z07iwIHY85ycx5ky5WaWLl0UxKkEIbWnloQ0eUfSKrjrrrsqL/Lt27fl5x9+WOVT/aRI\nVtD3D1lZPwdasG/fbwHIzJxAqwO7qlysO/U5kQ8/HMVR/BVnL7t5nxx28E+U8QmwDOgAfAG0BzKA\n7cC5QGfgMeA/gKviYrinTx+K//frysTTLnINs8v3VYnzkcGDyep4HBBLXkoKEk8tCWl0yeiWSXYM\n8M2n/BEjhjD7zjsr+/rHL13KY4sWVanTrVsbFhYUfNP1s+OTg16jrLwvFd04+/Y9QGxy6b8AJbQ4\nEPuu3S+BLcAgYPmH72KsYQ9wFLFVq6V0549kYuRwAONr1jEHaEkssZRFWlFYHrvAR6yYFn6gSgy9\ne/fmN7+7Oa4L6wRYtapKnQ4dO/L04qdreQZFDsHd0+4RC7tpKyws9AtGjvQLRo70wsLCpP88WTF2\nycnxueBzwbvk5Pi0adPqFde0adO819GdvNfRnXzatGmV+4wceYGPHHmBFxYWVqmTl5fnnbKyKmPo\nkJnpmZlHeWtO8Nac4G2J+FxwDx5zwU/u08ezsr5VWacNGVXq3ATeGvyOYDsbcxjkmfTwDDo5ZDu0\ndjA3sjxCxAeDXw9+F/jV4G0jEZ8JXhr3um2thcNch7mek9PF8/LyqrzX+Pc5bdq0g85t9fNV0/lv\nqN+1NA3BtbN219va7tDQD+Bs4H1iXbO3JqiTzPOWcg538T3cBSTRxSPZieOCkSOrXFxvBv+WWeXr\ndsrK8mnTplW58MVf0Kv/PC8vz9sGP5sL3hY8Ly/Pc3K6VF5cMzOPqlLnW8G/8TG0xQ76eTn4F+C/\nAc+xTM8hwy8F/zF4Fvg/gJ8C3g08AxwiHqG1Z9DBobNngV8KfgN4yyCJ/Fdw7FZBrPFxn9ynz0HJ\nacTgwVWS3ZH8HRzu99UYHwak6Uj7JEGsVb4B6AW0AN4BBtRQL7lnrpEd7j92TRff9pHIN5+WIwd/\nOv7e4MGVxxwxeHCNPz/UBbr6J/Rp06Yd9Lz6J/rqcX632gW7+qfntpZ50M9bkZHwE/1c8KMzc4L9\nY8WtaedzwfeBbwM/CXwy+ELwh8B7go8GHw9+NngXcAPPBM8OtjMwPwX8Z0FrYSR4Dvht4LPA24Bn\nZrarjLtdJOuw7/PEPn0OOlf6lC+ppi5JItXGJIYB6919E4CZzSc2lvd+qFHVUvzUx+ozYWrqI//R\nJZdUzpe/+qZJADiwk1g/94vAjeXl9AT2AB3Ky1lObECpDJgPrFq1in+peH1gINAqqPMOsO6997hq\n3z76AFnADfv2Mf2XU9nFXUAGS5acT05ZKb8lNog66Ze/JBO4NzjmxF/+kq8jrWgV9JmPXTqOSy45\nl4lx77umX1KZ9wX+GdhHmd/I53zF34L38VKsBtfyHnuB3wdlfye2FuBN4OuyfcCdRPg5cIBd7OI6\n4Erg6OD9zwROAloTG/w9HfgB0AlYAcwgQgb9MSKUZ21lUM8O3FhtYPrDY46hYG9sDODWYHbSofr+\nq+vbuzdPL15cpaxg4cLKsZKCkMZrROorpWY3mdlPgNHufk3w/FJgmLtPrFbPUynueEVFRYwdO46s\nfd0B2JO5Ecgk+0APACKsYw7lXAwUA7cRu6h/D/gUeBdo0bIle/fuJRvIIZYsegHfBrKJ9cNtAoYC\nBrwFnAr0IZZc3iLWHDuZ2EV0DbHJDMfjHA3sJ3bbhi20IJZGHNjDUcF7OEDsIg6xhFGRjA4E8USA\n3YBjxIZdy+L2dLKDkv2VZ8WCh2NBDC2I3TaiT/C+WgXvaz0wMjjqC0DvQYP429q1TAzqTA/2nRMc\ncSJUSV57MjfSjr3MCuaB3pKVxaTbbjto+ur4sWMrE/UtWVk8tmhRwov4kdw4r2DhQiUBSXlpf4O/\nppAkhgw5gw9XLWdOcIm8HiODCDdRxgZgAdCV2AWyN7FPzacC/0TsYvk2cE9GS3aWXUtr1gHwNeto\nS0nlMScS4UvOpzVfBq/6Gv+P3VU+GV9HLzLpF+x/Eq1b/4HI159WiSuTDOYQu5hWX8w1HLiW2Pyd\nA8AI4KfAxUA58ATwq4xsviibGZRGgNuBAo6iJ2Dstg9p7fu+aY1Q9YJekTDj477zmGMoCz7RX33T\nJFYuXVrlpnUFwKw+ffjy8y8r68R/6q9IAoebcVXbWVk1zaAKe1aXSG01hSmwW4Fj4573CMoOMnXq\n1Mrt3NxccnNzGzKuI7Zj0wfMYT/fBwqBbJy9lPEk8I9ALrE7cs4h1u3zc2KdMRUXwR1AhkWAIXzN\n7KB0OF9yMTdU3l5hM/Bjvq7cayAT+RAqk0gme+0z9nisAyon51by82/kzjv/nRuCFk4m65jDgcoj\nrIEqXUffgbh+AAAG/ElEQVSrgUlAxV/TemKtiPbB8zbAt7/dneL/vYvS0rbB6zzNlCmT4244l11l\nlTHAI//Qv3Ie/4gRl3DLnXdC/Cf6Rx6pemfTpUsPOsc1de1Uv0gf7qI9evToWl3Ya6qvxCCpLhqN\nEo1G63eQ2g5iNOSDWO9GxcB1FrHu9IE11Kv3AE5DqRg0/i14Hngf8LsPMVgbibQ87IyerKxveVZW\np0M+z8xsV3nMrKxvHTQo7V51KmlNg9vVB18PNdW0U1ZW5YypRDN2qg9szwW/YOTIKnWOZCqvBoBF\nkoN0n90Uew+cTay7fj2Qn6BOMs9bUhUWFla5mLbLzPQOmZm1mlVUcZz4OrV9fiRx1vbiW9vplsm6\nwGuap0hy1CVJpNSYxJFK5TEJSJ/+68ZYLZ0KK7JFJCbtB66PVKonCRGRVFSXJBFpqGBERCT9KUmI\niEhCShIiIpKQkoSIiCSkJCEiIgkpSYiISEJKEiIikpCShIiIJKQkISIiCSlJiIhIQkoSIiKSkJKE\niIgkpCQhIiIJKUmIiEhCShIiIpKQkoSIiCSkJCEiIgkpSYiISEL1ShJmdqGZvWtmZWY2pNrPJpvZ\nejNbZ2aj4sqHmNlqM/vAzO6JK88ys/nBPq+b2bH1iU1EROqvvi2JNcD5wNL4QjMbCFwEDATOAe4z\ns4rvVb0fuNLd+wH9zGx0UH4lsMPd+wL3ADPrGVvootFo2CEcEcWZPOkQIyjOZEuXOOuiXknC3Yvd\nfT1Q/Yu1zwXmu/sBd98IrAeGmVlXoI27rwjqPQqcF7dPQbD9FPDD+sSWCtLlD0dxJk86xAiKM9nS\nJc66aKgxie7A5rjnW4Oy7sCWuPItQVmVfdy9DPjCzI5uoPhEROQIZB6ugpm9CHSJLwIcmOLuzzVU\nYBzcOhERkcbm7vV+AK8AQ+Ke5wO3xj0vBE4DugLr4srHAffH1wm2M4Dth3g910MPPfTQo/aP2l7f\nD9uSqIX4T/6LgHlmdjexbqTjgeXu7ma208yGASuAy4A5cfvkAW8C/wdYkuiF3F2tDBGRRlCvJGFm\n5wH3Ah2BP5nZO+5+jruvNbMFwFpgP3C9B00AYAIwF8gGnnf3wqD8YeAxM1sPfEaslSEiIiGyb67d\nIiIiVaXtimszu93MtpjZyuBxdtgxVTCzs83s/WDB4K1hx5OImW00s7+a2SozWx52PBXM7GEzKzGz\n1XFl7c1ssZkVm1mRmbULM8YgppriTLm/SzPrYWZLzOw9M1tjZhOD8pQ6pzXEeWNQnjLn1Mxamtmb\nwf+Z98zs10F5qp3LRHHW+lymbUvCzG4HvnL3/wg7lnhmFgE+ILbO43+Jjb2Mc/f3Qw2sBmb2EXCK\nu38edizxzOwM4GvgUXc/KSibAXzm7jODxNve3fNTMM6U+7sM1id1dfd3zKw18DaxdUlXkELn9BBx\nXkwKnVMza+Xuu80sA3gNuBkYSwqdy0PEeRa1PJdp25IIpOIA9jBgvbtvcvf9wHxif+ipyEjBvwF3\nXwZUT1zxiy0L+GYRZmgSxAkp9nfp7tvc/Z1g+2tgHdCDFDunCeKsWEeVMufU3XcHmy2J/f/5nBQ7\nl5AwTqjluUy5C0Qt3WBm75jZf4XdvItTfSFh/ILBVOPAi2a2wsyuDjuYw+js7iUQu5gAnUOO51BS\n8e8SADP7NnAy8AbQJVXPaVycbwZFKXNOzSxiZquAbUDU3deSgucyQZxQy3OZ0knCzF602M0AKx5r\ngn9/DNwH9Hb3k4mdhJRoiqaZ77n7EOBHwISg+yRdpGo/acr+XQZdOE8Bk4JP6tXPYUqc0xriTKlz\n6u7l7j6YWGvsH80slxQ8l9Xi/L6ZjaAO5zKZ6ySSzt1HHmHVh4CGXP1dG1uB+DvY9gjKUo67/z34\n9xMzW0isq2xZuFElVGJmXdy9JOi73h52QDVx90/inqbM36WZZRK78D7m7s8GxSl3TmuKM1XPqbt/\naWbPA0NJwXNZIYjzz8BQd4+/GesRncuUbkkcSvCLqHAB8G5YsVSzAjjezHqZWRax9R6LQo7pIGbW\nKvjEhpkdBYwidc4hxPpNqy/QvDzYzgOerb5DSKrEmcJ/l78H1rr77LiyVDynB8WZSufUzDpWdNGY\nWQ4wElhFip3LBHG+U5dzmc6zmx4l1mdZDmwErq3oEwxbMK1sNrEk/LC7/ybkkA5iZscBC4k1izOB\neakSp5n9AcgFOgAlwO3AH4EngZ7AJuAid/8irBghYZxnkmJ/l2b2PeBVYrf2r7g9w78Cy4EFpMg5\nPUSc/0yKnFMzO5HYwHTFpI/H3P3fLXYz0lQ6l4nirPV1M22ThIiINLy07W4SEZGGpyQhIiIJKUmI\niEhCShIiIpKQkoSIiCSkJCEiIgkpSYiISEJKEiIiktD/B/7k7a/86AalAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from sklearn.linear_model import Ridge\n", "clf = Ridge(alpha=0.01) #Is the alpha that affects the weights!!\n", "clf.fit(x_5, y_n)\n", "print clf.coef_\n", "predict_y_n = clf.predict(x_p)\n", "plt.scatter(x, y, c ='blue')\n", "plt.scatter(x, y_n, c ='red')\n", "plt.plot(x_predict, predict_y_n, c ='black', label = '3rd Degree')" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[ 0.00715796 0.06315591 0.30537708 -0.02658066 0.00070794]]\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VPWd//HXJzeIcpVrC4p44ZJWBKqotZVUC4g/FxQf\ntdiq2Nq1WxGw0q2AXaEuaHGLF3a9IhXwoYVoV8EWAwgOSquCCKsFlNSKEAQsXqDcAiGf3x9zEoaY\nQGYyyZlk3s/HYx6c+c45Zz5zjOdzzvd2zN0REZH0lhF2ACIiEj4lAxERUTIQERElAxERQclARERQ\nMhARERJIBmaWYWZvm9mC4H1rM1tsZu+b2SIzaxmz7ngzKzKzDWY2MKa8r5m9Y2YbzeyB5PwUERFJ\nVCJ3BmOA9THvxwEvu3t3YBkwHsDM8oCrgZ7AYOBhM7Ngm0eAG929G9DNzAYlGL+IiCRBXMnAzDoD\nlwFPxBQPBWYHy7OBK4LlIcBcdy91901AEdDPzDoCzd19VbDenJhtREQkBPHeGdwP/DsQO2y5g7vv\nAHD37UD7oLwTsCVmva1BWSegOKa8OCgTEZGQ1DgZmNn/A3a4+1rAjrGq5rcQEWlgsuJY90JgiJld\nBuQCzc3sKWC7mXVw9x1BFdAnwfpbgZNjtu8clFVX/iVmpsQiIpIAdz/WRfuX1PjOwN0nuPsp7n4a\nMBxY5u7XAS8CNwSrjQDmB8sLgOFmlmNmXYEzgJVBVdIuM+sXNChfH7NNVd+b8q+JEyeGHkNjiFFx\nKs5UfzWUOBMRz51BdX4DFJjZj4GPiPYgwt3Xm1kB0Z5Hh4Cb/UiUI4FZQFNgobsXJiEOERFJUELJ\nwN2XA8uD5c+A71az3j3APVWUrwbOSuS7RUQk+TQCOQny8/PDDuG4GkKMoDiTTXEmV0OJMxGWaP1S\nfTAzT+X4RERSkZnhddWALCIijZeSgYiIKBmIiIiSgYiIoGQgIiIoGYiICEoGIiKCkoGIiKBkICIi\nKBmIiAhKBiIigpKBiIigZCAiIigZiIgIcSQDM2tiZm+a2RozW2dmdwflE82s2MzeDl6Xxmwz3syK\nzGyDmQ2MKe9rZu+Y2UYzeyC5P0lEpHbKysq44447OHz4cNih1Jt4noFcAnzH3fsAvYCLzezC4OP7\n3L1v8CoEMLOeRB+B2RMYDDwcPPMY4BHgRnfvBnQzs0FJ+j0iIrW2bt06nn32WTIzM8MOpd7EVU3k\n7vuCxSbBtp8H76t6iMJQYK67l7r7JqAI6GdmHYHm7r4qWG8OcEW8gYuI1JUVK1bwrW99K+ww6lVc\nycDMMsxsDbAdiLj7+uCjW8xsrZk9YWYtg7JOwJaYzbcGZZ2A4pjy4qBMRCQlpGMyyIpnZXcvA/qY\nWQtgsZn1Bx4G7nJ3N7PJwDTgJ8kKcNKkSRXL+fn5jfoZpCKSGlasWMHEiRPDDqPGIpEIkUikVvtI\n+BnIZvYfwD53nxZT1gV40d17mdk4wN19avBZITAR+Ah4xd17BuXDgf7u/rMqvkPPQBaRerV582bO\nPfdctm/fzpFmzoalTp+BbGZty6uAzCwXGACsDdoAyg0D/hosLwCGm1mOmXUFzgBWuvt2YJeZ9Qsa\nlK8H5scTtIhIXSmvImqoiSBR8VQTfQWYHZzAM4Cn3H2pmc0xs95AGbAJ+CmAu683swJgPXAIuDnm\nMn8kMAtoCiws74EkIhK2dGwvgFpUE9UHVROJSH3r1asXM2fO5Nxzzw07lIQlUk2kZCAiEvj88885\n5ZRT+Oyzz8jOzg47nITVaZuBiEhj95e//IV+/fo16ESQKCUDEZFAurYXgJKBiEiFFStW8O1vfzvs\nMEKhNgMREeDAgQO0bduWbdu20bx587DDqRW1GYiIJGj16tX06NGjwSeCRCkZiIgAr732Wtq2F4CS\ngYgIkN6Nx6A2AxERysrKaNu2LevXr6djx47H3yDFqc1ARCQB69at46STTmoUiSBRSgYikvYikUja\nT4+vZCAiaW/p0qVccsklYYcRKrUZiEhaO3z4MG3btmXDhg2NpppIbQYiInFas2YNX/3qVxtNIkiU\nkoGIpLVly5Zx8cUXhx1G6JQMRCStKRlExfPYyyZm9qaZrTGzdWZ2d1De2swWm9n7Zrao/NGYwWfj\nzazIzDaY2cCY8r5m9o6ZbTSzB5L7k0REaubgwYP85S9/SfueRBBHMnD3EuA77t4H6AVcbGYXAuOA\nl929O7AMGA9gZnnA1UBPYDDwsB15qOgjwI3u3g3oZmaDkvWDRERqauXKlXTr1o3WrVuHHUro4qom\ncvd9wWKTYNvPgaHA7KB8NnBFsDwEmOvupe6+CSgC+plZR6C5u68K1psTs42ISL1ZunSpqogCcSUD\nM8swszXAdiDi7uuBDu6+A8DdtwPtg9U7AVtiNt8alHUCimPKi4MyEZF6pfaCI7LiWdndy4A+ZtYC\nWGRm+UDlgQBJHRgwadKkiuX8/HzV7YlIUuzbt4/Vq1c3isnpIpEIkUikVvuIKxmUc/fdZrYQOAfY\nYWYd3H1HUAX0SbDaVuDkmM06B2XVlVcpNhmIiCTLn//8Z3r37k2zZs3CDqXWKl8o//rXv457H/H0\nJmpb3lPIzHKBAcAaYAFwQ7DaCGB+sLwAGG5mOWbWFTgDWBlUJe0ys35Bg/L1MduIiNSLZcuWpf0U\nFLHiuTP4CjA7OIFnAE+5+9KgDaHAzH4MfES0BxHuvt7MCoD1wCHg5pi5JUYCs4CmwEJ3L0zKrxER\nqaFly5Zx7733hh1GytDcRCKSdnbt2kXnzp3ZuXMnTZo0CTucpNPcRCIiNbB8+XLOP//8RpkIEqVk\nICJpZ9GiRQwYMCDsMFKKkoGIpBV356WXXmLw4MFhh5JSlAxEJK1s3LiRgwcP8vWvfz3sUFKKkoGI\npJWXXnqJSy+9lCNTpQkoGYhImlEVUdXUtVRE0sa+ffvo0KEDxcXFtGzZ8vgbNFDqWioicgyvvPIK\nffv2bdSJIFFKBiKSNlRFVD0lAxFJC+pSemxKBiKSFoqKijhw4AC9evUKO5SUpGQgImlBXUqPTclA\nRNJCYWGhqoiOQV1LRaTR279/P+3bt2fLli20atUq7HDqnLqWiohUIRKJ0Lt377RIBIlSMhCRRm/h\nwoWqIjqOeB572dnMlpnZOjN718xGBeUTzazYzN4OXpfGbDPezIrMbIOZDYwp72tm75jZRjN7ILk/\nSUTkCHdn/vz5DB06NOxQUlo8j70sBW5z97Vm1gxYbWZLgs/uc/f7Ylc2s55EH4HZk+hD7182szOD\nRoBHgBvdfZWZLTSzQe6+qPY/R0TkaG+//TZNmzYlLy8v7FBSWo3vDNx9u7uvDZb3ABuATsHHVTVU\nDAXmunupu28CioB+ZtYRaO7uq4L15gBXJBi/iMgxPf/881xxxRXqUnocCbUZmNmpQG/gzaDoFjNb\na2ZPmFn5pB+dgC0xm20NyjoBxTHlxRxJKiIiSfXCCy9w5ZVXhh1GyounmgiAoIroOWCMu+8xs4eB\nu9zdzWwyMA34SbICnDRpUsVyfn4++fn5ydq1iDRyRUVFfPrpp5x33nlhh1KnIpEIkUikVvuIa5yB\nmWUBfwRecvcHq/i8C/Ciu/cys3GAu/vU4LNCYCLwEfCKu/cMyocD/d39Z1XsT+MMRCRh9957L3//\n+9959NFHww6lXtXHOIPfAetjE0HQBlBuGPDXYHkBMNzMcsysK3AGsNLdtwO7zKyfRSvxrgfmxxmH\niMhxqYqo5mp8Z2BmFwKvAu8CHrwmAD8g2n5QBmwCfuruO4JtxgM3AoeIVistDsq/AcwCmgIL3X1M\nNd+pOwMRSci2bdvIy8tjx44d5OTkhB1OvUrkzkDTUYhIo/TYY4+xfPlynnnmmbBDqXeajkJEJFDe\npVRqRncGItLo7Nq1i5NPPpmtW7fSvHnzsMOpd7ozEBEhOhfRRRddlJaJIFFKBiLS6KiKKH6qJhKR\nRmXPnj107tyZoqIi2rVrF3Y4oVA1kYikvQULFvDNb34zbRNBopQMRKRR+f3vf88111wTdhgNjqqJ\nRKTR+Oyzz+jatSvFxcVp3XisaiIRSWt/+MMfGDhwYFongkQpGYhIo/H73/+eH/zgB2GH0SCpmkhE\nGoWPP/6Yr3/963z88cc0bdo07HBCpWoiEUlb8+bNY+jQoWmfCBKlZCAijYJ6EdWOkoGINHh/+9vf\n2Lx5MxdffHHYoTRYSgYi0uDNnTuX733ve2Rlxf0kXwnUOBmYWWczW2Zm68zsXTMbHZS3NrPFZva+\nmS0ys5Yx24w3syIz22BmA2PK+5rZO2a20cweSO5PEpF04u6qIkqCeO4MSoHb3P1rwAXASDPrAYwD\nXnb37sAyYDyAmeUBVwM9gcHAw8FjLgEeAW50925ANzMblJRfIyJpZ9WqVZSUlHDBBReEHUqDVuNk\n4O7b3X1tsLwH2AB0BoYCs4PVZgPlUwUOAea6e6m7bwKKgH7BM5Obu/uqYL05MduIiMRl5syZ/PjH\nP+bItaYkIqEKNjM7lehzj98AOpQ/89jdt5tZ+2C1TsDrMZttDcpKgeKY8uKgXEQkLnv37uXZZ5/l\n3XffDTuUBi/uBmQzawY8R/QB93uAyqPCNEpMROrFc889xze/+U06ddL1ZG3FdWdgZllEE8FT7j4/\nKN5hZh3cfUdQBfRJUL4VODlm885BWXXlVZo0aVLFcn5+Pvn5+fGELCKN2MyZM/n5z38edhihi0Qi\nRCKRWu0jrukozGwOsNPdb4spmwp85u5Tzex2oLW7jwsakJ8GziNaDbQEONPd3czeAEYDq4A/AdPd\nvbCK79N0FCJSpY0bN3LRRRexZcsWsrOzww4npSQyHUWN7wzM7ELgh8C7ZraGaHXQBGAqUGBmPwY+\nItqDCHdfb2YFwHrgEHBzzJl9JDALaAosrCoRiIgcy+9+9zuuu+46JYIk0UR1ItLglJaWcvLJJ/PK\nK6/Qo0ePsMNJOZqoTkTSwsKFCznttNOUCJJIyUBEGpyZM2dy4403hh1Go6JqIhFpULZt20ZeXh5b\ntmyhWbNmYYeTklRNJCKN3qOPPso111yjRJBkujMQkQajpKSELl26sGzZMvLy8sIOJ2XpzkBEGrWC\nggJ69eqlRFAHlAxEpEFwdx588EFGjx4ddiiNkpKBiDQIb7zxBl988QWXXXZZ2KE0SkoGItIgTJ8+\nnVtuuYWMDJ226oIakEUk5W3dupWzzjqLDz/8kJYtWx5/gzSnBmQRaZQeffRRfvjDHyoR1CHdGYhI\nSjtw4ABdunTh1VdfpXv37mGH0yDozkBEGp2nn36avn37KhHUMd0ZiEjKOnz4MD179mTGjBn0798/\n7HAaDN0ZiEij8uyzz9KuXTsuuuiisENp9OJ67KWISH0pKyvj7rvv5je/+Q1mcV3kSgJqfGdgZjPN\nbIeZvRNTNtHMis3s7eB1acxn482syMw2mNnAmPK+ZvaOmW00sweS91NEpDH54x//SFZWFoMHDw47\nlLQQTzXRk8CgKsrvc/e+wasQwMx6En38ZU9gMPCwHUntjwA3uns3oJuZVbVPEUlj7s6UKVOYMGGC\n7grqSY2TgbuvAD6v4qOq/ksNBea6e6m7bwKKgH5m1hFo7u6rgvXmAFfEF7KINHbLli1j9+7dDBs2\nLOxQ0kYyGpBvMbO1ZvaEmZWPCOkEbIlZZ2tQ1gkojikvDspERCpMmTKF8ePHa+qJelTbI/0wcJq7\n9wa2A9NqH5KIpLPXX3+dDz/8kGuuuSbsUNJKrXoTufs/Yt7OAF4MlrcCJ8d81jkoq668WpMmTapY\nzs/PJz8/P+F4RST13XnnnYwbN47s7OywQ2kwIpEIkUikVvuIa9CZmZ0KvOjuZwXvO7r79mD558C5\n7v4DM8sDngbOI1oNtAQ4093dzN4ARgOrgD8B08sbnqv4Pg06E0kjS5YsYeTIkaxbt07JoBYSGXRW\n4zsDM3sGyAfamNlmYCLwHTPrDZQBm4CfArj7ejMrANYDh4CbY87qI4FZQFNgYXWJQETSS1lZGePG\njWPKlClKBCHQdBQikhLmzZvHb3/7W1auXKnupLWUyJ2BkoGIhO7QoUP07NmTxx57jEsuuSTscBo8\nzU0kIg3SjBkzOO2005QIQqQ7AxEJ1Z49ezjzzDP505/+RN++fcMOp1HQnYGINDj3338/+fn5SgQh\n052BiIRm8+bN9OnTh1WrVnHaaaeFHU6joTsDEWlQxo4dy6hRo5QIUoCeZyAioVi8eDGrV69mzpw5\nYYci6M5AREJQUlLCqFGjmD59Orm5uWGHIygZiEgI7rvvPrp3787ll18edigSUAOyiNSrzZs307dv\nX1atWkXXrl3DDqdRUgOyiKS8W2+9ldGjRysRpBg1IItIvSkoKGD9+vU888wzYYcilaiaSETqxY4d\nOzj77LOZP38+5513XtjhNGqaqE5EUpK7c+WVV5KXl8fdd98ddjiNXp0+z0BEJFFPP/00f//735k3\nb17YoUg1dGcgInVq69at9OnTh0WLFtGnT5+ww0kLddqbyMxmmtkOM3snpqy1mS02s/fNbJGZtYz5\nbLyZFZnZBjMbGFPe18zeMbONZvZAPMGKSMPi7tx0002MHDlSiSDFxdO19ElgUKWyccDL7t4dWAaM\nBwiegXw10BMYDDxsRx5d9Ahwo7t3A7qZWeV9ikgj8eCDD/LJJ58wYcKEsEOR46hxMnD3FcDnlYqH\nArOD5dnAFcHyEGCuu5e6+yagCOhnZh2B5u6+KlhvTsw2ItKIvP7669xzzz0UFBTomcYNQG0HnbV3\n9x0A7r4daB+UdwK2xKy3NSjrBBTHlBcHZSLSiOzcuZPvf//7zJgxQ4PLGohk9yZKemvvpEmTKpbz\n8/PJz89P9leISBKVlZVx7bXXMnz4cIYMGRJ2OGkhEokQiURqtY+4ehOZWRfgRXfvFbzfAOS7+46g\nCugVd+9pZuMAd/epwXqFwETgo/J1gvLhQH93/1k136feRCINzJQpUygsLGTZsmWqHgpJfcxNZMGr\n3ALghmB5BDA/pny4meWYWVfgDGBlUJW0y8z6BQ3K18dsIyIN3OLFi3nooYeYO3euEkEDU+NqIjN7\nBsgH2pjZZqJX+r8BnjWzHxO96r8awN3Xm1kBsB44BNwcc4k/EpgFNAUWunthcn6KiITpr3/9K9de\ney3/+7//S6dOagpsaDToTERqbdu2bVxwwQXcc889XHPNNWGHk/Y0hbWI1Lu9e/cyZMgQfvKTnygR\nNGC6MxCRhB0+fJirrrqKVq1a8eSTT3JkbKmESRPViUi9cXdGjRrFrl27KCgoUCJo4JQMRCRu7s5t\nt93G6tWrWbJkCTk5OWGHJLWkZCAicXF3xo0bx/Lly1m6dCktWrQIOyRJAiUDEYnLxIkTeemll3jl\nlVdo3bp12OFIkigZiEiNuDt33XUXzz33HJFIhDZt2oQdkiSRkoGIHNfhw4cZM2YMr732GkuXLqV9\n+/bH30gaFCUDETmmkpISrrvuOj755BNeffVVWrZsefyNpMHRoDMRqdbu3bu57LLLKCsro7CwUImg\nEVMyEJEqbdq0iYsuuoju3bszb948mjZtGnZIUoeUDETkS15++WXOP/98brjhBh566CEyMzPDDknq\nmNoMRKSCu/Pb3/6W++67j7lz5+phUmlEyUBEANi1axc33XQTH3zwAW+++SannHJK2CFJPVI1kYjw\n6quvcvbZZ9O6dWtee+01JYI0pDsDkTRWUlLCnXfeyVNPPcXjjz/O5ZdfHnZIEpKk3BmY2SYz+z8z\nW2NmK4Oy1ma22MzeN7NFZtYyZv3xZlZkZhvMbGAyYhCR+KxZs4bzzz+f9957j7Vr1yoRpLlkVROV\nAfnu3sfd+wVl44CX3b07sAwYD2BmeUQfj9kTGAw8bJr7VqTe7Nq1izFjxjB48GDGjBnDCy+8oBHF\nkrRkYFXsaygwO1ieDVwRLA8B5rp7qbtvAoqAfohIwhYtWsRVAwdy1cCBLFq0qMp13J158+aRl5fH\nvn37WLduHTfccIOeQyBA8toMHFhiZoeBx9z9CaCDu+8AcPftZlZ+6dEJeD1m261BmYgkYNGiRYy4\n8kqm7t8PwIgVK5j9/PMMGjSoYp3XX3+d22+/nS+++IKCggIuvPDCsMKVFJWsZHChu28zs3bAYjN7\nn2iCiJXQ8ysnTZpUsZyfn69+zyKVPD5tGlP372dEecH+/Tw+bRqDBg1iw4YNTJgwgbfeeou77rqL\n66+/XgPIGqFIJEIkEqnVPpKSDNx9W/DvP8zsBaLVPjvMrIO77zCzjsAnwepbgZNjNu8clFUpNhmI\npKNFixbx+LRpANw0duxRV/zV2b1nDzfccAMLFy7kl7/8Jc888wy5ubl1HaqEpPKF8q9//eu491Hr\nZGBmJwAZ7r7HzE4EBgK/BhYANwBTgRHA/GCTBcDTZnY/0eqhM4CVtY1DpDGqSRXQTWPHMmLFCti/\nnw+AezMyyN2wgf6XXcb777+vB9BIzbh7rV5AV2AtsAZ4FxgXlJ8EvAy8DywGWsVsMx74G7ABGHiM\nfbtIY1NYWOjDBgzwYQMGeGFh4THXHTZggM8C9+A1C3zYgAFHrVNSUuLf//73vUlWtmdmZPjll1/u\ne/furcufICkuOHfGdS6v9Z2Bu38I9K6i/DPgu9Vscw9wT22/W6ShqcmVfk33M2XKAxQXb2Lnzq3s\n3XuAsrKfAN9g6dI7eO211+Lep6Q3jUAWqUfHauytyk1jx3Ld8uVw8CAAv8jOpn+LFgwefDnuTjYt\nOcRhYAYEe92/P4tp0x5XMpC4KBmIhGzHzp1cNTA6EP+msWMBmDbtcQD69+/LF57Lv9GBUnZTemgX\nL720mKbuPMhhcviM0WSwm3dDi18aB4tWL6UmM/NUjk/STyI9eypvH1tNdGtODpllZUwrLQVgbFYW\n/8xoxsGDU4B3gCeBTOB8YBjgnJR1O/eVHrm7mA2MIot/8gQAubm38/zzs3VnkMbMDHePazShZi0V\nqaHyE/mQJUsYsmQJI668strRvpW3GzjwKgYOvAqAUXfcwcST2jHxpHa0b9OGaaWlXAfkAReVllJ6\nsITobC5bifa8vpfojC63AC1ompvzpe84/fQuDBiwgAEDFigRSEJ0ZyBSQ1cNHBhNAsH72cCCAQP4\nw+LFFessWrSooopn7NibALjyyhHs3z8VgJycW4FsDh6cChTThEl8jTI2AR2AU4DldOEA7wFNgV+Q\nkfE7ysruB6JX/XfcMYoH77qL/wraEf49J4enFixQApAKidwZqM1A0trxqn1iP/90584vbf/pzp0V\nV/z9+/dlypT/rjjxr1gxgh49erB//7WcyGwOs5sDB7OArwBjgbaU0IIN/JO7OUxrYDRZHGAHzTgH\ngIM5W7nzzl+wfPkCAMaOjV71n3POORVxPZVAdZXIl8TbF7U+X2icgdShwsJC75Cb67OC/vsdcnOP\n6vdf+fNWOTneJiur4n2brCzPyWnlMMthlmdktAmW9zq84XCtZ2XlegZ4DvgZwb8wxGFbMHRgrJs1\n82Z8zZvxNc/KOvGo72iXk3PcsQgilZHAOANVE0lKqm1DbU1UVe0z7fTT2f35bgBat27BrR98cNTn\nN1sWGd4dgD18BPwn0IPoeMupRKfg+oIMTqSMU2nZ5G/8vGQ/44GcYB+jLZvdPgM4Uu2zfPnbABzc\n+SE/WrPmmFVRIsejaiJJWfGc3JM1MCte7wIfffAB04P3t372DyB6et8OrAdKvS0H+TawDngPmAic\nA5wFdOVE3uJBIIt/Mpp36dr5dE794ANim3z79P46OW2Prva5447oZ+VdTEXqXby3EvX5QtVEjUJV\n1TGTJ0+udkqGmkzBUBOTJ0/2Lie18y4ntfPJkyd7YWGhDxgwzAcMGOaFhYU+efJkbxHsfxZ4S/DJ\n4H8Avxs8DzwDPBe8OXgmOJztMN3hZYefHVXF08KyvhT3hX36HLMqqibHStVEEi/CmI5C5Hgqj7p9\nd/9+pt15J/eXlQHRK/9Rd9zB28uXA1U31Falcs+dt956ixn3PQjAmX3PZuXLL1dc5Y/+1a/Yl9WS\n0tJpwD9Yvnw47dq1YTdncyNbKGMfzgEeInqd34Noz/6/0YEMeuDkcJgWZGQso6ysBVBMTs7vaV52\ngGml6wC4jS/flXdo25b/eP75irui2ce5Kxo0aBCz41hfJFnUZiBxi7c+v3Ld/AXAv0HF+18Av8vI\nqEgOt+bkkA0VXSdvz81l9vPPH3Wyz/+XyygoKKzouZOVNZITSvcyHSgJ9vkjoAvRGRFfBdaTjWNE\n++6fQJMmmykp+Rfgh0A34N9pwfNHEggECST6ncer36/8O8rj1slc6pvaDKROxJ78+/bvz39PmRJX\nfX7sFMsARRkZEJwwAf4M3F9WdmS+noMHmXb66UwMGnJH3TaGt956i3t/9SumE63DHzV7Nvv5AZAN\nfAClB+kITAA+C3azHMgn+rDtXOBDurCPdZQ35eblzeC9915m//4BwCfk5v6FK68ewcQXFwLwy9vG\ncM4558TcfRy7fv8sIO/ss1nQti2gq3ppYOKtV6rPF2oziGu647r4zsmTJx9Vh90mI6NG9flV1c+X\n19+PGDHC2+XkVOzzpODf2H02J8vhHofbPDu7uZ+Y3cS/AX52UH9/IngGTR2udhjvTcj128E3g5eC\nj4Wj2gNagGdknFBRv5+T08oLCwu/FGe8x0n1+5KKSKDNIPQT/jGDa2DJINkn7jBONpW/s/LJ//wq\nTtzlv7f8pDp58mTPze1Q0f8+J6eV5+S0O+p9ZmYLP4FunsupDtmeQ4ZfCN4N3MAh0+F0h+865Htu\nRpaPBH8L/PPge3PJqDi5Z2Q0+dLJ/7vf/W61CShZ/ffDSNYix6NkEKK6OHEno1dNVSerylfDse/7\n9+lzzJN/VVfcI0aMOOrkHx18NdPhI4dXg5P6MIcbHS5xaOKQ5fBVh286tHU435uQ5035RrD+TD/y\ntbP89NPzjvreE8BbZ2YedXIfMWLEUb2Hkn0sRRqKRJJBaG0GZnYp8ADRyfJmuvvUsGJJhnjnqU9U\n7PQH5XOWv89HAAAJiElEQVTfVJ4LJ3b649g5bK5bvpwxd9551JQJy5cPp6zMaVraGYAMNvCjmO+7\nkGhDarkZwG6GMZJPKaOE/XTk2Wf/xP79p5HFOMo4QFnZHuAmorPtnAp8ARwEzgW+B+wERgL/Guw1\nOv9OSTD/TnT+nnEcPBh9cHtu7u089NBs3nrrLSYGDcjdKg0I4+BBFnz8MZs+LX/UtojEI5TeRGaW\nAWwELgE+BlYBw939vUrreRjx1dSUKVMqerdUNVr1yT59yGnbFTj+ibu699cNGVJxMj8yvfEDQOyk\nZ/9V5fsW9hOme+lRMU1o3oqP//kAR/ry9KQFHzCdQwD8DMgkk1EcZhcwE6OEM8liN04Jh9kTbHci\n0V45GWRl/Y3M0gNcTxltgP8hkz00BR6qUZyVe+lUdSwqJ9WaTBoXq/JANvX0kcYskd5EYVX/nA+8\nFPN+HHB7Fesl884pqSoPWDqhiuqTrKyWR9WTZ2W1PGoOmsr16LHvc3M7+OTJkz0np1XFNmbNg8/L\nazvOdxjrzRjgzRjg0OOoz0+khT8M/mFQ1/4L8BMzsh2+49l09Sw6eSaZfgb46UGjbA64keUZtPRM\n2ju0dview5xgoNWP3KzVUXH2Pv30L1XBnPaVrxxzwFdtGm7dE6uWU/2+pAsaSpsBcBXweMz7a4Hp\nVayX/KOUJF1OavelE2BTulR7YoYe3oLsipNXc8zhVod9DrsdvuHwQFDPXuQw2Vu06OQwwWGJw3yH\nMxwGeg7dPYczHJp7Nhn+LfBvUD5CtkOwXlsHHPA24KeAZ2LetOkJngN+Bfj1QRKbAP4++C7wJ8Fb\nWPZxk1TsibxyffxY8NZBw3NNT9SJ0MldpGqJJIOUH2cwadKkiuX8/Hzy8/NDi6UqpwB7iNaIH6AY\n2EX0HLyH6NCqW4L3e9lHtJb8MFCGE20ymQ5YsE4RMI1o3/kD7N27EygAXidaLXOAbBbzbaAJ0ced\nXEW0rq0F8CYwnX+QSXOgDSWZ+zDLpqS0MyVAZs5WzuvZ9aiBUu2Ah83oFlTH/TInh1/eeWfMlMlz\nAarta18udhzBrGDgVV23nwwaNEjVPCJAJBIhEonUbifxZo9kvIhWExXGvG+w1UT/E7yag2dmNnf4\nH4eHPDu7hWdnt3F4xOFRb06mPw5eAn4I/LYqqpUq94Ov3EWzZUbOcXv6tDI7qofN8a7iZxGdP6e2\nV9ixV+mVeyTNUs8dkXpFA6omyiQ6S0AXosNB1wI9q1ivLo5T0tRkIrRqB1oFJ+1jnciP1+2zcnVM\nTQaE1cfYBQ3GEglXIskgtLmJgq6lD3Kka+lvqljHw4qvtqrqvVJ5MrbY6prK8/VU90jFePZZXQ+b\n+nhWQH18h4hUrcH0JqrpixS/MziW4w1yOt5I38rrx25XXZWOrshFxD2xO4OUb0BurCpPVfzz/v25\nfcqUikbY23NzmT12bJXbVXeVremPRSRRmsK6jiQyyElVKyKSDIlUEykZ1CGd3EUkDEoGIiKSUDLI\nqKtgRESk4VAyEBERJQMREVEyEBERlAxERAQlAxERQclARERQMhAREZQMREQEJQMREUHJQEREqGUy\nMLOJZlZsZm8Hr0tjPhtvZkVmtsHMBsaU9zWzd8xso5k9UJvvFxGR5EjGncF97t43eBUCmFlP4Gqg\nJzAYeNjMyidNegS40d27Ad3MrMFP5VnrB1HXg4YQIyjOZFOcydVQ4kxEMpJBVTPjDQXmunupu28C\nioB+ZtYRaO7uq4L15gBXJCGGUDWEP5CGECMozmRTnMnVUOJMRDKSwS1mttbMnjCzlkFZJ2BLzDpb\ng7JOQHFMeXFQJiIiITpuMjCzJUEdf/nr3eDffwEeBk5z997AdmBaXQcsIiLJl7SH25hZF+BFd+9l\nZuOIPpB5avBZITAR+Ah4xd17BuXDgf7u/rNq9qkn24iIJCDeh9tk1ebLzKyju28P3g4D/hosLwCe\nNrP7iVYDnQGsdHc3s11m1g9YBVwPTK9u//H+GBERSUytkgFwr5n1BsqATcBPAdx9vZkVAOuBQ8DN\nMc+vHAnMApoCC8t7IImISHhS+hnIIiJSP1J6BPKxBrWlAjO71MzeCwbQ3R52PNUxs01m9n9mtsbM\nVoYdTzkzm2lmO8zsnZiy1ma22MzeN7NFMT3UQlNNnCn1t2lmnc1smZmtCzp5jA7KU+p4VhHnqKA8\n1Y5nEzN7M/h/Zp2Z3R2Up9rxrC7OuI9nSt8ZmNlE4J/ufl/YsVRmZhnARuAS4GOibSDD3f29UAOr\ngpn9HfiGu38ediyxzOxbwB5gjrv3CsqmAp+6+71Bgm3t7uNSMM6U+tsMxvB0dPe1ZtYMWE10vM+P\nSKHjeYw4v08KHU8AMzvB3feZWSbwZ2AsMIQUOp7HiPO7xHk8U/rOIJCqjcj9gCJ3/8jdDwFzif5R\npyIjBf9bu/sKoHKCGgrMDpZnkwKDEquJE1Lob9Pdt7v72mB5D7AB6EyKHc9q4iwfa5QyxxPA3fcF\ni02I/v/zOSl2PKHaOCHO45lyJ4gqVDWoLRVUHliXygPoHFhiZqvM7F/DDuY42rv7DoieOID2Icdz\nLCn5t2lmpwK9gTeADql6PGPifDMoSqnjaWYZZraG6BiqiLuvJwWPZzVxQpzHM/RkEOegtpS5hWxg\nLnT3vsBlwMig2qOhSNV6zJT82wyqXp4DxgRX3pWPX0oczyriTLnj6e5l7t6H6B3Wt80snxQ8npXi\nvMjM+pPA8axt19Jac/cBNVx1BvBiXcYSp63AKTHvOwdlKcfdtwX//sPMnidaxbUi3KiqtcPMOrj7\njqB++ZOwA6qKu/8j5m1K/G2aWRbRE+xT7j4/KE6541lVnKl4PMu5+24zWwicQwoez3JBnH8CznH3\n5TEf1eh4hn5ncCzBwS4XO6gtFawCzjCzLmaWAwwnOtgupZjZCcFVGGZ2IjCQ1DqOxtF1mwuAG4Ll\nEcD8yhuE5Kg4U/Rv83fAend/MKYsFY/nl+JMteNpZm3Lq1bMLBcYAKwhxY5nNXGuTeR4pnpvojlE\n6xQrBrWV19elgqC71oNEk+pMd/9NyCF9iZl1BZ4nejubBTydKnGa2TNAPtAG2EF0ypIXgGeBk4lO\nX3K1u38RVoxQbZzfIYX+Ns3sQuBV4F2i/60dmACsBApIkeN5jDh/QGodz7OINhCXd754yt1/a2Yn\nkVrHs7o44z53pnQyEBGR+pHS1UQiIlI/lAxERETJQERElAxERAQlAxERQclARERQMhAREZQMREQE\n+P8LK7iuOGpAeQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "clf = Ridge(alpha=100000) #Is the alpha that affects the weights!!\n", "clf.fit(x_5, y_n)\n", "print clf.coef_\n", "predict_y_n = clf.predict(x_p)\n", "plt.scatter(x, y, c ='blue')\n", "plt.scatter(x, y_n, c ='red')\n", "plt.plot(x_predict, predict_y_n, c ='black', label = '3rd Degree')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.11" } }, "nbformat": 4, "nbformat_minor": 0 }