{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "This short guide will help you get started with the anaconda python distribution on your PC or MAC. If you happen to run linux, please contact Professor Hicks directly.\n", "\n", "# Installing Python\n", "\n", "Both the PC and MAC installers can be downloaded at [http://continuum.io/downloads](http://continuum.io/downloads).\n", "\n", "Take care that you are downloading Python 3.x (which is the default) for your platform. The 64bit version (also the default) is what you want in 99% of all cases. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Once you have downloaded the installer, do the following:\n", "\n", "1. Navigate to your your browser's download folder \n", "2. Open the python executable (e.g. Anaconda3-4.3.1-Windows-x86_64.exe [PC] or Anaconda3-4.3.1-MacOSX-x86_64.pkg [MAC])\n", "3. Click through the install screens and choose to install \"Just Me\" and don't change the default install location. Leave any \"Advanced Options\" at their default values. \n", "\n", "\n", "# Using Python\n", "\n", "Now that python is installed and updated, we can use it. There are several ways to run python code, and I will show you a method allowing you to use your web browser as a \"front-end\" for python.\n", "\n", "We will be using the Ipython Notebook (now called Jupyter). The easiest way to open it is to locate and use the Anaconda Launcher. Hopefully, the install process will put this in your `Program Files` (PC) or `Applications` (Mac). If you see a new icon labeled Anaconda Launcher, click it and choose ipython notebook. Sometimes, the install process doesn't create the Anaconda icon. In this case, it is possible to locate it in the Anaconda folder in your local home directory.\n", "\n", "In case you can't find the Anaconda icon, you can also launch it from the terminal or command prompt window, type ``ipython notebook``. \n", "\n", "If you have made it this far, your web browser should pop up and you should see a window like this (note: the actual filenames and folders you will see will reflect the actual information you have on your computer):\n", "![](../site_pics/ipython_main.png)\n", "\n", "This listing allows you to see existing and running Ipython Notebooks. You won't have any, so let's create one. You can navigate to a folder of your choosing and create a new notebook by clicking on `New` and then `Python 3` [note the screenshot below is old and show `Python 2`] as shown below:\n", "\n", "![](../site_pics/ipython_main2.png)" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "## Testing the installation\n", "\n", "Having created a new notebook, you should see a window that looks like this:\n", "![](../site_pics/notebook1.png)\n", "\n", "You can re-name your new untitled notebook by clicking on ``Untitled`` and changing the name:\n", "![](../site_pics/notebook2.png)" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "Ipython Notebooks have execution blocks where you can enter and execute python code. These blocks look like this:\n", "![](../site_pics/execute.png)\n", "\n", "Copy and paste the following code into the first (and only at this point) execution block:\n", "```\n", "%matplotlib inline\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib pyplot as plt\n", "import statsmodels.formula.api as sm\n", "import sympy as sympy\n", "```\n", "\n", "You should see output like this after pressing the `shift enter` keys at the same time to execute the code:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import statsmodels.formula.api as sm\n", "import sympy as sympy" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "Next, we will make sure things are working. Paste this code into the next execution cell and press ``shift enter``:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
timeprice
count100.000000100.000000
mean50.5000009.866077
std29.0114922.946530
min1.0000002.643740
25%25.7500008.113189
50%50.5000009.724360
75%75.25000011.569507
max100.00000018.147543
\n", "
" ], "text/plain": [ " time price\n", "count 100.000000 100.000000\n", "mean 50.500000 9.866077\n", "std 29.011492 2.946530\n", "min 1.000000 2.643740\n", "25% 25.750000 8.113189\n", "50% 50.500000 9.724360\n", "75% 75.250000 11.569507\n", "max 100.000000 18.147543" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# create some toy data\n", "time=np.arange(1,101)\n", "price = np.random.normal(10,3,(np.size(time)))\n", "\n", "# load into pandas dataframe\n", "df = pd.DataFrame(zip(time,price),columns=['time','price'])\n", "\n", "# run summary statistics\n", "df.describe()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can also visualize our data easily using the pandas package:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAEACAYAAABVtcpZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEbVJREFUeJzt3XusZWV9xvHvI4MVBB3wMiCMHUq03lCwQqnacmrVkNYi\nrY2pbVO8xJrGChJbBUzlpEnToMXapNEmFSjaSrWoVKwXRstGrPWGAw53TaFKhUGQe4Pcfv1jrxlP\nDmfGfQ5n7bVn3u8n2cne61zeB2ad/ez1rluqCklSmx41dABJ0nAsAUlqmCUgSQ2zBCSpYZaAJDXM\nEpCkhvVWAknWJ7kwyRVJLk9yfLd8PskNSTZ1j6P7yiBJ2rH0dZ5Akv2A/arq0iR7AZcAxwKvBu6q\nqvf2MrAkaWJr+vrFVXUTcFP3/O4kVwEHdF9OX+NKkiY3lX0CSTYAhwFf7Ra9JcllSc5IsnYaGSRJ\nD9d7CXRTQecCJ1TV3cAHgIOAQ4EbgdP7ziBJWlpv+wQAkuwOfBr4bFW9b4mvbwDOr6pDFi33gkaS\ntAJVtazp9j6PDgpwBnDlwgJIsv+Cb/stYPNSP19VM/c49dRTB89gJjO1mMtMkz1Worcdw8CLgD8A\nvp1kU7fsFOA1SQ4FCrgOeFOPGSRJO9Dn0UFfZuktjc/2NaYkaXk8Y3gZ5ubmho7wMGaajJkmN4u5\nzNSfXncMr1SSmsVckjTLklCzsmNYkjT7LAFJapglIEkNswQkqWGWgCQ1zBKQpIZZApLUMEtAkhpm\nCUhSwywBSWqYJSBJDbMEJKlhloAkNcwSkKSGWQKS1DBLQJIaZglIUsMsAUlqmCUgSQ2zBCSpYWuG\nDiAtlizrPtmrpqoGGVcakiWgGTXtN+RhikcamtNBktQwS0CSGmYJSFLDLAFJapglIEkNswQkqWGW\ngCQ1zBKQpIZZApLUMEtAkhpmCUhSwywBSWpYbyWQZH2SC5NckeTyJMd3y/dNsjHJtUkuSLK2rwyS\npB1LX5fPTbIfsF9VXZpkL+AS4FjgdcAtVfXuJO8A9qmqkxb9bHlZ33aNLyU9/auIus5pZ5eEqlrW\nJXF72xKoqpuq6tLu+d3AVcABwDHA2d23nc24GCRJA5jKPoEkG4DDgK8B66pqS/elLcC6aWSQJD1c\n7zeV6aaCPg6cUFV3LbxrVFVVkiW3wefn57c9n5ubY25urt+gWtJQd/mS9NONRiNGo9Ej+h297RMA\nSLI78Gngs1X1vm7Z1cBcVd2UZH/gwqp6xqKfc5/AjBhqft59AtLyzdQ+gYzfPc4ArtxaAJ1PAcd1\nz48DzusrgyRpx/o8OujFwJeAb/OTj3UnA18HPgY8FbgeeHVV3b7oZ90SmBFuCUg7j5VsCfQ6HbRS\nlsDssASkncdMTQdJkmafJSBJDbMEJKlhloAkNcwSkKSGWQKS1DBLQJIaZglIUsMsAUlqmCUgSQ2z\nBCSpYZaAJDXMEpCkhvV+ZzFpZzHtu6h51VLNAktA2maab8retlOzwekgSWqYJSBJDbMEJKlhloAk\nNcwSkKSGWQKS1DBLQJIaZglIUsMsAUlqmCUgSQ2zBCSpYZaAJDXMEpCkhlkCktQwS0CSGmYJSFLD\nLAFJapglIEkNswQkqWGWgCQ1rNcSSHJmki1JNi9YNp/khiSbusfRfWaQJG1f31sCZwGL3+QLeG9V\nHdY9PtdzBknSdvRaAlV1MXDbEl9Kn+NKkiYz1D6BtyS5LMkZSdYOlEGSmjdECXwAOAg4FLgROH2A\nDJIkYM20B6yqm7c+T/JB4Pylvm9+fn7b87m5Oebm5vqOJkk7ldFoxGg0ekS/I1W1Omm2N0CyATi/\nqg7pXu9fVTd2z08EDq+q31v0M9V3Lk0mCeN9+VMdtYExg+u4VlsSqmpZ+1x73RJIcg5wFPDEJN8H\nTgXmkhzK+C/uOuBNfWaQJG1f71sCK+GWwOxwS6C/8VzHtdpWsiXgGcOS1DBLQJIa9lNLIMl+3fH8\nn+tePyvJG/qPJknq2yRbAv8IXAA8pXv9HeDEvgJJkqZnkhJ4YlV9FHgQoKruBx7oNZUkaSomKYG7\nkzxh64skRwJ39BdJkjQtk5wn8DbGZ/X+XJKvAE8CfqfXVJKkqZjoPIEka4CfZ7zlcHU3JdRfKM8T\nmBmeJ9DfeK7jWm29nCeQ5LHAycBbq2ozsCHJK1aYUZI0QybZJ3AWcB/wwu71D4C/7C2RJGlqJimB\ng6vqNMZFQFXd028kSdK0TFICP06yx9YXSQ4GftxfJEnStExydNA88DngwCQfAV4EvLbHTJKkKdlh\nCSR5FLAP8CrgyG7xCVX1w76DSZL691MPEU1ySVX9wpTybB3TQ0RnhIeI9jee67hWW1+Xkt6Y5E+T\nrE+y79bHCjNKkmbIJFsC17PER6SqOqinTG4JzBC3BPobz3Vcq20lWwLeWUw7ZAn0N57ruFZbL/cY\nTvIqHv7XcQewuapuXs5gkqTZMskhoq8Hfgm4kPHHpaOAbwEHJfmLqvpQj/kkST2apAR2B55ZVVsA\nkqwDPgz8IvAlwBKQpJ3UJEcHrd9aAJ2bu2W30l1KQpK0c5pkS+DCJP8OfIzxdNCrgFF3ddHb+wwn\nSerXJIeIPgr4bcaXiwD4T+DjfR6+49FBs8Ojg/obz3Vcq62Xo4Oq6qEk3wTuqKqNSfYE9gLuWmFO\nSdKMmOSmMn8E/Cvw992iA4Hz+gwlSZqOSXYMvxl4MXAnQFVdCzy5z1CSpOmY6H4CVbXt/gHd/Yad\nzJSkXcAkJXBRkncCeyZ5GeOpofP7jSVJmoZJjg7aDXgD8PJu0eeBD3p0UBs8Oqi/8VzHtdp6u4Bc\nkicDTOtaQZbA7LAE+hvPdVyrbVXvJ5Cx+SS3ANcA1yS5JcmpGb8zSJJ2cjvaJ3Ai4xPEDq+qfapq\nH+CIbtmJ0wgn7cqSTP0hLbbd6aAklwIvW3w/4SRPAjZW1aG9hXI6aGY4HbSrjDce07+rXdtq315y\nzVI3lO+WTXLNIUnSjNtRCdy/wq9JknYSO5oOehD4v+383B5VNcldyc4EfgO4uaoO6ZbtC3wU+Fng\neuDVVXX7op9zOmhGOB20q4w3HtO/q13bqk4HVdVuVbX3dh6TTgedBRy9aNlJjPcpPB34YvdakjSA\nSc4YXrGquhi4bdHiY4Czu+dnA8f2mUGStH29lsB2rFtwp7ItwLoBMkiSGKYEtukm/p2klKSBDHGo\n55Yk+1XVTUn2Z3zP4oeZn5/f9nxubo65ubnppJthnuwjaaHRaMRoNHpEv2Oiawc9ogGSDcD5C44O\nejdwa1WdluQkYG1VnbToZzw6aAkeqbMrjenRQVp9vV1AbqWSnAMcBTyR8fz/u4B/Y3zT+qfiIaLL\nYgnsSmNaAlp9M1cCK2UJLM0S2JXGtAS0+lb7shGSpF2cJSBJDbMEJKlhloAkNcwSkKSGWQKS1DBL\nQJIaZglIUsMsAUlqmCUgSQ2zBCSpYZaAJDXMEpCkhg1xUxlJAxnixkReuXS2WQJSU4a4RLdmmdNB\nktQwS0CSGmYJSFLDLAFJapglIEkNswQkqWGWgCQ1zBKQpIZZApLUMEtAkhpmCUhSwywBSWqYJSBJ\nDbMEJKlhloAkNWynv5/Afffdx0UXXTT1cdeuXcvhhx8+9XElaTVlFu/6k6QmzXXLLbewbt3+7L33\nXL+hFnjggdt42tP2YNOmi6c2Jmy9K9QQNwVxzJ1/vOHGnMX3mF1VEqpqWXfy2em3BAAe/ejHc8cd\nG6c44sU88MApUxxPkvrhPgFJapglIEkNswQkqWGD7RNIcj1wJ/AgcH9VHTFUFklq1ZA7hguYq6of\nDZhBkpo29HTQsg5lkiStriFLoIAvJPlmkjcOmEOSmjXkdNCLqurGJE8CNia5uqq2nX01Pz+/7Rvn\n5uaYm5ubfkJJmmGj0YjRaPSIfsdMnDGc5FTg7qo6vXu9rDOG169/Bvfee0ufERe5mOc85xQ2b/aM\nYcfcWcYbbsxZeI9pxUrOGB5kOijJnkn27p4/Fng5sHmILJLUsqGmg9YBnxx/smUN8M9VdcFAWSSp\nWYOUQFVdBxw6xNiSpJ8Y+hBRSdKALAFJapglIEkNswQkqWGWgCQ1zBKQpIZZApLUMEtAkhpmCUhS\nwywBSWqYJSBJDbMEJKlhloAkNcwSkKSGWQKS1DBLQJIaZglIUsMsAUlqmCUgSQ2zBCSpYZaAJDXM\nEpCkhq0ZOsDO6vLLv0ySoWNIWmSov8uqGmTcR8oSeESm/Y9u6UiT8W9zUk4HSVLDLAFJapglIEkN\nswQkqWGWgCQ1zBKQpIZZApLUMM8TkNQrT6qcbZaApJ554tYsczpIkhpmCUhSwwYpgSRHJ7k6yXeS\nvGOIDJKkAUogyW7A3wFHA88CXpPkmdPOsTKjoQMsYTR0gCWMhg6whNHQAZYwGjrAdoyGDrCE0dAB\nljAaOsCqGGJL4Ajgu1V1fVXdD/wL8MoBcqzAaOgASxgNHWAJo6EDLGE0dIAljIYOsB2joQMsYTR0\ngCWMhg6wKoYogQOA7y94fUO3TJI0ZUMcIrrqx4vdf/+dPO5xv7nav/Zh7r33Gh7zmEt48MFbueee\n3oeTpN5l2nfDSXIkMF9VR3evTwYeqqrTFnzPznmLHkkaWFUt60SJIUpgDXAN8GvAD4CvA6+pqqum\nGkSSNP3poKp6IMmfAJ8HdgPOsAAkaRhT3xKQJM2OmTpjOMn6JBcmuSLJ5UmOHzrTVkl2S7IpyflD\nZwFIsjbJuUmuSnJlt69lcElO7v79Nif5SJKfGSDDmUm2JNm8YNm+STYmuTbJBUnWzkCm93T/fpcl\n+USSxw+dacHX3pbkoST7zkKmJG/p/l9dnuS07f38NHMlOSLJ17v3hW8kOXzKmZZ8v1zuuj5TJQDc\nD5xYVc8GjgTePEMnkp0AXMn0r4a1PX8LfKaqngk8Fxh8Si3JBuCNwPOr6hDG032/O0CUsxifjLjQ\nScDGqno68MXu9dCZLgCeXVXPA64FTp6BTCRZD7wM+J8p54ElMiX5VeAY4LlV9Rzgr2chF/Bu4M+r\n6jDgXd3radre++Wy1vWZKoGquqmqLu2e3834je0pw6aCJAcCvw58kBm4RGH3ifGXq+pMGO9nqao7\nBo4FcCfjFXPP7gCAPYH/nXaIqroYuG3R4mOAs7vnZwPHDp2pqjZW1UPdy68BBw6dqfNe4O3TzLLV\ndjL9MfBX3cmlVNUPZyTXjcDWrbe1THld38775QEsc12fqRJYqPtUeRjjP46h/Q3wZ8BDP+0bp+Qg\n4IdJzkryrST/kGTPoUNV1Y+A04HvMT7y6/aq+sKwqbZZV1VbuudbgHVDhlnC64HPDB0iySuBG6rq\n20NnWeBpwK8k+WqSUZIXDB2ocxJwepLvAe9h+lty2yx6v1zWuj6TJZBkL+Bc4ISu4YbM8grg5qra\nxAxsBXTWAM8H3l9VzwfuYfrTGw+T5GDgrcAGxltweyX5/UFDLaHGR0PMyrQeSd4J3FdVHxk4x57A\nKcCpCxcPFGehNcA+VXUk4w9jHxs4z1ZnAMdX1VOBE4EzhwjRvV9+nPH75V0LvzbJuj5zJZBkd8b/\nQf9UVecNnQd4IXBMkuuAc4CXJPnQwJluYPxp7Rvd63MZl8LQXgB8papuraoHgE8w/v83C7Yk2Q8g\nyf7AzQPnASDJaxlPNc5CWR7MuMAv69b3A4FLkjx50FTj9f0TAN06/1CSJwwbCYAjquqT3fNzGV8X\nbaoWvF9+eMH75bLW9ZkqgYzvQ3cGcGVVvW/oPABVdUpVra+qgxjv5PyPqvrDgTPdBHw/ydO7RS8F\nrhgw0lZXA0cm2aP7t3wp453ps+BTwHHd8+OAwT9gJDma8SfbV1bVvUPnqarNVbWuqg7q1vcbGO/k\nH7owzwNeAtCt84+uqluHjQTAd5Mc1T1/CeOd+1Ozg/fL5a3rVTUzD+DFjOfdLwU2dY+jh861IN9R\nwKeGztFleR7wDeAyxp+SHj90pi7X2xkX0mbGO6V2HyDDOYz3SdzH+GKFrwP2Bb7A+A/1AmDtwJle\nD3yH8RE4W9f19w+U6cdb/z8t+vp/A/sOnQnYHfhwt05dAszNyDr1AsZz8JcC/wUcNuVMS75fLndd\n92QxSWrYTE0HSZKmyxKQpIZZApLUMEtAkhpmCUhSwywBSWqYJSBJDbMEJKlh/w/WqiJyDhDESwAA\nAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot histogram of price\n", "df.price.plot(kind='hist')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWwAAAEACAYAAACXqUyYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXt4XVWd97+/3NOkNE3vpXdugkIBFRDFRgHldXzBK4qK\niPMqvs+M4mV0AJ9XU8fncXScGZ0ZdGbQVlBg1IIgeAORFKTKvRRKa1taoGmTtE3StEmTJulZ7x+/\nLM4+6+z77Zx9zu/zPHlyrvvss8/e3/3d3/Vba5FSCoIgCEL5U1PqFRAEQRD8IYItCIKQEUSwBUEQ\nMoIItiAIQkYQwRYEQcgIItiCIAgZwVWwiWgxET1IRJuJ6Dki+szU4+1EdD8RbSOi+4ioLZ3VFQRB\nqF7IrQ6biOYDmK+U2khErQCeBPAuAFcDOKCU+hYR/T2AmUqp61JZY0EQhCrF1WErpXqVUhunbg8D\n2ALgeACXArh56mU3g0VcEARBSBDfGTYRLQNwFoBHAcxTSvVNPdUHYF7sayYIgiAU4Euwp+KQOwBc\nq5Q6bH1OcaYi/dsFQRASps7rBURUDxbrHyul7pp6uI+I5iuleoloAYB9Nu8TERcEQQiBUorsHveq\nEiEAPwTwvFLqO5anfgngqqnbVwG4y3zv1IfK39TfV7/61ZKvQzn9yfaQbSHbw/7PDS+H/UYAHwGw\niYiennrsegD/COBnRPTXAF4EcLnHcgRBEISIuAq2UuqPcHbhF8W/OoIgCIIT0tMxJTo6Okq9CmWF\nbI88si0Kke3hjGvHmUgLJlJJLVsQBKFSISKoMI2OgiAIQvkggi0IgpARRLAFQRAyggi2IAhCRhDB\nFgRByAgi2IIgCBlBBFsQBCEjiGALgiBkBBFsQRCEjCCCLQiCkBFEsAVBEDKCCLYgCEJGEMEWBEHI\nCCLYgiAIGcFTsIloDRH1EdGzlsfOIaLHiOhpInqciF6f7GoKgiAIfhz2WgCXGI99C8D/U0qdBeAr\nU/cFQRCEBPEUbKXUwwAGjYd7AMyYut0GYE/M65U4t94KfO97pV4LQRAE/3hNwuvEdQD+SETfBov+\nG+JbpXR44QWgt7fUayEIguCfsIL9QwCfUUr9gojeD2ANgIvNF3V2dr5yu6Ojo6zmapucBA4eLPVa\nCIJQ7XR1daGrq8vXa33N6UhEywDco5Q6fer+IaXUcVO3CcBBpdQM4z1lPafj9dcDmzYBv/pVqddE\nEAQhTxJzOu4golVTt98KYFvI5ZSMiQlgaKjUayEIguAfz0iEiG4HsArAbCLaDa4K+SSAG4moEcDo\n1P1MIYItCELW8BRspdQVDk+dG/O6pIoItiAIWaNqezqKYAuCkDWqWrAPHwZyuVKviSAIgj+qWrCV\nYtEWBEHIAlUt2IDEIoIgZIeqFezJSf4vgi0IQlaoWsHWDlt6OwqCkBWqWrBra8VhC4KQHapasNvb\nRbAFQcgOVS3Ys2eLYAuCkB1EsEWwBUHICCLYItiCIGQEEWwRbEEQMoIItgi2kDL33APcf3+p10LI\nIlUr2JOTIthCaejq4j9BCErYKcIyz8QEMGeOdJwR0md0FBgfL/VaCFmkqgVbHLZQCkZHZZRIIRyu\nkQgRrSGiPiJ61nj800S0hYieI6JvJruKyTAxAcyaJYItpM/YmOx3Qji8Muy1AC6xPkBEbwFwKYAz\nlFKvAfDthNYtUcRhC6VidBQ4dKjUayFkEVfBVko9DGDQePj/AviGUmpi6jX7E1q3RNFd00dGgGPH\nSr02QjUxNiaCLYQjTJXISQDeTER/JqIuInpd3CuVBhMTQGMj0NIikxgI6TI6Kld2QjjCNDrWAZip\nlDqPiF4P4GcAVti9sLOz85XbHR0d6OjoCPFxyTAxAdTXAzNm8MHT1lbqNRKqBYlEBCtdXV3o8lnn\nSUop9xcQLQNwj1Lq9Kn7vwHwj0qp9VP3dwA4VynVb7xPeS27VORyPLRqLgecfjpw223AGWeUeq2E\nauGMM4Bt2zgaEQQTIoJSiuyeCxOJ3AXgrVMLPhlAgynW5c7kJLtrInbWcnkqpMnoKHD0KP8JQhC8\nyvpuB7ABwMlEtJuIrgawBsCKqVK/2wF8NPnVjBcdhwAciUjnGSFNRkf5v7SdCEFxzbCVUlc4PHVl\nAuuSGqZgi8MW0mRsDGht5f1u9uxSr42QJapyLBERbKGUjI4C8+ZJw6MQHBFsEWwhRZRihz13rgi2\nEJyqFey6qTCo0gW7t7fUayBYmZgAampkWAQhHFUr2NXgsHM54IQTpHysnBgdBZqbgeOOE4ctBEcE\nu4IFe2gIOHJEhKGcGBsDmppEsIVwiGBXsGD3T1XHizCUD9phV/J+JyRHVY6HbRXsSu44MzDA/6Xe\nt3yQSKR82buXO9O1tvIYQzVlaGfLcJWSR/d0BCq744wIdvlhjUQq1ShkkaEhYMkS4KyzgAULgIaG\n8pzGreoddiVfmopglx/WSEQcdvlw6BAwfz7Q3c33P/hBoKentOtkR1U67GoTbBGG8mFsTCKRckSf\nSDVNTeVZXVX1gj19OldSVOIkBrrRURy2N7feCqxfn/znjI6yGFSyUcgipmA3N+fHfCknqlawdceZ\nmhpuZKhEtzMwwN9TBNube+8FHn00+c+RRsfyZHQUmDYtf18Eu4ywOmwg+27nppuAPXuKHx8YABYt\nEsH2Q29vOttJ6rDLkyNHJBIpOYcOAb/6VfHjlSbYN95o7w4HBoBly0QY/NDTk45gV3IddpanPpNI\npAx44gngH/6h+PFKE+zeXqCvr/jxgQFg6VJx2H5Iy2GbkUiZTsoUiptvBr785VKvRThEsMuA4WH7\nWT1Mwc5y55nJSWDfPnvB7u9nhy2C7Y52hmlGIo2N3EmjkmadGRrKVyZlDTPDzmQkQkRriKhvanYZ\n87kvEFGOiNqTW71ojIwA4+PFj1s7zgDZdtj79rFL27ev+Dlx2P7QIxqm6bCBbO93dmR5cmEzw86q\nw14L4BLzQSJaDOBiAC8lsVJx4ddhZ7m3oy7uNx12LsffaenS7B5EadHby243LYethaHSGh6PHMmu\nOaiISEQp9TCAQZun/gXAlxJZoxgJIthZdTo9PVxLbgr2oUM8HsLMmdk9iNKit5e7JaflsJua+Hal\ndU/PssOu2I4zRHQZgG6l1KYE1idWhoftI5FKE+yVK4sFe2AAaG9nMRfBdqenBzjppNJEIlkVODuy\n7rCzUIcdaCwRIpoG4AZwHPLKw06v7+zsfOV2R0cHOjo6gq1dREZG/DvsXbvSW6846ekBzjwTuOWW\nwsf7+0Ww/dLby4K9cWPyn6UbHYHKi0RGR7O7rx05wvNsatIU7K6uLnT5HGkq6OBPJwBYBuAZIgKA\nRQCeJKJzlFJFzV5WwS4Fbg67zvLNs+6wX/1qFgKrGAwM8DRUlSYKSdDbC5x2WvoOOyuRyFNP8bqe\neKL767I8WUYpIxHTzK5evdrxtYEiEaXUs0qpeUqp5Uqp5QC6AZxtJ9blQLVk2AsW8KSu1lhERyJN\nTVwVMzFRunUsd3p7ufxxctL+BB8nWYxEvvpV4Be/8H7dkSMsclnc1yqi0ZGIbgewAcDJRLSbiK42\nXlLWZf8jI1wtMTlZ+HglCva8efaCTSSxiBd6G6axnbIWiRw7Bjz0kL96cS1wWdzXKiLDVkpd4fH8\ninhXJ16Gh/n/+HhhBFJJHWe8BBtgITp0KH9fKKS3l8dC1oI9a1Zyn2VGIuXe0eSZZ3jf8SPYR47w\n/8OHs7evyVgiZYAWbHNns+s4k8U6bKVYpLVgWzvP9Pfnhee447LpetIgl+NtOG9eeg47S5FIVxeP\naOnXYdfUZHNfq4hIJOtYHbYVu0ik3A8cOwYG+DKuqcnbYWfxIEqDgQEeXrepKZ3tZNZhl/t+19UF\nnHeef4c9Z075fyc7RLDLgJER/m/ubKZgt7ZmcxIDHYcAzo2OgAi2GzoOAdIT7Kx0TT92DHj4YeDi\ni/077HnzsrmvHTliP5ZIuQ3OVdGCPTzMG97LYetJDLK2o1kF20+GXe7kcsB73pPuQWIV7DSioyw1\nOj7zDO9fS5b4y3N1LXM5fycnTIddV8e6UG4VLxUv2O3t3g4bKP+Dxw6/gp2VDPvoUS4fS7M9obc3\nvw3TdthR6rC3b89PAZcUXV1ARwefYLwctlJ5wc7CvmZiCjZQnrFIxQq2UhyJOAl2nVEfU+6Xp3aY\ngu3U6JiVSET/TnZDxSZFT09pI5GwJuHrXwd+/vP41ssOLdiNjd6CPTEB1NZmd+waJ8Eut0qRihXs\nsTGgoYFzKa9IBKgsh53LAYODfPAA2RNsPdxpGqSZYSvF37Gxke9H2edGRpJdV51fv/nN/gRbl8Vl\nJX4zMTNsgK8sxGGnxPAwj1Znt7PZCXY5OOwNG4D3vtf/662C3d7O6z8xwQfytGn575iVk5E+sabp\nsNMU7LEx3h9rpo66KLPOJN0NfNMm3i7z5/sTbN3xJCvxm4lEIiVmeJgbEhsbs+OwX3gB2LnT/+ut\ngl1byxHI/v2F+TUgDtsN6zZMejuZolBfz1eBusNJEJJ22DoOASrfYU9O8lWpqQki2CkyMsKC3dDg\n3XEGKA+H3dcXrOebVWyAfCySVcGuBoetK0Q0YY1C0kOZWgXbT6+/LDtsfSIlY9zRcuztWLGCHTQS\nKQeHvW9fsJZ/O8HWy7B2r86KYFd6hm132R22UiTpSOSRR4ALLuDbfh32tGnZ2des2OXXQHk67KDD\nq2YGHYk0NPiLRMrFYesxvHXDlBPDw3wZN316/jHtsBsaCh12OZyM/JC2wx4by5d+Auk4bFOww1aK\nJBmJjI3xsrUZqPRIxO5ECpSnYFe0w9YZdlYcthYqP7GIdtfWyzjd2zGrkcjRo3wZmpbD1mOI6EbA\nOLeTUsUHu7VbuqYcI5G+Pt6X9L5V6Y2OToItkUiKWDPsrDhsXUftJxYx4xCgMjLsJUvSc9jWOASI\ndzs98ADwwQ8WPhZ3JJKkYFtnXxGHXT5UrGBnMcPu6+OB9IM4bCtZF+yjR4HFi/nElcsl/3lJCvbe\nvfkZ7TV2jY5hIhHdqzCp/dVOsCul0fHd7wa6uwsfy1KG7SnYRLSGiPqI6FnLY/9ERFuI6BkiupOI\nZiS7msHxikTKradjLsclea96VTSHbdfoWA4nIz+Mj7NotrZyx5+kMbdhnII9OFj8HeJy2HpQoiQd\ntvVE5qdretqNjlu2AHfcEew9SgH33Vcs2JUWiawFcInx2H0AXq2UWglgG4Dr416xqARtdCy1qB08\nyFcECxbEH4m0tPBOmYZrjcLRo/x7zZuXTo5t57CHh+MZfGpwsPhKyUmwg+53R47wb5qUMPb2Bo9E\n9HdraWGRM2d5ipvf/x649dZg79m3z/7KpKIiEaXUwwAGjcfuV0rpw/9R8GS8ZcXISLo9HQcHgW9+\nM/z7dUPPrFnhIxGnRseaGnY/enzwckVXx8yfn06ObQp2XR3vF3EcpIODfBK2niTjikSOHMkPO+Bn\n2NOgmJGIvhp1E2HtsInYKCW9r/X0BN9uu3bx/4oWbB98HMCvY1hOrLg5bLuOM1Ed9rZtwH/+Z/j3\n79vHB8msWeEd9ty5wIED/GdO0ZSFHHt8nAU7LYdttw2DbKdcDrjtNvvnBgf5ees+FVckMjKSjx+S\nuCo0BRvwdtnWKbbs1uvHPwaeey6+dezpCb7dnATbKcPOaiTiCBF9GcC4Usphty0dQcv6ojrs4eFo\n8/Nph93eHl6w6+tZAHbtKp6XsNSRjx90JJKUw96zB1izJi/IpsMGggn2xo3ARz9qHzXp/NqaY8dV\nh60FJqkGvjCCbZ3E1m69fvIT7j0ZF3EKdpYcduiOM0T0MQDvAHCh02s6Oztfud3R0YEO3dc1BXRZ\n36FD/jLspiYeocxPpxU7hof5syYnixs0/aAPkiiRCMDLGBjIXzJrsuSw58xJxmH/7nfAl74E/N3f\nAR/+MPDii9EEe/163meGhoq3t1Wwly/n23HVYesM+9ix5ATb3C5ebtPLYff08AkzLsJGInPmFAt9\nqQW7q6sLXT7PZqEEm4guAfBFAKuUUo4/o1Ww0yZoWR9R3u3MmRP88/SBc/AgMHt28PfrSMSPwz56\n1Hl277lzgd272alayYJgWx32X/4S//IPHAA+9jHg2muB//ovYNGiaJHI+vX8v7/fXrDnzCk8+cYd\niSiVzFWT2egIBHPYdtuwt5dLHeOipyf4/rxrF7ByZTCHnUYkYprZ1atXO77WT1nf7QA2ADiFiHYT\n0ccB/DuAVgD3E9HTRPS9qCsdN0EjESBabKAbWcLGIkEaHbduBY4/Pt9Dz4oWfZMsCHbSGXZ/P59M\nFy/mCQCefLLY8frdTrkcjxe9aJH9CXZwEDjhhOJIJK5Gx6QikbGxwkZNjZ8M2ykSmZjgktW4BHti\ngrfrxETx1bMbO3cCZ54ZLMPOXCSilLrC5uE1CaxLrDg1OirlHFtEybGjCnaQRsdvfhP45Cftn3MS\n7Kxk2C0tyWXY/f3AiSe6v8avYD/3HG/nE090Fuy3vMXbYeshcYOgI5GGhvgFe9++wm7pGj8O2ykS\n0b9lXJFIXx9fvRw9yp/j54p2cpI///TTgV/9qnjd7ZZRjhl2xfZ0dCrrm5zksaPNHRKIJmr6wInq\nsHUk4lQLvHkzd3v+27+1f16LvknWHHYSgn3ggP22seJ3O61fD6xaZX+CHRvjfHnRokKHbSfYc+fy\nPhNkstckO6nYNTgC0Rx2by/34I3LYev2myAGq7ubRX7OnPKLRIJQsYLt5LCd4hCgtA5bHyjNzXxC\ncRrU/mtfA77whcJR+qzMn59dwdYZ9ty57Drj7uijIxE3/G6nhx5yFmw9PVt7u3ckUlfH3zdIBJRk\nWZ9dgyMQrayvpwc47TR+/8hI9HUMI9i7dgErVtibMreejuKwDX7yk2R6RTll2G6CHTXDbmiIFonM\nncu3nRoen32Wnd3f/I3zci6/3L4DTxYG5dEVOvX1fDC6RUOXXRZ81nCzy74dfgRbKX+CPXOmdyQC\nAAsXBosLdCSSRIZt1+AIeFeJuJX1aYFduDAel62XF+R43bWLq3Xs3hN0LJEzzwSef7748aNHgbvu\nCparm5jd5k1KKtjd3cCVV3qvZFByufyPYBeJJOGwDx/mkebCCPbICF9Ca9fs1PDY2cllaS0tzsua\nPp0vP03KfVAegHd0Xd3i1vA4Ngbcc0/w2CSuSGTrVt63lizxFmyvOmyAG5CDCnY5RiJWh21GIvPn\nxy/YQR22k2AHiUSU4vaL620G4/jmN4FrruF2jf/4D2d3nssBV19tH3s+84z79yipYN95J/8/eDDe\n5ep619raYJGI3Y/Z2ws8+KD3Zw4PhxdsfZDoXN1OBDZuBP70J+BTnwq+fCA7kYiugXdreNy+3X68\naTeUYvG0a5C14mc7rV/Ps4kDzoLd3s5/psM2IxGABTuIkFkjkXIRbLOsz4xEFiwI/j2dSEuw7SKR\nkRHWj02b+CpLs3Mn8G//Bjz+OLBuHXD//VwlZGdGh4eBH/3I/rfzaoAuqWCvW8dfPu5R8nQcAgSL\nROx2gF//mjtaeBHFYVvjEMA+Ern1Vq4Msbt080MWBFs3OgLuDnvrVv4fRLCHhnjbmfXpJn6iI93g\nCLBgHzhQ+LyTw04iEkkiw4670dEaicRRKRIlEtHHgdXdBuk4c/Ag/65f/zrwxS/ml3PttawTS5YA\n55wD3H03cMop9v0JdHuXnTiXrWD39nIm+5a3xO+wrYId1WHv2cOXKV6NJcPDwNKl0Ry2xi4S2bED\nePWrgy9bk5WyPi2obg57yxb+H0Sw/cQhgPeJTalCwZ49O1qjIxDceSYdiYRpdHQr6yunSKSujr+L\ntVE/SIY9NAS0tQFXXMFasm4d8Mtf8lXf5z9f+Fo9+qOJfkxPWGKlbAX7F78A3vEOdpZxC7Yu6QOi\nO+zubs6Xn3jC/TOjRCJ+HPYLL/AlVliq3WH7aXAEvLfTjh3cYWnFCr5fikbHUkQibo2OOp7S383J\nYZcqEhkd5d9h4cL8+lnfF2Q87IMHWbBraoB/+ifOsq+9FrjxRvvexW6CnSmHfccdwPvex18+yUik\noSFalciePSyUf/qT+2dGiUS8HLZSnJFVumD7ddhbt/LVTFDB9tPBwms76eoQt/YGLdjTp7N70zXW\nTsJQTg7bqUrEzWGPjfHvZjc3plKFDjtqJJLLscGZP99/t/6XXuLerbW1fN88zp1+l8ZG/u2OHcs/\nNjTEJwoAuPBC4KSTgHPP5dsmra32v4+bw7Z7zEpJBPvAAQ7n3/52/vJJRiKNjdHqsLu7+cTiJdhx\nRyJWEejr47O93lHCkAXB9uOwcznOBc86qzSRyGOPAeedl78/bRqvk3VdtGDX1LAh0fu3UyRSLhm2\nrpM2u6UD7oJtbXAECiORwUEWw+bmeCKRAwf4ezc0+O/Wv3NnfgAuwL9gExW7bO2wNevWATffbP+5\nTuOC6/0rMw777ruBt72Nf+SkHbZdJOI0mp6Tw37/+4ENG5x7HyrFn7loUfGg9X7wikSixiFA5WTY\n3d18oM6bV5pI5Mkngde+Nn+fqPgEqwUbKGx4dBKGtjYuN/V7Qk3KYe/bxz0B7caocRNsa0kfUBiJ\nWEeV1IIdZUYf6/L8RiK604zGFHpz/a2Ygm112EC+N7UdFZNhr1vHrhVIxmFbM+woPR3HxviHPess\n/uFeeMH+fWNjfBLQUyQFPYi8IpE4BDvO6a+Swo/D3roVOPXU4OM8xBGJjI9zh4kzzyx83E2wrQ2P\nTnXYRMFikaQybKf8GgjmsFtaWASPHePfUAusFrco83VaBdtvJKIbHDVW82Lm7ybmfmY6bDfcIpGG\nhow47MFB4JFHuMERSM9ha6EKkmHv3cs7R00N8IY3OMciw8P5Ti9m7a0f0nDYUaa/smZ4SWJ12Hpo\nUvOzt2zhiYqDCrbfSETPf2n3nZ97jp2aWVHg5bD193DrtBUkFtGRyLRpvM3i6insVCECuE/Ea1ZZ\n1NTk44CensJlRo1FTIft56rRTbDHx/nYcLrqthNsv9Gkm8NevrzYYR854v1bpi7Yf/wji58WOGvG\nFxdWwdYDPekD0O2g0T+kFvc9e9j5AO6Cffhw/vPCCLZXhh2HYAPhHNnBg9xzMg1nbu04U1fHYmc6\njq1bwwm230ikpobF0O5Ae/JJ4HWvK37crMW2i0S0i7MbdAwI5rCt8yfG6bKdGhwB/l2cqkTsHKpe\nL3OijaiVIn4ikYce4ghT4ybYbu4asI9E4nDYy5cX79v793uPxZ+6YL/0UuEQl1Gn5rJDT16gsV7O\nuTns+np2eLpGs7vbn2BbTxBBBXtignceq5joy2idhccl2GFy7L4+3g5erddxYI1EAPscO4pg+51Y\nwkkEn3iiML/WmCdY64w/en9wikM0QRy2jkTc1jUMYSMRuzpm3fBojUSA6JUifgT7+98HOjqAm27i\n+26C7ZZfA9Ectluj44oVxcdUWQr2yy9ziY0mCYetpwfT+BVsoHAn2LOHGxIBzrG3bXM+Y4aNRPbv\n5wPe2tBTX88HgN6pSumwtXPcvj3653thjUQA+xw7bIbtNxIBnLeT2eCosXae0UOragGzOmy7ChGN\nX+eZyxWKf7kItt1MOtphe0UiQY5/q2DryMFs5O/rA777XeDb3wY+8Qm+qrb+9qbDdus9bO5nQRy2\nVySyf3/hlWssgk1Ea4ioj4ietTzWTkT3E9E2IrqPiHx+BZ6+asmS/P2kHLZVsK0Nj34EW/+Y1kik\nsZEbmx5/vPg9USIRp4NENzwePszfx27+xqCEObi1EO3Y4f66XI6db5SJiE2HvXhxvpMMwAf28DD/\nJkFntPYbiQD228mpwREodNg6DtHRh75a8rr0thsA6rHHgFtuKXxMlwbqE3yc1T9xNToCeYftFYmM\nj7MZ8RuTWJdXW8ufa4piby9wwQXAo4+y3px4YmEUFSQSMQeACuqwnQze7Nm8bKv2xeWw1wK4xHjs\nOgD3K6VOBvDA1H1fmA5bV4nEmZGagh3EYVtbnru78w4bcI5FokQiZoOjRjc8anftlH0GIUmHvWUL\n10dHmYvRdNgf/Shf3moH9Ze/8PgMRMEctlLRBdupwRGwF2yNbnR0qsHW2EUF//M/3O3ZijUOcVrX\nsLg1OoZ12F6RyMMP8/bxa9rME4Cd4dPfo60NuPde4De/KV43/R4/GXZYh+0UiWi9mDOnMBbZv99e\nC6x4CrZS6mEAZiHOpQB0ufjNAN7ltRyN6bC1W4hzZgdrWR8Qj8MG3AU7bCTi5rCtgh0HYdxYfz8f\nIF6CrRt5du4Mt27HjrEwW1vrV61icdIHnK4QAZwFe2gIePe7Cx8bGeF9zO/AWXYi6BSHAO6CHcRh\nmy7zkUeK9yVdIeK2rmFxa3R0qxJxcth+IpF77+X/fiY2UKpYsM19enyc7+tRGevqik9CaWXYTr+N\nFmw9UYdG18G7ETbDnqeU0s1BfQAcfuZCJid5p9B9+jVx59heDtuphAcodtimYP/5z8VXA0lGInEK\ndpia9/5+7tnnFYls2MA7oFOtuhd6LGzrlQQRD6jzL//C93V+DTgLdl8fDyJvFYAg7hpwFmy7ChHA\nv8N2E4YFC1iM9NXE6Cjw1FP2gp2kww5TJeLU6NjXx++xbg/riUkpHtd85kznGZasDA2x2bKesEyH\nra9Y7Tr/aIJm2FGqRJwaHbVgmw478UZHpZQC4CvQ6OnhlTQdbtydZ+wy7CCNjocO8YFjnlz0bbN1\nN41IJA6WLgVefDHYew4cYMHWY1A7sWEDj2AWRbDteoxdfjkL9caN+QoRwFmw9YFvvSIIUiECxOuw\n/TY6NjWxkOgI6okneJ29BDuuDHt8nL+z03jhXhm2XSSyfTu7W+tJeP583uePHeNG/LEx4PWv9+ew\nTXcNFAu221WC9T1Byvr0fjYxkZ8o2g9asM3jxhqJWB22H8H2nDXdgT4imq+U6iWiBQBsi746Oztf\nud3R0YH6+o6C/FoTd+cZO4ftNxLRDnvfPl4vU0TmzOGD07pTWCMRc4Q2L/r67IdNtTrs97zH//Lc\nOOEE4Ke5cNGnAAAgAElEQVQ/DfYePdN4U1N+ZneTAwf4e1x6KfCVr4RbNzO/1jQ08ITD3/mOP8HW\nB/5f/pJvIAxSIQIUC7ZbgyMQTyQC5PPduXM5Dvmrv+Ic2/x+SThst27pQLiyvg0bigW2vp63z759\nHIe8850sVH4ctp1gmycstxze7j1+Gh31fqa7pfttT6qv56t58+pK64V22F1dXejq6sLTTxe3WZiE\nddi/BHDV1O2rANxl96LOzs5X/jo6Oorya03cDtvMsK07m1vHGb0uhw4VNzhqZs8uHrA+SiRiDkyj\nScJhn3BCcAes44QTT3TOsf/0Jx6x7KSTwjtsa6cZk2uu4fFnXnwxX8PvR7DN7+AXUwSfe45/I6dL\n57Y23mcmJ8NHIkBhXLBhAw+OdvRooVAmlWF7CV2YRsdt2+yXqb+nFuxp09J12EEzbB2JBOmWrrGL\nRUyH3dHRgc7OTjQ2duK66zpdl+enrO92ABsAnEJEu4noagD/COBiItoG4K1T9z0xK0Q0STvsII2O\n2mGbDY4auxlGokQi27YBJ59s/zm9vbxjL13qf3luaMEOUpGj3elJJzkL9oYNwPnn8/Y6eNCfWzJx\nikQA3qYf+hCf7PVrvCIRazlgVMF2i0MALi+bMYPF2hTs5ub89GRukQiQL+1TirfpG99YPAlCUpGI\nl9CFaXQ0K0Q0CxcCmzfzdn3rW/Njj3jhR7D9OGxdcqjHEXHLsK2RiDnwkx/sarGdMmyneNSKnyqR\nK5RSC5VSDUqpxUqptUqpAaXURUqpk5VSb1NK+fLHaTnsKGV9YR22GYn4EUV9qWy3g82axY1Oxx/v\nvr5BaGtzHnTGCZ3/nniic8OjFuyaGu7GHqZSxCkS0Vx/fWHc4uawV6wodthRMmwvwQbynWdMwSbi\n+3v3+o9Etm/n/XfhwmLBTjIScRMLr0ZHu67pgLNgr1nD82JOmxbNYZsnLD8Ou6EhH1UEiUTCOmzr\n73PsWP5qy5phj42xaTnuOPflpdrTMQmHrRRP5qsF0rpBNHE6bLspoawniOZmdlxWx9DXZ+9Ot29n\n52qXibW3xxuHaILEIkrxyae93dlhT0ywoJ17bvDlW3Fz2ACfPD/ykfx9N8HWvVL1PhElw965k8sK\nzznH/T06xzYFG8gLth+HvXcv59fnn8+PmVdsSUUiXg1eQTvOaOFxikTWrwf+9//m+0Ectrk8O4ft\nJdh6/Q4dChaJxOGwdVxbU1PosPfvZ23xysdTFezdu50FO6zDfuop4L3v5QMUyF8yWhtPwjpsJ8G2\ni0S0owCKXdGNNwJf/nLxsrZvt49DgLzAJCHYXiV6mqEh3mEbGpwFe+NGXqY+QMMKtpfDNnET7OOP\n531A58FhI5G77uIKmc9/Pn9CcsJNsNvbgzlsHYcAxY3YSZX1eXXaCNPoCDg7bIAbVQF+rx/B3r69\n+ArdLsP2ikSAvGAHqRKJw2FbtcLqsP1UiAAlcNhOkUhYh/2zn/HlzQMP8H0zDgGClfVZHXaYRkeg\n2BU99hhP5GvilF/rZQClddjWKEFHImbUo+MQzYoVyThsk6Ymfo85joR2oNYZq8NEIg89BHzuc9ww\ndu213s7Hj8P22+ho3abmvmRGInFl2F6dNsKU9QH2gr10KXD22fnjq6XFOxLZupVNlD6RWT/HrBIJ\n4rCDjCUSpNOMxmx0tOqT1pJcrgwF+8gRXlm7lQrrsJViwf7854Hf/54fsxPsIGV92mGHbXQECg8y\npViwX3yx2EW4CXZbG4tEKQXbGiW0tbFImqPnmYKdlsMmss9V9SWnVbCDRiKnngr89V9z1OMVhWj8\nOGw/kciOHXwlevrp+femFYl4Oezxcfu2GTeHbed2L7oof7wC/hz22rU8VIHZ6S2Kwx4aCjaWSJBO\nMxozErGau4YGvn3woL8GRyBFwdZxiJ1TCeuwH3uMd6TPfhbo6uL82izpA8L1dHRrdLTLsM1IRB9k\nO3bwc6edxuVhVtwEu7aWD/xycdiAfSxiJ9hhGh2DOmzAPhbRDtR02EEEe/Zs4N//3bkTiR1eDru/\n39thz5nD++k55+T3UTvBTioScXN4NTW8TnpCYStOZX1E9iJUU1O4jbwc9sQED4J19dXFz1m1Y2yM\n18VuTkoTvxm2GYmEcdhmJGI1dzrHLjuH7ZRfA+Ed9k9/CnzgA3zZtWAB8PTTzpFIEIfd18c7m12L\nbdBI5LHH+ABcubIwFlEq3+joxNq19p1qohBUsK1CZ1aK7N7N29U6X97y5Rx9BZ2lJqjDBuwFO45I\nJAyzZrGLnpgoNgxaQLwEu6aG92PrCTCqYD/4YH68Djf8jGPhVCni1Oj44IP+Kpy8HPZvf8v7mO40\nZX6OjkR0HOKnY4vfDNvsOBO1Dts0dzrHLjvBdsqvgXAOO5fjOOQDH+D7F17IObZTJOK340xLCy/b\nLg4BigU7lyu+THUS7E2b8q/p6+P1cnMDl17qfjUQhgUL+AD348rMKMF02LffznW01gOkqYl3vN27\ng62XW8cZJ5wctlWwdccT60GSBLNm8cnMOrSqRjt1r0gE4BPeqlWF73XLsFtb+TGniZ/vu6+4t6Qd\nfkaKc8qx7SIRoPB7uOHlsNesAT7+cfvnrNrhN7/W70sjwzZPqFXrsDds4IPjtNP4vlWwTYcTpNGx\npobPvnZxCMA/2MhI/tLQrirFepA9+igL9hlnFDpstzgkSWpq2K34iS1MZ2rt7bh3L/CtbwFf+1rx\n+8Lk2HrwpyCYA/MAecFevpzXce9e/j3iGJ7WDatgm/h12ADw61/zSVDjlWHX1vJynQRvYMC7Kki/\n12uMDCfB9tPt3g03h93Xx3Hn5ZfbP28VbL/5NRCsSiRKhu3W6AjkR+wrO8H2cthBBVvHIZpVq7iL\ndH9/tEZHgH9MJ4et8zedY5txCJA/yMbHgWef5U4X2mHrRptSCTbgX1DtHLY++L/4Re4ybhfphBHs\nOB32tGn8Gy9bxqMrJh2HAPwZ+/bZC3YQh93SUnhy8YpEAPdYpL/fe2hcP3EIENxh+8XNYf/kJ8C7\n3uV8hdTYyMfU0aPBHHaQOuyoDtup0RHIj4ldto2Odkyf7m/GYM2xY8C6dYVn3bY2znvvv989EvEj\n2DNmODtsoLDh0S6C0QfZpk0sXq2t/J6WFj5xAe412EnjV1DtHPaOHdzp4Y9/BG64wXn5QRse42p0\ntDrQU07hK7EgDY5h0Z8R1WGbeEUigHtpX38/v99tuAQ/cQhg3z09lwv321lxcthKucchAJ/ctMv2\n0y1dk2aG7RWJZM5h19Tk+/f74aGH+IcxBe/CC7mBIkqjI+DusIHCHNvaLV2jDzKdX2ussUgWHLbZ\n6KhL+z72MR6j2ukSOkwtdphGR3M2EKCwSqgSBduMRABvh93U5O6yozhsr9ng/eDksF9+mb/7m97k\n/n4t2H66pWv81mFbI5G46rDNRseyy7CVcnfYQLDu6T/+MQ8GZHLhhf7K+rwEe+FC9+oNq2C7OWyd\nX2uslSLbtrl/RpKEjUQAXucTT3Qf8tVr+ddcU/x83GV9AAv2M8+kE4k0NvJ+FzUSMZkxg8VYV90E\njUQGBjiScxNsv2JhVyXiFSn4wclhDw6yA/U6GWjxTdJhK8WvjxqJ2Dns7m7+/n7ceyqCPTDA7smt\npd5vjj08DPziF8CVVxY/d/75fFBE6ekIcD5+0UXOz/sVbDuHvWkTH3w7d+aHCk2bsJEIAHz1q8AP\nfuB+EHmNCvi73xVPpBB3WR/Agn3sWDoOG+DPidth19bycaPNjF0k4uWwzz3XW7D9RCJODjtKfg04\nO+zDh70HQwKiOWy/GfbwMGtL0IHYvCKROXN4rHU/44gAKQm2l7sG/DvsO+/kSyS7M2lTE3dd9Wp0\n9CqVq6lx33jW3o5OkUhvL3/v17wm/7h22Lt38w8VdUcPy9Kl3JPTrhOExmni2re/3Xu41/Z23oZm\nByOAM889e4qHnIzTYVsFGyi9YNfX8z4Z1olaYxG7SMQpwz5yhH/HlSuTi0SiNjgCeRE098dDh4IJ\ndhIOu7aW9WL//uDuGvBudJw7l42qnxMmkJJgO43SZ8Vvad/atZyhOnHjjcD73lf4WFCH7YVXo6M+\noM46q/DkcMopfPnz9NOly68B3h4LFwIvveT8Gj1xbViRcXLx+/dz47LpqOJw2EoVCvbs2SygaUQi\ngLNgA1z+6Nf9mZiC7ddh6xOu21jmgP9IxK7RMWpJn8ZuiNVDh/zVz2vxDeqwDx5ko+AVVTU18bKD\nNjgC3g5bmwk/2x9I0WE7NThq/HSe2bWLu3e/853OrznllGJHZXXYXh1n/OAViegxkM1xKOrquLfW\nnXeWVrAB71gkau9Ap0qRPXv4fxIOe3yct73+fYl4e/s9GKLy3vfy/IR2fO5zwU9IGi3YuiLDFBgn\nwTaHxnWKqKJEInE4bMB+iNUgkUhvL28bvy74uOPyndfcJuwFeD/r6QnnsL0aHevqWK8SF2wiup6I\nNhPRs0R0GxE5Hm5xOexbbuGJXoMe2EEbHb3wqhIB+ECxGzho5Uqet63cBTvogEkmTpUi3d38PwmH\nbRcX3HYbcPHFwZYblk99Kj9oU5xowdZ5qxnXHXecu8OeNYvfYw6poAkSiZiNjuXgsGfM4EZ8v93S\n9Xv8Npg2N4d32C0tvI10T1Q7gzdnTsKCTUTLAHwCwNlKqdMB1AL4oNPrN2/2zj29HHYuB/zoR+5x\niBNBy/q88HLYANcov+1txY+vXMk7YqkqRDRe42JHddjLl9tHLk4OO46OM3YVQsuWhXe25YJVsO3c\nrFNJrBZsIvdYJOsO223+SCcaG3m/8CPYOhIJ47B1rKhPRnZ6MXdu8g77EIAJANOIqA7ANAB77F74\nm9+wYF96qfsCvRz2Qw/xjnnWWcFXNm6HbZ0l20mwr7zSPs884wz+X+4OO+gIdybLlhVXggAs2HZu\nKo5IxE6wKwE/gu3msAFnwVYqWFmfnWAn6bD9RiLaYQfhuOP8nWx0JBLGYQOFsYhdz+jEBVspNQDg\nnwG8DGAvgINKqd+brxsaAj75SS4D87q08XLY2l2HKdBPotHRKxJxYuVK/sGWLYu2DlFZutR9gKao\nkYiTYHd388nKzmHHEYmUqvImSbRg25X0Ad4ZNlA4rIAVv+OIAMmV9enPNx12kEhk375gDhtgwQ4S\niYRx2EDh72Nn8L7yleJCCSdCjQVHRCcA+CyAZQCGAPyciD6slLrV+rq3vKUT8+YBDz8M1NZ2oKOj\nw3GZbg47l+NBcewGGvJD0LFEvND51/i4s8N2YtYszvTjmlg3LDNnul/RRI1ElizhE0IuV9ios2cP\nNwyLw/bPzJlcWWSX0QPuGbYWsZNOAu65p/g1fuMQwL5KJK6TpJ3DDtLoCIRz2H5Mgo5EvHpcOqEd\n9rFjvP3M7dXf34U77ujytaywg3e+DsAGpVQ/ABDRnQDOB1Ag2AMDndi0yf9GdxKQzZv5ea9KEyfi\njkSI8rFIUIcN+BtgPWm8Iqj+/mg5e1NTfqYV67gse/YAb3gDsGVL4evjcNiVKthekUhbm/1YIf39\n+fHUnSIRv3EI4N41PSpODttvhg2Ec9h+rtibmzlyieKwh4fz+6f5mR0dhWZ29erVjssKm2FvBXAe\nETUTEQG4CMDz5otuusnfBgfcO8488EDhkJNBibvREcjHIkEddrmgIyincZQPHIhev2wXi3R3J+uw\nKzUSGRx0/n7HH59vzLUyMFCcYZulfX4rRADnrulJOewgkQiQbIbd1xctw9Zj0EfVirAZ9jMAbgHw\nBAA9LP9/m68LUk7l5rD/8Idogh1kijC/aIedVcGuq+Od1cySNVEbHYFiwdZjYhx/fHIZdqU7bLvv\nt2ABn2DNnoL9/fkMe+ZM3r779hW+JkgkkrbDTiMS8VslMjkZ3mHrSCQOrQhdh62U+pZS6tVKqdOV\nUlcppVw6Onvj5LAnJ7lCxCX+9sTa6BhHxxkg77DDRCLlgluOHbXRESgWbD2xsZ4lxUpYh211fNUa\nidTWsljt3Vv4uN0Ub2YsEtRhJ5lhh210jBKJ+G10BMI7bB2JlFSw40Y7bPOS7amnuNNN2G69QPyN\njkD2IxHAPceOYx5EN8FOwmFXqmDPnOleJQJwO4HulKQxBdsuxw6SYTs1OsblsMOW9U2fnj9pBSGo\nYEdx2IcPF/dyDEPZCHZjI1+mm4P5/OEPPGxqFLTDVkoE20pbG2ejdiQRieiZ6O0Ozjg6zlRqWV9j\nI//19TmfkBYvLhTsXI5/W+vM706CHTUSScJhK+X/6rWmhqto/LaXaWbM8B+JAOKwi7BzfFHza4DP\nvrW1vLMReY8d4IfZs/PjF8ThMEqBk8MeG+PvFXXnCuKwpazPnfZ2FmQ3h22tqz90iLeF1ZzYCXa5\nRCLmSXx0lNfdr7kKMyTAFVfw2OxexOWwS9bomBRm55mjR3mexje/OfqyGxt5h4ir/nnWLO563dqa\n/ASvSeEk2DoOifq9dC22HnxfC7Y+OK3xVxwzzlS7YFsdtrXBURM1Ekl6LBGrw/bb4BiF5cvzZY9u\nNDez4Qu7b5VFo2MSmALy5z8Dp54a/lLESkMDb7C4BHv2bHaPWY1DAHfBjmMM6aYmXk5PD9/v7mbB\nbmjgk4FuVwDEYXsxcyaf/Jy+n51gm7/hq17FIyhaxx2JGokk5bD9NjimQVMTHythDUzFRiKmw44j\nDtE0NsYv2Lt2lc9OFQYnwY6jBltjjUX27Ml3ojEPUOma7o6XwzYzbGsNtqa1FVi1Kt/jUal4IpEk\nHLbfBsc0aG4OH4cAFdroCBQLSJyC3dAQbyQyezYvL8sO26msLy6HDRQLtp7c2Myxwzjs+npuXJuc\n5PuV7LDb21nE/GbYTr/h+94H/PznfFtPUuF3m9lViezfH8/J3TyBpxGJ+KW5OdpVfsVGIrq0b/du\n4O67ueX3jW+MZ9naYcfRaQbI76RZFmynKpE4arA1WrAnJni5uvTKKti5HOfcQX8bokKXXemCDTgL\n9vz5LJ6684xdhg0Al13GRujQoWD5NVDssLVDj1Jyq7Fz2OVy9RrVYevBnyqu0XH2bODTnwZe9zrg\n+98Hvv3t+A7AuB329OksMOWyU4UhzUikp4ezUi3KVkc1Pp7PtYNiFexKj0QA5+Ohvp7Ft7eX7zs5\n7LY24IILgHvvDRaHAMWNjsPDwRy6G3YZdrk47Ne8BnjXu8K/v2Id9g03cOVFby/w29/yDB5xEXeV\nCBGLWtYdtp1gx+WagLxgW/NroNBhh8mvNeKw8yxenI9F7DJszfvfz7FIkAZHoNhh9/XFt5+UokrE\nL6eeCnzmM+HfX7GNjq2tPDlsEmVycQs2ULmCHfRS2Y1ly7hx1ppfA8UOO2h+rRHBzmOtFHFrh7js\nMh5QbefOaJFInIJdzlUiUanYRsckibusD2DBzvJO5eawgzgvN5YuZRF56aVCwRaHHQyvSAQoFmy7\nDBvgxuY3vQlYuzaYYJuNjn198e0n5VwlEpWKjUSSJO6yPqByHXacgt3YyNvp8ccLIxGrowrTLV1T\nbRm2X4ftFokAHIts3Fg+kYgeyEsP95vlQdVM9Gh/g4Mi2L6Ju9ERyL7DnjGDT2LmmNhxCjbAscgj\njzg77DgiEaWqw2H7zbC9SjMvu4wbgMslEjErfirJYRPx/t7TI4LtmyQc9t//PXDVVfEtL21qangH\nsnZWOnaMnUBcZX0AC7bu5agxHXbUSOToURaguMo2yw09S1GQSMTtN2xvBz70IZ5Mwi/19ewU9Qk+\nTsEGCicxKKdGxzhobeWr2aiCXaG7dzFJNDqWeiLdONCxiBaEgQF23nEKn95OSTrsSo5DABbqb3zD\n/cSmBXtigvd1L8G7+eZg60CUH/myuTneaiKgcBKDSmp0BPLfpWSNjkTURkTriGgLET1PROdFW5Vk\nSaLRsRIwc+y44xDAXrDjdNhjY5UdhwAsltdd515BtXAhu979+/kEHMeolCbWWCRJh11JkQiQHyQu\najf+KD7quwB+rZR6HxHVASjrw6WxkXfkSr1kDktagt3WViiocTvsShdsP9TXc7vK5s3xRlpWrJUi\ncVaJAIUOu5IaHQH+LnGM7BnqHExEMwBcoJRaAwBKqUmllMMUuuVBEo2OlUAagn3GGVyVYCXuDFsE\nm1m0CHjmmeQEWxx2OFpb46koC3vRtBzAfiJaS0RPEdFNRFTWCWISGXYlYA4AlYRgz5sH/LcxRbNk\n2MmQlmCPjvLvFmWMDRMzwxbBLiZsQFAH4GwAf6uUepyIvgPgOgBfsb6os7PzldsdHR3oiDKTbkTE\nYdtjDgAVZy9HN8RhJ8OiRUBXF3D22cksX48nouOQOHsla4edy2V/JEyT6dOdI56uri50dXX5Wk5Y\nwe4G0K2Uenzq/jqwYBdgFexSk0RZXyVgF4msXJn855o9HcM67KYmLmETwWYWLwa2bAEuvjiZ5WuH\nfeRIvHEIkHfYIyP5WV4qBTeHbZrZ1atXOy4nVCSilOoFsJuITp566CIAm8MsKy0kErEnjQzbjrjH\nEjlyRAQbYIc9OZl8o2PcJX1Avnt6pcUhQL7RMSpRaiY+DeBWImoA8AKAq6OvTnJIWZ89pRLsuMcS\nGRmRDBvId/9POsOOu8ERyJ/EK61CBCh9hg2l1DMAXh99FdKhsZGdnAh2IXaCnXaGLWV98aEF22ng\np6hYBTvuE3slO+w5c+IZY75qqpK1gxPBLsQU7KBjJIclCYctgp3vnJSkw9aNjitWxLvslhbuaVtp\nvRwB4MMfBj7wgejLqaqxRADpOGNiLesbH2cRjWOWei8aGzlrnZiQsr44aWjgqCKLkYh22JU2jgjA\nutPUFH05VSfY4rALsZb16ZK+JLo0m+gRzEZGxGHHzc03A6edlsyydaNjkhl2JUYicVE1gi2RiD3W\nSCStBkeNPkAlw46Xt789uf08LYddaZFIXFSNYIvDtqe1lQ+Sycn0Ghytnz08HI/DlrK+dEijSkQc\ntjNVI9jisO2pqeHuxUND6TU4avQBGseMM1LWlw6NjeyADx+OvxLFWiUiDtueqhFscdjO6Fgk7UhE\nO2yJRLJDYyOPuT17dvxtHdY6bHHY9lSNYIvDdkZXipQqw5ZGx+zQ2Ai8/HL8cQhQ2XXYcVE1gi0O\n2xldKZJVhz02Jhl2WjQ1JSfY1gxbIhF7RLCFgkgkzUbHuBy2nnFGMuzkScNhSyTiTNUIthYE6ThT\njBbstBsd43DYNTX8mw4OisNOg8ZGbqBO2mGLYNtTNYItDtuZUjU6xuGwAXbZBw+KYKeBPo6SzrAl\nErGnagRbGh2dyXKVCJCf2FQikeTRv1MS+0ldHY+BfeCAOGwnqkawxWE709YG7NkDKJWuS21pid5x\nBmDBbmysrAHvyxU9HkYSDhvgk+7goDhsJ0SwBcycCWzbxg2OcU755IUeSySqw25qkjgkLZKMRAD+\nHYnk93SiagRbIhFn2tpYsNOMQ4B4uqYD7LDlAE+HpAV72jR212kMQJZFIm0WIqoloqeJ6J64Vigp\ndDWBCHYxbW2cG6Yt2HF0TQdYsCW/TofGRnbASZV/asEW7Ila5HYtgOcBZGITNzSIYNuhx78ulcOO\no9FRHHY6NDbyWNtJlce2tPD+INgT2mET0SIA7wDwAwApJp/haWwUwbZDC3aanWaAeMv6RLDTYdas\n+GeasSIO250o58l/BfBFAJkpwBGHbU8lOOyJifjWS3DmhBOAP/85ueXLidedUIJNRO8EsE8p9TQR\ndTi9rrOz85XbHR0d6OhwfGkqNDZKT0c7Wlp4u5Qyw47qsCcn41svwZ0kK4mmTau+8syuri50dXX5\nem1Y+TofwKVE9A4ATQCOI6JblFIftb7IKtjlwGtfm/5lfxYgYpddCod9+DCLbZQrn+ZmIJeLb72E\n0tHSEu3knUVMM7t69WrH14YSbKXUDQBuAAAiWgXg70yxLkfuvLPUa1C+zJyZXKmWEy0t3EmioSGa\na9M9HYXsM20ad+AS7IkrIJBNnHHuvht41avS/czmZj44o+TXejlpdvgRkkN3nBHsiSzYSqn1ANbH\nsC5CCTn11PQ/s6aGHVXUS+DXvEZKwSqFa64RwXaDVELXH0Skklq2UDnMn8+NTHv2lHpNBKE8ICIo\npWxPW1IzIZSUlhbJLAXBL9JjXygpra3VVxUgCGERhy2UFOkoIQj+EcEWSkprq/RSFAS/SCQilJSW\nluhlfYJQLYhgCyVFMmxB8I8ItlBSxGELgn9EsIWSIg5bEPwjgi2UFHHYguAfEWyhpIjDFgT/iGAL\nJWX6dBltTxD8ImOJCCVlaIjHxF60qNRrIgjlgdtYIiLYgiAIZYSbYEskIgiCkBGizJq+mIgeJKLN\nRPQcEX0mzhUTBEEQConisCcAfE4p9WoA5wH4GyIqwTD42cDvJJvVgmyPPLItCpHt4UxowVZK9Sql\nNk7dHgawBcDCuFas0pCdsBDZHnlkWxQi28OZWDJsIloG4CwAj8axPEEQBKGYyIJNRK0A1gG4dspp\nC4IgCAkQqayPiOoB3AvgN0qp7xjPSU2fIAhCCGKvwyYiAnAzgH6l1OcirJsgCILggyiC/SYADwHY\nBEAv5Hql1G9jWjdBEATBQmI9HQVBEIR4SaSnIxFdQkRbiWg7Ef19Ep9Rrjh1KCKidiK6n4i2EdF9\nRNRW6nVNEyKqJaKnieieqftVuz2IqI2I1hHRFiJ6nojOrdbtQUTXTx0rzxLRbUTUWK3bwg+xCzYR\n1QL4DwCXADgNwBVV1qHGqUPRdQDuV0qdDOCBqfvVxLUAnkc+Pqvm7fFdAL9WSp0K4AwAW1GF22Oq\nHPgTAM5WSp0OoBbAB1GF28IvSTjscwDsUEq9qJSaAPA/AC5L4HPKEocORccDuBTcSIup/+8qzRqm\nDxEtAvAOAD8AoFu/q3J7ENEMABcopdYAgFJqUik1hOrcHofABmcaEdUBmAZgL6pzW/giCcE+HsBu\ny0WjytoAAAH0SURBVP3uqceqDqND0TylVN/UU30A5pVotUrBvwL4IoCc5bFq3R7LAewnorVE9BQR\n3URELajC7aGUGgDwzwBeBgv1QaXU/ajCbeGXJARbWjHxSoeiO8Adig5bn5sad7YqthMRvRPAPqXU\n08i76wKqaXsAqANwNoDvKaXOBjAC45K/WrYHEZ0A4LMAloGHtWgloo9YX1Mt28IvSQj2HgCLLfcX\ng1121TDVoegOAD9WSt019XAfEc2fen4BgH2lWr+UOR/ApUS0C8DtAN5KRD9G9W6PbgDdSqnHp+6v\nAwt4bxVuj9cB2KCU6ldKTQK4E8AbUJ3bwhdJCPYTAE4iomVE1ADgAwB+mcDnlCVTHYp+COB5o/fn\nLwFcNXX7KgB3me+tRJRSNyilFiulloMblP6glLoS1bs9egHsJqKTpx66CMBmAPeg+rbHVgDnEVHz\n1HFzEbhhuhq3hS8SqcMmov8F4DvgVt8fKqW+EfuHlClOHYoAPAbgZwCWAHgRwOVKqYOlWMdSQUSr\nAHxBKXUpEbWjSrcHEa0EN8A2AHgBwNXgY6XqtgcRfQksyjkATwH4PwCmowq3hR+k44wgCEJGkCnC\nBEEQMoIItiAIQkYQwRYEQcgIItiCIAgZQQRbEAQhI4hgC4IgZAQRbEEQhIwggi0IgpAR/j852X/Q\nd4JSlgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# or a trace plot through time:\n", "df.price.plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next we will test the statsmodels library:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " OLS Regression Results \n", "==============================================================================\n", "Dep. Variable: price R-squared: 0.002\n", "Model: OLS Adj. R-squared: -0.008\n", "Method: Least Squares F-statistic: 0.2010\n", "Date: Tue, 25 Aug 2015 Prob (F-statistic): 0.655\n", "Time: 07:35:12 Log-Likelihood: -249.35\n", "No. Observations: 100 AIC: 502.7\n", "Df Residuals: 98 BIC: 507.9\n", "Df Model: 1 \n", "Covariance Type: nonrobust \n", "==============================================================================\n", " coef std err t P>|t| [95.0% Conf. Int.]\n", "------------------------------------------------------------------------------\n", "Intercept 9.6340 0.596 16.160 0.000 8.451 10.817\n", "time 0.0046 0.010 0.448 0.655 -0.016 0.025\n", "==============================================================================\n", "Omnibus: 0.653 Durbin-Watson: 2.125\n", "Prob(Omnibus): 0.721 Jarque-Bera (JB): 0.250\n", "Skew: 0.036 Prob(JB): 0.883\n", "Kurtosis: 3.234 Cond. No. 117.\n", "==============================================================================\n", "\n", "Warnings:\n", "[1] Standard Errors assume that the covariance matrix of the errors is correctly specified.\n" ] } ], "source": [ "import statsmodels.formula.api as smf\n", "formula = \"price ~ time\"\n", "results = smf.ols(formula,df).fit()\n", "print(results.summary())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Markdown Text and Notebook Annotations\n", "\n", "Changing the execution cell type to markdown here:\n", "\n", "![](../site_pics/markdown1.png)\n", "\n", "allows you to enter markdown compliant passages and get nicely formatted text when you execute the cell by pressing ``shift enter``. For example, executing the following markdown code in a markdown execution cell:\n", "\n", "```\n", "This is some markdown code. I can type sentences and paragraphs here. Or, if I want a numbered list, here it is:\n", "\n", "1. Dogs\n", "2. Cats\n", "3. Elephants\n", "\n", "I can include links in my markdown passages:\n", "![](https://whyevolutionistrue.files.wordpress.com/2012/08/cute-kitten-saying-hello.jpg)\n", "\n", "And I can include latex equations inline: $\\frac{\\alpha}{\\beta}$. Or I can put them on their own line\n", "$$\n", "W_j = \\sum_{i\\ne j} W_{ij} \\epsilon_i\n", "$$\n", "```\n", "\n", "results in this output:\n", "\n", "This is some markdown code. I can type sentences and paragraphs here. Or, if I want a numbered list, here it is:\n", "\n", "1. Dogs\n", "2. Cats\n", "3. Elephants\n", "\n", "I can include links in my markdown passages:\n", "![](../site_pics/cute-kitten-saying-hello.jpg)\n", "\n", "And I can include latex equations inline: $\\frac{\\alpha}{\\beta}$. Or I can put them on their own line\n", "$$\n", "W_j = \\sum_{i\\ne j} W_{ij} \\epsilon_i\n", "$$" ] } ], "metadata": { "hide_input": false, "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.6.6" } }, "nbformat": 4, "nbformat_minor": 1 }