{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import sklearn" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Linear regression with one variable" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "data1 = pd.read_csv('ex1data1.txt', header=None, names=['x', 'y'])" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "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", "
xy
06.110117.5920
15.52779.1302
28.518613.6620
37.003211.8540
45.85986.8233
\n", "
" ], "text/plain": [ " x y\n", "0 6.1101 17.5920\n", "1 5.5277 9.1302\n", "2 8.5186 13.6620\n", "3 7.0032 11.8540\n", "4 5.8598 6.8233" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data1.head()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeMAAAFwCAYAAACLhRPVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHRdJREFUeJzt3X2sZHd93/H3tyxBCq60ICrbRbav/4n6YIMtCEgsyNtW\ntpysRYkUrXCa1BApclQLQ1TtAtko3rZZ0ezKyGlLaJTYFsXg1kqKE9sKshN5A1kUnoQfFpzSSl43\ntPbSNGxr8wdB1bd/nDM7M/fO3Hvm6fzOmXm/pNG9c+7cmTM7O/czv6fvLzITSZJUzt8ofQKSJG06\nw1iSpMIMY0mSCjOMJUkqzDCWJKkww1iSpMIahXFEXBERT0bENyPibETcWR8/HhHfiYhv1JebV3u6\nkiStn2iyzjgiLgMuy8ynIuIS4OvAe4DDwMuZ+fHVnqYkSetrX5MbZeZLwEv1969ExHPAG+sfx4rO\nTZKkjTDzmHFEbAHXA39WH/pARDwdEfdGxP4lnpskSRthpjCuu6h/F/hgZr4CfBK4GrgOeBG4e+ln\nKEnSmms0ZgwQEa8GHgX+MDPvmfDzLeCRzLx223GLX0uSNkpmzjSE22jMOCICuBf41mgQR8Tlmfli\nffWngGeXcVLqjog4npnHS5+HZudr12++fv01TyO0URgDB4CfBZ6JiG/Ux34ZuDUirgMSeB64fdYT\nkCRp0zWdTf2nTB5f/sPlno4kSZvHClzay+nSJ6C5nS59AlrI6dInoPY0nsA19wNEpGPGkqRNMU/u\n2TKWJKkww1iStDwRh9heACpiPxGHCp1RLxjGkqRlOgOcuBjI1dcT9XFN4ZixJGm5hgF8CjgCHCPz\nQtmTas88uWcYS5KWr6rK+DxwNZnnip5Ly5zAJUkqr2oZH6Hau+DIjjFk7WAYS5KWZ9hFfaxuER9j\ndAxZE9lNLUlanmrW9JmxMeIqiA+Q+Vix82qRY8aSJBXmmLEkST1kGEuSVJhhLElSYYaxJEmFGcaS\nJBVmGEuSVJhhLElSYYaxJEmFGcaS1JR79WpFDGNJas69erUSlsOUpFls+F692pu1qSWpDRu8V6/2\nZm1qSVo19+rVChjGktSUe/VqReymlqSm3KtXDThmLElSYY4ZS5LUQ4axJEmFGcaS1HVW/lp7hrEk\ndZ+Vv9acE7gkqQ+s/NUbzqaWpHVm5a9ecDa1JK0rK3+tNcNYkrrOyl9rz25qSeo6K3/1imPGkiQV\n5pixJEk9ZBhLUh9Y+GOtGcaS1A8W/lhjjhlLUl9Y+KMXnMAlSevOwh+d5wQuSVpnFv5YW4axJPWB\nhT/Wmt3UktQHFv7oDceMJUkqzDFjSZJ6qFEYR8QVEfFkRHwzIs5GxJ318ddHxBMR8e2IeDwcu5Ak\naWaNuqkj4jLgssx8KiIuAb4OvAd4P/CXmXkyIj4MvC4zP7Ltd+2mliRtjJV1U2fmS5n5VP39K8Bz\nwBuBdwOfqm/2KaqAliRJM5h5zDiqBefXA18GLs3M8/WPzgOXLu3MJEnaEDOFcd1F/XvABzPz5dGf\nZdXfvdqp2ZIkraF9TW8YEa+mCuJPZ+bD9eHzEXFZZr4UEZcD353yu8dHrp7OzNNznq8kSZ0SEQeB\ngwvdR8MJXEE1Jvy/M/OXRo6frI/9ekR8BNjvBC5J0iZbWdGPiHgn8AXgGYZd0R8FvgI8BFwJnAMO\n57YdRAxjSdImsQKXJEmFWYFLkqQeMowlSSrMMJYkqTDDWJKkwgxjSZIKM4wlSSrMMJYkqTDDWJKk\nwgxjSZIKM4wlSSrMMJYkqTDDWJKkwgxjSZIKM4wlSSrMMJYkqTDDWJKkwgxjSdpUEYeI2L/t2H4i\nDhU6o41lGEvS5joDnLgYyNXXE/VxtSgyc7UPEJGZGSt9EEnSfIYBfAo4Ahwj80LZk+q3eXLPMJak\nTRexBTwPXE3muaLnsgbmyb3+dlM71iFJi6v+jh4BrgaO7Pi7qlb0N4wd65CkxQz/bh6rW8THGP27\nqtb0u5vasQ5Jml/Vk3hm7O9m9Xf1AJmPFTuvntvMMWPHOiRJHbJZY8bgWIckaS30N4wd65AkrYn+\nhjEcYHSMuPp6rD4uSavhSg6tQP/HjCWpTeO9chd2XNfG28wJXJLUNldyaBeGsSS1xZUcmmLzZlNL\nUgmu5NCSGcaSNAtXcmgF7KaWpFlYtUp7cMxYkqTCHDOWJKmHDGNJkgozjCVJKswwliSpMMNYkqTC\nDGNJkgozjCVJKswwliSpMMNYkqTCDGOpC9ywXtpohrHUDWcY3WxguBnBmZInJakd1qaWusIN66W1\n4EYRUt+5Yb3UeyvbKCIi7ouI8xHx7Mix4xHxnYj4Rn25edYTljTCDeuljdV0zPh+YHvYJvDxzLy+\nvnx+uacmbRA3rJc2WqMwzswvAt+b8CO7n6XlOMDoGHH19Vh9XNKaW3Q29Qci4umIuDf8BC/NL/Ox\nHZO1Mi+Q+VihM5LUokXC+JNUY1vXAS8Cdy/ljCRJ2jD75v3FzPzu4PuI+B3gkWm3jYjjI1dPZ+bp\neR9XKqYqwHFmrAVb9QgdsAUrba6IOAgcXOg+mi5timrJxSOZeW19/fLMfLH+/peAH8/Mn5nwey5t\n0noYn2R1Ycd1SWKF64wj4kHgBuANwHngLqpPAddRzap+Hrg9M88v46SkzrIwh6Q9WPRDaoOFOSTt\nYmVFPyTVLMwhaQUMY6kpC3NIWhG7qaWmnE0tqQHHjCVJKswxY0mSesgwliSpMMNYkqTCDGNJkgoz\njLsm4tCOpTIR++uZvJKkNWQYd88ZRteuDte2nil5UhvDD0OSCjCMu2a4qfyJuuyiGxG0yw9Dklrn\nOuOusv5xOW4GIWkBrjNeF9Y/LqsK3lNUH4ZOGcSSVs0w7hrrH5fnhyFJLbObumusf1zW+IehCzuu\nS9IerE0tLcoPQ5IWZBhLklSYE7gkSeohw1iSpMIMY0mSCjOMJUkqzDCWJKkww1iSpMIMY0mSCjOM\nJUkqzDBeN+7Hu5l83aVeM4y7aLE/rO7Hu5l83aUeM4xLmxS8cBZ4YK4/rFVN5cFOT1u4ycFm8HWX\nes3a1AOlNgiYvkvQSeAo825wX/1Bfh64ut6KUZvA110qztrUiynTzTe9RfMC825w7368m8nXXeot\nw3igZDdf9RjjwTvvH9bxlvY5hs/JP8zrzNdd6jW7qbcr0c03/EM66JIedFHPvsG9+/FuJl93qTPc\nz3hRO0Nx9S3jyWPGDwB31F3Vo7fzD6skdZxhvIjpE6lWG8i2aCRprRjGizAUJUlLYBhLklSYS5sk\nSeohw1iSpMIMY0mSCjOMJUkqzDCWJKkww1iSpMIMY0mSCjOMJUkqzDCWJKkww1iSpMIMY0mSCmsU\nxhFxX0Scj4hnR469PiKeiIhvR8Tj4SbmWgcRh9j+fzlif72RiCStRNOW8f3AzduOfQR4IjN/DPjj\n+rrUd2eAExcDebiV5pmSJyVpvTXetSkitoBHMvPa+vqfAzdk5vmIuAw4nZl/Z8LvuWuT+mUYwKeA\nI6x6T2tJa2We3Nu3wONdmpnn6+/PA5cucF9Sd2ReIOIU8DxwtUEsadUWCeOLMjMjYmoTOyKOj1w9\nnZmnl/G40kpULeMjwNXAESJsGUuaKiIOAgcXuo8Fu6kPZuZLEXE58KTd1Oq9YRf1sbqFPH5dkvYw\nT+4tsrTpD4Db6u9vAx5e4L6krjjAaPBWX4/VxyVpJRq1jCPiQeAG4A1U48O/Cvw+8BBwJXAOOJwT\nWg62jCVJm2Se3GvcTT0vw1iStEna7qaWJElLYBhLklRYN8LYEoSSpA3WjTC2BKH6yg+SkpagG2E8\nXD5ygmo9s+s61Rd+kJS0sG7Npq6CeFCC8NzqzkpaImtZSxrR79nUO0sQuiWj+qEK3kEt61MGsaRZ\ndSOMx0sOnmPYZW0gq/v8IClpQd3opq4mu5wZa1FUf9AOkPnYSk9QWoS1rCVtYwUuqW1+kJS0jWEs\nSVJh/Z7AJUnShjKMJUkqzDCWJKkww1iSpMIMY/WHdaAlrSnDWH1iHWhJa8kwVjNdaJW6oYikNWUY\nq6lutEqtAy1pDRnGaqYrrVLrQEtaQ1bg2mTzlHIsuc2ldaAl9YAVuDSr2bqey7dKDzAavMPW+oGW\nz0OSlsqW8aYbBvApqqCd3Mq0VSpJjbhRRF+V3vmnSddz6XOUpJ5Yz27qLiypWb1yM5Wbdj1nPraj\nBZx5wSCWpMV1P4y7sqRmlUrNVB7vaj43cg7LHwvejA9VkjSXfnRTNx3X7Lu2Zyq32fXsmLOkDbHe\nY8Yll9S0YRM+cGzCc5S08dZ5zPgqxsc1r1qr7s1ldBf3oRvY6lmSNFH3wxjOAo8CJ+ugOllfP1vy\npHYJv+NzhOIy1s+uZmx99HkOvh99PrMEfvl1ypLUSX0I42uAW4CjdVf10fr6Nbv+1upbitPC7/4p\nx6eH4jJmKq9uEtjo8zxD1bI9BZyZKfDbnCwmSX2TmSu9VA+xhPuCrYRM2Gp4+/0Jn0jYP/H6cs5p\ncJ9bUx5ra+JjwqEJx/YnHGr932n25/nb9WXyc5t+H6t7zl68ePHSocs8udePCVzzTvxpY8LQtIll\nu004W9XM4lU+39HnU1nfyXSStIB1ncA1f/fmqicMTRsD3WtsdBVdyqvsBh5/Psfqi+O+krQk3Q/j\nRSY3rXLC0PTwu2rK8UmBvMwPCqvZRGH0ecLoOY5+oDCQJWkB/eimnu+BV1tkYnrBjA8B90w4Pl5I\noy9rbkef5+D7SvV8rE8tSWPWu+jH7A/c3Y0NrEYlSWvLMJ58At0L5S6ekyRpKQzjySdgK1SS1BrD\nePpJ9GN8VpLUe91f2lSqVnJXayL3oZ60JGnl2gvjkvsQd7cm8vrv1SxJ2lM73dSDECzRPdz1MWO7\n0CVprXR3zLj6dvWlEyfPUj4MQOZDI8e6NXN53fdqlqQN0uUx47a6h7d3+x4GbgQev3iL3YK4xBju\nqrvQ12lcep2eiySNaCeM29sy7wDVfseDms83Aj8C3AQ0GZNtdwy3nW0F12lcep2eiyRdtHA3dUSc\nA/4v8P+AH2bm27b9fNhcX3X38LAl/G+Bp4E3A/8cuBT4RZqMybY5hjtL8Y9FCoWs07j0Oj0XSWup\nyJhxRDwPvCUz/6rRSa0ykKv7/g3grcDtwG8BXwPupgrnZmOyXRzDXXQiWhef07zW6blIWjslx4x3\nf9B2uxX/GngG+GL9FeADNB2TnWUMt80xzEW2Xezu0q7ZrdNzkaTaMsI4gT+KiK9FxC9Muc3y9u3d\n3QHg14Afbjv+RKMx2fHtAv8+w/Hn4YeJ8aBtdwxznuIl7YxLt2OdnoskjcrMhS7A5fXXvwU8Bbxr\n288zYSsHXxd8vF0vsD/ht+vLVsLnEz6VsH/bbQ5N+f1DF29b3e4TCVddPF5d3z/hMT9RP97Ony//\n+c32WKPPqcm/QZcv6/RcvHjxsraXKlpn+52lrjOOiLuAVzLz7pFj+XPw1SfhS/8A3vFV+JfPZT66\ntAcdP4HBBK4jDMdVT1G1jB/a/Zcn3l+zyUJtjGF2vXiJJG2oiDgIHBw5dFe2OYErIn4UeFVmvhwR\nr6Vaz/svMvPxkdtkwutaCZBVbE24V9C2NbvXbRclqRdan00dEVcDn6uv7gM+k5kf2/Wk9iq60aXA\n2Stod2utVuPX3XkukqRWdLccZtOT6lJXbJNz2e3Dw2ByVxeeiySpNf0P4+oXulHUYXrQfgi4p1GL\ntyvPRZLUmn6G8eTQexOzFOloS3WuZ4GjDFu87wPeC7x3pMU8DGcLVEjSRunyRhG72b5W9yrgM1Sl\nLLtW1OEMVRAP1h+/iarFe/7iLUbXGs9aoMKNECRpI5UP4/HKUm8CHgVuIfMZulbUYXiuR4EHqVrv\ntwIfZHthk8pogYo/AU6NPZfSRUQkSZ1Qvpt6eMMtqu7cN9dBPDhefgby9q70YTf6UWCLKnz3M9od\nvfN3xtc8T5vQ5TizJPVaX7upt9cbvn2s9Zh5YWVB3LxbeNhirbrRHwT+E/BjVF3Wp6gCedgdnfnY\nWIhW3x8Bbti1NOg8JS/nf16SpA5oN4wnBULZesPNuoWH3dOngPuodoL6z1R1sH+Fas/kJ4ALVN3R\nk8+/SdAuZyMEu7slqU9aq9E5vbZz2XrDs9R7hvfnsNb24Pfel3B47PlNO/+9Hmv7v9G0f7NlPy8v\nXrx48bK0C6VrU09Sl8OsWnltjX/OWsmryfKjSWO5lWbju4sWEZmnq95lVZLUui6PGc8//jmf5t20\nTbqFp3WlV5qO7x5gNHiHXd8HLt5i+zjz4HbzBbH7/kpST7QVxuOBsOoJRuPLpbaYNlmq+Xj1tCC9\niaaBt8yg3Yv7/kpSr7S3tGl8EwXY2WX7AHAHmS+M/PKqd1yav1u4S3W0d55btzbckKQN0v1ymONl\nIrePwZ5kvMzkYuG2jPW6u28EwdSfGXiStLG6H8Y7f7jFeKGMQYB+DXgHcGRkclPzoFtWq7XLrV9J\nUid1eQLXTpMmGA3X4d4HXLLttqeA1064n53jz9VY7p/sOlmqiaZjz5IkLaBMy3hyi/MB4KPALwK/\nxbD28xNUrWQYbSnvfl/LDU2XCEmSGupPN/XksdhrgUeAG8h8ob7+Baqaz58G7mywhnf59ZytFS1J\nmkF/uqknLfOBK4EbgKN1S/RO4I8b3t/i9ZwncYmQJKkF5WtTD1QB/QLDUL0E+B7VmPIP2L794Pb7\nXU2Bi70LdUiStKAy64z37m7+FvDPgJ+su6zHtx+c/DvOeJYkFdfdbuomM5HHi4KcA36Sqst6MMv6\nCPD9Cb9p61WS1GttbRQBe81EXqwalhWnJEmd0N2W8Ty1mwdjy6NhOn3M2f17JUm91U4YzzcTuXnA\nWpxDktRj5WpTN/vl2db4WpxDklRYf4p+zHYHWzQJWItzSJI6oMtjxkOz7FvcdP2wxTkkST3W1tKm\n0XHfB4CzO36+PaBnC1iXN0mSequtpU2/ybD7+KtUIXlkpEjHeEGPKpgHOzQ9Xt/+LPD2+tj3qSZy\nuXRJktQpXe6mHtaNhocvHqvGg0/V1x8fuf0ZqjrVX6ZqHf934FHgufr4WVy6JElaE22F8XDcd3j9\nNVQB/Rq2b4047GY+SrWd4meptlb8bH39KE7QkiStiX0tPc6gbnQ17gsn9/yN6vaDFvW7gD+tvz5N\nNbPaIJYkrYW2in6MTqw6SdW6/QG77cg0nEn9ZuCTwDvrr29muTszSZJUVIn9jN8O/AVV1/Q5hl3X\nN128xXAm9UngduBngH9ff729Pr5zZnXEoYmh3nQplSRJBZQI4+8zOkY8eUemaqkSXFN/vRK4pf46\nenz70iVrVEuSeqf7FbjG72zv3ZmsxCVJKqjLS5uWZe+WbxW8w6VUBrEkqeP6FcZNdmdqWkJzGsed\nJUktKx/GVfgdHgvAKvwOTwzA3Vq+y6lR7bizJKlV5cO4CrkbGSxvGpbHvJHtAVgF91WMt3yvGgnt\nxWtUuzeyJKll3ZjANQzg19RHfsD2qlzV7a6iKot5C5kv7Li+3BPfwr2RJUkz6u8Erip0TwA/V19O\nXAzi8fHaa6iWOB2tw/Joff2a+rbLGe9ddNxZkqQZlG0ZD5YqVU4BlwBvodoY4qfr45MmaW0xqdU6\nPmZ8Ycf1Zie8+H1IkjbWPC3jtrZQfN2EiVYfAu4HfgX4EeCvgX8D3AP8beAp4BW2d1fvtY540XXG\nTdYyS5I0RZfD+DfZ2dI8SdXN/FXgvcC/oip5eUt9+XfAz5N5f31H+6lKZt4w4b6atZwlSVqxLo8Z\nT5qd/EJ9/MeBX6PalemfAP8HuA74NPCOkRnWJy7e126zpR3vlST1THtjxtPHed9EtS3iu4Db6qOD\nzSMGM6wnz67e+WCTxnsfAO4Ym21tt7MkaUWKtIwj4uaI+POI+K8R8eEpNzpM1YIdtlary/uAzwB3\nUAXx9sIaX6KaXf2lhuO+k9YZ3wF8onNFPKz0JUmqLRTGEfEqqrHdm4G/B9waEX93wk1vrL8OupVP\nAb9BNVZ8C/DZ+ucHqLqsb6pv8w6qAH9ro+7mzMd2hHbVIv5ZulfEw0pfkiRg8Zbx24D/lpnnMvOH\nwH8E/vGE2x2pLyeAQai+CLyXzBdGtlEE+IcMw3uw5/E8ZS2Hurh5hJW+JEm1fQv+/huBvxi5/h3g\n7TtuNSzgMQjEnbOcqzHeE/XPfx743Fh3c8Rgotbs47w7J3V1I/Sq5zX6b1L+nCRJrVu0Zdx89tde\ns5zHf/7WnY+UF+aacLWczSNWw5nfkiQWD+P/AVwxcv0KqtbxmP0R//qfwuP74eWALbYH4moDc/HN\nI1ahyx8SJEmNRcTBiDg+uMx1H4ssbYqIfcB/Af4R8D+BrwC3ZuZzI7fJrCZpTa9qtYlVrzbxOUvS\nBihSgSsifoKqhOWrgHsz82OLnpQkSX3V3XKYhrEkaUN0uRymJEmaop0wttqUJElTtdUyttqUJElT\ntLlRxGL7DEuS1APdn8DlPsOSpDXX7QlcVpuSJGmitiZwWW1KkqQp2moZVyUp4QAR+8dKUjqrWpK0\n4doJ4+E+w8NZ1aPXnVUtSdpg7Rb9WPUevq5nliT1UPsVuKrgHezhe2rJy5tczyxJ6p32w3iVs6pX\n3fKWJGkF2l5nPDqr+sKO68t70C1czyxJKqDb64wr1azqQfCOzqpeFtczS5J6Zr22UGyr5S1J0hTd\nL4e5atWs6TNjwVsF8gEyH2vlHCRJG80wliSpsD6MGUuSpG0MY0mSCjOMJUkqzDCWJKkww1iSpMIM\nY0mSCms3jN1BSZKkHdoLY3dQkiRpon2tPEq1ccMRLEspSdIO7VTgqr51ByVJ0trrcgUud1CSJGmK\ndsK4ahEfA04YyJIkjWt3owh3UJIkrTl3bZIkqbAujxlLkqQpDGNJkgozjCVJKswwliSpMMNYkqTC\nDGNJkgozjCVJKswwliSpMMNYkqTCDGNJkgozjCVJKswwliSpMMNYkqTCDGNJkgozjCVJKswwliSp\nsLnDOCKOR8R3IuIb9eXmZZ6YJEmbYpGWcQIfz8zr68vnl3VS6o6IOFj6HDQfX7t+8/XbLIt2U8dS\nzkJddrD0CWhuB0ufgBZysPQJqD2LhvEHIuLpiLg3IvYv5YwkSdowu4ZxRDwREc9OuLwb+CRwNXAd\n8CJwdwvnK0nS2onMXPxOIraARzLz2gk/W/wBJEnqkcycaRh337wPFBGXZ+aL9dWfAp5dxglJkrRp\n5g5j4Ncj4jqqWdXPA7cv55QkSdosS+mmliRJ81tpBa6IOBcRz9RFQb6yysfSYiLivog4HxHPjhx7\nfT2J79sR8bgz5rtryutnYZ4eiIgrIuLJiPhmRJyNiDvr477/emCX12+m999KW8YR8Tzwlsz8q5U9\niJYiIt4FvAL8h8FEvIg4CfxlZp6MiA8Dr8vMj5Q8T0025fW7C3g5Mz9e9OS0q4i4DLgsM5+KiEuA\nrwPvAd6P77/O2+X1O8wM7782alM7gasHMvOLwPe2HX438Kn6+09R/QdTB015/cD3X+dl5kuZ+VT9\n/SvAc8Ab8f3XC7u8fjDD+2/VYZzAH0XE1yLiF1b8WFq+SzPzfP39eeDSkiejuViYp0fqZaLXA1/G\n91/vjLx+f1Yfavz+W3UYH8jM64GfAO6ou9LUQ1mNZzjbr18szNMjdRfn7wEfzMyXR3/m+6/76tfv\nd6lev1eY8f230jAerEPOzP8FfA542yofT0t3vh4PISIuB75b+Hw0g8z8btaA38H3X2dFxKupgvjT\nmflwfdj3X0+MvH4PDF6/Wd9/KwvjiPjRiPib9fevBW5iSmEQddYfALfV398GPLzLbdUx9R/wgamF\neVRWRARwL/CtzLxn5Ee+/3pg2us36/tvZbOpI+JqqtYwVMVFPpOZH1vJg2lhEfEgcAPwBqrxqV8F\nfh94CLgSOAcczswLpc5R0014/e6i2vVnrDDPyBikOiIi3gl8AXiGYVf0R4Gv4Puv86a8fr8M3MoM\n7z+LfkiSVFgbS5skSdIuDGNJkgozjCVJKswwliSpMMNYkqTCDGNJkgozjCVJKswwliSpsP8PFSgx\n0vcaHPkAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(8, 6))\n", "plt.plot(data1['x'], data1['y'], 'rx')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from sklearn import linear_model" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 1.19303364]\n", "-3.89578087831\n" ] } ], "source": [ "clf = linear_model.LinearRegression()\n", "clf.fit(data1[['x']].values, data1['y'].values)\n", "\n", "print clf.coef_\n", "print clf.intercept_" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeMAAAFwCAYAAACLhRPVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+05Hdd3/Hne1lAIDabNJikIeTuqXoE+ZFUSZQFs/5I\nGs6mFFsag0URrY2CiBY2gLGSqjnqXsPBtoCeGjDyI5qjgiapNKGHKC5VwCaQIBR7zEYoyUYiK0nk\nl+HdP77fufOduTN35+f3x8zzcc6cvfO5c2c+k5t7X/fz4/v+RGYiSZKas6vpDkiStO4MY0mSGmYY\nS5LUMMNYkqSGGcaSJDXMMJYkqWEThXFEnBkR742Ij0bEnRHx42X7lRHxqYi4rbxdtNzuSpK0emKS\n64wj4jTgtMy8PSJOAP4ceB5wCfBAZr5uud2UJGl17Z7kQZl5L3Bv+fGDEfEx4Izy07GkvkmStBam\nXjOOiA3gHOBPy6aXRcSHI+KaiNizwL5JkrQWpgrjcor6d4CXZ+aDwJuAvcDZwD3A1QvvoSRJK26i\nNWOAiHgkcCPwh5n5+hGf3wBuyMynDrVb/FqStFYyc6ol3InWjCMigGuAv6gGcUScnpn3lHe/G7hj\nEZ1Se0TElZl5ZdP90PT83nWb37/ummUQOlEYA/uAFwIfiYjbyrafAl4QEWcDCdwFXDZtByRJWneT\n7qb+E0avL//hYrsjSdL6sQKXjufWpjugmd3adAc0l1ub7oDqM/EGrplfICJdM5YkrYtZcs+RsSRJ\nDTOMJUmLE3GA4QJQEXuIONBQjzrBMJYkLdJh4KqtQC7+vaps1xiuGUuSFqsfwJvAQeAKMo8126n6\nzJJ7hrEkafGKqox3AXvJPNJoX2rmBi5JUvOKkfFBirMLDm5bQ9Y2hrEkaXH6U9RXlCPiK6iuIWsk\np6klSYtT7Jo+PLBGXATxPjJvaqxfNXLNWJKkhrlmLElSBxnGkiQ1zDCWJKlhhrEkSQ0zjCVJaphh\nLElSwwxjSZIaZhhLktQww1iSJuVZvVoSw1iSJudZvVoKy2FK0jTW/KxeHZ+1qSWpDmt8Vu8kIjgB\neA/wo5nc1nR/6mZtaklaNs/q3VEEvww8AJwH/F3D3ekMw1iSJuVZvWNF8G0RJPAK4NcyiUz+qul+\ndYXT1JI0Kc/q3SaCE4GjwKOBLwGnZrLWa+iuGUuSahPBG4CXlHe/PZNbG+xOa7hmLElaugi+s5yS\nfgnw+nJK+taGu9Vpu5vugCSpGyI4Gbi/vPs54AmZPNBgl1aGI2NJaruGK39FEBFcQz+In5XJiQbx\n4hjGktR+jVX+iuA5wFeAHwR+sZyStuLYgrmBS5K6oObKXxE8HrivvHsv8LWZPLSs11sl7qaWpFVW\nQ+WvCAJ4B3Bp2XRuJh9cxmutKndTS9KqqqHyVwTPo5iSvhS4spySNohr4MhYktpusPLXsW335356\nzoatGtJHgCdn8vl5n3ddOU0tSatoSZW/ItgFPFxpOieT22fupwDDWJI0oQiuB/5NefezmZzcZH9W\nySy5Z9EPSVojEZwH/Gml6cRMPtdUf1RwA5ckdcGchT8i2F2WsOwF8feUG7QM4hYwjCWpG2Yu/BHB\nu4Evl3ePlCF8/bI6qum5ZixJXTFl4Y8IzoeBAxwel8nfL7WPcgOXJK28CQp/RPAo4IuVpudmcsPS\n+ybAoh+StNomKPwRwfvpB/GHyylpg7jlHBlLUhccp/BHBP8ceHflK74qc2B0rJo4TS1Jq2pM4Y9P\nc/r+M/j0OyuPvCCT99TeP20xjCVpjUTwUeDJ5d33ZfJtTfZHBYt+SNIaKA90qI6GH5W5demSOmii\nDVwRcWZEvDciPhoRd0bEj5ftJ0fELRHxiYi4OZZwiogkqRDBnrJwRy+In11u0DKIO26iaeqIOA04\nLTNvj4gTgD8Hnge8GPhMZh6KiFcBJ2Xmq4e+1mlqSZpTGcJbMvH3akst7dKmzLw3M28vP34Q+Bhw\nBvBc4NryYddSBLQkaUEi+IOhIH6UQbx6pr7OOIoLzs8B/gw4NTOPlp86Cpy6sJ5J0hqL4MwyhP9F\n2fSjTkmvrqk2cJVT1L8LvDwzH4jo/3GWmRkRy92aLUlrwCnp9TNxGEfEIymC+K2Z+a6y+WhEnJaZ\n90bE6cB9Y772ysrdWzPz1hn7K0krq6ye9a2Vpt2ZPNxUfzSZiNgP7J/rOSbcwBUUa8L3Z+ZPVtoP\nlW2/FBGvBva4gUuSphPB1wGfqDS9MJO3N9UfzWdpRT8i4lnAHwMfga3pk9cAHwCuB54IHAEuyaET\nRAxjSRrPKenVYwUuSeqICP4S+NpK065M3HezAjy1SZJaLoKnlaPhXhA/t9wlbRCvMcthSlJNnJLW\nOI6MJWnJIvi7oSDeZRCryjCWpCWJ4FvKEP5HZdN3OiWtUZymlqQlcEpa0zCMJWmBDGHNwmlqSVqA\nCL5rKIjPM4g1KUfGkjSHCAL4SqXpM5k8vqn+qJsMY0makVPSWhSnqSVpShG8eiiIn2oQax6OjCVp\nQhHsgoFTlD6eyZOa6o9Wh2EsSRNwSlrL5DS1JO0ggquHgvgZBrEWzZGxJI0QwSOBL1XbDGEti2Es\nSUOcklbdnKaWpFIEbx0K4icZxKqDI2NJay+CxwB/X21bixCOOAAcJvNYpW0PsI/Mmxrr1xpyZCxp\nrZUj4a0gLk9VWv0gLhwGrioDuBfEV5XtqpFhLGktRXDz0JT0WWsUwoViRHwFRSBvUATxFQMjZdUi\nMpd7rGZEZGau1//gklorghOBatjcm8npTfWnFYogvgvYS+aRRvuyAmbJve6OjCMObE2t9Nv2lGsg\nkrRNORLeCuJySnrdg3gPcBDYCxzc9ntVtehuGLvWIWlCEdwxNCX9NWs3JT1K//fmFeWIuDdlbSDX\nrNvT1P3/kTYp/rJzrUPSlghOBe6tNN2WyT9rqj+t427qpZgl97odxsULbOBah6QhFu5QU9ZrzRhc\n65C0TQT3DwXxiQax2q67Yexah6SKCDbKED65bPrv5QatzzXZL2kS3Z2mdq1DUqnWKWl/9+g41nPN\nWNLaGg5h4DGZfGHJL1qdlTu27b7W3vqtGUtaSxF841AQv6Wckl5uEINVq7QUjowldUprdkl7JYfG\ncGQsaWVFkENB/MgGg9grObRQhrGkVovg3KEQPlROSf9DQx3ySg4tnNPUklqrNVPSVe6m1nG4m1rS\nShixS3pX5rY2qZVcM5bUaRFcMBTEB8spaYNYK2130x2QJGjplLRUE8NYUqOckpacppbUkAiePxTE\n/84paa0rR8aSaueUtDTIMJZUG0NYGs1paklLF8EPDwXxvzKIpT5HxpKWJoIAvlJtM4Sl7QxjSUvh\nlLQ0OaeppTaIOLCttnHEnrL0YqdE8KqhIP4ug1jameUwpTZYgQPrI9gFPFxtM4S1jqxNLXVZP4A3\nKY7n61IQOyUtlaxNLXVZEbybFAfWb3YhiCP45aEgPtcglqY3URhHxJsj4mhE3FFpuzIiPhURt5W3\ni5bXTWkNdOjA+gh2lyH8il5bWT3rgw12S+qsSUfGbwGGwzaB12XmOeXt3YvtmrRGOnRgfRnCX+7d\nL0PY0bA0h4nCODPfB3x2xKf8AZQWYx/VNeLi3yvK9laI4NqhKeknG8LSYsx7nfHLIuL7gQ8Br8gO\nrHFJrZR504i2Y8D29ppF8Bjg76tthrC0WPNs4HoTxdrW2cA9wNUL6ZGk1ihHwltB7JS0tBwzj4wz\n877exxHx68AN4x4bEVdW7t6ambfO+rpSY4oCHIcHdjkXa7r7Ro5sOyyC/wFcWGnayOTupvojtVlE\n7Af2z/Uck15nHBEbwA2Z+dTy/umZeU/58U8Cz8jM7x3xdV5nrNWwAoU5jieCE4Hqe7kvk1Ob6o/U\nRUsr+hER1wHnA6cAR4HXUvwVcDbFruq7gMsy8+giOiW1VocLcxyPhTukxbACl1SHYpboLmBveRlS\np0XwYeBplaZTM7lv3OMl7cwKXNKydagwx/FE8DXlaLgXxLeXG7QMYqlmjoylSa3QmrFT0tLyOE0t\nLdMK7KaO4D7g8ZWmEzP5XFP9kVaRYSxppAjOAo5Umt6dyXMa6o600mbJvXkrcElqOaekpfZzA5e0\noiLIoSB+jEEstZNhLK2YCJ48FMK/Ue6S/kJjnZK0I6eppRXilLTUTY6M2ybiwLZrVyP2lDt5pZEi\neGAoiB9pEEvdYRi3z2Gqh8r3r2U93GSn1kbH/hiK4KllCJ9QNr2ynJL+hyb7JWk6XtrURitc/7j1\nOlTYwylpqZ28zniVrFj9405p+R9DwyEM7Mrc1iapIdamXhUrVP+4k4rg3aT4Y2izLUEcwbcMBfEP\nlVPSBrHUcYZx2wxOix4BrqC6hqzla+EfQ2UI/6/e/TKE39xglyQtkNPUbbMC9Y87rWVrxk5JS93j\nmrE0r5b8MRTBBcDNlabnZ/K7db2+pNkZxtIKcJe01G0eFCF1mCEsrS83cEkNi+CHhoL4QoNYWi+O\njKWGRBDAV6pthrC0ngxjqQFOSUuqcppaqlEElw8F8QUGsSRHxlINItgFPFxtM4Ql9RjG0pI5JS3p\neJymlpYkgkNDQXyeQSxpFMN41XTsPN5VFMHuMoQP9trKWtIfWOKL+n2XOswwbqP5frEepnqwRL+2\n8uFFd1PblSH85d79MoTrGA37fZc6zDBu2qjghTuBt830i7Woqdw76WmDBg85WCcR3Dw0Jf2NtU5J\n+32XOs3a1D1NHRAw/pSgQ8DlzHrAffEL+S5gb3kUo5YggscCD1XbGl0X9vsuNW6W3HNk3NfMNN/4\nEc3dzHrAfQvP411F5Uh4K4hrnJIe1yG/71JHGcY9TU7zFa8xGLyz/mIdHGkfof+e/MW8IBF8dGhK\n+usa3yXt913qNKephzUxzdf/Rdqbku5NUU9/wH1LzuNdRRHsAT5bbWs8hHv8vkut4XnG89oeissf\nGY9eM34b8NJyqrr6OH+xNsTCHZImZRjPY/xGquUGsiOaVovgc8BXV5rOyOTTTfVHUvsZxvMwFFUR\nwekwELoPZXJCU/2R1B2GsbQATklLmoeXNklziCCHgvjxBrGkOhjGWnsRfMNQCH++vGb4M411StJa\n8QhFrTWnpCW1gSNjraURU9InGMSSmmIYa61EcO5QCH+8nJJ+aOwXSdKSOU2tteGUtKS2cmSslTdi\nSvrRBrGkNjGMtbIiuHAohN9TTkl/qbFOSdIITlNrJTklLalLDGOtlOEQBh6RyVca6YwkTchpaq2E\nCC4dCuLfKKekDWJJrTdRGEfEmyPiaETcUWk7OSJuiYhPRMTN4SHmakgZwtf17pch/OIZn+wAw/8v\nR+wpDxKRpKWYdGT8FuCiobZXA7dk5tcD/7O8L9VmxC7pXQtYGz4MXLUVyP2jNA/P+bySNNZEYZyZ\n7wM+O9T8XODa8uNrgectsF/SWBH82FAI/2I5Gp7/CLLiCM0rKAJ5gzrOtJa09ubZwHVqZh4tPz4K\nnLqA/kg7qmWXdOYxIjaBu4C9BrGkZVvIburMzIgYOyqJiCsrd2/NzFsX8bpaH7VeqlRMTR8E9gIH\niXBkLGmsiNgP7J/rOTInm9mLYsruhsx8ann/48D+zLw3Ik4H3puZ3zDi66Y+ZFnqieCtwAsrTVdn\n8solvmBvjfiKcoQ8eF+SjmOW3JtnZPwHwIuAXyr/fdcczyUNiCBg8LKkmgp37KMavEUgX1G231TD\n60taQxONjCPiOuB84BSK9eGfAX4fuB54InAEuCRHjBwcGWtaVs+S1GWz5N7E09SzMow1qQj+CPi2\nStMrM7m6qf5I0izqnqaWFiKCXcDD1TZHw5LWiWGsRjklLUltqU1tCcK1E8FfDwXxCw1iSeuqHWFs\nCcK1EcGjyhA+s9dWVs96e4Pdmp1/SEpagHaEsSUI10IZwl/s3S9DuOujYf+QlDS3du2mLoK4V4Lw\nyPJ6pTpF8BngH1eaviOT9zbVn4XrB/AmReUu/5CU1li3d1NbgnDlRHAC8EC1bQVGwttZy1rSnNox\nTT1YcvAI/Slrz0juqHJKeiuIV2RKerTtf0j6/62kqbRjmrrY7HJ4YERR/ELbR6YlCDtk+FIl4JxM\nbm+kM3WwlrWkIVbgUmMiOAX4m2rbyo6Eq/xDUtIQw1iNsHCHJPXNknvtWDNWJ0WQQ0H8Tw1iSZpe\ne3ZTqzMiOIvipK4thrAkzc4w1lSckpakxXOaWhMZMSV9mkEsSYthGGtHETxj1Gg4k6MNdMY60JJW\nkmGsscoQ/kDvfgsKd1gHWtJKMoy1zYgp6ZOSuLjxUakHikhaUYaxtkRw0Zgp6WO0ZVRaBG+vDvSm\nQSxpFRjGArampP+wd3/blHRbRqXWgZa0ggzjdRZxYHhK+rf4ntOSuHjk45selXqgiKQVZRivqQj+\nfZA3VtuSOOl7uP5nGDf13PyodB/V0Xh/tL6v5n5I0kJZm3oNbVsXJt5IMeI9yLipZ08nkqSJeFBE\nV9V08s+I4w13Z/JwuQZ8F7C3nP5trI+S1HWreVDEehR6WOpO5QgODQXxg+UGrYcnnnrOvGnbCDjz\nmEEsSfNrfxi35ZKaZVriTuUyhA/2X4rI5KvLT9a3IWo9/qiSpJl0Y5q6Hxo7r2t23STTxRM/1bYp\n6V2ZQ211Tj275ixpTaz2mvECg6qVFvQHRwTvB7610nQkk72L6eSc1uWPKklrbZXXjM9icF3zrJWa\n3lzEdHF5zTCVIE7ipCR+bMG9nV3T1ylLUku1P4zhTuBG4FAZVIfK+3c22akd1kCvnGFtdK7rZyPI\n6jXDmUQSJ7GItfXq++x9XH0/06z7Nn+dsiS1UhfC+CnAxcDl5VT15eX9p+z4VcvfMDRuY9lbxrSP\nD8UZdypH8LfVteFdPPyXSbxxwZvAqu/zMMXIdhM4PNVmOqtnSdJ4mbnUW/ESC3gu2EjIhI0JH78n\n4Q0Je0beX0yfes+5Mea1Nka+JhwY0bYn4cDk/10zq7eZ/ztN/z7/W3kb/d7GP8fc79mbN2/eunCb\nJfe6sYFr1o0/dWwYGrexbKcNZ3PsLB51qtKI5138+62+n8LqbqaTpDms6gau2ac3l71haNwa6PHW\nRme4rnjEGcPvGhPEi58GHnw/V5Q3130laUHaPzKe51rY5Y4Ux41uD1Gsax9/1DvB5VoR7AIerrYN\nhHD/gcu5Zrja/8Jm+W+vkIjXCktSxWpfZzz9Cy+3yMT48PsJ4PXHDcUJ/lDYcUq6LtX32fu4ULwf\n61NL0gDDePCF23uwwXH+UBhRPet1mbyi9n5KkqZmGI/uQPtCeUyfPskTvv2JfPL3qg9tZDQsSZqZ\nYTy6A52oidyKKWlJ0txWczf1vJZ4ItIijNgl/fMGsSStl3pHxk1OD7fsoIkITgPuqbZlEo1PoUuS\n5tLukXGT5xC3rCZyORLeCuIkTqoE8Wqd1SxJOq56wrjJ6eEW1UQeMSX98q0DHVo4hS5Jqkc909TF\nh8ufHh69S/kSADKvr7TVOhUcwTcydMrUUPWsDVo0hS5Jml2bp6nrmh4ePknpEuAC4OatR+wUxEs4\n6akcCW8FcVEPfFsZy+VNoS//9Kr6rNJ7kaSKesK4vunhfRTlKHvTvhcAjwIuBCZZtx53LOLUa7gj\npqSfv22XdD1T6At7Ty2wSu9FkrbMPU0dEUeAz1HUT/5yZp479Pn6dlP3R8L/Bfgw8HTgFcCpwI8w\nSX3qOetZRwyNxNnhmuFpCpK0tUZ33VbpvUhaSY0U/YiIu4Bvysy/nahTywzk4rl/Bfhm4DLg14AP\nAVdThPNka7IzruEutXDHvMVLVmldepXei6SV0+Sa8c4vWu+04peAjwDvK/8FeBmTrslOs4ZbrmEO\nT0n/HD/9HxZeuGOe4iUtu7RrLqv0XiSpJzPnugF/BdxGMQL94RGfz4Q3JGyU/+6Z9zXH3uBAwlkJ\nv5nF616XcG3CJeXn9+zYh+rn+8/Vf3yvvXz8Sdz/H4uX6d9qeI8b5YttTPj4PSPew3L7uLz3vjrv\nxZs3byt7K6J1uq9ZxDT16Zl5T0Q8HrgFeFlmvq/y+cxiFLP8acVilNQ7b/cq4FeBo8DL6Y0gj7eb\nun9cYPV84qfQ2zw05mSlJIqR2jLXMGdZL23jQRmzWqX3ImllNX5QRES8FngwM6+udur74IPvhfd/\nOzzzg/CzH8u8cWEvOtiB3gaug5VA3QRuoXqd8eTPty38gvzs0KOelMQXWPYfGx058EKS1k1E7Af2\nV5peW2sYR8RjgUdk5gMR8TiKXcT/KTNvrjwmE06qJUCWMXIqNwudwt9cdz+nvKD6qRwsYbnc3b2O\nCiWpE2ofGUfEXuCd5d3dwNsz8xd27NSk08STPH7ZyqAN8iXV5uxtztpptFpc89ye9yJJqkXj09Qj\nX2CaTrVpKjZiz/CU9BHOuuYs/vqVlfXn8X88DK0xO60sSeuh+2FcfEHjRR0i+N/AOdW2ypT0TwCv\nn2jE24L3IkmqVzfDePTo8mlMU6RjgXYs3FH09U7gcvoj3h8ALgUuHblj2wIVkrRW2nxQxE6G6w2f\nBbydopRlbUUdRtSSfkxuL9xxmCKIe/Wvn0Yx4j1aeaJ+YZNpC1R4EIIkraXmw3iwstTTgBuBi8n8\nCDUcLjEihCmuv+YLO/T1cuA6itH7C4CXM1wZq1A9BOKPgM2B97I9aD0IQZLWUPPT1P0HblBM5z69\nDOJe+9J2IE9cS3p4Kr0/jX45sEERvnuoTkdv/5rBa57HbehynVmSOq2r09TD9YYvGxg9Zh5bdBCP\nGA3vGrhcafu0cH/EWkyjXwf8NvD1FFPWmxSB3J+OzrxpIESLjw8C5+9YW7q4v0kR7JszBbHT3ZLU\nKfWG8ahAqOdM3/KlRkxJE29M4sShvgxOC/enpzeBN1PU4f494OeBn6Y4M/kW4BjFdPTo/k8StIs5\nCMHpbknqktoKZo8r6l8cyDDcNnAgw2L6kYMHOgy+1mQHWcCLtw5p6H/dDyRcMnTIxOj+H++1FnkQ\nwjTvy5s3b968LexGEwdFHE/loIj61j8r67UTrQtPcvnRqLXcwmTru5MUNFl0BTIvq5Kk2rX3OuPi\nw/oCYUT1LIAkTtoWlpNsmNq57OXgxq3xfaq31KcbwSSpEW3ewDW4/rnkDUbDQVyuC+8UxMdbr95H\nNcz6a8gXMun67vCGrt7zLDeIl74OL0maX32XNg2OJmH7SPNtwEvJvLvyxVONHMecMTx+1DrPaLVN\ndbS3961dB25I0hpp7zR1r1ODZSKHp1EPMVhmcuJwGw5hKKek552m3fkgCMZ+zsCTpLXV/jDe/skN\nBgtl9AL4Q8AzgYOVzU3bgi6CRwD/UH3KHDxjeL5Ra5tHv5KkVmrzmvF2o66n7V+H+2bghKHHbgKP\n6zeRVIK42Bm+tVP6QuCPRqzx9ka0kxks1bmBQSxJWoJmRsajR5xvA14D/Ajwa/RrP99CMUoGOLh9\nl3Q+lOy6lmWOXr1ESJI0oS6NjEftTn4NcANwiKI29aXAcyhGyY9+Ay957bZd0klkxgksc/S6mIpY\nkiSN1aaDIqpnBfdqPZ8E/Ovh/VkzF+6YlmvGkqQptX9kvNO1xMV1uHfTr918QpADQfxVfP73xgTx\nskav464vnm7tWZKkHdQXxpMcVlCG6m2c/ZogL61+Kolf/zyP/e0dnnfxBS7qLNQhSVpb7alNXYZq\nkC+pNm9Vzhp3Da8FLiRJLdLe64yLD3dcyx1RuOOVSVzDZNWwDGRJUiu0ec147FpuBE8YdbJSEh8v\n79xUefC4NWfP75UkdVY9YTxmLbcM4U/2HzZQuGPygLU4hySpwxqpTR3kjUMPe0YmHxrxxdMdA2hx\nDklSw9q7Zlx2KoInUBkJw5hrhgefYIPJzgv2/F5JUuPavGZ8vCnpcV802fXDnt8rSeqwWsK4ukHr\nU5xxcxIbQw/YvjFruoC1OIckqbNqmaaG5GTu/537OeU+4IMUIXmwcv3wJnALmdeXX3SA/glNN5eP\nvxM4r2x7iGIjl5cuSZJapbXT1EnsvZ9Tnk8Ruu8qmzfL9eDN8v7NlS85DJwP/BnF6PivgRuBj5Xt\nd+KlS5KkFVHXmnF/3bd//9EUG7MeTW+U3NOfZr6c4jjFd1AcrfiO8v7luEFLkrQidtf0Or1Tj4p1\nXzh03K8oHt87NOLZwJ+U/36YYme1QSxJWgl1Ff2obqw6RDG6/SLFaPmLFFPWgxuz+jupnw68CXhW\n+e/T8VxhSdIKqfcIxcJ5FJc4HSx3Sfemri/cekR/J/Uh4DLge4FfLf+9rGzfvrM64sDIUB93bKMk\nSS3QRBg/RHWNuPj3YNneU1yqBE8p/30icHH5b7V9+NIla1RLkjqn3nKY8z/Z8U9nshKXJKlBrb20\naYGOP/Itgre38WvTIJYktV23wniS05kmLaE5juvOkqSaNR/GRfhdMhCARfhdMjIAdxr5LqZGtevO\nkqRaNR/GRchdQO/ypn55zAsYDsAiuM9icOR7ViW0569R7dnIkqSatWMDVz+AH122fJHhqlzF486i\nKIt5MZl3b7u/2I5v4NnIkqQpdXcDVxG6VwHfV96u2griwfXap1Bc4nR5GZaXl/efUj52Meu98647\nS5I0hWZHxr1LlQqbwAnAN1EcDPH8sn3UJq0NRo1aB9eMj227P1mH538OSdLammVkXE8Yw0kjNlr9\nBPAW4KeBRwFfAv4z8HrgnwC3Aw8yPF19vOuI573OeJJrmSVJGqPNYfxGto80D1FMM38QuBT4OYqS\nlxeXt/8K/CCZbymfaA9FyczzRzzXZCNnSZKWrM1rxqN2J99dtj8D+HmKU5n+LfB3wNnAW4FnVnZY\nX7X1XDvtlna9V5LUMfWtGY9f530axbGIzwZeVLb2Do/o7bAevbt6+4uNWu99G/DSgd3WTjtLkpak\nkZFxRFwUER+PiL+MiFeNedAlFCPY/mi1uP0A8HbgpRRBPFxY4/0Uu6vfP+G676jrjF8KvKF1RTys\n9CVJKs0VxhHxCIq13YuAJwMviIgnjXjoBeW/vWnlTeBXKNaKLwbeUX5+H8WU9YXlY55JEeDfPNF0\nc+ZN20JaYdj8AAAF7ElEQVS7GBG/kPYV8bDSlyQJmH9kfC7wfzPzSGZ+Gfgt4F+OeNzB8nYV0AvV\ne4BLyby7cowiwHfQD+/emcezlLXsa+PhEVb6kiSVds/59WcAn6zc/xRw3rZH9Qt49AJx+y7nYo33\nqvLzPwi8c2C6OaK3UWv6dd7tm7raEXrF+6r+N2m+T5Kk2s07Mp5899fxdjkPfv6bt79SHptpw9Vi\nDo9YDnd+S5KYP4z/H3Bm5f6ZFKPjAXsifvH74eY98EDABsOBuNzAnP/wiGVo8x8JkqSJRcT+iLiy\nd5vpOea5tCkidgP/B/hO4NPAB4AXZObHKo/JLDZpja9qtY5Vr9bxPUvSGmikAldEPIeihOUjgGsy\n8xfm7ZQkSV3V3nKYhrEkaU20uRymJEkao54wttqUJElj1TUyttqUJElj1HlQxHznDEuS1AHt38Dl\nOcOSpBXX7g1cVpuSJGmkujZwWW1KkqQx6hoZFyUpYR8RewZKUrqrWpK05uoJ4/45w/1d1dX77qqW\nJK2xeot+LPsMX69nliR1UP0VuIrg7Z3hu7ngy5u8nlmS1Dn1h/Eyd1Uve+QtSdIS1H2dcXVX9bFt\n9xf3oht4PbMkqQHtvs64UOyq7gVvdVf1ong9sySpY1brCMW6Rt6SJI3R/nKYy1bsmj48ELxFIO8j\n86Za+iBJWmuGsSRJDevCmrEkSRpiGEuS1DDDWJKkhhnGkiQ1zDCWJKlhhrEkSQ2rN4w9QUmSpG3q\nC2NPUJIkaaTdtbxKcXDDQSxLKUnSNvVU4Co+9AQlSdLKa3MFLk9QkiRpjHrCuBgRXwFcZSBLkjSo\n3oMiPEFJkrTiPLVJkqSGtXnNWJIkjWEYS5LUMMNYkqSGGcaSJDXMMJYkqWGGsSRJDTOMJUlqmGEs\nSVLDDGNJkhpmGEuS1DDDWJKkhhnGkiQ1zDCWJKlhhrEkSQ0zjCVJaphhLElSw2YO44i4MiI+FRG3\nlbeLFtkxSZLWxTwj4wRel5nnlLd3L6pTao+I2N90HzQbv3fd5vdvvcw7TR0L6YXabH/THdDM9jfd\nAc1lf9MdUH3mDeOXRcSHI+KaiNizkB5JkrRmdgzjiLglIu4YcXsu8CZgL3A2cA9wdQ39lSRp5URm\nzv8kERvADZn51BGfm/8FJEnqkMycahl396wvFBGnZ+Y95d3vBu5YRIckSVo3M4cx8EsRcTbFruq7\ngMsW0yVJktbLQqapJUnS7JZagSsijkTER8qiIB9Y5mtpPhHx5og4GhF3VNpOLjfxfSIibnbHfHuN\n+f5ZmKcDIuLMiHhvRHw0Iu6MiB8v2/3564Advn9T/fwtdWQcEXcB35SZf7u0F9FCRMSzgQeB3+xt\nxIuIQ8BnMvNQRLwKOCkzX91kPzXamO/fa4EHMvN1jXZOO4qI04DTMvP2iDgB+HPgecCL8eev9Xb4\n/l3CFD9/ddSmdgNXB2Tm+4DPDjU/F7i2/Phaiv/B1EJjvn/gz1/rZea9mXl7+fGDwMeAM/DnrxN2\n+P7BFD9/yw7jBN4TER+KiB9e8mtp8U7NzKPlx0eBU5vsjGZiYZ4OKS8TPQf4M/z565zK9+9Py6aJ\nf/6WHcb7MvMc4DnAS8upNHVQFusZ7vbrFgvzdEg5xfm7wMsz84Hq5/z5a7/y+/c7FN+/B5ny52+p\nYdy7Djkz/wZ4J3DuMl9PC3e0XA8hIk4H7mu4P5pCZt6XJeDX8eevtSLikRRB/NbMfFfZ7M9fR1S+\nf2/rff+m/flbWhhHxGMj4qvLjx8HXMiYwiBqrT8AXlR+/CLgXTs8Vi1T/gLvGVuYR82KiACuAf4i\nM19f+ZQ/fx0w7vs37c/f0nZTR8ReitEwFMVF3p6Zv7CUF9PcIuI64HzgFIr1qZ8Bfh+4HngicAS4\nJDOPNdVHjTfi+/dailN/BgrzVNYg1RIR8Szgj4GP0J+Kfg3wAfz5a70x37+fAl7AFD9/Fv2QJKlh\ndVzaJEmSdmAYS5LUMMNYkqSGGcaSJDXMMJYkqWGGsSRJDTOMJUlqmGEsSVLD/j8xXPUYvalKagAA\nAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(8, 6))\n", "\n", "plt.plot(data1['x'], data1['y'], 'rx')\n", "plt.plot(data1['x'], clf.predict(data1[['x']]))\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "collapsed": false }, "source": [ "#Linear regression with multiple variables" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": true }, "outputs": [], "source": [ "data2 = pd.read_csv('ex1data2.txt', header=None, names=['size', 'number of bedrooms', 'price'], dtype=float)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "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", "
sizenumber of bedroomsprice
021043399900
116003329900
224003369000
314162232000
430004539900
\n", "
" ], "text/plain": [ " size number of bedrooms price\n", "0 2104 3 399900\n", "1 1600 3 329900\n", "2 2400 3 369000\n", "3 1416 2 232000\n", "4 3000 4 539900" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data2.head()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from sklearn import preprocessing" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 2000.68085106 3.17021277] [ 7.86202619e+02 7.52842809e-01]\n" ] } ], "source": [ "scaler = preprocessing.StandardScaler().fit(data2[['size', 'number of bedrooms']].values)\n", "features = scaler.transform(data2[['size', 'number of bedrooms']].values)\n", "print scaler.mean_, scaler.std_ " ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 109447.79646964 -6578.35485416]\n", "340412.659574\n" ] } ], "source": [ "clf2 = linear_model.LinearRegression()\n", "clf2.fit(features, data2['price'].values)\n", "\n", "print clf2.coef_\n", "print clf2.intercept_" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "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.9" } }, "nbformat": 4, "nbformat_minor": 0 }