part 08.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$$m = \\frac{\\overline{X} \\cdot \\overline{Y} - \\overline{XY}}{(\\overline{X})^2 - \\overline{X^2}}$$"
      ],
      "text/plain": [
       "<IPython.core.display.Math object>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import Math\n",
    "Math(r'm = \\frac{\\overline{X} \\cdot \\overline{Y} - \\overline{XY}}{(\\overline{X})^2 - \\overline{X^2}}')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from statistics import mean\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEvhJREFUeJzt3X+M3PV95/Hni7OplqTmj+vKPYWyTnKxfPzhEjjxQyTn\n8SVEMajmH06kuhM6hBoHOBGdojRRlRN7UtXeP9f8UIKILz6unLg9H1FIqYiqpFfPobQVdTAuTQCF\nBoYDmqKpCKQtqxNp3vfHDmjZrD0z6/HOzofnQ1r5O9/v+zvft8fj1378nhlvqgpJ0uw7Z9oNSJIm\nw0CXpEYY6JLUCANdkhphoEtSIwx0SWrE0EBPsjvJo0lODH59Jcnta2r2JXl5UHMiyWfOXsuSpPVs\nG1ZQVd8H3guQ5BzgeeD+dUofqqqDk21PkjSqcUcuHwR+UFXPrXMsE+hHkrRB4wb6DcDSKY5dmeRk\nkgeTXHSGfUmSxpRRP/qfZDvwV8BFVdVfc+ztwE+r6tUkB4DPV9XuiXcrSTqlcQL9IHBrVX14hNpn\ngEur6qU1+/2PYyRpA6pq6Fh7nJHLr3KKcUuSnau2L2PlG8VL69VW1cx+3XHHHVPvwf6n38dbsf9Z\n7r2F/kc19F0ug5A+j5UXRD+6at+hlXyuw8D1SW4BXgOWWZm1S5I20UiBXlWvAvNr9n151faXgC9N\ntjVJ0jj8pOgYOp3OtFs4I/Y/XbPc/yz3DrPf/6hGflF0IhdLajOvJ0ktSEJN+EVRSdIWZqBLUiMM\ndElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCX\npEYY6JLUCANdkhoxNNCT7E7yaJITg19fSXL7OnVfSPJUkpNJLj477UqSTmXoD4muqu8D7wVIcg7w\nPHD/6pokB4B3V9V7klwO3AVcMfl2JUmnMu7I5YPAD6rquTX7rwPuAaiqh4Hzk+ycQH+SdEb6/T7H\njx+n3+9Pu5WzbtxAvwFYWmf/O4DVIf/CYJ8kTc3S0lEWFvZw9dUfY2FhD0tLR6fd0lk1cqAn2Q4c\nBO47e+1I0mT0+31uvvlWlpeP8corj7C8fIybb7616ZX60Bn6KgeAR6pqvUfjBeCXVt2+YLDvZywu\nLr6x3el06HQ6Y7QgSaPp9Xqce+4ulpf3DvbsZfv2BXq9HvPz81PtbZhut0u32x37vFTVaIXJEvAH\nVfW76xy7Britqq5NcgXwuar6mRdFk9So15OkM9Hv91lY2MPy8jFgL/AYc3P7efbZJ7d8oK+VhKrK\nsLqRRi5JzmPlBdGvrdp3KMlHAarqG8AzSf4S+DJw64a6lqQJmZ+f58iRO5mb28+OHZcwN7efI0fu\nnLkwH8fIK/SJXMwVuqRN1u/36fV67Nq1a2bDfNQVuoEuSVvcREcukqStz0CXpEYY6JLUCANdkhph\noEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6\nJDXCQJekRoz6Q6LPT3JfkieSfC/J5WuO70vycpITg6/PnJ12JUmnsm3Eus8D36iqf5VkG3DeOjUP\nVdXBybUmSRrH0EBPsgN4f1X9W4Cq+gnw4/VKJ9uaJGkco4xc3gn8TZK7B+OUw0nm1qm7MsnJJA8m\nuWjCfUqShhhl5LINuAS4raq+k+RzwKeBO1bVPAJcWFWvJjkAfB3Yvd6dLS4uvrHd6XTodDob61yS\nGtXtdul2u2Ofl6o6fUGyE/jTqnrX4Pb7gE9V1a+c5pxngEur6qU1+2vY9SRJb5aEqho61h46cqmq\nF4Hnkry+4v4A8Piai+1ctX0ZK98o3hTmkqSza9R3udwO3JtkO/A0cFOSQ0BV1WHg+iS3AK8By8AN\nZ6VbSdIpDR25TPRijlwkaWwTG7lIkmaDgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMM\ndElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREjBXqS85Pc\nl+SJJN9Lcvk6NV9I8lSSk0kunnyrkqTTGXWF/nngG1X1z4BfBp5YfTDJAeDdVfUe4BBw10S71Fte\nv9/n+PHj9Pv9abfyluNjPzuGBnqSHcD7q+pugKr6SVX9eE3ZdcA9g+MPA+cn2TnpZvXWtLR0lIWF\nPVx99cdYWNjD0tLRabf0luFjP1tSVacvSH4ZOAw8zsrq/DvAx6tqeVXN7wO/XVV/Mrj9h8CvV9WJ\nNfdVw64nrdbv91lY2MPy8jFgL/AYc3P7efbZJ5mfn592e03zsd86klBVGVa3bYT72gZcAtxWVd9J\n8jng08AdG2lscXHxje1Op0On09nI3egtotfrce65u1he3jvYs5ft2xfo9XqGylnmYz893W6Xbrc7\n9nmjrNB3An9aVe8a3H4f8Kmq+pVVNXcBx6rq6OD2k8C+qnpxzX25QtdYXCVOj4/91jHqCn3oDH0Q\nys8l2T3Y9QFWxi+rPQDcOLjwFcDLa8Nc2oj5+XmOHLmTubn97NhxCXNz+zly5E4DZRP42M+eoSt0\neGOO/hVgO/A0cBPwEaCq6vCg5ovAh4G/B25aOz8f1LhC14b0+316vR67du0yUDaZj/30jbpCHynQ\nJ8VAl6TxTWzkIkmaDQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCX\npEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGbBulKEkPeAX4KfBaVV225vg+\n4PdY+QHSAF+rqt+cYJ+SpCFGCnRWgrxTVT86Tc1DVXVwAj1JkjZg1JFLRqgd+hOpJUlnz6iBXsC3\nkhxP8munqLkyyckkDya5aEL9SZJGNOrI5aqq+mGSeVaC/Ymq+vaq448AF1bVq0kOAF8Hdq93R4uL\ni29sdzodOp3OhhqXpFZ1u1263e7Y56WqxjshuQP426r6ndPUPANcWlUvrdlf415Pkt7qklBVQ8fa\nQ0cuSc5L8vbB9tuADwHfXVOzc9X2Zax8o3hTmEuSzq5RRi47gfuT1KD+3qr6ZpJDQFXVYeD6JLcA\nrwHLwA1nrWNJ0rrGHrmc0cUcuUjS2CY2cpEkzQYDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXC\nQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSI0YK\n9CS9JH+e5NEkf3aKmi8keSrJySQXT7ZNSdIw20as+ynQqaofrXcwyQHg3VX1niSXA3cBV0yoxy2h\n3+/T6/XYtWsX8/Pz025HM8bnjzbDqCOXDKm9DrgHoKoeBs5PsvMMe9sylpaOsrCwh6uv/hgLC3tY\nWjo67ZY0Q3z+aLOkqoYXJU8DLwP/AByuqv+y5vjvA79dVX8yuP2HwK9X1Yk1dTXK9baSfr/PwsIe\nlpePAXuBx5ib28+zzz7pSktD+fzRJCShqjKsbtSRy1VV9cMk88C3kjxRVd/eSGOLi4tvbHc6HTqd\nzkbuZtP0ej3OPXcXy8t7B3v2sn37Ar1ez7+QGsrnjzai2+3S7XbHPm+kFfqbTkjuAP62qn5n1b67\ngGNVdXRw+0lgX1W9uOZcV+h6S/H5o0kYdYU+dIae5Lwkbx9svw34EPDdNWUPADcOaq4AXl4b5rNq\nfn6eI0fuZG5uPzt2XMLc3H6OHLnTv4waic8fbaahK/Qk7wTuB4qVEc29VfWfkhwCqqoOD+q+CHwY\n+HvgprXz80HNzK3QX+e7FHQmfP7oTIy6Qh975HImZjnQJWlaJjZykSTNBgNdkhphoEtSIwx0SWqE\ngS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjo\nktQIA12SGmGgS1IjRg70JOckOZHkgXWO7Uvy8uD4iSSfmWybkqRhto1R+3HgcWDHKY4/VFUHz7wl\nSdJGjLRCT3IBcA3wldOVTaQjSdKGjDpy+SzwSaBOU3NlkpNJHkxy0Zm3Jkkax9CRS5JrgRer6mSS\nDuuvxB8BLqyqV5McAL4O7F7v/hYXF9/Y7nQ6dDqd8buWpIZ1u1263e7Y56XqdItuSPJbwL8BfgLM\nAT8PfK2qbjzNOc8Al1bVS2v217DrSZLeLAlVNXSsPTTQ19zpPuATa1/8TLKzql4cbF8G/K+q2rXO\n+Qa6JI1p1EAf510uay9wCKiqOgxcn+QW4DVgGbhho/crSdqYsVboZ3wxV+iSNLZRV+h+UlSSGmGg\nS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrok\nNcJAl6RGGOiS1AgDXZIaYaBLUiNGDvQk5yQ5keSBUxz/QpKnkpxMcvHkWpQkjWKcFfrHgcfXO5Dk\nAPDuqnoPcAi4awK9SZLGMFKgJ7kAuAb4yilKrgPuAaiqh4Hzk+ycSIeaiH6/z/Hjx+n3+9NuRdJZ\nMuoK/bPAJ4E6xfF3AM+tuv3CYJ+2gKWloyws7OHqqz/GwsIelpaOTrslSWfB0EBPci3wYlWdBDL4\n0ozo9/vcfPOtLC8f45VXHmF5+Rg333yrK3WpQdtGqLkKOJjkGmAO+Pkk91TVjatqXgB+adXtCwb7\nfsbi4uIb251Oh06nM2bLGkev1+Pcc3exvLx3sGcv27cv0Ov1mJ+fn2pvktbX7Xbpdrtjn5eqU01R\n1ilO9gGfqKqDa/ZfA9xWVdcmuQL4XFVdsc75Nc71dOb6/T4LC3tYXj4G7AUeY25uP88++6SBLs2I\nJFTV0OnIht+HnuRQko8CVNU3gGeS/CXwZeDWjd6vJmt+fp4jR+5kbm4/O3Zcwtzcfo4cudMwlxo0\n1gr9jC/mCn1q+v0+vV6PXbt2GebSjBl1hW6gS9IWd9ZHLpKkrcVAl6RGGOiS1AgDXZIaYaBLUiMM\ndElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCX\npEYMDfQkP5fk4SSPJvlekt9ap2ZfkpeTnBh8febstCtJOpWhgV5V/w/YX1XvBfYC/zLJVeuUPlRV\nlwy+fnPSjW4F3W532i2cEfufrlnuf5Z7h9nvf1QjjVyq6tXB5s8NzvnROmVDf4DprJv1J4X9T9cs\n9z/LvcPs9z+qkQI9yTlJHgX+GuhW1ePrlF2Z5GSSB5NcNNEuJUlDjbpC/+lg5HIB8C+S7FtT8ghw\nYVVdDHwR+Ppk25QkDZOqGu+E5D8Ar1bVfz5NzTPApVX10pr9411MkgRAVQ0da28bVpDkF4DXquqV\nJHPA1cB/XFOzs6peHGxfxso3ipfW3tcoDUmSNmZooAP/BPjdJGFlRPPfq+p/JzkEVFUdBq5Pcgvw\nGrAM3HDWOpYkrWvskYskaWvatE+KJvlwkieTfD/JpzbrupOQ5EiSF5M8Nu1expXkgiR/NPhQ2F8k\nuX3aPY1jlA+2zYLBO8VOJHlg2r2MK0kvyZ8P/gz+bNr9jCvJ+UnuS/LE4Dl0+bR7GlWS3YPH/cTg\n11dO93d4U1boSc4Bvg98APgr4Djwkap68qxffAKSvA/4O+Ceqto77X7GkeQXgV+sqpNJ3s7KO5Ku\nm5XHHiDJeVX1apJ/BPwx8Imq+uNp9zWOJP8euBTYUVUHp93POJI8zcqbHNb7/MmWl+S/Af+nqu5O\nsg04r6p+POW2xjbI0eeBy6vqufVqNmuFfhnwVFU9W1WvAf8TuG6Trn3GqurbrP9hqi2vqv66qk4O\ntv8OeAJ4x3S7Gs+IH2zbspJcAFwDfGXavWzQ66+fzZwkO4D3V9XdAFX1k1kM84EPAj84VZjD5v0h\nvQNY3cTzzFiotCDJLuBi4OHpdjKeET/YtpV9FvgkMKsvWBXwrSTHk/zatJsZ0zuBv0ly92BscXjw\nbr1ZdAOwdLqCmfyuq/ENxi1fBT4+WKnPjBE+2LZlJbkWeHHwr6Qwm/9FxlVVdQkr/8q4bTCCnBXb\ngEuALw1+D68Cn55uS+NLsh04CNx3urrNCvQXgAtX3b5gsE+bYDA3/Corbzn9vWn3s1GDfyo/CPzz\nafcyhquAg4M59BKwP8k9U+5pLFX1w8GvfeB+Vkaos+J54Lmq+s7g9ldZCfhZcwB4ZPBncEqbFejH\ngX+aZCHJucBHgFl7tX9WV1cA/xV4vKo+P+1GxpXkF5KcP9h+/YNtJ6fb1eiq6jeq6sKqehcrz/s/\nqqobp93XqJKcN/jXHUneBnwI+O50uxrd4AOPzyXZPdj1AWDWRnYAv8qQcQuM9sGiM1ZV/5Dk3wHf\nZOWbyJGqemIzrj0JSf4H0AH+cZL/C9zx+ossW93gvzr+18BfDObQBfxGVf3BdDsb2bofbJtyT28l\nO4H7B/9txzbg3qr65pR7GtftwL2DscXTwE1T7mcsSc5j5QXRjw6t9YNFktQGXxSVpEYY6JLUCANd\nkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNeL/A9oOHBRRi679AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f7b7fb64278>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "xs = np.array([1,2,3,4,5,6], dtype=np.float64)\n",
    "ys = np.array([5,4,6,5,6,7], dtype=np.float64)\n",
    "\n",
    "plt.scatter(xs, ys)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.42857142857142866"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def best_fit_slope(xs, ys):\n",
    "    m = (\n",
    "         ((mean(xs) * mean(ys)) - mean(xs*ys))\n",
    "         / \n",
    "         ((mean(xs))**2 - mean(xs**2))\n",
    "        )\n",
    "    return m\n",
    "\n",
    "m = best_fit_slope(xs, ys)\n",
    "m"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}