aula-16-minimos-quadrados.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ajuste de Curvas por Mínimos Quadrados\n",
    "\n",
    "## Regressão Linear\n",
    "\n",
    "O exemplo mais simples de aproximação por mínimos quadrados é ajustar uma reta a um conjunto de pares de observações $(x_1, y_1), (x_2, y_2), \\dots, (x_n, y_n)$. A expressão matemática do ajuste por uma reta é\n",
    "\n",
    "$$ y = a_0 + a_1 x + e \\qquad (1) $$\n",
    "\n",
    "onde $a_0$ e $a_1$ são coeficientes representando a intersecção com o eixo $y$ e a inclinação, respectivamente, e $e$ é o erro ou resíduo entre o modelo e a observação, o qual pode ser representado, depois de se reorganizar a equação acima, por\n",
    "\n",
    "$$ e = y − a_0 − a_1 x $$\n",
    "\n",
    "Portanto, o erro ou resíduo é a discrepância entre o valor verdadeiro de $y$ e o valor aproximado, $a_0 + a_1 x$, previsto pela equação linear.\n",
    "\n",
    "### Critério para um Melhor Ajuste\n",
    "\n",
    "Uma estratégia para ajustar uma melhor reta pelos dados é minimizar a soma dos quadrados dos resíduos entre o $y$ medido e o $y$ calculado com o modelo linear\n",
    "\n",
    "$$ S_r = \\sum _{i = 1}^{n} e_i^2 = \\sum _{i = 1}^{n} (y_{i,medido} - y_{i,modelo})^2 = \\sum _{i = 1}^{n} (y_i - a_0 - a_1 x_i)^2 \\qquad (2) $$\n",
    "\n",
    "Esse critério tem diversas vantagens, incluindo o fato de que ele fornece uma única reta para um dado conjunto de dados.\n",
    "\n",
    "### Ajuste por Mínimos Quadrados por uma Reta\n",
    "\n",
    "Para determinar os valores de $a_0$ e $a_1$, a Equação (2) é derivada com relação a cada coeficiente\n",
    "\n",
    "$$\n",
    "\\dfrac{\\partial S_r}{\\partial a_0} = - 2 \\sum (y_i - a_0 - a_1 x_i) \\\\\n",
    "\\dfrac{\\partial S_r}{\\partial a_1} = - 2 \\sum [(y_i - a_0 - a_1 x_i) x_i]\n",
    "$$\n",
    "\n",
    "Igualando essas derivadas a zero será obtido um $S_r$ mínimo\n",
    "\n",
    "$$\n",
    "0 = \\sum y_i - \\sum a_0 - \\sum a_1 x_i \\\\\n",
    "0 = \\sum y_i x_i - \\sum a_0 x_i - \\sum a_1 x_i^2\n",
    "$$\n",
    "\n",
    "Agora, pode-se expressar essas equações como um conjunto de duas equações lineares simultâneas em duas variáveis ($a_0$ e $a_1$)\n",
    "\n",
    "$$\n",
    "n a_0 + \\left( \\sum x_i \\right) a_1 = \\sum y_i \\\\\n",
    "\\left( \\sum x_i \\right) a_0 + \\left( \\sum x_i^2 \\right) a_1 = \\sum x_i y_i\n",
    "$$\n",
    "\n",
    "Essas são as chamadas equações normais. Elas podem ser resolvidas simultaneamente para obter-se\n",
    "\n",
    "$$\n",
    "a_1 = \\dfrac{n \\sum x_i y_i - \\sum x_i \\sum y_i}{n \\sum x_i^2 - (\\sum x_i)^2} \\qquad (3) \\\\\n",
    "a_0 = \\overline{y} - a_1 \\overline{x} \\qquad (4)\n",
    "$$\n",
    "\n",
    "onde $\\overline{y}$ e $\\overline{x}$ são as médias de $y$ e $x$, respectivamente.\n",
    "\n",
    "### Quantificação do Erro da Regressão Linear\n",
    "\n",
    "Um “desvio padrão” para a reta de regressão pode ser determinado por\n",
    "\n",
    "$$ s_{y/x} = \\sqrt{\\dfrac{S_r}{n-2}} $$\n",
    "\n",
    "onde $s_{y/x}$ é chamado erro padrão da estimativa. O subscrito “$y/x$” indica que o erro é para um valor previsto de $y$ correspondente a um valor particular de $x$. Além disso, observe que agora estamos dividindo por $n − 2$ porque duas estimativas provenientes dos dados —$a_0$ e $a_1$ — foram usadas para calcular $S_r$; portanto, perdemos dois graus de liberdade.  \n",
    "\n",
    "Exatamente como no caso do desvio padrão, o erro padrão da estimativa quantifica a dispersão dos dados. Entretanto, erro padrão da estimativa quantifica a dispersão em torno da reta de regressão, em contraste com o desvio padrão original que quantificava a dispersão em torno da média.\n",
    "\n",
    "A fim de quantificar o \"quão bom\" é um ajuste de reta, dois conceitos utilizados são os coeficiente de determinação ($r^2$) e o coeficiente de correlação ($r = \\sqrt{r^2}$).\n",
    "\n",
    "$$ r^2 = \\dfrac{S_t - S_r}{S_t} $$\n",
    "\n",
    "$$ S_t = \\sum _{i = 1}^n (y_i − \\overline{y})^2 $$\n",
    "\n",
    "Para um ajuste perfeito, $S_r = 0$ e $r = r_2 = 1$, significando que a reta explica 100\\% da variação dos dados. Para $r = r_2 = 0$, $S_r = St$ e o ajuste não representa nenhuma melhora. Uma formulação alternativa para $r$ que é mais conveniente para implementação computacional é\n",
    "\n",
    "$$ r = \\dfrac{n \\sum x_i y_i - (\\sum x_i)(\\sum y_i)}{\\sqrt{n \\sum x_i^2 - (\\sum x_i)^2} \\sqrt{n \\sum y_i^2 - (\\sum y_i)^2}} \\qquad (5) $$\n",
    "\n",
    "## Regressão Polinomial\n",
    "\n",
    "O procedimento dos mínimos quadrados pode ser prontamente estendido para ajustar dados por um polinômio de grau mais alto. Por exemplo, suponha que se queira ajustar um polinômio de segundo grau ou quadrático\n",
    "\n",
    "$$ y = a_0 + a_1 x + a_2 x_2 + e \\qquad (6) $$\n",
    "\n",
    "Nesse caso, a soma dos quadrados dos resíduos é \n",
    "\n",
    "$$ S_r = \\sum _{i = 1}^n (y_i − a_0 − a_1 x_i − a_2 x_i^2)^2 \\qquad (7) $$\n",
    "\n",
    "Seguindo um procedimento análogo ao anterior, toma-se a derivada da Equação (7) com relação a cada um dos coeficientes desconhecidos do polinômio, como em\n",
    "\n",
    "$$\n",
    "\\dfrac{\\partial S_r}{\\partial a_0} = - 2 \\sum (y_i - a_0 - a_1 x_i - a_2 x_i^2) \\\\\n",
    "\\dfrac{\\partial S_r}{\\partial a_1} = - 2 \\sum [(y_i - a_0 - a_1 x_i - a_2 x_i^2) x_i] \\\\\n",
    "\\dfrac{\\partial S_r}{\\partial a_2} = - 2 \\sum [(y_i - a_0 - a_1 x_i - a_2 x_i^2) x_i^2]\n",
    "$$\n",
    "\n",
    "Essas equações podem ser igualadas a zero e reorganizadas para determinar o seguinte conjunto de equações normais\n",
    "\n",
    "$$\n",
    "n a_0 + \\left( \\sum x_i \\right) a_1 + \\left( \\sum x_i^2 \\right) a_2 = \\sum y_i \\\\\n",
    "\\left( \\sum x_i \\right) a_0 + \\left( \\sum x_i^2 \\right) a_1 + \\left( \\sum x_i^3 \\right) a_2 = \\sum x_i y_i \\\\\n",
    "\\left( \\sum x_i^2 \\right) a_0 + \\left( \\sum x_i^3 \\right) a_1 + \\left( \\sum x_i^4 \\right) a_2 = \\sum x_i^2 y_i\n",
    "$$\n",
    "\n",
    "Nesse caso, vê-se que o problema de determinar o polinômio de segundo grau por mínimos quadrados é equivalente a resolver um sistema de três equações lineares simultâneas. Na forma matricial, temos\n",
    "\n",
    "$$\n",
    "\\begin{bmatrix}\n",
    "n & \\sum x_i & \\sum x_i^2 \\\\\n",
    "\\sum x_i & \\sum x_i^2 & \\sum x_i^3 \\\\\n",
    "\\sum x_i^2 & \\sum x_i^3 & \\sum x_i^4\n",
    "\\end{bmatrix}\n",
    "\\begin{bmatrix}\n",
    "a_0 \\\\\n",
    "a_1 \\\\\n",
    "a_2\n",
    "\\end{bmatrix} =\n",
    "\\begin{bmatrix}\n",
    "\\sum y_i \\\\\n",
    "\\sum x_i y_i \\\\\n",
    "\\sum x_i^2 y_i\n",
    "\\end{bmatrix} \\qquad (8)\n",
    "$$\n",
    "\n",
    "O caso bidimensional pode ser facilmente estendido para um polinômio de grau $m$.\n",
    "\n",
    "$$ y = a_0 + a_1 x + a_2 x^2 + \\dots + a_m x^m + e $$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Nota\n",
    "\n",
    "O presente texto (conteúdo teórico) é um resumo baseado no livro Métodos Numéricos para Engenharia (Chapra e Canale)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Motivação: Evolução da População Paraibana\n",
    "\n",
    "Este exemplo lê um arquivo CSV, plota o gráfico de dispersão, compara modelos e formata os dados para visualização."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "População obtida por interpolação linear em: 1975 = 1225580; 1982 = 1552017; \n"
     ]
    },
    {
     "data": {
      "text/html": [
       "\n",
       "    <div class=\"bk-root\">\n",
       "        <a href=\"https://bokeh.org\" target=\"_blank\" class=\"bk-logo bk-logo-small bk-logo-notebook\"></a>\n",
       "        <span id=\"1056\">Loading BokehJS ...</span>\n",
       "    </div>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": [
       "\n",
       "(function(root) {\n",
       "  function now() {\n",
       "    return new Date();\n",
       "  }\n",
       "\n",
       "  var force = true;\n",
       "\n",
       "  if (typeof root._bokeh_onload_callbacks === \"undefined\" || force === true) {\n",
       "    root._bokeh_onload_callbacks = [];\n",
       "    root._bokeh_is_loading = undefined;\n",
       "  }\n",
       "\n",
       "  var JS_MIME_TYPE = 'application/javascript';\n",
       "  var HTML_MIME_TYPE = 'text/html';\n",
       "  var EXEC_MIME_TYPE = 'application/vnd.bokehjs_exec.v0+json';\n",
       "  var CLASS_NAME = 'output_bokeh rendered_html';\n",
       "\n",
       "  /**\n",
       "   * Render data to the DOM node\n",
       "   */\n",
       "  function render(props, node) {\n",
       "    var script = document.createElement(\"script\");\n",
       "    node.appendChild(script);\n",
       "  }\n",
       "\n",
       "  /**\n",
       "   * Handle when an output is cleared or removed\n",
       "   */\n",
       "  function handleClearOutput(event, handle) {\n",
       "    var cell = handle.cell;\n",
       "\n",
       "    var id = cell.output_area._bokeh_element_id;\n",
       "    var server_id = cell.output_area._bokeh_server_id;\n",
       "    // Clean up Bokeh references\n",
       "    if (id != null && id in Bokeh.index) {\n",
       "      Bokeh.index[id].model.document.clear();\n",
       "      delete Bokeh.index[id];\n",
       "    }\n",
       "\n",
       "    if (server_id !== undefined) {\n",
       "      // Clean up Bokeh references\n",
       "      var cmd = \"from bokeh.io.state import curstate; print(curstate().uuid_to_server['\" + server_id + \"'].get_sessions()[0].document.roots[0]._id)\";\n",
       "      cell.notebook.kernel.execute(cmd, {\n",
       "        iopub: {\n",
       "          output: function(msg) {\n",
       "            var id = msg.content.text.trim();\n",
       "            if (id in Bokeh.index) {\n",
       "              Bokeh.index[id].model.document.clear();\n",
       "              delete Bokeh.index[id];\n",
       "            }\n",
       "          }\n",
       "        }\n",
       "      });\n",
       "      // Destroy server and session\n",
       "      var cmd = \"import bokeh.io.notebook as ion; ion.destroy_server('\" + server_id + \"')\";\n",
       "      cell.notebook.kernel.execute(cmd);\n",
       "    }\n",
       "  }\n",
       "\n",
       "  /**\n",
       "   * Handle when a new output is added\n",
       "   */\n",
       "  function handleAddOutput(event, handle) {\n",
       "    var output_area = handle.output_area;\n",
       "    var output = handle.output;\n",
       "\n",
       "    // limit handleAddOutput to display_data with EXEC_MIME_TYPE content only\n",
       "    if ((output.output_type != \"display_data\") || (!output.data.hasOwnProperty(EXEC_MIME_TYPE))) {\n",
       "      return\n",
       "    }\n",
       "\n",
       "    var toinsert = output_area.element.find(\".\" + CLASS_NAME.split(' ')[0]);\n",
       "\n",
       "    if (output.metadata[EXEC_MIME_TYPE][\"id\"] !== undefined) {\n",
       "      toinsert[toinsert.length - 1].firstChild.textContent = output.data[JS_MIME_TYPE];\n",
       "      // store reference to embed id on output_area\n",
       "      output_area._bokeh_element_id = output.metadata[EXEC_MIME_TYPE][\"id\"];\n",
       "    }\n",
       "    if (output.metadata[EXEC_MIME_TYPE][\"server_id\"] !== undefined) {\n",
       "      var bk_div = document.createElement(\"div\");\n",
       "      bk_div.innerHTML = output.data[HTML_MIME_TYPE];\n",
       "      var script_attrs = bk_div.children[0].attributes;\n",
       "      for (var i = 0; i < script_attrs.length; i++) {\n",
       "        toinsert[toinsert.length - 1].firstChild.setAttribute(script_attrs[i].name, script_attrs[i].value);\n",
       "        toinsert[toinsert.length - 1].firstChild.textContent = bk_div.children[0].textContent\n",
       "      }\n",
       "      // store reference to server id on output_area\n",
       "      output_area._bokeh_server_id = output.metadata[EXEC_MIME_TYPE][\"server_id\"];\n",
       "    }\n",
       "  }\n",
       "\n",
       "  function register_renderer(events, OutputArea) {\n",
       "\n",
       "    function append_mime(data, metadata, element) {\n",
       "      // create a DOM node to render to\n",
       "      var toinsert = this.create_output_subarea(\n",
       "        metadata,\n",
       "        CLASS_NAME,\n",
       "        EXEC_MIME_TYPE\n",
       "      );\n",
       "      this.keyboard_manager.register_events(toinsert);\n",
       "      // Render to node\n",
       "      var props = {data: data, metadata: metadata[EXEC_MIME_TYPE]};\n",
       "      render(props, toinsert[toinsert.length - 1]);\n",
       "      element.append(toinsert);\n",
       "      return toinsert\n",
       "    }\n",
       "\n",
       "    /* Handle when an output is cleared or removed */\n",
       "    events.on('clear_output.CodeCell', handleClearOutput);\n",
       "    events.on('delete.Cell', handleClearOutput);\n",
       "\n",
       "    /* Handle when a new output is added */\n",
       "    events.on('output_added.OutputArea', handleAddOutput);\n",
       "\n",
       "    /**\n",
       "     * Register the mime type and append_mime function with output_area\n",
       "     */\n",
       "    OutputArea.prototype.register_mime_type(EXEC_MIME_TYPE, append_mime, {\n",
       "      /* Is output safe? */\n",
       "      safe: true,\n",
       "      /* Index of renderer in `output_area.display_order` */\n",
       "      index: 0\n",
       "    });\n",
       "  }\n",
       "\n",
       "  // register the mime type if in Jupyter Notebook environment and previously unregistered\n",
       "  if (root.Jupyter !== undefined) {\n",
       "    var events = require('base/js/events');\n",
       "    var OutputArea = require('notebook/js/outputarea').OutputArea;\n",
       "\n",
       "    if (OutputArea.prototype.mime_types().indexOf(EXEC_MIME_TYPE) == -1) {\n",
       "      register_renderer(events, OutputArea);\n",
       "    }\n",
       "  }\n",
       "\n",
       "  \n",
       "  if (typeof (root._bokeh_timeout) === \"undefined\" || force === true) {\n",
       "    root._bokeh_timeout = Date.now() + 5000;\n",
       "    root._bokeh_failed_load = false;\n",
       "  }\n",
       "\n",
       "  var NB_LOAD_WARNING = {'data': {'text/html':\n",
       "     \"<div style='background-color: #fdd'>\\n\"+\n",
       "     \"<p>\\n\"+\n",
       "     \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n",
       "     \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n",
       "     \"</p>\\n\"+\n",
       "     \"<ul>\\n\"+\n",
       "     \"<li>re-rerun `output_notebook()` to attempt to load from CDN again, or</li>\\n\"+\n",
       "     \"<li>use INLINE resources instead, as so:</li>\\n\"+\n",
       "     \"</ul>\\n\"+\n",
       "     \"<code>\\n\"+\n",
       "     \"from bokeh.resources import INLINE\\n\"+\n",
       "     \"output_notebook(resources=INLINE)\\n\"+\n",
       "     \"</code>\\n\"+\n",
       "     \"</div>\"}};\n",
       "\n",
       "  function display_loaded() {\n",
       "    var el = document.getElementById(\"1056\");\n",
       "    if (el != null) {\n",
       "      el.textContent = \"BokehJS is loading...\";\n",
       "    }\n",
       "    if (root.Bokeh !== undefined) {\n",
       "      if (el != null) {\n",
       "        el.textContent = \"BokehJS \" + root.Bokeh.version + \" successfully loaded.\";\n",
       "      }\n",
       "    } else if (Date.now() < root._bokeh_timeout) {\n",
       "      setTimeout(display_loaded, 100)\n",
       "    }\n",
       "  }\n",
       "\n",
       "\n",
       "  function run_callbacks() {\n",
       "    try {\n",
       "      root._bokeh_onload_callbacks.forEach(function(callback) {\n",
       "        if (callback != null)\n",
       "          callback();\n",
       "      });\n",
       "    } finally {\n",
       "      delete root._bokeh_onload_callbacks\n",
       "    }\n",
       "    console.debug(\"Bokeh: all callbacks have finished\");\n",
       "  }\n",
       "\n",
       "  function load_libs(css_urls, js_urls, callback) {\n",
       "    if (css_urls == null) css_urls = [];\n",
       "    if (js_urls == null) js_urls = [];\n",
       "\n",
       "    root._bokeh_onload_callbacks.push(callback);\n",
       "    if (root._bokeh_is_loading > 0) {\n",
       "      console.debug(\"Bokeh: BokehJS is being loaded, scheduling callback at\", now());\n",
       "      return null;\n",
       "    }\n",
       "    if (js_urls == null || js_urls.length === 0) {\n",
       "      run_callbacks();\n",
       "      return null;\n",
       "    }\n",
       "    console.debug(\"Bokeh: BokehJS not loaded, scheduling load and callback at\", now());\n",
       "    root._bokeh_is_loading = css_urls.length + js_urls.length;\n",
       "\n",
       "    function on_load() {\n",
       "      root._bokeh_is_loading--;\n",
       "      if (root._bokeh_is_loading === 0) {\n",
       "        console.debug(\"Bokeh: all BokehJS libraries/stylesheets loaded\");\n",
       "        run_callbacks()\n",
       "      }\n",
       "    }\n",
       "\n",
       "    function on_error() {\n",
       "      console.error(\"failed to load \" + url);\n",
       "    }\n",
       "\n",
       "    for (var i = 0; i < css_urls.length; i++) {\n",
       "      var url = css_urls[i];\n",
       "      const element = document.createElement(\"link\");\n",
       "      element.onload = on_load;\n",
       "      element.onerror = on_error;\n",
       "      element.rel = \"stylesheet\";\n",
       "      element.type = \"text/css\";\n",
       "      element.href = url;\n",
       "      console.debug(\"Bokeh: injecting link tag for BokehJS stylesheet: \", url);\n",
       "      document.body.appendChild(element);\n",
       "    }\n",
       "\n",
       "    const hashes = {\"https://cdn.bokeh.org/bokeh/release/bokeh-2.0.2.min.js\": \"ufR9RFnRs6lniiaFvtJziE0YeidtAgBRH6ux2oUItHw5WTvE1zuk9uzhUU/FJXDp\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-2.0.2.min.js\": \"8QM/PGWBT+IssZuRcDcjzwIh1mkOmJSoNMmyYDZbCfXJg3Ap1lEvdVgFuSAwhb/J\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-2.0.2.min.js\": \"Jm8cH3Rg0P6UeZhVY5cLy1WzKajUT9KImCY+76hEqrcJt59/d8GPvFHjCkYgnSIn\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-2.0.2.min.js\": \"Ozhzj+SI7ywm74aOI/UajcWz+C0NjsPunEVyVIrxzYkB+jA+2tUw8x5xJCbVtK5I\"};\n",
       "\n",
       "    for (var i = 0; i < js_urls.length; i++) {\n",
       "      var url = js_urls[i];\n",
       "      var element = document.createElement('script');\n",
       "      element.onload = on_load;\n",
       "      element.onerror = on_error;\n",
       "      element.async = false;\n",
       "      element.src = url;\n",
       "      if (url in hashes) {\n",
       "        element.crossOrigin = \"anonymous\";\n",
       "        element.integrity = \"sha384-\" + hashes[url];\n",
       "      }\n",
       "      console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n",
       "      document.head.appendChild(element);\n",
       "    }\n",
       "  };\n",
       "\n",
       "  function inject_raw_css(css) {\n",
       "    const element = document.createElement(\"style\");\n",
       "    element.appendChild(document.createTextNode(css));\n",
       "    document.body.appendChild(element);\n",
       "  }\n",
       "\n",
       "  \n",
       "  var js_urls = [\"https://cdn.bokeh.org/bokeh/release/bokeh-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-2.0.2.min.js\"];\n",
       "  var css_urls = [];\n",
       "  \n",
       "\n",
       "  var inline_js = [\n",
       "    function(Bokeh) {\n",
       "      Bokeh.set_log_level(\"info\");\n",
       "    },\n",
       "    function(Bokeh) {\n",
       "    \n",
       "    \n",
       "    }\n",
       "  ];\n",
       "\n",
       "  function run_inline_js() {\n",
       "    \n",
       "    if (root.Bokeh !== undefined || force === true) {\n",
       "      \n",
       "    for (var i = 0; i < inline_js.length; i++) {\n",
       "      inline_js[i].call(root, root.Bokeh);\n",
       "    }\n",
       "    if (force === true) {\n",
       "        display_loaded();\n",
       "      }} else if (Date.now() < root._bokeh_timeout) {\n",
       "      setTimeout(run_inline_js, 100);\n",
       "    } else if (!root._bokeh_failed_load) {\n",
       "      console.log(\"Bokeh: BokehJS failed to load within specified timeout.\");\n",
       "      root._bokeh_failed_load = true;\n",
       "    } else if (force !== true) {\n",
       "      var cell = $(document.getElementById(\"1056\")).parents('.cell').data().cell;\n",
       "      cell.output_area.append_execute_result(NB_LOAD_WARNING)\n",
       "    }\n",
       "\n",
       "  }\n",
       "\n",
       "  if (root._bokeh_is_loading === 0) {\n",
       "    console.debug(\"Bokeh: BokehJS loaded, going straight to plotting\");\n",
       "    run_inline_js();\n",
       "  } else {\n",
       "    load_libs(css_urls, js_urls, function() {\n",
       "      console.debug(\"Bokeh: BokehJS plotting callback run at\", now());\n",
       "      run_inline_js();\n",
       "    });\n",
       "  }\n",
       "}(window));"
      ],
      "application/vnd.bokehjs_load.v0+json": "\n(function(root) {\n  function now() {\n    return new Date();\n  }\n\n  var force = true;\n\n  if (typeof root._bokeh_onload_callbacks === \"undefined\" || force === true) {\n    root._bokeh_onload_callbacks = [];\n    root._bokeh_is_loading = undefined;\n  }\n\n  \n\n  \n  if (typeof (root._bokeh_timeout) === \"undefined\" || force === true) {\n    root._bokeh_timeout = Date.now() + 5000;\n    root._bokeh_failed_load = false;\n  }\n\n  var NB_LOAD_WARNING = {'data': {'text/html':\n     \"<div style='background-color: #fdd'>\\n\"+\n     \"<p>\\n\"+\n     \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n     \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n     \"</p>\\n\"+\n     \"<ul>\\n\"+\n     \"<li>re-rerun `output_notebook()` to attempt to load from CDN again, or</li>\\n\"+\n     \"<li>use INLINE resources instead, as so:</li>\\n\"+\n     \"</ul>\\n\"+\n     \"<code>\\n\"+\n     \"from bokeh.resources import INLINE\\n\"+\n     \"output_notebook(resources=INLINE)\\n\"+\n     \"</code>\\n\"+\n     \"</div>\"}};\n\n  function display_loaded() {\n    var el = document.getElementById(\"1056\");\n    if (el != null) {\n      el.textContent = \"BokehJS is loading...\";\n    }\n    if (root.Bokeh !== undefined) {\n      if (el != null) {\n        el.textContent = \"BokehJS \" + root.Bokeh.version + \" successfully loaded.\";\n      }\n    } else if (Date.now() < root._bokeh_timeout) {\n      setTimeout(display_loaded, 100)\n    }\n  }\n\n\n  function run_callbacks() {\n    try {\n      root._bokeh_onload_callbacks.forEach(function(callback) {\n        if (callback != null)\n          callback();\n      });\n    } finally {\n      delete root._bokeh_onload_callbacks\n    }\n    console.debug(\"Bokeh: all callbacks have finished\");\n  }\n\n  function load_libs(css_urls, js_urls, callback) {\n    if (css_urls == null) css_urls = [];\n    if (js_urls == null) js_urls = [];\n\n    root._bokeh_onload_callbacks.push(callback);\n    if (root._bokeh_is_loading > 0) {\n      console.debug(\"Bokeh: BokehJS is being loaded, scheduling callback at\", now());\n      return null;\n    }\n    if (js_urls == null || js_urls.length === 0) {\n      run_callbacks();\n      return null;\n    }\n    console.debug(\"Bokeh: BokehJS not loaded, scheduling load and callback at\", now());\n    root._bokeh_is_loading = css_urls.length + js_urls.length;\n\n    function on_load() {\n      root._bokeh_is_loading--;\n      if (root._bokeh_is_loading === 0) {\n        console.debug(\"Bokeh: all BokehJS libraries/stylesheets loaded\");\n        run_callbacks()\n      }\n    }\n\n    function on_error() {\n      console.error(\"failed to load \" + url);\n    }\n\n    for (var i = 0; i < css_urls.length; i++) {\n      var url = css_urls[i];\n      const element = document.createElement(\"link\");\n      element.onload = on_load;\n      element.onerror = on_error;\n      element.rel = \"stylesheet\";\n      element.type = \"text/css\";\n      element.href = url;\n      console.debug(\"Bokeh: injecting link tag for BokehJS stylesheet: \", url);\n      document.body.appendChild(element);\n    }\n\n    const hashes = {\"https://cdn.bokeh.org/bokeh/release/bokeh-2.0.2.min.js\": \"ufR9RFnRs6lniiaFvtJziE0YeidtAgBRH6ux2oUItHw5WTvE1zuk9uzhUU/FJXDp\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-2.0.2.min.js\": \"8QM/PGWBT+IssZuRcDcjzwIh1mkOmJSoNMmyYDZbCfXJg3Ap1lEvdVgFuSAwhb/J\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-2.0.2.min.js\": \"Jm8cH3Rg0P6UeZhVY5cLy1WzKajUT9KImCY+76hEqrcJt59/d8GPvFHjCkYgnSIn\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-2.0.2.min.js\": \"Ozhzj+SI7ywm74aOI/UajcWz+C0NjsPunEVyVIrxzYkB+jA+2tUw8x5xJCbVtK5I\"};\n\n    for (var i = 0; i < js_urls.length; i++) {\n      var url = js_urls[i];\n      var element = document.createElement('script');\n      element.onload = on_load;\n      element.onerror = on_error;\n      element.async = false;\n      element.src = url;\n      if (url in hashes) {\n        element.crossOrigin = \"anonymous\";\n        element.integrity = \"sha384-\" + hashes[url];\n      }\n      console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n      document.head.appendChild(element);\n    }\n  };\n\n  function inject_raw_css(css) {\n    const element = document.createElement(\"style\");\n    element.appendChild(document.createTextNode(css));\n    document.body.appendChild(element);\n  }\n\n  \n  var js_urls = [\"https://cdn.bokeh.org/bokeh/release/bokeh-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-2.0.2.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-2.0.2.min.js\"];\n  var css_urls = [];\n  \n\n  var inline_js = [\n    function(Bokeh) {\n      Bokeh.set_log_level(\"info\");\n    },\n    function(Bokeh) {\n    \n    \n    }\n  ];\n\n  function run_inline_js() {\n    \n    if (root.Bokeh !== undefined || force === true) {\n      \n    for (var i = 0; i < inline_js.length; i++) {\n      inline_js[i].call(root, root.Bokeh);\n    }\n    if (force === true) {\n        display_loaded();\n      }} else if (Date.now() < root._bokeh_timeout) {\n      setTimeout(run_inline_js, 100);\n    } else if (!root._bokeh_failed_load) {\n      console.log(\"Bokeh: BokehJS failed to load within specified timeout.\");\n      root._bokeh_failed_load = true;\n    } else if (force !== true) {\n      var cell = $(document.getElementById(\"1056\")).parents('.cell').data().cell;\n      cell.output_area.append_execute_result(NB_LOAD_WARNING)\n    }\n\n  }\n\n  if (root._bokeh_is_loading === 0) {\n    console.debug(\"Bokeh: BokehJS loaded, going straight to plotting\");\n    run_inline_js();\n  } else {\n    load_libs(css_urls, js_urls, function() {\n      console.debug(\"Bokeh: BokehJS plotting callback run at\", now());\n      run_inline_js();\n    });\n  }\n}(window));"
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "\n",
       "\n",
       "\n",
       "\n",
       "  <div class=\"bk-root\" id=\"e04e5d14-a67a-47de-b6ef-e35e363416a1\" data-root-id=\"1001\"></div>\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": [
       "(function(root) {\n",
       "  function embed_document(root) {\n",
       "    \n",
       "  var docs_json = {\"b5e34aba-00b6-48d6-8802-a395bc9ea058\":{\"roots\":{\"references\":[{\"attributes\":{\"background_fill_alpha\":0.02,\"background_fill_color\":\"blue\",\"below\":[{\"id\":\"1010\"}],\"center\":[{\"id\":\"1013\"},{\"id\":\"1017\"}],\"left\":[{\"id\":\"1014\"}],\"plot_height\":400,\"plot_width\":400,\"renderers\":[{\"id\":\"1035\"},{\"id\":\"1044\"},{\"id\":\"1049\"},{\"id\":\"1054\"}],\"title\":{\"id\":\"1057\"},\"toolbar\":{\"id\":\"1025\"},\"x_range\":{\"id\":\"1037\"},\"x_scale\":{\"id\":\"1006\"},\"y_range\":{\"id\":\"1038\"},\"y_scale\":{\"id\":\"1008\"}},\"id\":\"1001\",\"subtype\":\"Figure\",\"type\":\"Plot\"},{\"attributes\":{\"end\":3402044.0,\"start\":502156.0},\"id\":\"1038\",\"type\":\"Range1d\"},{\"attributes\":{\"data\":{\"x\":[1970,1980,1991,2000,2010],\"y\":{\"__ndarray__\":\"gL6KXnZTLkHAefpofXw2QUCwQXW5ij5BQFcy7gCRQkGA5AwLYjpGQQ==\",\"dtype\":\"float64\",\"shape\":[5]}},\"selected\":{\"id\":\"1068\"},\"selection_policy\":{\"id\":\"1069\"}},\"id\":\"1051\",\"type\":\"ColumnDataSource\"},{\"attributes\":{},\"id\":\"1019\",\"type\":\"WheelZoomTool\"},{\"attributes\":{\"overlay\":{\"id\":\"1024\"}},\"id\":\"1020\",\"type\":\"BoxZoomTool\"},{\"attributes\":{},\"id\":\"1021\",\"type\":\"SaveTool\"},{\"attributes\":{},\"id\":\"1066\",\"type\":\"Selection\"},{\"attributes\":{\"source\":{\"id\":\"1041\"}},\"id\":\"1045\",\"type\":\"CDSView\"},{\"attributes\":{\"data\":{\"x\":[1970,1980,1991,2000,2010],\"y\":[1002156,1449004,2015576,2447212,2902044]},\"selected\":{\"id\":\"1064\"},\"selection_policy\":{\"id\":\"1065\"}},\"id\":\"1041\",\"type\":\"ColumnDataSource\"},{\"attributes\":{\"fill_alpha\":{\"value\":0.1},\"fill_color\":{\"value\":\"black\"},\"line_alpha\":{\"value\":0.1},\"line_width\":{\"value\":4},\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1048\",\"type\":\"Square\"},{\"attributes\":{},\"id\":\"1022\",\"type\":\"ResetTool\"},{\"attributes\":{},\"id\":\"1067\",\"type\":\"UnionRenderers\"},{\"attributes\":{\"source\":{\"id\":\"1046\"}},\"id\":\"1050\",\"type\":\"CDSView\"},{\"attributes\":{\"data_source\":{\"id\":\"1051\"},\"glyph\":{\"id\":\"1052\"},\"hover_glyph\":null,\"muted_glyph\":null,\"nonselection_glyph\":{\"id\":\"1053\"},\"selection_glyph\":null,\"view\":{\"id\":\"1055\"}},\"id\":\"1054\",\"type\":\"GlyphRenderer\"},{\"attributes\":{\"line_color\":\"red\",\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1052\",\"type\":\"Line\"},{\"attributes\":{},\"id\":\"1068\",\"type\":\"Selection\"},{\"attributes\":{\"fill_color\":{\"value\":\"blue\"},\"line_color\":{\"value\":\"#1f77b4\"},\"radius\":{\"units\":\"data\",\"value\":1},\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1033\",\"type\":\"Circle\"},{\"attributes\":{},\"id\":\"1023\",\"type\":\"HelpTool\"},{\"attributes\":{},\"id\":\"1008\",\"type\":\"LinearScale\"},{\"attributes\":{},\"id\":\"1069\",\"type\":\"UnionRenderers\"},{\"attributes\":{\"axis_label\":\"Ano\",\"formatter\":{\"id\":\"1060\"},\"ticker\":{\"id\":\"1011\"}},\"id\":\"1010\",\"type\":\"LinearAxis\"},{\"attributes\":{},\"id\":\"1064\",\"type\":\"Selection\"},{\"attributes\":{},\"id\":\"1015\",\"type\":\"BasicTicker\"},{\"attributes\":{},\"id\":\"1063\",\"type\":\"UnionRenderers\"},{\"attributes\":{},\"id\":\"1006\",\"type\":\"LinearScale\"},{\"attributes\":{},\"id\":\"1011\",\"type\":\"BasicTicker\"},{\"attributes\":{},\"id\":\"1062\",\"type\":\"Selection\"},{\"attributes\":{\"data_source\":{\"id\":\"1046\"},\"glyph\":{\"id\":\"1047\"},\"hover_glyph\":null,\"muted_glyph\":null,\"nonselection_glyph\":{\"id\":\"1048\"},\"selection_glyph\":null,\"view\":{\"id\":\"1050\"}},\"id\":\"1049\",\"type\":\"GlyphRenderer\"},{\"attributes\":{\"axis\":{\"id\":\"1010\"},\"ticker\":null},\"id\":\"1013\",\"type\":\"Grid\"},{\"attributes\":{\"data\":{\"x\":[1975,1982],\"y\":{\"__ndarray__\":\"AAAAAGyzMkF00UUXka43QQ==\",\"dtype\":\"float64\",\"shape\":[2]}},\"selected\":{\"id\":\"1066\"},\"selection_policy\":{\"id\":\"1067\"}},\"id\":\"1046\",\"type\":\"ColumnDataSource\"},{\"attributes\":{},\"id\":\"1065\",\"type\":\"UnionRenderers\"},{\"attributes\":{\"axis\":{\"id\":\"1014\"},\"dimension\":1,\"ticker\":null},\"id\":\"1017\",\"type\":\"Grid\"},{\"attributes\":{\"axis_label\":\"Pop. urbana - PB\",\"formatter\":{\"id\":\"1039\"},\"ticker\":{\"id\":\"1015\"}},\"id\":\"1014\",\"type\":\"LinearAxis\"},{\"attributes\":{\"source\":{\"id\":\"1051\"}},\"id\":\"1055\",\"type\":\"CDSView\"},{\"attributes\":{\"fill_color\":{\"value\":\"black\"},\"line_width\":{\"value\":4},\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1047\",\"type\":\"Square\"},{\"attributes\":{},\"id\":\"1018\",\"type\":\"PanTool\"},{\"attributes\":{\"line_color\":\"green\",\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1042\",\"type\":\"Line\"},{\"attributes\":{\"fill_alpha\":{\"value\":0.1},\"fill_color\":{\"value\":\"blue\"},\"line_alpha\":{\"value\":0.1},\"line_color\":{\"value\":\"#1f77b4\"},\"radius\":{\"units\":\"data\",\"value\":1},\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1034\",\"type\":\"Circle\"},{\"attributes\":{\"source\":{\"id\":\"1032\"}},\"id\":\"1036\",\"type\":\"CDSView\"},{\"attributes\":{\"line_alpha\":0.1,\"line_color\":\"green\",\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1043\",\"type\":\"Line\"},{\"attributes\":{\"text\":\"\"},\"id\":\"1057\",\"type\":\"Title\"},{\"attributes\":{\"line_alpha\":0.1,\"line_color\":\"red\",\"x\":{\"field\":\"x\"},\"y\":{\"field\":\"y\"}},\"id\":\"1053\",\"type\":\"Line\"},{\"attributes\":{\"format\":\"%1.2e\"},\"id\":\"1039\",\"type\":\"PrintfTickFormatter\"},{\"attributes\":{\"data_source\":{\"id\":\"1032\"},\"glyph\":{\"id\":\"1033\"},\"hover_glyph\":null,\"muted_glyph\":null,\"nonselection_glyph\":{\"id\":\"1034\"},\"selection_glyph\":null,\"view\":{\"id\":\"1036\"}},\"id\":\"1035\",\"type\":\"GlyphRenderer\"},{\"attributes\":{\"data_source\":{\"id\":\"1041\"},\"glyph\":{\"id\":\"1042\"},\"hover_glyph\":null,\"muted_glyph\":null,\"nonselection_glyph\":{\"id\":\"1043\"},\"selection_glyph\":null,\"view\":{\"id\":\"1045\"}},\"id\":\"1044\",\"type\":\"GlyphRenderer\"},{\"attributes\":{},\"id\":\"1060\",\"type\":\"BasicTickFormatter\"},{\"attributes\":{\"bottom_units\":\"screen\",\"fill_alpha\":0.5,\"fill_color\":\"lightgrey\",\"left_units\":\"screen\",\"level\":\"overlay\",\"line_alpha\":1.0,\"line_color\":\"black\",\"line_dash\":[4,4],\"line_width\":2,\"render_mode\":\"css\",\"right_units\":\"screen\",\"top_units\":\"screen\"},\"id\":\"1024\",\"type\":\"BoxAnnotation\"},{\"attributes\":{\"active_drag\":\"auto\",\"active_inspect\":\"auto\",\"active_multi\":null,\"active_scroll\":\"auto\",\"active_tap\":\"auto\",\"tools\":[{\"id\":\"1018\"},{\"id\":\"1019\"},{\"id\":\"1020\"},{\"id\":\"1021\"},{\"id\":\"1022\"},{\"id\":\"1023\"}]},\"id\":\"1025\",\"type\":\"Toolbar\"},{\"attributes\":{\"end\":2020,\"start\":1960},\"id\":\"1037\",\"type\":\"Range1d\"},{\"attributes\":{\"data\":{\"x\":[1970,1980,1991,2000,2010],\"y\":[1002156,1449004,2015576,2447212,2902044]},\"selected\":{\"id\":\"1062\"},\"selection_policy\":{\"id\":\"1063\"}},\"id\":\"1032\",\"type\":\"ColumnDataSource\"}],\"root_ids\":[\"1001\"]},\"title\":\"Bokeh Application\",\"version\":\"2.0.2\"}};\n",
       "  var render_items = [{\"docid\":\"b5e34aba-00b6-48d6-8802-a395bc9ea058\",\"root_ids\":[\"1001\"],\"roots\":{\"1001\":\"e04e5d14-a67a-47de-b6ef-e35e363416a1\"}}];\n",
       "  root.Bokeh.embed.embed_items_notebook(docs_json, render_items);\n",
       "\n",
       "  }\n",
       "  if (root.Bokeh !== undefined) {\n",
       "    embed_document(root);\n",
       "  } else {\n",
       "    var attempts = 0;\n",
       "    var timer = setInterval(function(root) {\n",
       "      if (root.Bokeh !== undefined) {\n",
       "        clearInterval(timer);\n",
       "        embed_document(root);\n",
       "      } else {\n",
       "        attempts++;\n",
       "        if (attempts > 100) {\n",
       "          clearInterval(timer);\n",
       "          console.log(\"Bokeh: ERROR: Unable to run BokehJS code because BokehJS library is missing\");\n",
       "        }\n",
       "      }\n",
       "    }, 10, root)\n",
       "  }\n",
       "})(window);"
      ],
      "application/vnd.bokehjs_exec.v0+json": ""
     },
     "metadata": {
      "application/vnd.bokehjs_exec.v0+json": {
       "id": "1001"
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "from bokeh.io import output_notebook\n",
    "from bokeh.plotting import figure, show\n",
    "from bokeh.models import PrintfTickFormatter, Range1d\n",
    "\n",
    "from scipy import stats\n",
    "from scipy import interpolate\n",
    "\n",
    "# tabela de dados\n",
    "df = pd.read_csv(\"file-pop-pb.csv\")\n",
    "x = df['ano']\n",
    "y = df['pop. urbana']\n",
    "\n",
    "# grafico de dispersao\n",
    "f = figure(plot_width=400, plot_height=400)\n",
    "f.scatter(x,y,fill_color=\"blue\",radius=1,alpha=1)\n",
    "\n",
    "# formatacao\n",
    "f.background_fill_color = \"blue\"\n",
    "f.background_fill_alpha = 0.02\n",
    "f.x_range=Range1d(1960, 2020)\n",
    "e = 5.e5\n",
    "f.y_range=Range1d(y.min() - e,y.max() + e)\n",
    "f.xaxis.axis_label = \"Ano\"\n",
    "f.yaxis.axis_label = \"Pop. urbana - PB\"\n",
    "f.yaxis[0].formatter = PrintfTickFormatter(format=\"%1.2e\")\n",
    "\n",
    "# interpolacao\n",
    "\n",
    "f.line(x,y,color='green')\n",
    "\n",
    "# estimando o valor da populacao em 1975 e 1982 por interpolacao\n",
    "p = interpolate.interp1d(x, y)\n",
    "xx = [1975,1982]\n",
    "yy = p(xx)\n",
    "print(\"População obtida por interpolação linear em: 1975 = {:d}; 1982 = {:d}; \".format(int(yy[0]),int(yy[1])))\n",
    "f.square(xx,yy,color='black',line_width=4)\n",
    "\n",
    "# ajuste por minimos quadrados\n",
    "\n",
    "slope, intercept, r_value, p_value, std_err = stats.linregress(x,y)\n",
    "px = intercept + slope*x\n",
    "f.line(x,px,color='red')\n",
    "\n",
    "output_notebook()\n",
    "show(f)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pop. urbana</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ano</th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1970</th>\n",
       "      <td>1002156</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980</th>\n",
       "      <td>1449004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1991</th>\n",
       "      <td>2015576</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000</th>\n",
       "      <td>2447212</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2010</th>\n",
       "      <td>2902044</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      pop. urbana\n",
       "ano              \n",
       "1970      1002156\n",
       "1980      1449004\n",
       "1991      2015576\n",
       "2000      2447212\n",
       "2010      2902044"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.set_index('ano')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Exemplo\n",
    "\n",
    "Vamos resolver um problema de regressão linear por meio das equações normais passo a passo. Em primeiro lugar, importemos os módulos de computação numérica e de plotagem. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# importação de módulos \n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import sympy as sy\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Consideremos a simples tabela abaixo de um experimento fictício que busca correlacionar densidade e diâmetro médio de grãos em alimentos.\n",
    "\n",
    "densidade | diâmetro médio de grão\n",
    "-|-\n",
    "1.0|0.5\n",
    "2.1|2.5\n",
    "3.3|2.0\n",
    "4.5|4.2\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Vamos escrever os dados como _arrays_."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# tabela de dados\n",
    "x = np.array([1,2,3,4]) # densidade\n",
    "y = np.array([0.5,2.5,2.0,4.0]) # diâmetro"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Agora, calculamos os coeficientes linear $\\alpha_0$ e angular $\\alpha_1$ pelas fórmulas das equações normais vistas em aula."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "m = np.size(x)\n",
    "alpha1 = (m*np.dot(x,y) - np.sum(x)*np.sum(y))/(m*np.dot(x,x)-np.sum(x)**2)\n",
    "alpha0 = np.mean(y) - alpha1*np.mean(x)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Podemos agora escrever a equação da reta de regressão usando o _array_ `x` como abscissa. Este será o nosso _modelo de ajuste_."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "y2 = alpha0 + alpha1*x"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Enfim, plotamos o gráfico de dispersão dos valores _medidos_ juntamente com o _modelo de ajuste_ da seguinte forma:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXhV1bn48e+bEOZRQECGBARBEgKEgAwyKIoD1gmuMpQrKBe14nBt7a/Uq6UqbfXWoTjUiwNQDCoCKqVY1MogFdEwKYJYigwRNGEKhATJ8P7+WIeTgYSckJzsc07ez/OcJ3vtvXLOu7Phzcraa68lqooxxpjwF+V1AMYYY6qGJXRjjIkQltCNMSZCWEI3xpgIYQndGGMiRC2vPrhFixYaFxfn1ccbY0xYWr9+/QFVbVnaMc8SelxcHKmpqV59vDHGhCUR2V3WMetyMcaYCGEJ3RhjIoQldGOMiRCe9aGXJjc3l7S0NE6cOOF1KCaI6tatS7t27YiJifE6FGMiSkgl9LS0NBo1akRcXBwi4nU4JghUlYMHD5KWlkbHjh29DseYiBJSXS4nTpygefPmlswjmIjQvHlz+yvMmCAIqYQOWDKvAewaGxMcIZfQjTHGnJ2AE7qIRIvIRhFZWsqxOiLypojsEJF1IhJXlUGGq7i4OA4cOFDpOmWZM2cOU6dOPavvLWrfvn2MHj36rGPYt29fpWMwJlKlpEBcHERFua8pKcH7rIq00O8FtpVx7DbgsKp2Bp4GHq9sYKb6nHfeeSxcuPCsvtcSujFlS0mBKVNg925QdV+nTAleUg8ooYtIO2Ak8HIZVa4D5vq2FwLDJQw7Snft2kW3bt2YPHkyCQkJjB8/ng8//JBBgwbRpUsXPvvsMwAOHTrE9ddfT2JiIv379+eLL74A4ODBg4wYMYLevXtz++23U3Q1qNdee41+/frRq1cvbr/9dvLz80/7/KeeeoqEhAQSEhJ45plnSo1x9uzZXHDBBQwdOpR//vOf/v0ZGRmMGjWKvn370rdv32LHip7f4MGDSUpKIikpiU8++cS/PyEhATi91X/NNdewcuVK8vPzmThxIgkJCfTo0YOnn36ahQsXkpqayvjx4+nVqxc5OTmsX7+eoUOH0qdPH6644gr2799f0ctgTMR48EHIzi6+Lzvb7Q8KVS33hUvSfYBhwNJSjm8B2hUp/xtoUUq9KUAqkNqhQwctaevWrcV3DB2qOnu22z550pXnzXPl48dd+Y03XPnIEVdetMiVMzJceckSV96//7TPK+nbb7/V6Oho/eKLLzQ/P1+TkpJ00qRJWlBQoO+8845ed911qqo6depUnT59uqqq/uMf/9CePXuqqurdd9+tv/3tb1VVdenSpQpoRkaGbt26Va+55ho9efKkqqreeeedOnfuXFVVjY2N1YyMDE1NTdWEhATNysrSY8eOaffu3XXDhg3F4tu3b5+2b99e09PT9ccff9SBAwfqXXfdpaqqY8eO1Y8//lhVVXfv3q3dunU77fyOHz+uOTk5qqr6zTffaJ8+ffznHR8fr6qqs2fP9r+nqurIkSN1xYoVmpqaqpdddpl//+HDh1VVdejQofr555+rqurJkyd1wIABmp6erqqqb7zxhk6aNKnUn/Vp19qYCCSi6trmxV8iZ/+eQKqWkavLHYcuItcA6aq6XkSGlVWttN8VpfzymAXMAkhOTg7JxUw7duxIjx49AIiPj2f48OGICD169GDXrl0ArFmzhkWLFgFw6aWXcvDgQTIzM1m9ejWLFy8GYOTIkTRr1gyAf/zjH6xfv56+ffsCkJOTw7nnnlvsc9esWcMNN9xAgwYNALjxxhv5+OOP6d27t7/OunXrGDZsGC1buonWbr75Zr755hsAPvzwQ7Zu3eqve/ToUY4dO0ajRo38+3Jzc5k6dSqbNm0iOjra/72B6NSpEzt37uTuu+9m5MiRjBgx4rQ627dvZ8uWLVx++eUA5Ofn06ZNm4A/w5hI06GD62YpbX8wBPJg0SDgWhG5GqgLNBaR11T1p0XqpAHtgTQRqQU0AQ5VOrqVKwu3Y2KKl+vXL15u0qR4uUWL4uXWrQP6yDp16vi3o6Ki/OWoqCjy8vIAinWlnHKqh6m0niZV5ZZbbuH3v/99mZ9b2nuWpqyerIKCAtauXUu9evXK/N6nn36aVq1asXnzZgoKCqhbt+5pdWrVqkVBQYG/fGq8eLNmzdi8eTPLly/n+eefZ8GCBbz66qunnUN8fDxr164N6FyMiXQzZrg+86LdLvXru/3BUG4fuqpOU9V2qhoHjAE+KpHMAZYAt/i2R/vqhGQLvCoMGTKEFN9djZUrV9KiRQsaN25cbP97773H4cOHARg+fDgLFy4kPT0dcH3wu0v82h4yZAjvvPMO2dnZHD9+nLfffpvBgwcXq3PRRRexcuVKDh48SG5uLm+99Zb/2IgRI3juuef85U2bNp0Wd2ZmJm3atCEqKop58+aV2o8fFxfHpk2bKCgoYO/evf77BgcOHKCgoIBRo0bx6KOPsmHDBgAaNWrEsWPHAOjatSsZGRn+hJ6bm8tXX30V6I/VmIgzfjzMmgWxsSDivs6a5fYHw1k/+i8ij+D6cpYArwDzRGQHrmU+poriC0nTp09n0qRJJCYmUr9+febOdfeDf/Ob3zB27FiSkpIYOnQoHXx/V3Xv3p3HHnuMESNGUFBQQExMDM8//zyxsbH+90xKSmLixIn069cPgMmTJxfrbgFo06YN06dPZ8CAAbRp04akpCR/Up45cyZ33XUXiYmJ5OXlMWTIEF588cVi3/+zn/2MUaNG8dZbb3HJJZf4u3egsOU/aNAgf7dTQkICSUlJAHz33XdMmjTJ33o/9dfGxIkTueOOO6hXrx5r165l4cKF3HPPPWRmZpKXl8d9991HfHx81fzgjQlD48cHL4GXJF41pJOTk7XkAhfbtm3jwgsv9CSemmz9+vXcf//9rFq1qto+0661MWdHRNaranJpx+xJ0RouNTWVsWPHcu+993odijGmkkJqtkVT/ZKTkys02sUYE7qshW6MMRHCEroxxkQIS+jGGBMhLKEbY0yEsIQeRMOGDePU0Myrr76aI0eOlFn3hx9+4JJLLuHKK6/kscceq64QjTERxEa5VJNly5ad8XirVq1YsWJFNUVjjKkW+/dDy5ZQq3pSbVi30IMxcXwgU+geP36cW2+9lb59+9K7d2/effddwE26NWbMGBITE7n55pvJycnxv2/RhSz+8pe/kJiYSM+ePZk8eTIAL730En379qVnz56MGjWKbN/kD7t372b48OEkJiYyfPhw9uzZU/mTNMYE3/bt0LkzzJ5dfZ9Z1jSMwX6dmrq1qIpMqfraa6r16xefkrJ+fbe/MgKZQnfatGk6zzeN7+HDh7VLly6alZWlTz75pH+62M2bN2t0dLR/atlT0+Ru2bJFu3XrpgcOHFBV1YMHD6qq+suqqg8++KDOnDlTVVWvueYanTNnjqqqvvLKK/4pfMOdTZ9rIlJuruqpf9sFBaoPPaS6Y0eVfgRnmD43bFvowZw4/tRcJlFRUaVOofv+++/zhz/8gV69ejFs2DBOnDjBnj17WL16NT/9qZu3LDExkcTExNPe+6OPPmLUqFE0b94cgHPOOQeALVu2MHjwYHr06EFKSop/Uqu1a9cybtw4ACZMmMCaNWsqf4LGmOC4/XYYNgyystxsXI88AuefX20fH7Z96GX1PFRFj0R5U+hGR0ezaNEiunbtetr3lrdQk5Yxd87EiRN555136NmzJ3PmzGFl0al/K/D+xphqtnEjdOrkpvC+6y4YORKKTHxXncK2hV7WBPHBmji+qCuuuIJnn33Wn5w3btwIFJ9Wd8uWLf6l6Yo6NZXuoUNuuvhTX48dO0abNm3Izc31vwfAwIEDeeONNwBISUnh4osvDt6JGWMqZvduSE6Gp5925aQkuPFG1zr3QNgm9Bkz3ETxRQVz4viiHnroIXJzc0lMTCQhIYGHHnoIgDvvvJOsrCwSExN54okn/FPhFhUfH8+0adMYMmQIbdu25Ve/+hUAjz76KBdddBGXX3453bp189efOXMms2fPJjExkXnz5vGnP/0p+CdojCnb0aOwfLnbjo2F116D++7zNqZTyupcD/arsjdFVd0N0NhYtz5fbGzlb4hWt8mTJ2tBQYHXYXjCboqasPWzn6nWrataZCBDdaIya4qGsuqcOL6qDRkyhKNHj1JQUEB0dLTX4RhjyqIKS5ZAz55ufPS0aTBpEvgGNoSSsO1yCXerV6/2L9ZsjAlh6ekwdiw8+6wrt2vn+s1DUMgldI3cpUiNj11jE/J27SpM4K1awapV8PjjnoYUiHITuojUFZHPRGSziHwlIr8tpc5EEckQkU2+1+SzCaZu3bocPHjQ/sNHMFXl4MGD1K1b1+tQjCnbnDnwq1/Bvn2u3LdvtT2+XxmBRPgjcKmqZolIDLBGRN5T1U9L1HtTVadWJph27dqRlpZGRkZGZd7GhLi6devSrl07r8MwptDJk/DCC3DRRTBgAPziFzB5Mpx3nteRVUi5Cd13VzXLV4zxvYLShI6JiaFjx47BeGtjjClbbi488QSMG+cSesOG7hVmAupDF5FoEdkEpAMfqOq6UqqNEpEvRGShiLQv432miEiqiKRaK9wY46l16+COO9wolgYNYMMG+OMfvY6qUgJK6Kqar6q9gHZAPxFJKFHlr0CcqiYCHwJzy3ifWaqarKrJLVu2rEzcxhhTOdu2ueGIu3e7cuvW3sZTBSo0ykVVjwArgStL7D+oqj/6ii8BfaokOmOMqSo5Oe6JzjffdOUJE+Bf/3JjyyNEIKNcWopIU992PeAy4OsSddoUKV4LbKvKII0xptLq1IFPPoGtW105OtqzSbSCJZBRLm2AuSISjfsFsEBVl4rII7hHUJcA94jItUAecAiYGKyAjTEmYMuXuwmeli+HevVgzRqoXdvrqIImkFEuXwC9S9n/cJHtacC0qg3NGGPOUkGBW8qsTh03N/n+/W6K2whO5hCCT4oaY8xZO3kSrrsOfv97Vx42DNavd8m8BrCEbowJfydPuq+1a0OzZsXHkNegRWEsoRtjwtvbb0PHjvD99648Zw7ce6+nIXnFEroxJvzk57u+cYCEBOjfv7CVXoNZQjfGhJf8fBg4EP77v125SxdYtKh61p8McaE/fZgxxoAbqdKmjRs/PmqUW/7NFGMtdGNM6Fu40LXAN2925V/+Em6+2duYQpAldGNMaMrKgr173fbw4XD//W61IFMm63IxxoQeVddP3qoVfPCBG4oYBisGec0SujEmNKjC2rVuPnIRmD7d9ZmbgFmXizEmNLz9NgwaBH//uyvfeKNL7iZgltCNMd5JS4PUVLf9k5/Ayy+7/nJzVqzLxRjjDVU370p+PmzcCDExcNttXkcV1qyFboypPrm5MHu2e6pTBF58Ed55p0bNtxJM1kI3xlSflSvh1lvdwhI33QR9+3odUUSxFroxJrhSU2HxYrd92WWwejX8x394G1OEsha6MSa4HnwQ9uyB6693i04MHux1RBHLWujGmKp15IhL4ocOufLLL8O6dS6Zm6Cyn7Axpmrt3QtPPOHW8QRo3x4aN/Y2phqi3IQuInVF5DMR2SwiX4nIb0upU0dE3hSRHSKyTkTighGsMSZELV7skjhAjx7w7bcwdqy3MdVAgbTQfwQuVdWeQC/gShHpX6LObcBhVe0MPA3YpAvG1CTLl8Mbb7hhiWCTaHmk3ISujm9pEGJ8Ly1R7Tpgrm97ITBcxAaWGhOxvvsOxoyBb75x5SefhM8/dw8HGc8E1IcuItEisglIBz5Q1XUlqrQF9gKoah6QCTQv5X2miEiqiKRmZGRULnJjjHdq1YJVqwrnJ2/Y0C08YTwVUEJX1XxV7QW0A/qJSEKJKqW1xku24lHVWaqarKrJLVu2rHi0xhjv/PnPMHGi227VCnbtsvHkIaZCo1xU9QiwEriyxKE0oD2AiNQCmgCHqiA+Y4yXtEi7LDMT0tPhxAlXrlPHm5hMmQIZ5dJSRJr6tusBlwFfl6i2BLjFtz0a+EhVT2uhG2PCyO7dcNFF8NFHrvzLX8KyZVC3rrdxmTIF8qRoG2CuiETjfgEsUNWlIvIIkKqqS4BXgHkisgPXMh8TtIiNMcGVm+tubrZqBbVrQ06O228PBoU88aohnZycrKmn5kE2xoSGGTPc7Ieffmo3OUOUiKxX1eTSjtmvXGNqupwcNyc5QJcu0K9fYT+5CSuW0I2pyfbsgQsugLm+x0huugmef95Nb2vCjiV0Y2qi9HT3tX17GDkSunXzNh5TJSyhG1PTTJsGPXtCVlbhqkEDB3odlakCNh+6MTXB99+7pzkbNnTreDZvbo/pRyBroRsT6X74wd3sPDUbYv/+8Itf2INBEcgSujGRKD8f1q93261aueGI//mf3sZkgs4SujGR6MEH4eKLYf9+V77nHujc2duYTNBZH7oxkWLzZtc33q4d3Hkn9OkDrVt7HZWpRtZCNyYSHDniRqo8+qgrx8a6mRBtWYIaxRK6MeHq6FFYsMBtN20KCxfCH/7gbUzGU5bQjQlXM2e6VYN27nTlq66CZs28jcl4yhK6MeFCFZYsgQ0bXPnee2HdOujUydu4TMiwhG5MuMjOhilT4JlnXLlRI+jb19uYTEixhG5MKNu9G6ZPd63zBg3cYhOvvOJ1VCZEWUI3JpR9+CE8/jhs3erK3bvbI/umTJbQjQkleXnuZue777ryxInwr39BfLynYZnwYAndmFAi4rpUlixx5eho96CQMQEIZJHo9iKyQkS2ichXInJvKXWGiUimiGzyvR4OTrgmVKWkQFycW3YyLs6VTYA2bICbb3arBEVHw4oV8PLLXkdlwlAgj/7nAT9X1Q0i0ghYLyIfqOrWEvU+VtVrqj5EE+pSUtzgi+xsV96925UBxo/3Lq6wcfgwrF4N27e7ecrPOcfriEyYKreFrqr7VXWDb/sYsA1oG+zATPh48MHCZH5Kdrbbb0qRnw/33w9PPeXKw4e7h4N69vQ2LhP2KtSHLiJxQG9gXSmHB4jIZhF5T0RKvYMjIlNEJFVEUjMyMiocrAlNe/ZUbH+Npeq+RkfDt99CWlrhsXr1vInJRJSAE7qINAQWAfep6tEShzcAsaraE3gWeKe091DVWaqarKrJLVu2PNuYTYjp0KFi+2uk1auhV6/CtTwXLixsoRtTRQJK6CISg0vmKaq6uORxVT2qqlm+7WVAjIi0qNJITciaMQPq1y++r359t7/Gy8tzX1u1gtq1CxN6dLR3MZmIFcgoFwFeAbapaqlNChFp7auHiPTzve/BqgzUhK7x42HWLDdjq4j7OmtWDb8hqupGrkyd6spdu8Jnn0FCgrdxmYgWyCiXQcAE4EsR2eTb92ugA4CqvgiMBu4UkTwgBxijeqrD0NQE48fX8AR+Sna2+/NExK0Q1KhR4TGbm9wEmXiVd5OTkzU1NdWTzzYmKD78EG66CT7+2J7sNEEjIutVNbm0Y/akqDGVUVAAhw657d694coroW5db2MyNZatKWrM2VKFESOgVi34+9/dep7z53sdlanBLKEbU1HffuvmNxCBCRPc7Ieq1kduPGddLsZUxEcfuZud773nyrfcAuPGWTI3IcESujHlOX4cvv7abV98MTz0EPTr521MxpTCulyMKc8117gHgr780j0cNH261xEZUypL6MaU5qOPYPBg1z/+m9+4r1H2B60JbfYv1JiS1qxxMyDOm+fKw4bBoEGehmRMICyhGwPw3Xfwj3+47UGD4PXX7dFXE3asy8UYcCtybN7shiTGxMCYMV5HZEyFWQvd1Ex5eW4GsSNHXPmpp9wj+zEx3sZlTCVYC93UTNu2wR13wMmTbkbErl29jsiYSrOEbmqO9evdgsz/9V/Qowekprr5V4yJENblYmqOF16ARx6BEydcOSnJnvA0EcUSuolcR4/CtGnw73+78uOPw5YtNhuiiViW0E3kOn4cnnsOli935RYtoEkTb2MyJoisD91Elr/+FVascKNW2rSBnTvBFiQ3NYS10E1k2bzZtciPHXNlS+amBrGEbsJbRoZbjHn1ald+4AGX1Iuu5WlMDVFuQheR9iKyQkS2ichXInJvKXVERGaKyA4R+UJEkoITrjE+p9bCbdDAJfAdO1y5Th23glANlpLi1t+IinJfU1K8jshUl0D+5ecBP1fVDSLSCFgvIh+o6tYida4CuvheFwF/9n01purNnu2Welu+HOrXdyNXangSPyUlxc1ikJ3tyrt3uzLY1DQ1QbktdFXdr6obfNvHgG1A2xLVrgP+os6nQFMRaVPl0ZqaS9UtyAxuTvJatQof27dk7vfgg4XJ/JTsbLffRL4K9aGLSBzQG1hX4lBbYG+RchqnJ31EZIqIpIpIakZGRsUiNTVXRgYMGFA4ne24cW4JuHPO8TauELRnT8X2m8gScEIXkYbAIuA+VT1a8nAp36Kn7VCdparJqprc0kYfmPKceqKzeXM3BLFhQ1e2pzvL1KFDxfabyBJQQheRGFwyT1HVxaVUSQPaFym3A/ZVPjxTY82cCRde6PoLoqLg7bdh1Civowp5M2a42wpF1a/v9pvIF8goFwFeAbap6lNlVFsC/KdvtEt/IFNV91dhnKYmyMlxL3DzrFx5pZsN0QRs/Hg3K3BsrPtDJjbWle2GaM0gqqf1jBSvIHIx8DHwJeC7K8WvgQ4AqvqiL+k/B1wJZAOTVDX1TO+bnJysqalnrGJqkiNHoGdPuPVWt4anMaZUIrJeVZNLO1bu8ABVXUPpfeRF6yhw19mFZ2q0PXtcB2/TpnDLLW79TmPMWbEnRY13/vd/XT/5d9+58iOPwNCh3sZkTBizAbymev3wg/vaqpW7yanqRrEYYyrNWuim+mRnQ0KCm6McoFMn+OUvbX5yY6qItdBNcOXnu8WXhw1z4+eefhouslkhjAkGa6Gb4PrTn+CSS+Crr1z5pz+FLl28jcmYCGUtdFP1vvgCoqMhPh5uu82NYune3euojIl41kI3VevkSbjiCvj1r125SRMYPdoe1zemGlhCN5V37Jh7HFHVzYS4cCHMmeN1VMbUOJbQTeUtWAC33w6ff+7KgwZBs2bexmRMDWR96KbiVOFvf3OrA11+uXvCs1cv6NPH68iMqdEsoZuKKyiAX/3K3ey8/HK3wIQlc2M8Z10uJjB798L998OPP7oRLEuXwrvveh2VMaYIS+gmMNu2wQsvFPaTx8VBTIynIRljirMuF1M6VXjuOddPPmUKjBjhVhxu1crryIwxZbAWuimdiFu38/33C/dZMjcmpFlCN4W2bIGRI+HAAVdesADeesvbmIwxAbOEbgpFRcHmzfD1167csKE94WlMGLE+9JpM1U1le/IkPPWUm2/l22/tZqcxYSqQRaJfFZF0EdlSxvFhIpIpIpt8r4erPkxTpQp8S8OKwPHj7nVqbVlL5saErUC6XObgFn8+k49VtZfv9UjlwzJBk5oK3bq5YYgAM2fC//2fda0YEwHKTeiquho4VA2xmGD68Uf3NTYWWrd2rXKwRG5MBKmqm6IDRGSziLwnIvFlVRKRKSKSKiKpGRkZVfTRply33QY33ui2W7aE1ashOdnbmIwxVa4qEvoGIFZVewLPAu+UVVFVZ6lqsqomt2zZsgo+2pQpM7NwOykJ+vcv7Ds3xkSkSid0VT2qqlm+7WVAjIi0qHRk5ux9/rmbOOuDD1z5rrvgoYfcsERjTMSq9P9wEWkt4jpiRaSf7z0PVvZ9TQUVFMC+fW47MRFuuskldWNMjVHuOHQReR0YBrQQkTTgN0AMgKq+CIwG7hSRPCAHGKN6agycqTZjxriRK5s2uflXXnrJ64iMMdWs3ISuqmPLOf4c8FyVRWQCt307nH++m4984kQ4dMhGrRhTg1mnarjatAni4+HVV1356qvhpz+1fnJjajD73x9Ojh93DwYB9OwJTzwBN9zgbUzGmJBhc7mEk0mTYNUq2LUL6tVzKwgZY4yPtdBD3QcfFI4p/5//gcWLXTI3xpgSLKGHsm++cSsFPf+8KycmwqBB3sZkjAlZltBDzb59sGiR277gAvjb3+DnP/c2JmNMWLCEHmoefhhuvRWOHXPlq69248qNMaYcltC9lpcHs2a5hSUAHnkENm6ERo28jcsYE3YsoXstPR3uuw/mznXl886DTp28jckYE5Zs2KIXNm6E996DX//aJfANG6BrV6+jMsaEOWuhe+Hdd+Hpp+Ggbw6zbt3skX1jTKVZQq8OWVmuNf7JJ678wAOwYwc0b+5tXMaYiGJdLtVBBP7yF2jYEAYOhAYNvI7IGBOBLKEHy/Ll7kbna6+5BL51KzRu7HVUxpgIZl0uwfL99+7m5/ffu7Ilc2NMkFlCrypHj8LYsTB/vitPmABffulGsRhjTDWwLpfKUnV95A0bQloaHDjg9kdF2dzkxphqZRmnMhYsgN69ISfHJe9Vq+Cee7yOyhhTQ1lCryhVyM11261audfhw65sLXJjjIfKzUAi8qqIpIvIljKOi4jMFJEdIvKFiCRVfZghIjvbTV/7+OOuPHSoG81i/eTGmBAQSJNyDnDlGY5fBXTxvaYAf658WCHm1MyH9etDr14214oxJiSVm9BVdTVw6AxVrgP+os6nQFMRaVNVAXpu7lyIjS0cfvjCCzBunLcxGWNMKaqi07ctsLdIOc237zQiMkVEUkUkNSMjowo+OkhOnCjsFx80CG66CaKjvY3JGGPKURUJvbRZpbS0iqo6S1WTVTW5ZcuWVfDRQXDypFvq7YEHXLlzZ3jxRQjVeI0xxqcqxqGnAe2LlNsB+6rgfavX9u1uCtvatd3Qw+7dvY7ImGrz3Q/f8fp7r7Pru100qNeAwX0GM3LISKSMWUAPHDnAslXL2L5rO5lZmTRp2ITk+GRGDh1J7Zja/nq3//b207533MhxDE0eGrRz2Ze+j7+u/Ct79u/hwJEDXDP0Gn4y7Cflft+hzEPMXzaf7d9uJ6ZWDP169GP0iNHUii5Mk9/s+oa33n+Lfen7aNqoKSMGjQjquVRUVST0JcBUEXkDuAjIVNX9VfC+1efVVyirSaUAAA35SURBVGHyZPeofs+eMHWq1xEZU21yfszhmXnP0CW2C9P+axo/HPiBOe/OoU5MHS4feHmp3/P9ge8p0ALGjRzHuc3P5fuM73lt6Wtk5WQx4ScTitWd8JMJ9Ligh79cr069CsX35JwnGdBrAAN7DQyo/snckzRv2pzeF/bm3Y/eDeh7CgoKeHb+szSs15AHJj1AVnYWc96Zg6oy9uqxABw4fIBn5z/LwF4DufWGW9mxZwfzl82nUf1GJHUPjcF9gQxbfB1YC3QVkTQRuU1E7hCRO3xVlgE7gR3AS8DPghZtVUpPh5073fYNN7ihiBdc4G1MpkZbu3kt9z9xP7l5ucX2v7L4FZ5//fmgfe5nX3zGydyTTLx+Im3PbUtS9ySuGHQFH376Iaql9p6S0DmBiddPJL5zPC2btaTHBT24avBVbNy28bS69erWo0nDJv5X0RZ8MMS1jWP0iNH069Ev4M/a+u+t7E/fz6QbJtGhTQe6n9+dGy+/kTUb1pDzYw4Aq1JX0bRRU8ZePZY2LdswuM9gBvQcwPtr3wfg2PFjPPDkAyxdtdT/vmk/pHHXY3exfuv6qj/RUpTbQlfVseUcV+CuKouoOuTnQ//+0KWLG0ferFlhn7kxHunTvQ8L/r6Azds3kxyfDEDOiRw2btvI5FGTAfjX7n/xbMqzZ3yfKwdfydWDrw74c3em7aRzbOdiyS/+/HiWrFjCwSMHadGsRUDvc+LHE9SvW/+0/Qv+voCUpSm0aNaCQb0HMbjPYKIktB7C25m2k9YtW3NOk3P8++LPjycvP489+/bQtWNXdqbt5MJOFxb7vvjz41m7eS35+fk0atCIiddN5PnXn+fCThfSvnV7Xl70Mn179KVP9z7Vch41Zy6X/Hy37NvIkW7EynPPwfnnex2VMX61Y2rTr0c//rnxn/6E/tmXn1GvTj1/l0XsebH8zx3/c8b3aVCvYvPtZ2Zl0qxxs2L7Gjds7D8WSEI/lHmI9z95n6sGX1Vs/7XDrqVrx67UqV2Hr3d+zcL3F5KVncXIISMrFGOwZWZl0rhB8RlRG9ZvSJREkZmVCcDRrKOnJfRGDRtRUFBAVnYWTRo1Ib5zPEP7DuXVxa/SJa4LeXl5jLlqTLWdR81J6K+/7mZA/OgjuOQSuDrwFowx1eXipIuZMWsGh48eplnjZvxz0z/p37M/0VFu2GztmNqce865Z/3+01+YzqEj7rGSzrGduWe8m3tISgxWO9XVUtZN0aKOZh3lT6/9ie7nd+ey/pcVOzZyaGHibt+6PQVawLKPl50xoacsTWHdF+v85ZN5J9mZtpM3lr1ReB53TS/Wmq4SZZxq0Z9ByZ+Tfzxfkd03XnYjX+34ik83f8r/u/X/Ubd23aqN8wwiO6F/+aWb1nbQILj5Zjcj4rBhXkdlTJnat25PhzYd+GTTJ/Tq1ovd+3Zz6w23+o9Xtsvl7nF3k1+QD0BMrRgAmjRsQubxzGL1jh13T0eXbLWWlJmVydNzn+a8c89j0g2Tyv0F0LFtR078eIKjWUf9fwWU9JNLflLsZuwri18h6cIkel/Y27+vSaMmZ/ycimrSsAn/3vvvYvuysrMo0AL/z6Bxw8b+1vopx44fIyoqiob1Gvr3HTxykMNHDyMIGYcz6NiuY5XGeiaRm9BV3ROdDRrAp59CTAxcf73XURlTrouTLub9f75PVnYW57c/n9YtWvuPVbbLpXnT09ex7dSuE4s/XExuXq4/yW/duZWmjZqWWv+UzGOZPDn3Sc479zwmj5rs/yviTPZ+v5eYWjHUq1v2SJfGDRoX+0VSu1ZtGjVoVKm/TMrTqV0nlq1e5v/LCGDbzm3Uiq5Fh/M6+Ots3r652Pdt3bmV2PNiifY9eJifn88ri1+hZ9eexLWNY/7f5tO5Q+eq/2uiDKF1Z6KysrLgySfdw0EibrGJZcu8jsqYCumX0I+jWUdZnbqaQb0HFTt2qsvlTK+K9qGfGg0y5505fJf+HRu2bWD5muVc1v8yf4t747aNPPzcwxw+6p6gPnLsCH+c80eaNGzCTVfcRFZ2FplZmWRmZVJQUADA5u2b+Xj9x3yX/h0ZhzJYs2ENS1YuYXDSYP8vjmDIy89j7/d72fv9XnLzcjmadZS93+8l/VC6v86Kz1bw8HMP+8vdz+9Om3PbMPvt2ezZv4dtO7ex6INFXJx0sX+Y5dDkoRw+epg3//4m+zP2s2bDGtZuWsuIASP87/Puinc5dvwY40aOY/hFw+nYriOvvv0qBVoQtPMtKrJa6GvWwC9+Ad26uZufPXqU/z3GhJi6derSJ74P67eu998cDaZ6detx34T7mL9sPr+b9Tvq16vPZQMu47IBhf3hOT/m8MPBH/zdNVv/vZX0Q+mkH0pn2jPTir3fjHtn0KJpC6Kjo1mVuoq33n8LVaVFsxZcO+xahvUbFtTzOXLsCI/932P+csb6DFavX80FsRfw84k/B1x3yg8Hf/DXiYqK4u5xdzP/b/N54tUnqB1Tm74JfRk9YrS/TotmLbh73N0sWL6A1amradKoCTdfdbN/DPo3u77hg7UfcN+E+/yjfSZeP5FH//woy9csP+2GcTBIWeNMgy05OVlTU1Mr9yaqbuTKkSOue0UVvvoKEhKqJkhjPDIzZSbNGjVjwrUTyq9sahQRWa+qpf6mD+8uFxF45hmYObNwKThL5iaMHc85TupXqWz991Yu7X+p1+GYMBP+XS7z5rkHgwIYXmVMqJvxfzM4nnOc6y+9nrbnljppqTFlCv+E3qqV1xEYU2V+d9/vvA7BhLHw7nIxxhjjZwndGGMihCV0Y4yJEJbQjTEmQlhCN8aYCGEJ3RhjIoQldGOMiRCW0I0xJkIElNBF5EoR2S4iO0TkV6UcnygiGSKyyfeaXPWhFpeSAnFxEBXlvqakBPsTjTEmtJX7pKiIRAPPA5cDacDnIrJEVbeWqPqmqk4NQoynSUmBKVMgO9uVd+92ZYDx46sjAmOMCT2BtND7ATtUdaeqngTeAK4Lblhn9uCDhcn8lOxst98YY2qqQBJ6W2BvkXKab19Jo0TkCxFZKCLtS3sjEZkiIqkikpqRkXEW4Tp79lRsvzHG1ASBJPTSpjEsOYn6X4E4VU0EPgTmlvZGqjpLVZNVNblly5YVi7SIDh0qtt8YY2qCQBJ6GlC0xd0O2Fe0gqoeVNUffcWXgD5VE17pZsyA+vWL76tf3+03xpiaKpCE/jnQRUQ6ikhtYAywpGgFEWlTpHgtsK3qQjzd+PEwaxbExrpp0GNjXdluiBpjarJyR7moap6ITAWWA9HAq6r6lYg8AqSq6hLgHhG5FsgDDgETgxgz4JK3JXBjjCkU3muKGmNMDRO5a4oaY4zxs4RujDERwhK6McZECEvoxhgTISyhG2NMhLCEbowxEcISujHGRAjPxqGLSAawuwreqgVwoArex2t2HqElEs4jEs4B7DxKilXVUifD8iyhVxURSS1rkH04sfMILZFwHpFwDmDnURHW5WKMMRHCEroxxkSISEjos7wOoIrYeYSWSDiPSDgHsPMIWNj3oRtjjHEioYVujDEGS+jGGBMxwiahi8irIpIuIlvKOC4iMlNEdvgWq06q7hjLE8A5DBORTBHZ5Hs9XN0xBkJE2ovIChHZJiJfici9pdQJ6esR4DmE/PUQkboi8pmIbPadx29LqVNHRN70XYt1IhJX/ZGeWYDnMVFEMopcj8lexFoeEYkWkY0isrSUY8G9FqoaFi9gCJAEbCnj+NXAe7hFrfsD67yO+SzOYRiw1Os4AziPNkCSb7sR8A3QPZyuR4DnEPLXw/fzbejbjgHWAf1L1PkZ8KJvewzwptdxn+V5TASe8zrWAM7lfmB+af92gn0twqaFrqqrccvbleU64C/qfAo0LbHWqecCOIewoKr7VXWDb/sYbg3ZtiWqhfT1CPAcQp7v55vlK8b4XiVHOlwHzPVtLwSGi4hUU4gBCfA8Qp6ItANGAi+XUSWo1yJsEnoA2gJ7i5TTCMP/oMAA35+d74lIvNfBlMf3J2NvXIuqqLC5Hmc4BwiD6+H7E38TkA58oKplXgtVzQMygebVG2X5AjgPgFG+LryFItK+mkMMxDPAL4GCMo4H9VpEUkIv7bdcuP2G34Cbp6En8CzwjsfxnJGINAQWAfep6tGSh0v5lpC7HuWcQ1hcD1XNV9VeQDugn4gklKgSFtcigPP4KxCnqonAhxS2dEOCiFwDpKvq+jNVK2VflV2LSEroaUDR39jtgH0exXJWVPXoqT87VXUZECMiLTwOq1QiEoNLhCmquriUKiF/Pco7h3C6HgCqegRYCVxZ4pD/WohILaAJIdz1V9Z5qOpBVf3RV3wJ6FPNoZVnEHCtiOwC3gAuFZHXStQJ6rWIpIS+BPhP3+iK/kCmqu73OqiKEJHWp/rTRKQf7voc9Daq0/lifAXYpqpPlVEtpK9HIOcQDtdDRFqKSFPfdj3gMuDrEtWWALf4tkcDH6nvrlyoCOQ8StyDuRZ33yNkqOo0VW2nqnG4G54fqepPS1QL6rWoVVVvFGwi8jpu1EELEUkDfoO7cYKqvggsw42s2AFkA5O8ibRsAZzDaOBOEckDcoAxofYfz2cQMAH40tfnCfBroAOEzfUI5BzC4Xq0AeaKSDTuF84CVV0qIo8Aqaq6BPeLa56I7MC1Bsd4F26ZAjmPe0TkWiAPdx4TPYu2AqrzWtij/8YYEyEiqcvFGGNqNEvoxhgTISyhG2NMhLCEbowxEcISujHGRAhL6MYYEyEsoRtjTIT4/9El4nXb0Kf9AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mod = plt.plot(x,y2,'r:'); # modelo\n",
    "med = plt.scatter(x,y,c='b'); # medição\n",
    "plt.legend({'modelo de ajuste':mod, 'medição':med}); # legenda\n",
    "\n",
    "# esta linha adiciona a equação de ajuste ao gráfico na posição (x,y) = (2.8,0.8)\n",
    "# com fonte tamanho 14 e cor RGB = [0.4,0.5,0.4].\n",
    "plt.annotate('y= {0:.2f} + {1:.2f}x'.format(alpha0,alpha1),(2.8,0.8),fontsize=14,c=[0.4,0.5,0.4]);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Na prática, podemos calcular regressão linear usando o módulo `scipy.stats`. Vide _Code session 7_."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Execução de Código Implementado"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "from IPython.core.display import HTML\n",
    "\n",
    "def css_styling():\n",
    "    styles = open(\"styles/custom.css\", \"r\").read()\n",
    "    return HTML(styles)\n",
    "css_styling();"
   ]
  }
 ],
 "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.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}