{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import statsmodels.api as sm\n", "from matplotlib import pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEKCAYAAADTgGjXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAGiZJREFUeJzt3X+UXHV9//HnK0HBVBSBCCnJZiOk2kAVZUsB+0ORSqqUYFskui1YbHOsyI/WlkLTL1Vxj6j9dRCpXbUV7FZMD1pSVChEoVQP4gYDIUEkQjaQBghafrTbw6+8v3/cz5LJsjN79+7cmbkzr8c5c+69n3tn5n3nwr7zuffzQxGBmZlZEXPaHYCZmVWXk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZW2F7tDqBsBx54YPT397c7DDOzSlm/fv2jETF/uuO6Pon09/czOjra7jDMzCpF0lie43w7y8zMCnMSMTOzwpxEzMysMCcRMzMrzEnEzMwKcxIxM+sxIyPQ3w9z5mTLkZHin9X1TXzNzGy3kRFYtQrGx7PtsbFsG2BwcOaf55qImVkPWb16dwKZMD6elRfhJGJm1kO2bZtZ+XScRMzMekhf38zKp+MkYmbWQ4aGYN68PcvmzcvKi3ASMTPrIYODMDwMixeDlC2Hh4s9VAe3zjIz6zmDg8WTxmSuiZiZWWFOImZmVpiTiJmZFeYkYmZmhbU9iUiaK+n7kq5N2/tLukHSvWn5ippjL5S0RdI9kk5sX9RmZgYdkESAc4G7a7YvANZFxFJgXdpG0jJgJXA4sBy4XNLcFsdqZmY12ppEJC0E3g58rqZ4BXBFWr8COKWm/KqIeCoi7ge2AEe3KlYzM3uhdtdE/hY4H9hVU3ZQROxI6w8BB6X1Q4AHao57MJWZmVmbtC2JSDoJeCQi1tc7JiICiAKfvUrSqKTRnTt3ziZMMzNroJ01kTcCJ0vaClwFHC/pn4CHJS0ASMtH0vHbgUU171+Yyl4gIoYjYiAiBubPn19W/GZmPa9tSSQiLoyIhRHRT/bA/JsR8dvAWuCMdNgZwDVpfS2wUtLekpYAS4HbWhy2mZnV6MSxsy4B1kh6LzAGvBMgIjZJWgNsBp4FzoqI59oXppmZKXvs0L0GBgZidHS03WGYmVWKpPURMTDdce1unWVmZhXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZWgpER6O+HOXOy5chIuyMqRyf2EzEzq7SREVi1CsbHs+2xsWwbmje3eadwTcTMrMlWr96dQCaMj2fl3cZJxMysybZtm1l5lTmJmJk1WV/fzMqrzEnEzKzJhoZg3rw9y+bNy8q7jZOImVmTDQ7C8DAsXgxSthwe7r6H6uDWWWZmpRgc7M6kMZlrImbWU3ql/0aruCZiZj2jl/pvtIprImbWM3qp/0arOImYWc/opf4breIkYmY9o5f6b7SKk4iZ9Yxe6r/RKk4iZtYzeqn/Rqu4dZaZ9ZRe6b/RKq6JmJlZYU4iZtaR3CmwGnw7y8w6jjsFVodrImbWcdwpsDqcRMys47hTYHU4iZhZx3GnwOpwEjGzjuNOgdXhJGJmHcedAqvDrbPMrCO5U2A1uCZiZmaFOYmYmVlhTiJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVljbkoikRZK+JWmzpE2Szk3l+0u6QdK9afmKmvdcKGmLpHskndiu2M3MLDNtEpF0qKS90/qbJJ0jab8mfPezwAcjYhlwDHCWpGXABcC6iFgKrEvbpH0rgcOB5cDlkuY2IQ4zMysoT03kauA5SYcBw8Ai4J9n+8URsSMibk/rTwJ3A4cAK4Ar0mFXAKek9RXAVRHxVETcD2wBjp5tHGZmVlyeJLIrIp4F3gF8KiL+BFjQzCAk9QOvB74LHBQRO9Kuh4CD0vohwAM1b3swlZmZWZvkSSLPSHoXcAZwbSp7UbMCkPRSstrOeRHxRO2+iAggCnzmKkmjkkZ37tzZpEjNzGyyPEnkd4FjgaGIuF/SEuCLzfhySS8iSyAjEfGVVPywpAVp/wLgkVS+nexW2oSFqewFImI4IgYiYmD+/PnNCNXMzKYwbRKJiM3AnwITzy/uj4iPz/aLJQn4PHB3RPx1za61ZLUe0vKamvKVkvZOiWwpcNts4zAzs+LytM76dWADcF3aPlLS2iZ89xuB3wGOl7Qhvd4GXAL8qqR7gRPSNhGxCVgDbE6xnBURzzUhDjMzKyjPfCIfImsFdRNARGyQ9KrZfnFE/CegOrvfUuc9Q4DnNjMz6xC5HqxHxOOTynaVEYyZmVVLnprIJknvBuZKWgqcA3yn3LDMzKwK8tREzibrJf4U8CXgCeC8MoMyM7NqmLYmEhHjwOr0MjMze17dJCLp32jQ0S8iTi4lIjMzq4xGNZG/bFkUZmZWSXWTSETc3MpAzMysehrdzloTEe+UtJEpbmtFxGtLjczMzDpeo9tZ56blSa0IxMzMqqduE9+a4djfHxFjtS/g/a0Jz8zMOlmefiK/OkXZrzU7EDMzq55Gz0T+gKzG8SpJd9bs2hf4dtmBmZlZ52v0TOSfgW8AHyPNc548GRE/KTUqMzOrhEZNfB8HHgfeJWku2TS1ewEvlfTSiNjWohjNzKxDTTvsiaQPkA0H/zC7R+8NwE18zcx6XJ5RfM8DXh0RPy47GDMzq5Y8rbMeILutZWZmtoc8NZH7gJskfY1sOHgAJs2LbmZmPShPEtmWXi9OLzMzMyDffCIfbkUgZmZWPXlaZ80Hzieb3XCfifKIOL7EuMzMrALyPFgfAX4ALAE+DGwFvldiTGZmVhF5ksgBEfF54JmIuDkizgRcCzEzs1wP1p9Jyx2S3g78F7B/eSGZmVlV5EkiH5X0cuCDwKeAlwF/WGpUZmZWCXlaZ12bVh8H3lxuOGZmViV5Wmf9I1NPj3tmKRGZmVll5Hmwfi3wtfRaR3Y763/KDMrMqmtkBPr7Yc6cbDky0u6IrEx5bmddXbst6UvAf5YWkZlV1sgIrFoF4+PZ9thYtg0wONi+uKw8eWoiky0FXtnsQMys+lav3p1AJoyPZ+XWnfI8E3mS7JmI0vIh4E9LjsvMKmhbnanq6pVb9eW5nbVvKwIxs+rr68tuYU1Vbt2p4e0sSS+R9HuS/jq93i3JI/ma2ZSGhmDevD3L5s3Lyq071U0ikn4O2Az8Etl4WVuBE4FvS9pP0kdbEaCZVcfgIAwPw+LFIGXL4WE/VO9mjW5nXQqsiogbagslnQDcBWwqMzAza72Rkewh+LZt2S2ooaGZJ4DBQSeNXtIoiSyYnEAAIuJGSc8A7ygvLDNrNTfPtSIaPROZI2nvyYWS9iEb0Xd8iveYWUW5ea4V0SiJXAlcLWnxRIGkfmAN8MVyw6pP0nJJ90jaIumCdsVh1m3cPNeKqJtEIuKjwHXALZIelfQocDNwQ0Rc3KoAa0maC3wa+DVgGfAuScvaEYtZt6nXDNfNc62Rhk18I+KyiOgjm9VwSUQsjohPtSa0KR0NbImI+yLiaeAqYEUb4zHrGm6ea0XkGvYkIp6MiCfLDiaHQ4AHarYfTGV7kLRK0qik0Z07d7YsOLN2m83gh26ea0UUGTur40XEcEQMRMTA/Pnz2x2O2YwUTQQTravGxiBid+uqmSaSrVth165s6QRi06laEtkOLKrZXpjKzLrCbBKBW1dZO0ybRCTNk/T/JH02bS+VdFL5oU3pe8BSSUvS8CsrgbVtisWs6WaTCNy6ytohT03kH4GngGPT9nagLUOeRMSzwAeA64G7gTUR4Z7z1jVmkwjcusraIU8SOTQiPgE8A5A6GarUqBqIiK9HxM9ExKER4XYj1lVmkwjcusraIU8SeVrSS0jzrEs6lKxmYmZ1FH04PptE4NZV1g7TzicC/AVZp8NFkkaANwLvKTMosyqbzRhUE/uLDoLowQ+t1RQR0x8kHQAcQ3Yb69aIeLTswJplYGAgRkdH2x2G9ZD+/qknZlq8OGs2a1YFktZHxMB0x9WtiUh6w6SiHWnZJ6kvIm6fTYBm3cqtpKyXNLqd9VcN9gVwfJNjMesKniLWekndJBIRb25lIGbdYmhoz2ci4FZS1r3ydDbcR9IfSfqKpKslnZfmFDGzKbiVlPWSaR+sS1oDPAn8Uyp6N7BfRJxacmxN4QfrZmYzN+sH6zWOiIjaOTu+JWlz8dDMzKxb5OlseLukYyY2JP0C4H/am5lZrprIUcB3JE00UOwD7pG0EYiIeG1p0ZmZWUfLUxNZTjaz4a+k15JUdhLw6+WFZtY+s5ncyayXTFsTiYgxSa8gm8djr5pydza0rjSbYUvMek2eJr4XA3cCl5J1QPwr4C9LjstsVmZTk/DkTmb55Xkm8k6y4eCfLjsYs2aYbU3Cw5aY5ZfnmchdwH5lB2LWLLOtSXhyJ7P88iSRjwHfl3S9pLUTr7IDMytqtjUJT+5kll+e21lXAB8HNgK7yg3HbPZmOwDibOf0MOsleZLIeERcWnokZk3SjAEQPbmTWT55ksgtkj4GrKVmWlw38bVO5ZqEWevkSSKvT8tjaso8n4h1NNckzFojT2dDzytiZmZTylMTQdLbgcOB5+cRiYiPlBWUmZlVQ54e658BTgPOBgScCiwuOS4zM6uAPP1EjouI04H/jogPA8cCP1NuWGZmVgV5ksj/peW4pJ8GngEWlBeSmZlVRZ5nItdK2g/4JHA7Wcusz5YalZmZVcK0NZGIuDgiHouIq8mehbwmIi4qPzTrVZ7Lw6w66iYRST8v6eCa7dOBNcDFkvZvRXDWeyZG4B0bg4jdI/A6kZh1pkY1kb8HngaQ9MvAJcCVwOPAcPmhWdU0owbhuTzMqqXRM5G5EfGTtH4aMJxuaV0taUP5oVmVNGs2QM/lYVYtjWoicyVNJJm3AN+s2Zerk6L1jmbVIDyXh1m1NEoiXwJulnQNWTPfWwAkHUZ2S8vsec2qQXguD7NqqZtEImII+CDwBeAXIyJq3nN2+aFZlTSrBjE4CMPDsHgxSNlyeNiDKZp1qoa3pSLi1inKflheOFZVzZjDY4JH4DWrjjw91s2m5RqEWW9qSxKR9ElJP5B0p6Svph7xE/sulLRF0j2STqwpP0rSxrTvUklqR+xW3+AgbN0Ku3ZlSycQs+7XrprIDcAREfFa4IfAhQCSlgEryYadXw5cLmlues/fAb8PLE2v5a0O2szM9tSWJBIR/x4Rz6bNW4GFaX0FcFVEPBUR9wNbgKMlLQBeFhG3pgf8VwKntDzwLuKhRcysGTqhv8eZwJfT+iFkSWXCg6nsmbQ+udwKaFbHQDOz0moikm6UdNcUrxU1x6wGngWa+u9gSaskjUoa3blzZzM/uit4aBEza5bSaiIRcUKj/ZLeA5wEvKWmD8p2YFHNYQtT2XZ23/KqLa/33cOk8b0GBgai3nG9ykOLmFmztKt11nLgfODkiKj9N/FaYKWkvSUtIXuAfltE7ACekHRMapV1OnBNywPvEh5axMyapV2tsy4D9gVukLQhzeNORGwiG25+M3AdcFZEPJfe837gc2QP238EfKPlUXcJDy1iZs3SlgfrEXFYg31DwAv+nEXEKHBEmXH1iomH56tXZ7ew+vqyBOKH6mY2U53QOsvawEOLmFkzeNgTMzMrzEnEzMwKcxLpIO5FbmZV42ciHcK9yM2silwT6RDuRW5mVeQk0iHci9zMqshJpEO4F7mZVZGTSIdwL3IzqyInkQ7h6WXNrIrcOquDuBe5mVWNayJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZWmJOImZkV5iRiZmaFOYmYmVlhTiJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZWmJOImZkV5iRiZmaFOYmYmVlhTiJmZlaYk4iZmRXW1iQi6YOSQtKBNWUXStoi6R5JJ9aUHyVpY9p3qSS1J2ozM5vQtiQiaRHwVmBbTdkyYCVwOLAcuFzS3LT774DfB5am1/KWBmxmZi/QzprI3wDnA1FTtgK4KiKeioj7gS3A0ZIWAC+LiFsjIoArgVNaHrGZme2hLUlE0gpge0TcMWnXIcADNdsPprJD0vrk8hkbGYH+fpgzJ1uOjBT5FDMzA9irrA+WdCNw8BS7VgN/RnYrq6zvXgWsAujr63u+fGQEVq2C8fFse2ws2wYYHCwrGjOz7lVaTSQiToiIIya/gPuAJcAdkrYCC4HbJR0MbAcW1XzMwlS2Pa1PLq/33cMRMRARA/Pnz3++fPXq3Qlkwvh4Vm5mZjPX8ttZEbExIl4ZEf0R0U92a+oNEfEQsBZYKWlvSUvIHqDfFhE7gCckHZNaZZ0OXDPT7962bWblZmbWWEf1E4mITcAaYDNwHXBWRDyXdr8f+BzZw/YfAd+Y6efX3NnKVW5mZo2V9kwkr1Qbqd0eAoamOG4UOGI23zU0tOczEYB587JyMzObuY6qiZRtcBCGh2HxYpCy5fCwH6qbmRXV9ppIqw0OOmmYmTVLT9VEzMysuZxEzMysMCcRMzMrzEnEzMwKcxIxM7PClA2K270k7QTGZvi2A4FHSwinHbrpXKC7zqebzgW663x8LrA4IuZPd1DXJ5EiJI1GxEC742iGbjoX6K7z6aZzge46H59Lfr6dZWZmhTmJmJlZYU4iUxtudwBN1E3nAt11Pt10LtBd5+NzycnPRMzMrDDXRMzMrDAnEUDSqZI2SdolqW4rBklbJW2UtEHSaCtjzGsG57Jc0j2Stki6oJUxzoSk/SXdIOnetHxFneM69tpM91src2naf6ekN7QjzjxynMubJD2ersMGSRe1I848JP2DpEck3VVnf2WuC+Q6n3KuTUT0/Av4WeDVwE3AQIPjtgIHtjve2Z4LMJdsYq9XAS8G7gCWtTv2OrF+ArggrV8AfLxK1ybPbw28jWySNQHHAN9td9yzOJc3Ade2O9ac5/PLwBuAu+rsr8R1mcH5lHJtXBMBIuLuiLin3XE0Q85zORrYEhH3RcTTwFXAivKjK2QFcEVavwI4pY2xFJHnt14BXBmZW4H9JC1odaA5VOm/m2lFxH8AP2lwSFWuC5DrfErhJDIzAdwoab2kVe0OZhYOAR6o2X4wlXWigyJiR1p/CDioznGdem3y/NZVuR554zwu3f75hqTDWxNaKapyXWai6demZyalknQjcPAUu1ZHxDU5P+YXI2K7pFcCN0j6Qcr+LdWkc+kYjc6ndiMiQlK95oQdcW2M24G+iPgfSW8D/hVY2uaYLFPKtemZJBIRJzThM7an5SOSvkpWvW/5H6omnMt2YFHN9sJU1haNzkfSw5IWRMSOdCvhkTqf0RHXZgp5fuuOuh4NTBtnRDxRs/51SZdLOjAiqjgOVVWuSy5lXRvfzspJ0k9J2ndiHXgrMGUriAr4HrBU0hJJLwZWAmvbHFM9a4Ez0voZwAtqWh1+bfL81muB01NroGOAx2tu4XWSac9F0sGSlNaPJvsb8+OWR9ocVbkuuZR2bdrdoqATXsA7yO53PgU8DFyfyn8a+HpafxVZa5Q7gE1kt47aHnuRc0nbbwN+SNbapiPPJcV5ALAOuBe4Edi/atdmqt8aeB/wvrQu4NNp/0YatBBs9yvHuXwgXYM7gFuB49odc4Nz+RKwA3gm/T/z3qpel5znU8q1cY91MzMrzLezzMysMCcRMzMrzEnEzMwKcxIxM7PCnETMzKwwJxGrBEkH1Iw++pCk7Wn9MUmbWxzLkanH78T2yUVHQk6jDx84RfnLJV2ZRpD9kaSReiMYz0ajc5H0IUl/3OzvtO7iJGKVEBE/jogjI+JI4DPA36T1I4Fdzf4+SY1GcziSrL/ERGxrI+KSJofweeC+iDgsIg4FtgBfaPJ3QGvOxbqYk4h1g7mSPqtsHpV/l/QSAEmHSrouDcp4i6TXpPJ+Sd9MA9Gtk9SXyr8g6TOSvgt8IvWE/wdJt0n6vqQVqaf2R4DTUk3oNEnvkXRZ+oyDJH1V0h3pdVwq/9cUx6bpBoiUdBhwFHBxTfFHgNdJerWyeSGurTn+MknvSesXSfqepLskDdf0UL5J0sfTufxQ0i9Ndy6TYqr3W56avusOSZ0wzIy1mJOIdYOlwKcj4nDgMeA3U/kwcHZEHAX8MXB5Kv8UcEVEvBYYAS6t+ayFZD15/4hsAMhvRsTRwJuBTwIvAi4CvpxqRl+eFMulwM0R8TqyuR02pfIzUxwDwDmSDmhwPsuADRHx3ERBWv8+2XwxjVwWET8fEUcALwFOqtm3VzqX84C/iGw490bnUqveb3kRcGI635Onic26UM8MwGhd7f6I2JDW1wP9kl4KHAf8S/rHOMDeaXks8Btp/YtkE19N+JeaP95vBU6ueS6wD9A3TSzHA6fD83/4H0/l50h6R1pfRJb4yhhT6s2SzgfmAfuTJbF/S/u+kpbrgf68HzjNb/lt4AuS1tR8vvUQJxHrBk/VrD9H9i/wOcBj6bnJTPxvzbqA34xJk3xJ+oWZfKCkNwEnAMdGxLikm8gSUj2bgSMlzYmIXekz5gCvIw3nzZ53EfZJx+xDVkMYiIgHJH1o0vdM/E7PMbP/9+v+lhHxvvR7vB1YL+moiKjqgItWgG9nWVeKbNjr+yWdCs/Pl/26tPs7ZCPQAgwCt9T5mOuBs2ueK7w+lT8J7FvnPeuAP0jHz5X0cuDlwH+nBPIasqlWG8W+hezW1Z/XFP85sC4itgFjwDJJe0vaD3hLOmYiYTyaag+/1eh7cpzLRDx1f0tJh0bEdyPiImAnew6dbj3AScS62SDwXkkTo/tOTOV6NvC7ku4Efgc4t877LyZ7BnKnpE3sftD9LbI/4hsknTbpPeeS3VLaSHbbaBlwHbCXpLuBS8hGUJ3OmWTDrv9I0k6yxPM+gIh4AFhDNtz9GrKEQ0Q8Bnw2lV9PNnT7dBqdS616v+UnJW2UdBdZcr4jx3daF/EovmYdTtKrga8B50TE19sdj1ktJxEzMyvMt7PMzKwwJxEzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCvv/212WbtntAFMAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data = sm.datasets.longley.load()\n", "data.exog = sm.add_constant(data.exog)\n", "model = sm.OLS(data.endog, data.exog)\n", "mod_fit = model.fit()\n", "res = mod_fit.resid # residuals\n", "probplot = sm.ProbPlot(res)\n", "probplot.qqplot()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEKCAYAAADTgGjXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAGvJJREFUeJzt3X+UXWV97/H3J0HBVBSBCCnJZCKk2sBVlCkFtK0CLSlSgrcXiU4LFtssK/Lj1l4KN71UxVmCtPYuROodtRXsVEwXWlJUEKJwqS7ACQZCgkiATCA3QNDyo3fu4le+94/9DDkZ5pzZs3PO2fuc83mtddbe+9l7zv7OI843z97PD0UEZmZmRcwqOwAzM+tcTiJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYXuUHUCr7b///tHf3192GGZmHWXt2rVPRsTc6a7r+iTS39/P6Oho2WGYmXUUSWN5rvPjLDMzK8xJxMzMCnMSMTOzwpxEzMysMCcRMzMrzEnEzKzHjIxAfz/MmpVtR0aKf1fXd/E1M7OdRkZgxQoYH8+Ox8ayY4DBwZl/n1siZmY9ZOXKnQlkwvh4Vl6Ek4iZWQ/ZsmVm5dNxEjEz6yF9fTMrn46TiJlZDxkagjlzdi2bMycrL8JJxMyshwwOwvAwLFwIUrYdHi72Uh3cO8vMrOcMDhZPGpO5JWJmZoU5iZiZWWFOImZmVpiTiJmZFVZ6EpE0W9JPJF2fjveVdJOkB9L2DTXXXihpk6T7JZ1QXtRmZgYVSCLAucB9NccXAGsiYjGwJh0jaQmwHDgUWApcKWl2m2M1M7MapSYRSfOB9wJfrileBlyV9q8CTqkpvyYinouIh4FNwJHtitXMzF6p7JbI/wTOB3bUlB0QEdvS/mPAAWn/IOCRmuseTWVmZlaS0pKIpJOAJyJibb1rIiKAKPDdKySNShrdvn377oRpZmYNlNkSeSdwsqTNwDXAsZL+EXhc0jyAtH0iXb8VWFDz8/NT2StExHBEDETEwNy5c1sVv5lZzystiUTEhRExPyL6yV6Yfz8i/gBYDZyRLjsDuC7trwaWS9pT0iJgMXBnm8M2M7MaVZw76xJglaQPA2PA+wEiYoOkVcBG4EXgrIh4qbwwzcxM2WuH7jUwMBCjo6Nlh2Fm1lEkrY2IgemuK7t3lpmZdTAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRM7MZGhmB/n6YNSvbjoyUHVF5qjhOxMysskZGYMUKGB/PjsfGsmNo3rrlncQtETOzGVi5cmcCmTA+npX3IicRM7MZ2LJlZuXdzknEzGwG+vpmVt7tnETMzGZgaAjmzNm1bM6crLwXOYmYmc3A4CAMD8PChSBl2+Hh3nypDu6dZWY2Y4ODvZs0JnNLxMw6nsdtlMctETPraB63US63RMyso3ncRrmcRMyso3ncRrmcRMyso3ncRrmcRMyso3ncRrmcRMyso3ncRrncO8vMOp7HbZTHLREzMyvMScTM2s6DA7uHH2eZWVt5cGB3cUvEzNrKgwO7i5OImbWVBwd2FycRM2srDw7sLk4iZtZWHhzYXZxEzKytPDiwu7h3lpm1nQcHdg+3RMzMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRMzMrzEnEzMwKKy2JSFog6QeSNkraIOncVL6vpJskPZC2b6j5mQslbZJ0v6QTyordzMwy0yYRSQdL2jPtv1vSOZL2acK9XwQ+HhFLgKOAsyQtAS4A1kTEYmBNOiadWw4cCiwFrpQ0uwlxmJlZQXlaItcCL0k6BBgGFgD/tLs3johtEXFX2n8WuA84CFgGXJUuuwo4Je0vA66JiOci4mFgE3Dk7sZhZmbF5UkiOyLiReB9wOcj4r8B85oZhKR+4O3AHcABEbEtnXoMOCDtHwQ8UvNjj6YyMzMrSZ4k8oKkDwBnANenslc1KwBJryVr7ZwXEc/UnouIAKLAd66QNCppdPv27U2K1MzMJsuTRP4IOBoYioiHJS0CvtaMm0t6FVkCGYmIb6bixyXNS+fnAU+k8q1kj9ImzE9lrxARwxExEBEDc+fObUaoZmY2hWmTSERsBP4CmHh/8XBEXLq7N5Yk4CvAfRHxuZpTq8laPaTtdTXlyyXtmRLZYuDO3Y3DzMyKy9M76/eAdcAN6fhwSaubcO93An8IHCtpXfqcCFwC/LakB4Dj0zERsQFYBWxMsZwVES81IQ4zMysoz3oinyDrBXULQESsk/Sm3b1xRPwboDqnj6vzM0OA1z8zM6uIXC/WI+LpSWU7WhGMmZl1ljwtkQ2SPgjMlrQYOAf4UWvDMjOzTpCnJXI22Sjx54CvA88A57UyKDMz6wzTtkQiYhxYmT5mZmYvq5tEJP0rDQb6RcTJLYnIzMw6RqOWyF+3LQozM+tIdZNIRNzazkDMzKzzNHqctSoi3i9pPVM81oqIt7Y0MjMzq7xGj7POTduT2hGImZl1nrpdfGumY/9oRIzVfoCPtic8MzOrsjzjRH57irLfbXYgZmbWeRq9E/lTshbHmyTdU3Nqb+CHrQ7MzMyqr9E7kX8Cvgt8hrTOefJsRPyipVGZmVlHaNTF92ngaeADkmaTLVO7B/BaSa+NiC1titHMzCpq2mlPJH2MbDr4x9k5e28A7uJrZtbj8sziex7w5oj4eauDMTOzzpKnd9YjZI+1zMzMdpGnJfIQcIukb5NNBw/ApHXRzcysB+VJIlvS59XpY2ZmBuRbT+ST7QjEzMw6T57eWXOB88lWN9xrojwijm1hXGZm1gHyvFgfAX4KLAI+CWwGftzCmMzMrEPkSSL7RcRXgBci4taIOBNwK8TMzHK9WH8hbbdJei/wf4B9WxeSmZl1ijxJ5NOSXg98HPg88Drgv7Y0KjMz6wh5emddn3afBt7T2nDMzKyT5Omd9Q9MvTzumS2JyMzMOkaeF+vXA99OnzVkj7P+o5VBmVl1jYxAfz/MmpVtR0bKjsjKlOdx1rW1x5K+DvxbyyIys8oaGYEVK2B8PDseG8uOAQYHy4vLypOnJTLZYuCNzQ7EzKpv5cqdCWTC+HhWbr0pzzuRZ8neiShtHwP+osVxmVkFbamzFF29cut+eR5n7d2OQMys+vr6skdYU5Vbb2r4OEvSayT9saTPpc8HJXkmX7MeNTQEc+bsWjZnTlZuvaluEpH0n4CNwG+QzZe1GTgB+KGkfSR9uh0Bmll1DA7C8DAsXAhSth0e9kv1XtbocdblwIqIuKm2UNLxwL3AhlYGZmatMzKSvQzfsiV7FDU0lD8RDA46adhOjZLIvMkJBCAibpb0AvC+1oVlZq3ibrrWTI3eicyStOfkQkl7kc3oOz7Fz5hZxbmbrjVToyRyNXCtpIUTBZL6gVXA11obVn2Slkq6X9ImSReUFYdZp3I3XWumukkkIj4N3ADcJulJSU8CtwI3RcTF7QqwlqTZwBeA3wWWAB+QtKSMWMw6Vb3uuO6ma0U07OIbEVdERB/ZqoaLImJhRHy+PaFN6UhgU0Q8FBHPA9cAy0qMx6zjuJuuNVOuaU8i4tmIeLbVweRwEPBIzfGjqWwXklZIGpU0un379rYFZ9ZuRSZDdDdda6Yic2dVXkQMR8RARAzMnTu37HDMcplpQpjoZTU2BhE7e1nlTSSbN8OOHdnWCcSK6rQkshVYUHM8P5WZdbQiCcG9rKwKpk0ikuZI+h+SvpSOF0s6qfWhTenHwGJJi9L0K8uB1SXFYtY0RRKCe1lZFeRpifwD8BxwdDreCpQy5UlEvAh8DLgRuA9YFREeOW8dr0hCcC8rq4I8SeTgiPgs8AJAGmSolkbVQER8JyJ+JSIOjgj3J7GuUCQhuJeVVUGeJPK8pNeQ1lmXdDBZy8TMGpjJi/IiCcG9rKwKpl1PBPgrskGHCySNAO8EPtTKoMw63Uznp5oom+mkiJ4M0cqmiJj+Imk/4Ciyx1i3R8STrQ6sWQYGBmJ0dLTsMKzH9PdPvXjTwoVZl1qzqpO0NiIGpruubktE0jsmFW1L2z5JfRFx1+4EaNbN3HPKekWjx1l/0+BcAMc2ORazruFlZK1X1E0iEfGedgZi1k2GhnZ9JwLuOWXdKc9gw70k/Zmkb0q6VtJ5aU0RM6vDPaesV0z7Yl3SKuBZ4B9T0QeBfSLi1BbH1hR+sW5mNnO7/WK9xmERUbtmxw8kbSwempmZdYs8gw3vknTUxIGkXwf8T3szM8vVEjkC+JGkic6JfcD9ktYDERFvbVl0ZmZWaXlaIkvJVjb8rfRZlMpOAn6vdaGZVUeRxZ/MesG0LZGIGJP0BrJ1PPaoKfdgQ+sJM53CxKyX5OniezFwD3A52QDEvwH+usVxmbVEkRaFF38yqy/PO5H3k00H/3yrgzFrpaItCk9hYlZfnnci9wL7tDoQs1Yr2qLw4k9m9eVJIp8BfiLpRkmrJz6tDsys2Yq2KLz4k1l9eR5nXQVcCqwHdrQ2HLPWKTopYtG1Psx6QZ4kMh4Rl7c8ErMW251JEb34k9nU8iSR2yR9BlhNzbK47uJrncYtCrPmy5NE3p62R9WUeT0R60huUZg1V57Bhl5XxMzMppSnJYKk9wKHAi+vIxIRn2pVUGZm1hnyjFj/InAacDYg4FRgYYvjMjOzDpBnnMgxEXE68O8R8UngaOBXWhuWmZl1gjxJ5P+l7bikXwZeAOa1LiQzM+sUed6JXC9pH+Ay4C6ynllfamlUZmbWEaZtiUTExRHxVERcS/Yu5C0RcVHrQzOrz+t7mFVD3SQi6dckHVhzfDqwCrhY0r7tCM5sKhOz8Y6NQcTO2XidSMzar1FL5H8BzwNI+k3gEuBq4GlguPWhWbfa3VaE1/cwq45G70RmR8Qv0v5pwHB6pHWtpHWtD826UTNWCfT6HmbV0aglMlvSRJI5Dvh+zblcgxTNJmtGK8Lre5hVR6Mk8nXgVknXkXXzvQ1A0iFkj7TMZqwZrQiv72FWHXWTSEQMAR8Hvgq8KyKi5mfObn1o1o2a0YoYHIThYVi4EKRsOzzsiRXNytDwsVRE3D5F2c9aF451u91Z06OWZ+M1q4Y8I9bNmsatCLPuUkoSkXSZpJ9KukfSt9KI+IlzF0raJOl+SSfUlB8haX06d7kklRG77b7BQdi8GXbsyLZOIGadq6yWyE3AYRHxVuBnwIUAkpYAy8mmnV8KXClpdvqZvwP+BFicPkvbHbSZme2qlCQSEd+LiBfT4e3A/LS/DLgmIp6LiIeBTcCRkuYBr4uI29ML/quBU9oeeA/y9CJm1kgVxnucCXwj7R9EllQmPJrKXkj7k8uthZoxMNDMulvLWiKSbpZ07xSfZTXXrAReBJr671tJKySNShrdvn17M7+6p3h6ETObTstaIhFxfKPzkj4EnAQcVzMGZSuwoOay+alsKzsfedWW17v3MGl+r4GBgah3nTXm6UXMbDpl9c5aCpwPnBwRtf/WXQ0sl7SnpEVkL9DvjIhtwDOSjkq9sk4Hrmt74D3G04uY2XTK6p11BbA3cJOkdWkddyJiA9l08xuBG4CzIuKl9DMfBb5M9rL9QeC7bY+6x3h6ETObTikv1iPikAbnhoBX/JmKiFHgsFbGZbuaeHm+cmX2CKuvL0sgfqluZhOq0DvLKszTi5hZI572xMzMCnMSMTOzwpxEOoxHkJtZlfidSAfxCHIzqxq3RDqIR5CbWdU4iXQQjyA3s6pxEukgHkFuZlXjJNJBPILczKrGSaSDeGlZM6sa987qMB5BbmZV4paImZkV5iRiZmaFOYmYmVlhTiJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZWmJOImZkV5iRiZmaFOYmYmVlhTiJmZlaYk4iZmRXmJGJmZoU5iZiZWWFOImZmVpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZWmJOImZkV5iRiZmaFlZpEJH1cUkjav6bsQkmbJN0v6YSa8iMkrU/nLpekcqI2M7MJpSURSQuA3wG21JQtAZYDhwJLgSslzU6n/w74E2Bx+ixta8BmZvYKZbZE/hY4H4iasmXANRHxXEQ8DGwCjpQ0D3hdRNweEQFcDZzS9ojNzGwXpSQRScuArRFx96RTBwGP1Bw/msoOSvuTy2dsZAT6+2HWrGw7MlLkW8zMDGCPVn2xpJuBA6c4tRL472SPslp17xXACoC+vr6Xy0dGYMUKGB/PjsfGsmOAwcFWRWNm1r1a1hKJiOMj4rDJH+AhYBFwt6TNwHzgLkkHAluBBTVfMz+VbU37k8vr3Xs4IgYiYmDu3Lkvl69cuTOBTBgfz8rNzGzm2v44KyLWR8QbI6I/IvrJHk29IyIeA1YDyyXtKWkR2Qv0OyNiG/CMpKNSr6zTgetmeu8tW2ZWbmZmjVVqnEhEbABWARuBG4CzIuKldPqjwJfJXrY/CHx3pt9f82QrV7mZmTXWsncieaXWSO3xEDA0xXWjwGG7c6+hoV3fiQDMmZOVm5nZzFWqJdJqg4MwPAwLF4KUbYeH/VLdzKyo0lsi7TY46KRhZtYsPdUSMTOz5nISMTOzwpxEzMysMCcRMzMrzEnEzMwKUzYpbveStB0Ym1S8P/BkCeHMRNVjrHp8UP0Yqx4fVD/GqscH1Y+xXnwLI2LuFOW76PokMhVJoxExUHYcjVQ9xqrHB9WPserxQfVjrHp8UP0Ydzc+P84yM7PCnETMzKywXk0iw2UHkEPVY6x6fFD9GKseH1Q/xqrHB9WPcbfi68l3ImZm1hy92hIxM7Mm6IkkIukyST+VdI+kb0nap851SyXdL2mTpAvaHOOpkjZI2iGpbk8JSZslrZe0TtJoBeMrsw73lXSTpAfS9g11rmtrHU5XJ8pcns7fI+kdrY5phvG9W9LTqb7WSbqozfH9vaQnJN1b53yp9ZczxrLrcIGkH0jamP5/fO4U1xSrx4jo+g/Zeu57pP1LgUunuGY22WJXbwJeDdwNLGljjL8KvBm4BRhocN1mYP8S6nDa+CpQh58FLkj7F0z1v3O76zBPnQAnki2yJuAo4I421lme+N4NXN/u/+Zq7v+bwDuAe+ucL63+ZhBj2XU4j2wFWYC9gZ8167/DnmiJRMT3IuLFdHg7u67XPuFIYFNEPBQRzwPXAMvaGON9EXF/u+43UznjK7UO072uSvtXAae08d715KmTZcDVkbkd2EfSvArFV6qI+N/ALxpcUmb9AbliLFVEbIuIu9L+s8B9wEGTLitUjz2RRCY5k6mX1j0IeKTm+FFeWclVEMDNktZKWlF2MJOUXYcHRMS2tP8YcECd69pZh3nqpMx6y3vvY9Ijju9KOrQ9oeVW9n93eVWiDiX1A28H7ph0qlA9ds2iVJJuBg6c4tTKiLguXbMSeBEYaWdsE/LEmMO7ImKrpDcCN0n6afpXUFXia6lGMdYeRERIqtf1sGV12KXuAvoi4j8knQj8C7C45Jg6TSXqUNJrgWuB8yLimWZ8Z9ckkYg4vtF5SR8CTgKOi/QAcJKtwIKa4/mprGmmizHnd2xN2yckfYvscURT/gA2Ib5S61DS45LmRcS21Ax/os53tKwOp5CnTlpebw1Me+/aPzYR8R1JV0raPyKqMh9UmfWXSxXqUNKryBLISER8c4pLCtVjTzzOkrQUOB84OSLG61z2Y2CxpEWSXg0sB1a3K8Y8JP2SpL0n9sk6DEzZG6QkZdfhauCMtH8G8IrWUwl1mKdOVgOnp94xRwFP1zyWa7Vp45N0oCSl/SPJ/m78vE3x5VFm/eVSdh2me38FuC8iPlfnsmL1WFZvgXZ+gE1kz/rWpc8XU/kvA9+pue5Esl4LD5I9wmlnjO8jewb5HPA4cOPkGMl60NydPhvaGWOe+CpQh/sBa4AHgJuBfatQh1PVCfAR4CNpX8AX0vn1NOidV1J8H0t1dTdZx5Rj2hzf14FtwAvpv8EPV6n+csZYdh2+i+xd4D01fwdPbEY9esS6mZkV1hOPs8zMrDWcRMzMrDAnETMzK8xJxMzMCnMSMTOzwpxErCNI2q9mBtTHJG1N+09J2tjmWA5Po44njk9WwRmLlc0ovP8U5a+XdHWaUfVBSSOqMyvx7mj0u0j6hKQ/b/Y9rbs4iVhHiIifR8ThEXE48EXgb9P+4cCOZt9PUqPZHA4n62M/EdvqiLikySF8BXgoIg6JiIPJxjp9tcn3gPb8LtbFnESsG8yW9KW0TsL3JL0GQNLBkm5IEy3eJuktqbxf0vfTZHhrJPWl8q9K+qKkO4DPptHtfy/pTkk/kbQsjer+FHBaagmdJulDkq5I33GAsjVr7k6fY1L5v6Q4NmiaSR8lHQIcAVxcU/wp4G2S3qxsbYrra66/Ik3rg6SLJP1Y0r2ShmtGSd8i6dL0u/xM0m9M97tMiqleXZ6a7nW3JM8/1oOcRKwbLAa+EBGHAk8Bv5/Kh4GzI+II4M+BK1P554GrIuKtZJNxXl7zXfPJRhP/Gdmkjt+PiCOB9wCXAa8CLgK+kVpG35gUy+XArRHxNrL1JTak8jNTHAPAOZL2a/D7LAHWRcRLEwVp/ydk67o0ckVE/FpEHAa8hmy+uAl7pN/lPOCvIpv6vdHvUqteXV4EnJB+35Onic26UNdMwGg97eGIWJf21wL9ymYrPQb45/SPcYA90/Zo4D+n/a+RLWY14Z9r/nj/DnByzXuBvYC+aWI5FjgdXv7D/3QqP0fS+9L+ArLE14q5k94j6XxgDrAvWRL713RuYtK9tUB/3i+cpi5/CHxV0qqa77ce4iRi3eC5mv2XyP4FPgt4Kr03mYn/W7Mv4Pdj0mJckn59Jl8o6d3A8cDRETEu6RayhFTPRuBwSbMiYkf6jlnA20hTirPrU4S90jV7kbUQBiLiEUmfmHSfiXp6iZn9f79uXUbER1J9vBdYK+mIiKjS5IzWYn6cZV0psqm3H5Z0Kry8fvTb0ukfkc1WCzAI3Fbna24Ezq55r/D2VP4s2RKjU1kD/Gm6frak1wOvB/49JZC3kC092ij2TWSPrv6ypvgvgTURsQUYA5ZI2lPSPsBx6ZqJhPFkaj38l0b3yfG7TMRTty4lHRwRd0TERcB2dp1K3HqAk4h1s0Hgw5ImZuydWPb1bOCPJN0D/CFwbp2fv5jsHcg9kjaw80X3D8j+iK+TdNqknzmX7JHSerLHRkuAG4A9JN0HXEI2i+t0ziSbov1BSdvJEs9HACLiEWAV2RT2q8gSDhHxFPClVH4j2TTv02n0u9SqV5eXSVov6V6y5Hx3jntaF/EsvmYVJ+nNwLeBcyLiO2XHY1bLScTMzArz4ywzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwv4/zHGcRSZeLawAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import scipy.stats as stats\n", "probplot = sm.ProbPlot(res, stats.t, distargs=(4,))\n", "fig = probplot.qqplot()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEKCAYAAADTgGjXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAGoBJREFUeJzt3X2UXXV97/H3J8GCuaIIpJibZDIBojZwFWVKwYdeQSpBkeBtkei0YPE2y0p5qPRSuOmlKs4Stdq7ENE7aAvYUUwXWlJUKESlVBfiBAMhoZFAMoE0QNDycO90hYd87x/7N2RnmHPmzM45Z+9z5vNa66y992+fh+/Z6yTf+e3fkyICMzOzImaUHYCZmXUuJxEzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMyssH3KDqDVDj744Ojt7S07DDOzjrJmzZonImL2ZM/r+iTS29vL8PBw2WGYmXUUSSONPM+3s8zMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRM7NpZmgIenthxoxsOzRU/L26vouvmZntNjQEy5fD6Gh2PDKSHQP090/9/VwTMTObRlas2J1AxoyOZuVFOImYmU0jW7dOrXwyTiJmZtNIT8/UyifjJGJmNo0MDMCsWXuWzZqVlRfhJGJmNo3098PgICxYAFK2HRws1qgO7p1lZjbt9PcXTxrjuSZiZmaFOYmYmVlhTiJmZlaYk4iZmRVWehKRNFPSzyXdlI4PlHSrpAfS9tW5514iaZOkjZJOKi9qMzODCiQR4Hzg/tzxxcDqiFgErE7HSFoMLAOOAJYAV0ma2eZYzcwsp9QkImke8B7gq7nipcC1af9a4LRc+fURsTMiNgObgGPaFauZmb1U2TWR/w1cBOzKlR0SEdvT/qPAIWl/LvBw7nmPpDIzMytJaUlE0inA4xGxptZzIiKAKPDeyyUNSxresWPH3oRpZmZ1lFkTeStwqqQtwPXACZL+DnhM0hyAtH08PX8bMD/3+nmp7CUiYjAi+iKib/bs2a2K38xs2istiUTEJRExLyJ6yRrMfxARvw+sAs5KTzsLuDHtrwKWSdpX0kJgEXBXm8M2M7OcKs6ddTmwUtKHgRHg/QARsV7SSmAD8DxwTkS8UF6YZmamrNmhe/X19cXw8HDZYZiZdRRJayKib7Lnld07y8zMOpiTiJmZFeYkYmZmhTmJmJlZYU4iZmZWmJOImdkUDQ1Bby/MmJFth4bKjqg8VRwnYmZWWUNDsHw5jI5mxyMj2TE0b93yTuKaiJnZFKxYsTuBjBkdzcqnIycRM7Mp2Lp1auXdzknEzGwKenqmVt7tnETMzKZgYABmzdqzbNasrHw6chIxM5uC/n4YHIQFC0DKtoOD07NRHdw7y8xsyvr7p2/SGM81ETPreB63UR7XRMyso3ncRrlcEzGzjuZxG+VyEjGzjuZxG+VyEjGzjuZxG+VyEjGzjuZxG+VyEjGzjuZxG+Vy7ywz63get1Ee10TMzKwwJxEzazsPDuwevp1lZm3lwYHdxTURM2srDw7sLk4iZtZWHhzYXZxEzKytPDiwuziJmFlbeXBgd3ESMbO28uDA7uLeWWbWdh4c2D1cEzEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRMzMrrLQkImm+pB9K2iBpvaTzU/mBkm6V9EDavjr3mkskbZK0UdJJZcVuZmaZSZOIpMMk7Zv23yHpPEkHNOGznwcujIjFwLHAOZIWAxcDqyNiEbA6HZPOLQOOAJYAV0ma2YQ4zMysoEZqIjcAL0g6HBgE5gPf2NsPjojtEXF32n8GuB+YCywFrk1PuxY4Le0vBa6PiJ0RsRnYBByzt3GYmVlxjSSRXRHxPPA+4IsR8T+AOc0MQlIv8Cbgp8AhEbE9nXoUOCTtzwUezr3skVRmZmYlaSSJPCfpA8BZwE2p7GXNCkDSK8hqOxdExNP5cxERQBR4z+WShiUN79ixo0mRmpnZeI0kkT8EjgMGImKzpIXA15vx4ZJeRpZAhiLi26n4MUlz0vk5wOOpfBvZrbQx81LZS0TEYET0RUTf7NmzmxGqmZlNYNIkEhEbgD8HxtovNkfEZ/b2gyUJ+Bpwf0R8IXdqFVmth7S9MVe+TNK+KZEtAu7a2zjMzKy4RnpnvRdYC9ycjo+StKoJn/1W4A+AEyStTY93A5cDvyPpAeDEdExErAdWAhtSLOdExAtNiMPMzApqZD2Rj5P1gvoRQESslXTo3n5wRPwLoBqn31njNQOA1z8zM6uIhhrWI+KpcWW7WhGMmZl1lkZqIuslfRCYKWkRcB7wk9aGZWZmnaCRmsi5ZKPEdwLfBJ4GLmhlUGZm1hkmrYlExCiwIj3MzMxeVDOJSPpH6gz0i4hTWxKRmZl1jHo1kb9qWxRmZtaRaiaRiLi9nYGYmVnnqXc7a2VEvF/SOia4rRURb2hpZGZmVnn1bmedn7antCMQMzPrPDW7+OamY/9oRIzkH8BH2xOemZlVWSPjRH5ngrKTmx2ImZl1nnptIn9MVuM4VNK9uVP7Az9udWBmZlZ99dpEvgF8H/g0aZ3z5JmI+FVLozIzs45Qr4vvU8BTwAckzSRbpnYf4BWSXhERW9sUo5mZVdSk055I+hOy6eAfY/fsvQG4i6+Z2TTXyCy+FwCvi4hftjoYMzPrLI30znqY7LaWmZnZHhqpiTwE/EjSd8mmgwdg3LroZmY2DTWSRLamx6+lh5mZGdDYeiKfaEcgZmbWeRrpnTUbuIhsdcP9xsoj4oQWxmVmZh2gkYb1IeBfgYXAJ4AtwM9aGJOZmXWIRpLIQRHxNeC5iLg9Is4GXAsxM7OGGtafS9vtkt4D/BtwYOtCMjOzTtFIEvmUpFcBFwJfBF4J/GlLozIzs47QSO+sm9LuU8DxrQ3HzMw6SSO9s/6WiZfHPbslEZmZWcdopGH9JuC76bGa7HbW/21lUGZWXUND0NsLM2Zk26GhsiOyMjVyO+uG/LGkbwL/0rKIzKyyhoZg+XIYHc2OR0ayY4D+/vLisvI0UhMZbxHw680OxMyqb8WK3QlkzOhoVm7TUyNtIs+QtYkobR8F/rzFcZlZBW2tsRRdrXLrfo3cztq/HYGYWfX19GS3sCYqt+mp7u0sSS+X9N8lfSE9PijJM/maTVMDAzBr1p5ls2Zl5TY91Uwikv4LsAF4O9l8WVuAk4AfSzpA0qfaEaCZVUd/PwwOwoIFIGXbwUE3qk9n9W5nXQEsj4hb84WSTgTuA9a3MjAza52hoawxfOvW7FbUwEDjiaC/30nDdquXROaMTyAAEXGbpOeA97UuLDNrFXfTtWaq1yYyQ9K+4wsl7Uc2o+/oBK8xs4pzN11rpnpJ5DrgBkkLxgok9QIrga+3NqzaJC2RtFHSJkkXlxWHWadyN11rpppJJCI+BdwM3CHpCUlPALcDt0bEZe0KME/STOBLwMnAYuADkhaXEYtZp6rVHdfddK2Iul18I+LKiOghW9VwYUQsiIgvtie0CR0DbIqIhyLiWeB6YGmJ8Zh1HHfTtWZqaNqTiHgmIp5pdTANmAs8nDt+JJXtQdJyScOShnfs2NG24MzarchkiO6ma81UZO6syouIwYjoi4i+2bNnlx2OWUOmmhDGelmNjEDE7l5WjSaSLVtg165s6wRiRXVaEtkGzM8dz0tlZh2tSEJwLyurgkmTiKRZkv6XpKvT8SJJp7Q+tAn9DFgkaWGafmUZsKqkWMyapkhCcC8rq4JGaiJ/C+wEjkvH24BSpjyJiOeBPwFuAe4HVkaER85bxyuSENzLyqqgkSRyWER8FngOIA0yVEujqiMivhcRr42IwyLC/UmsKxRJCO5lZVXQSBJ5VtLLSeusSzqMrGZiZnVMpaG8SEJwLyurgknXEwH+kmzQ4XxJQ8BbgQ+1MiizTjfV+anGyqY6KaInQ7SyKSImf5J0EHAs2W2sOyPiiVYH1ix9fX0xPDxcdhg2zfT2Trx404IFWZdas6qTtCYi+iZ7Xs2aiKQ3jyvanrY9knoi4u69CdCsm7nnlE0X9W5nfb7OuQBOaHIsZl3Dy8jadFEziUTE8e0MxKybDAzs2SYC7jll3amRwYb7SfqYpG9LukHSBWlNETOrwT2nbLqYtGFd0krgGeDvUtEHgQMi4vQWx9YUblg3M5u6vW5YzzkyIvJrdvxQ0obioZmZWbdoZLDh3ZKOHTuQ9FuA/7Q3M7OGaiJHAz+RNNY5sQfYKGkdEBHxhpZFZ2ZmldZITWQJ2cqG/zU9FqayU4D3ti40s+oosviT2XQwaU0kIkYkvZpsHY99cuUebGjTwlSnMDGbThrp4nsZcC9wBdkAxM8Df9XiuMxaokiNwos/mdXWSJvI+8mmg3+21cGYtVLRGoWnMDGrrZE2kfuAA1odiFmrFa1RePEns9oaSSKfBn4u6RZJq8YerQ7MrNmK1ii8+JNZbY3czroW+AywDtjV2nDMWqfopIhF1/owmw4aSSKjEXFFyyMxa7G9mRTRiz+ZTayRJHKHpE8Dq8gti+suvtZpXKMwa75Gksib0vbYXJnXE7GO5BqFWXM1MtjQ64qYmdmEGqmJIOk9wBHAi+uIRMQnWxWUmZl1hkZGrH8FOAM4FxBwOrCgxXGZmVkHaGScyFsi4kzg3yPiE8BxwGtbG5aZmXWCRpLIf6TtqKT/DDwHzGldSGZm1ikaaRO5SdIBwOeAu8l6Zl3d0qjMzKwjTFoTiYjLIuLJiLiBrC3k9RFxaetDM6vN63uYVUPNJCLpNyW9Jnd8JrASuEzSge0IzmwiY7PxjoxAxO7ZeJ1IzNqvXk3k/wDPAkj6beBy4DrgKWCw9aFZt9rbWoTX9zCrjnptIjMj4ldp/wxgMN3SukHS2taHZt2oGasEen0Ps+qoVxOZKWksybwT+EHuXEODFM3Ga0Ytwut7mFVHvSTyTeB2STeSdfO9A0DS4WS3tMymrBm1CK/vYVYdNZNIRAwAFwLXAG+LiMi95tzWh2bdqBm1iP5+GByEBQtAyraDg55Y0awMdW9LRcSdE5T9onXhWLfbmzU98jwbr1k1NDJi3axpXIsw6y6lJBFJn5P0r5LulfSdNCJ+7NwlkjZJ2ijppFz50ZLWpXNXSFIZsdve6++HLVtg165s6wRi1rnKqoncChwZEW8AfgFcAiBpMbCMbNr5JcBVkmam13wZ+CNgUXosaXfQZma2p1KSSET8U0Q8nw7vBOal/aXA9RGxMyI2A5uAYyTNAV4ZEXemBv7rgNPaHvg05OlFzKyeKoz3OBv4VtqfS5ZUxjySyp5L++PLrYWaMTDQzLpby2oikm6TdN8Ej6W556wAngea+vetpOWShiUN79ixo5lvPa14ehEzm0zLaiIRcWK985I+BJwCvDM3BmUbMD/3tHmpbBu7b3nly2t99iBpfq++vr6o9Tyrz9OLmNlkyuqdtQS4CDg1IvJ/664ClknaV9JCsgb0uyJiO/C0pGNTr6wzgRvbHvg04+lFzGwyZfXOuhLYH7hV0tq0jjsRsZ5suvkNwM3AORHxQnrNR4GvkjW2Pwh8v+1RTzOeXsTMJlNKw3pEHF7n3ADwkv+mImIYOLKVcdmexhrPV6zIbmH19GQJxI3qZjamCr2zrMI8vYiZ1eNpT8zMrDAnETMzK8xJpMN4BLmZVYnbRDqIR5CbWdW4JtJBPILczKrGSaSDeAS5mVWNk0gH8QhyM6saJ5EO4hHkZlY1TiIdxEvLmlnVuHdWh/EIcjOrEtdEzMysMCcRMzMrzEnEzMwKcxIxM7PCnETMzKwwJxEzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRMzMrzEnEzMwKcxIxM7PCnETMzKwwJxEzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwpxEzMysMCcRMzMrrNQkIulCSSHp4FzZJZI2Sdoo6aRc+dGS1qVzV0hSOVGbmdmY0pKIpPnAu4CtubLFwDLgCGAJcJWkmen0l4E/Ahalx5K2BmxmZi9RZk3kr4GLgMiVLQWuj4idEbEZ2AQcI2kO8MqIuDMiArgOOK3tEZuZ2R5KSSKSlgLbIuKecafmAg/njh9JZXPT/vjyKRsagt5emDEj2w4NFXkXMzMD2KdVbyzpNuA1E5xaAfxPsltZrfrs5cBygJ6enhfLh4Zg+XIYHc2OR0ayY4D+/lZFY2bWvVpWE4mIEyPiyPEP4CFgIXCPpC3APOBuSa8BtgHzc28zL5VtS/vjy2t99mBE9EVE3+zZs18sX7FidwIZMzqalZuZ2dS1/XZWRKyLiF+PiN6I6CW7NfXmiHgUWAUsk7SvpIVkDeh3RcR24GlJx6ZeWWcCN071s7dunVq5mZnVV6lxIhGxHlgJbABuBs6JiBfS6Y8CXyVrbH8Q+P5U3z93Z6uhcjMzq69lbSKNSrWR/PEAMDDB84aBI/fmswYG9mwTAZg1Kys3M7Opq1RNpNX6+2FwEBYsACnbDg66Ud3MrKjSayLt1t/vpGFm1izTqiZiZmbN5SRiZmaFOYmYmVlhTiJmZlaYk4iZmRWmbFLc7iVpBzDS5Lc9GHiiye/ZbFWPserxQfVjrHp8UP0Yqx4flBfjgoiYPdmTuj6JtIKk4YjoKzuOeqoeY9Xjg+rHWPX4oPoxVj0+qH6Mvp1lZmaFOYmYmVlhTiLFDJYdQAOqHmPV44Pqx1j1+KD6MVY9Pqh4jG4TMTOzwlwTMTOzwpxEpkDS6ZLWS9olqS9X3ivpPyStTY+vVCm+dO4SSZskbZR0UhnxjSfp45K25a7bu8uOCUDSknSdNkm6uOx4JiJpi6R16boNlx0PgKS/kfS4pPtyZQdKulXSA2n76orFV5nfoKT5kn4oaUP6d3x+Kq/MNZyIk8jU3Af8N+CfJzj3YEQclR4faXNcYyaMT9JiYBlwBLAEuErSzPaHN6G/zl2375UdTLouXwJOBhYDH0jXr4qOT9etKt0/ryH7feVdDKyOiEXA6nRclmt4aXxQnd/g88CFEbEYOBY4J/32qnQNX8JJZAoi4v6I2Fh2HLXUiW8pcH1E7IyIzWSrQx7T3ug6xjHApoh4KCKeBa4nu342iYj4Z+BX44qXAtem/WuB09oaVE6N+CojIrZHxN1p/xngfmAuFbqGE3ESaZ6FqTp8u6S3lx3MOHOBh3PHj6SyKjhX0r3pVkMVqulVvlZ5AdwmaY2k5WUHU8chEbE97T8KHFJmMDVU7TeIpF7gTcBPqfg1dBIZR9Jtku6b4FHvr9HtQE9EHAV8DPiGpFdWKL7STBLvl4FDgaPIruHnSw22s7wt/d5OJrvt8dtlBzSZyLqCVq07aOV+g5JeAdwAXBART+fPVfEaTruVDScTEScWeM1OYGfaXyPpQeC1QNMbPIvEB2wD5ueO56Wylms0XklXAze1OJxGlHatpiIitqXt45K+Q3YbbqK2urI9JmlORGyXNAd4vOyA8iLisbH9KvwGJb2MLIEMRcS3U3Glr6FrIk0gafZYQ7WkQ4FFwEPlRrWHVcAySftKWkgW310lx0T6BzHmfWQdA8r2M2CRpIWSfo2sQ8KqkmPag6T/JGn/sX3gXVTj2k1kFXBW2j8LuLHEWF6iSr9BSQK+BtwfEV/Inar0NSQi/GjwQfYje4Ss1vEYcEsq/11gPbAWuBt4b5XiS+dWAA8CG4GTy76WKaavA+uAe8n+ocwpO6YU17uBX6TrtaLseCaI71DgnvRYX5UYgW+S3RJ6Lv0OPwwcRNaj6AHgNuDAisVXmd8g8DayW1X3pv9L1qbfYmWu4UQPj1g3M7PCfDvLzMwKcxIxM7PCnETMzKwwJxEzMyvMScTMzApzErGOIOmg3Eyrj+ZmXn1S0oY2x3JUfrZXSacWne03zcZ78ATlr5J0XZpJ+EFJQ62YkqPed0kz3P5Zsz/TuouTiHWEiPhlpJlWga+QZl4lm65iV7M/T1K92RyOIuu/Pxbbqoi4vMkhfA14KCIOj4jDyCbNvKbJnwHt+S7WxZxErBvMlHR1WoPhnyS9HEDSYZJuTpMU3iHp9am8V9IP0qR7qyX1pPJrJH1F0k+Bz6aR4X8j6S5JP5e0NI1i/yRwRqoJnSHpQ5KuTO9xiKTvSLonPd6Syv8hxbF+sgkTJR0OHA1cliv+JPBGSa+T9A5JN+Wef6WkD6X9SyX9LM1PNphGQSPpR5I+k77LLyS9fbLvMi6mWtfy9PRZ90iq4rQr1mJOItYNFgFfiogjgCfJZhCAbG3qcyPiaODPgKtS+ReBayPiDcAQcEXuveYBb4mIj5GN8v9BRBwDHA98DngZcCnwrVQz+ta4WK4Abo+INwJvJhtRDnB2iqMPOE/SQXW+z2JgbUS8MFaQ9n8O/MYk1+LKiPjNiDgSeDlwSu7cPum7XAD8ZWRT3df7Lnm1ruWlwEnp+546SWzWhTwBo3WDzRGxNu2vAXrTTKhvAf4+/TEOsG/aHke2eBdk0158Nvdef5/7z/tdwKm5doH9gJ5JYjkBOBNe/I//qVR+nqT3pf35ZInvl419vSk5XtJFwCzgQLIk9o/p3NiEfmuA3kbfcJJr+WPgGkkrc+9v04iTiHWDnbn9F8j+Ap8BPJnaTabi/+X2BfxujFvoS9JvTeUNJb0DOBE4LiJGJf2ILCHVsgE4StKMiNiV3mMG8Eayudl62PMuwn7pOfuR1RD6IuJhSR8f9zlj1+kFpvZvv+a1jIiPpOvxHmCNpKMjohXJ0SrKt7OsK0W2DsNmSadDNkOqpDem0z8hm50XoB+4o8bb3EK2YNFYu8KbUvkzwP41XrMa+OP0/JmSXgW8Cvj3lEBeT7b0ab3YN5HduvqLXPFfkC2RuhUYARYrm5X5AOCd6TljCeOJVHv4vXqf08B3GYun5rWUdFhE/DQiLgV2sOc0+jYNOIlYN+sHPixpbLbbsYW7zgX+UNK9wB8A59d4/WVkbSD3SlrP7obuH5L9J75W0hnjXnM+2S2ldWS3jRYDNwP7SLofuBy4s4HYzyabkv5BSTvIEs9HACLiYWAl2bTlK8kSDhHxJHB1Kr+FbFr7ydT7Lnm1ruXnJK2TdB9Zcr6ngc+0LuJZfM0qTtLrgO8C50XE98qOxyzPScTMzArz7SwzMyvMScTMzApzEjEzs8KcRMzMrDAnETMzK8xJxMzMCnMSMTOzwv4/lEXAaP2HfqYAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "probplot = sm.ProbPlot(res, stats.t, distargs=(4,), loc=3, scale=10)\n", "fig = probplot.qqplot()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEKCAYAAAAfGVI8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYlOXVx/HvEWxgQ0RjVHbV2AAFZVE0ryJijIkFTUzUrCUBJIgVSyxgixIRNNgQXIqIIAYVA2JBxCCKgiy92aKAiMjakLhKPe8f94Ou65bZMvPM7Pw+1zXXzDw7M8+ZhZ0zdzu3uTsiIpK9too7ABERiZcSgYhIllMiEBHJckoEIiJZTolARCTLKRGIiGQ5JQIRkSynRCAikuWUCEREslz9uANIxG677ea5ublxhyEikv7Wr4elS2HtWmbBZ+7epLKnZEQiyM3NpbCwMO4wRETS16ZNMGAA3HgjmMGDD2KXXroskaeqa0hEJNMtWQLHHQdXXAHHHgsLF8IllyT8dCUCEZFMtWED9O4NrVrB22/DiBHw/POQk1Oll0laIjCzYWa22swWljjWysymm9lcMys0syOTdX4RkTpt9mxo0wZ69YKOHWHxYjj//NAtVEXJbBEMB04udawvcJu7twJuju6LiEiivv0Wrr8ejjwSPv0Uxo6FMWNgjz2q/ZJJGyx296lmllv6MLBTdHtnYGWyzi8iUudMnQpdusB770HnztCvHzRqVOOXTfUYwZVAPzP7CLgbuKG8B5pZ16j7qLCoqChlAYqIJGLUKMjNha22CtejRiXxZF9/HQZ/27UL4wKTJsGQIbWSBCD1ieBioIe77wP0AIaW90B3L3D3PHfPa9Kk0mmwIiIpM2oUdO0Ky5aBe7ju2jVJyeCFF6BFCxg4EK68MswIOvHEWj1FqhPBhcDY6PaTgAaLRSTj9OwJxcU/PlZcHI7Xms8/hwsugN/+FnbYAaZNg/79oWHDWjxJkOpEsBJoF90+AXgvxecXEamx5curdrxK3MPg7yGHwOjRcNNNMGcOHH10Lbx42ZI2WGxmo4Hjgd3MbAVwC3ARcJ+Z1Qe+A7om6/wiIsnStGnoDirreI2sXAndu8O4cdC6dRgLaNmyhi9auWTOGjq3nB+1TtY5RURSoXfvMCZQsnuoQYNwvFrcYdgwuPpqWLcO+vaFHj2gfmqqAGllsYhIFeXnQ0FBWMBrFq4LCsLxKvvggzD426VL+PY/fz5ce23KkgBkSNE5EZF0k59fzQ/+LTZtggceCCPM9eqFWUFdu4b5qCmmRCAikmqLFoUFYTNmhFlBgwbBPvvEFo66hkREUmX9erj9djj8cHj/fRg5EiZMiDUJgFoEIiKpMXNmaAUsWADnnAP33Qe77x53VIBaBCIiyVVcHAZ/27YNi8TGjQvrA9IkCYBaBCIiyTNlClx0UegGuuiiUCRu553jjuon1CIQEalta9ZAt27Qvj1s3gyTJ4f5pWmYBECJQESkdj33HDRvDoMHw1VXhTGBE06IO6oKKRGIiNSGoqKwsODUU2GXXeCNN+Cee8KS4zSnRCAiUhPuYfC3WTN48km49dawjeRRR8UdWcI0WCwiUl0rVsDFF4e1AEceCUOHhr0DMoxaBCIiVbV5cxj8bd48DATfc0/oCsrAJABqEYiIVM2WqaBTpoRZQYMHw/77xx1VjahFICKSiE2bwjf/ww4LYwAFBaE1kOFJANQiEBGp3MKF0KlTKBNx2mmhUuhee8UdVa1Ri0BEpDzr1oVZQEccAR9+GGYHjRtXp5IAqEUgIlK2GTNCkbhFi8L6gHvvhd12izuqpFCLQESkpG++CSuCjz46lIqYMCGUi66jSQDUIhAR+cErr4QZQR98EGoF3XUX7LRT3FElnVoEIpK2Ro2C3Nywe2NubrifFF99FRJAhw7hZFOmhAHhLEgCoBaBiKSpUaPCFr7FxeH+smXhPtRwr+DSxo8Pq4NXrQr7Btx6a0bUB6pNSWsRmNkwM1ttZgtLHb/MzN42s0Vm1jdZ5xeRzNaz5w9JYIvi4nC8VqxeHXYK69gRGjcOg8N9+2ZdEoDkdg0NB04uecDM2gMdgZbu3hy4O4nnF5EMtnx51Y4nzD0M/h5yCDzzTNhDuLAQ8vJq+MKZK2mJwN2nAl+UOnwx0Mfd10WPWZ2s84tIZmvatGrHE/LRR6FM9Pnnw4EHwpw50KsXbLNNDV4086V6sPhA4Fgzm2Fmr5pZmxSfX0QyRO/eP+2ladAgHK+yzZvD4G/z5mEg+N574fXXQ+loSXkiqA/sCrQFrgXGmJmV9UAz62pmhWZWWFRUlMoYRSQN5OeHcj45OWAWrgsKqjFQ/N57oThc9+6hVPSCBXDFFVCvXlLizkSpTgQrgLEevAVsBspcpeHuBe6e5+55TZo0SWmQIpIe8vNh6dLwhX7p0iomgY0bw+DvYYfBvHlhr4BJk2C//ZIUbeZKdSL4N9AewMwOBLYBPktxDCJSC1I2x7865s2Dtm3huuvg5JNh8eJQNK7sDoisl8zpo6OBN4GDzGyFmXUGhgH7RVNKnwAudHdPVgwikhxb5vgvWxYm4WyZ4x97Mli3Dm66KcwA+ugjGDMGxo6Fn/885sDSm2XC53BeXp4XFhbGHYaIRHJzw4d/aTk5oQsnFm++GYrELVkCF1wA//xnWB+QxcxslrtXOi9WJSZEpMqSNse/Ov73P7jySvjlL8Pt55+HRx/N+iRQFUoEIlJlSZnjXx2TJsGhh8J994VZQYsWwW9+k+IgMp8SgYhUWa3O8a+OL78M3UAnnRQWg02dCg8+CDvumKIA6hYlAhGpslqb418dzzwTFoI9+ihcf32YIXTssSk4cd2l6qMiUi35+Sn64N9i1Sq47DJ46ilo1Qqeey5sISk1phaBiKQ3dxgxIrQCxo8P/U9vvaUkUIuUCESySFovAivLsmVh8PfCC0O10Hnz4MYbYeut446sTlEiEMkSabsIrCybN8OAAdCiRSgOd//98NprcPDBcUdWJykRiGSJpG/0UlveeQfatYNLL4VjjoGFC8PYwFb6uEoW/WZFskRaLQIry4YN0KcPtGwZ1gMMHw4vvhj6sCSplAhEskTaLAIry5w5cNRRcMMNYeOYxYvDuICKxKWEEoFIloh9EVhZvvsuDP62aQMrV4apoU89BT/7WYxBZR8lApEsEesisLJMmxbWA9x5Z9g6cvFi+P3vYwomu2lBmUgWSfkisLKsXRtaAQMGhH6piRNDqQiJjVoEIpI6EyeGKaEDBoSZQAsXKgmkASUCEUm+L76AP/857BbWoEFYE3DffbDDDnFHJigRiEiyPfVUWBU8cmRYtDBnTtg7QNKGxghEJDk++SQsChs7Fg4/PHQLtWoVd1RSBrUIRKR2ucMjj4Qicc89FxaJvfWWkkAaU4tARGrP0qWhgNGkSWGPgMGD4aCD4o5KKqEWgYjU3KZNoTBcixZhE/kBA2DKFCWBDKEWgYjUzJIlYdvIN98Ms4IefjhN6lZIoiptEZjZ/ma2bXT7eDO73Mx2SeB5w8xstZktLONnV5uZm9lu1QtbRGK3YUOoT9GqVagYOmIEPP+8kkAGSqRr6Glgk5n9AigA9gEeT+B5w4GTSx80s32Ak4B0qXkoIlU1axbk5UGvXnDGGaE8xPnnq0hchkokEWx2943AmcAD7n4tsGdlT3L3qcAXZfyoP/A3wKsSqIikgW+/DRvGH3UUrF4dNpL/179gjz3ijkxqIJExgg1mdi5wIXBadKxa+8SZWUfgY3efZ/rmIJJZpk6FLl3gvffCmEC/ftCoUdxRSS1IpEXwF+BooLe7f2hm+wKPVfVEZtYAuBG4OcHHdzWzQjMrLCoqqurpRKS2fP01XHJJ2DVs40Z4+WUYMkRJoA6pNBG4+2LgOmB2dP9Dd7+rGufaH9gXmGdmS4G9gdlmVmbhcXcvcPc8d89r0qRJNU4nIjX2/PNhSujAgXDllbBgAXToEHdUUssSmTV0GjAXeDG638rMxlf1RO6+wN13d/dcd88FVgBHuPuqqr6WiCTZZ5+Fwd9TToEdd4Q33oD+/aFhw7gjkyRIpGvoVuBI4CsAd58L7FfZk8xsNPAmcJCZrTCzzjWIU0RSwR3GjAnlIZ54Am66CWbPhrZt445MkiihwWJ3X1NqcHdzZU9y93Mr+XluAucWkVRZuRK6d4dx46B16zAWcNhhcUclKZBIi2CRmf0JqGdmB5jZA8AbSY5LRFLFHYYODa2AiRPDbKDp05UEskgiieAyoDmwDhgNfA1cmcygRCRFPvgATjwxTAtt1SoMBl9zDdRX9ZlsUum/trsXAz2ji4jUBVuKxPXqBfXqwaBBcNFFsJXqUGajchOBmT1LBat/3f30pEQkIsm1aFFYEDZjRpgVNGgQ7L133FFJjCpqEdydsihEJPnWrw+bxNxxB+y0E4waBeeeq/pAUn4icPdXUxmIiCTRzJmhFbBgQfjwv+8+0EJNiVTUNTTG3f9oZgsoo4vI3TWlQCTdFRfDLbfAP/8Je+4J48fDaadV/jzJKhV1DV0RXZ+aikBEpJZNmRIGgN9/P2wf2bcv7Lxz3FFJGip3ioC7fxLd7O7uy0pegO6pCU9EqmzNGujWDdq3D2sEXnkl7BqmJCDlSGSu2K/KOPab2g5ERGrBhAnQvHnYNP7qq2H+/JAQRCpQ0RjBxYRv/vuZ2fwSP9oRmJbswESkCoqK4IorYPToUC107Fg48si4o5IMUdEYwePAC8CdwPUljq9197J2HhORVHMPxeEuvzx0Cd16K9xwA2yzTdyRSQapaProGmANcK6Z1QP2iB6/g5nt4O7ac1gkTitWwMUXh+6gI48M9YJatIg7KslAlZaYMLNLCaWoP+WHqqMOaPqoSBw2bw47hF17LWzYEKaGXn55KBUhUg2JVJa6EjjI3T9PdjAiUon33w9TQqdMCYPAgwfD/vvHHZVkuERmDX1E6CISkbhs3Ah33w2HHho2ihk8GCZPVhKQWpFIi+ADYIqZPUcoRQ2Au/8zaVGJyA8WLAjlIWbODKuCBw6EvfaKOyqpQxJJBMujyzbRRURSYd06+Mc/wqVRozA76I9/VJE4qXWJ7EdwWyoCEZESZswIrYBFi+C888LG8bvtFndUUkclMmuoCfA3wi5l22057u4nJDEukez0zTdhw/h77w3dPxMmhD0DRJIokcHiUcDbwL7AbcBSYGYSYxLJTq+8EvYJ7t8/1ApatEhJQFIikUTQ2N2HAhvc/VV37wSoNSBSW776KkwJ7dAhbBU5ZQo89FDYPEYkBRJJBBui60/M7BQzOxzYNYkxiWSPceOgWTMYNgz+9rdQJK5du7ijkiyTSCK4w8x2Bq4GrgGGAD0qe5KZDTOz1Wa2sMSxfmb2tpnNN7NnzGyXakcukslWr4ZzzoEzzgg7hc2YAXfdBdtvH3dkkoUqTQTuPsHd17j7Qndv7+6t3X18Aq89HDi51LFJQItod7N3gRuqHLFIJnOHkSPhkEPgmWfg9tuhsBDy8uKOTLJYIrOGHqHsrSo7VfQ8d59qZrmljr1U4u504KyEohSpCz76KAwCP/88tG0bisQ1axZ3VCIJdQ1NAJ6LLpOBnYD/1cK5OxHKXJfJzLqaWaGZFRYVFdXC6URisnlzWA3cvHkYCL73Xnj99WolgVGjIDc3jCnn5ob7IjWVyIKyp0veN7PRwOs1OamZ9QQ2EqamlnfeAqAAIC8v7yctEpGM8O670KULvPYanHgiFBTAvvtW66VGjQpbDxcXh/vLloX7APn5tRSvZKVEWgSlHQDsXt0TmtmfgVOBfHfXB7zUTRs3hs3iW7YMtYKGDYOXXqp2EgDo2fOHJLBFcXE4LlITiYwRrCWMEVh0vQq4rjonM7OTCauU27l7cWWPF8lI8+ZBp06hSuiZZ8KAAbDnnjV+2eXlbAVV3nGRRCXSNbRjdV446kI6HtjNzFYAtxBmCW0LTLJQOGu6u3erzuuLpJ116+COO6BPH9h1V3jySfj972utSFzTpqE7qKzjIjVRYSIws+2BfGDLqFYh8JS7r6/shd393DIOD61yhCKZ4M03Q5G4JUvgggvCrmGNG9fqKXr3/vEYAUCDBuG4SE2UO0ZgZocCi4FjCfWFlgK/BqaZ2S5mdkcqAhRJa//7H1x5Jfzyl6Fg3AsvwKOP1noSgDAgXFAAOTmhkZGTE+5roFhqqqIWwf1AV3efVPKgmZ0ILAQWJTMwkbQ3aVL4ir50KVxyCdx5J+yYWE/qqFFhkHf58tC107t3Yh/o+fn64JfaV1Ei2LN0EgBw95fNbANwZvLCEkljX34JV18NjzwCBx4IU6fCsccm/HRNA5V0U9H00a3MbNvSB81sO0IlUs36kezzzDNhIdiIEXDDDWGGUBWSAGgaqKSfihLBCOBpM8vZciAqGTEGeCy5YYmkmVWr4A9/gN/9Dn72M3jrrbCF5HbbVf7cUjQNVNJNuYnA3e8AXgReM7PPzOwz4FVgkrvfnqoARWLlHr79N2sGzz4bPvzfeguOOKLaL1nedE9NA5W4VLiy2N0fdPemhN3J9nX3HHd/IDWhicRs2TL4zW/gwgtDtdC5c0N30NZb1+hle/cO0z5L0jRQiVNCJSbcfa27r012MCJpYfNmePDBUCTu9dfhgQdCraCDDy7z4VUtBKdpoJJuKl1ZLJJNxvd7hz1v6kybddN4dbtf82Wfhznj0pxyH1/dGUCaBirppDpF50Tqng0bmHv2nZz0t5bsv24xFzKc4797gfwbcyr8hq8ZQFIXVJoIzKyBmd1kZoOj+weY2anJD00kRebMgaOOotWYG3mW02jGYkZwIWCVfqhrBpDUBYm0CB4B1gFHR/c/BlReQjLfd9/BjTdCmzawciW/52n+yJN8ys9+9LCKPtQ1A0jqgkQSwf7u3hfYABAtJKudcooicXn9dWjVKpSFuOACWLKEWTm/K/OhFX2oawaQ1AWJJIL1URVSBzCz/QktBJGMUHJWT/Oma3nnpEvDauB162DixLBpTKNG1fpQ1wwgqRPcvcIL8CvCQrIiwtaSS4HjK3tebV5at27tItUxcqR7gwbu4H4SL/pSmvomzJecdLn72rVlPj4nx90sXI8cmfKQRWoNUOgJfMaaJ7BbpJk1BtoSuoSmu/tnScpLZcrLy/PCwsJUnlLqiNxc+HrZF/SnBxcygiUcTGeGsjLnGJYujTs6keQys1nunlfZ48pdR2BmpdfQfxJdNzWzpu4+uyYBiqRCm2VP8SCXsCtfcAc9uYNerGM7TLN6RL5X0YKyeyr4mQMn1HIsIrXnk0/g0kt5krHM4gh+zUTm0er7H2tWj8gPyk0E7t4+lYGI1Ap3GD4crroKvv2WOWf3of34q1n77Q//1TWrR+THEllQtp2ZXWVmY83saTO7MtqTQCS9fPghnHQSdOoEhx4K8+dz+BPXMXBwfc3qEalApYPFZjYGWAuMjA79CdjF3f+Q5Ni+p8FiqdCmTTBgQKgMutVW0Lcv/PWv4bZIFqvxYHEJLdy9WYn7/zGzxdUPTaQWLVkCnTvDm2+GktGDBmkAQKSKEvnKNNvM2m65Y2ZHAZV+PTezYWa22swWlji2q5lNMrP3outG1Qtbst6GDaGjv1UreOcdeOwxeO45JQGRakgkEbQG3jCzpWa2FHgTaGNmC8xsfgXPGw6cXOrY9cBkdz8AmBzdF6maWbMgLw969YIzzgitgvPOC4MAIlJliXQNlf4wT4i7T432OC6pI3B8dPtRYApwXXVeX7LQt9/Crbey+e57KLLd6cYzzJlxBr0nafBXpCYqTQTuvizqwtmn5OOruaBsD3ffsjBtFbBHNV5DstHUqdClC7z3Ho/W60KPTf1Ywy6Q4EYwIlK+RKaP3g7MB+4nLDK7B7i7pieO6mCUO2XJzLqaWaGZFRYVFdX0dJLGKtzq8euvoXt3aNcONm7kT7u/TKdNg0MSiGgjGJGaSaRr6I+EUtTra+F8n5rZnu7+iZntCawu74HuXgAUQJg+WgvnljRU4VaPjZ6Hbt1gxQro0QNuv50ndmxY5utoIxiR6ktksHghlPj6VTPjgQuj2xcC42rpdSVDlbXV4/bFn7H9X8+HU06BHXeEN96Af/4TGjbURjAiSZBIIrgTmGNmE81s/JZLZU8ys9GEGUYHmdkKM+sM9AF+ZWbvASdG9yWL/fibvPMHxrCYZpz2zRNw880weza0/X72sjaCEUmCRLqGHgXuAhYAmxN9YXc/t5wfdUj0NaTua9o0dAftyUoeojtnMI6Z5HHBni/z4m2H/eTxWwaEe/YMSaRp05AENFAsUn2JJIJid78/6ZFIVup9hzOt81D+sf4atmUdV3M3g7e/goH9yv+vmZ+vD36R2pRIInjNzO4k9O9/v0Wl9iOQGvvgA/IfuYj89a8wfdt2nL9uCBtyfsFAfcMXSalEEsHh0XXbEse0H4FU36ZNcP/9oX+nfn14+GHadunCeyoSJxKLRBaUaV8CqT2LFoUicTNmhFlBgwbB3nvHHZVIVkukRYCZnQI0B77fh8Dd/56soKQOWr8e+vSBO+6AnXeGxx+Hc85RfSCRNFBpIjCzQUADoD0wBDgLeCvJcUldMnNm2Cxm4UL405/g3nuhSZO4oxKRSCKdsse4+wXAl+5+G3A0cGByw5I6obgYrrkmrAP48ksYPz4sJVYSEEkriXQNfRtdF5vZz4HPgT2TF5LUCVOmhCJx//1v2C3srrtCl5CIpJ1EWgQTzGwXoB8wG1gKPJ7MoCSDrVkTPvjbR3MMXnklDAgrCYikrUoTgbvf7u5fufvTQA5wsLvfnPzQJONMmADNm8OQISz+7TUcvH4+W3Vo/9OKoiKSVspNBGbWxsx+VuL+BcAY4HYz2zUVwUmGKCoKg8CnnQaNGvHiLW/SZko/3vmoAe4/VBRVMhBJTxW1CB4G1gOY2XGEAnEjgDVE5aElO5S7X4A7jB4NzZrBU0/BbbfBrFl0G3bkTyqKas8AkfRV0WBxPXf/Irp9NlAQdQ89bWZzkx+apIPy9gvY/vMV/G7SxaE76KijYOjQ0C1E+XsDaM8AkfRUUYugnpltSRQdgFdK/CyhhWiS+UrvF2Bs5rzih/nVlc1g8uSwT8C0ad8nASh/bwDtGSCSnipKBKOBV81sHGEK6WsAZvYLQveQZIGS3+L3530m04GH6cZb3iYsEOvRA+rV+9FztGeASGYpNxG4e2/gamA48H/RHsNbnnNZ8kOTdNC0KdRjI1dzNws4lCOYTRcG07npy7DffmU+Jz8fCgogJydUkMjJCfdVUVQkPVXYxePu08s49m7ywpF0M6DbAvbo2Zm8zTMZx+l05yG+arAXBf+o+HnaM0Akc6jur5Rt3Tq45RZOuekIWuywlEt2+xdn8m+2ztlL3+5F6hgN+spPTZ8eSkUvXgznncd2997LgMaNGRB3XCKSFGoRyA+++QauugqOOQa+/hqeew4eewwaN447MhFJIiWCLFLuwjAIU0EPPRT694du3cIGMr/9bUyRikgqqWsoS5S3MGzrb77ijzOvhSFD4IAD4NVX4bjj4g1WRFJKiSBLlF4YBnBi8Tjadb8YWA3XXQe33ALbbx9LfCISn1i6hsysh5ktMrOFZjbazLar/FlSEyUXhjVhNaM5h3GcwSebdg/7B/fpoyQgkqVSngjMbC/gciDP3VsA9YBzUh1HtgnlHZx8RrKEQziTZ+jJHZzVdCa0bh13eCISo7i6huoD25vZBsJ+yCtjiiNr9O+xnAZXdePXm1/gDY6mM0NZ3uCQSheGiUjdl/IWgbt/DNwNLAc+Ada4+0ulH2dmXc2s0MwKi4qKUh1m3bF5MwwcyJm9mtNh61e5tdF9HMdrfJtziBaGiQgQT9dQI6AjsC/wc6ChmZ1X+nHuXuDuee6e10SbnVfPu+/C8cdD9+7Qti31lyzk1i8uZ6PXY+lSJQERCeIYLD4R+NDdi9x9AzAWOCaGOOqujRuhb19o2RIWLIBhw+Cll2DffeOOTETSUBxjBMuBtmbWgFDeugNQGEMcddO8edCpE8yeDWeeCQMGwJ57xh2ViKSxOMYIZgBPAbOBBVEM2voyARWuDP7uO+jVC/Ly4OOPw9aRY8cqCYhIpWKZNeTutwC3xHHuTFXeymCA/H3fCEXi3n4bLrww7Bq2667xBSsiGUUrizNEWSuDrfh/rL+4J/zvAdhnH3jxRfj1r+MJUEQylorOZYjSG7+fyCQWcCh/WXs/XHJJ2DZSSUBEqkGJIENs2fh9F75kKJ2YxEmsY1vO2uM1eOAB2HHHeAMUkYylRJAheveGc7YZy2KacQEj+Ac3cMz2cznznv+LOzQRyXAaI8gEq1aR/8yl5K9/mkVbt+KUDc/zRc7hPNBbi8JEpOaUCNKZO4wYAT16hJHif/yD5tdcw+ytt447MhGpQ5QI0tWyZfDXv8LEifDLX4aNYw4+OO6oRKQO0hhButm8GR58EJo3h2nTwu2pU5UERCRp1CJIJ++8ExaGTZsWpoI+/DDk5MQdlYjUcWoRpIMNG+DOO0ORuMWL4dFH4YUXlAREJCXUIojbnDmhSNzcuXDWWaEraI894o5KRLKIWgRx+e47uOEGaNMGVq2Cp5+GJ59UEhCRlFOLIA6vvx7GAt59F/7yF7jnHmjUKO6oRCRLqUWQSmvXwqWXwrHHwvr1YbOYYcOUBEQkVkoEqTJxIrRoAQ89BFdcEXYO+9Wv4o5KRESJIOk+/zzsEXDyydCwYZgaeu+9sMMOcUcmIgIoESSPe9glrFkzePzxsHvYnDlw9NFxRyYi8iMaLE6GTz4JewQ88wy0bh3GAlq2jDsqEZEyqUVQm9zhkUdCK+CFF+Cuu2D6dCUBEUlrahHUlg8/DJsIv/wyHHccDB4MBx4Yd1QiIpVSi6CmNm2C++8PM4JmzICBA+E//1ESEJGMoRZBTSxeDF26wJtvwm9+E4rE7bNP3FGJiFRJLC0CM9vFzJ4ys7fNbImZZdZUmg0b4I474PDDw+rgkSPhueeUBEQkI8XVIrgPeNHdzzKzbYAGMcVRdbNmhSJx8+fD2WeHbqHdd487KhGRakt5i8DMdgaOA4YCuPt6d/8q1XFU2bffwnXXwZHOxv3xAAAK9UlEQVRHQlER/Pvf8MQTSgIikvHi6BraFygCHjGzOWY2xMwaxhBH4qZODVNA+/YNxeIWL4aOHeOOSkSkVsSRCOoDRwAD3f1w4Bvg+tIPMrOuZlZoZoVFRUWpjjH4+mvo3h3atQuzgyZPhoIC2GWXeOIREUmCOBLBCmCFu8+I7j9FSAw/4u4F7p7n7nlNmjRJaYAAPP982Df44YfhqqvCmMAJJ6Q+DhGRJEt5InD3VcBHZnZQdKgDsDjVcZTrs8/gvPPglFNgp53gjTfCfgEN07v3SkSkuuKaNXQZMCqaMfQB8JeY4viBO4wZA5ddBl9+CbfcEnYQ23bbuCMTEUmqWBKBu88F8uI4d5lWroSLL4bx48PWkZMnw6GHxh2ViEhKZHeJCXcYMiQUiZs0Ce6+O6wSVhIQkSySvSUmPvgALroIXnkFjj8+FIn7xS/ijkpEJOWyr0WwaRP07x+KxBUWhllBkycrCYhI1squFsHChWFB2Ftvwamnhkqhe+8dd1QiIrHKjhbB+vVw221wxBGhS+jxx8PAcJKSwKhRkJsLW20VrkeNSsppRERqRd1vEcycGYrELVwIf/pT2Dg+iQvURo0K+9MUF4f7y5aF+wD5+Uk7rYhItdXdFkFxMVxzDbRtG9YFPPts+JRO8irlnj1/SAIlQ+nZM6mnFRGptrrZIpgyJWwY89//wl//GvYO3nnnlJx6+fKqHRcRiVvdahGsWRM++Nu3D/f/8x8YNChlSQCgadOqHRcRiVvdSQTPPhsWhg0ZErqE5s8P6wNSrHdvaFBqm50GDcJxEZF0lPmJoKgoDAKffjo0bgzTp0O/fj/9NE6R/PxQqTonB8zCdUGBBopFJH1l7hiBO4weDZdfHvYN+Pvfww5i22wTd2Tk5+uDX0QyR2YmghUrQpG4CRPgqKNg6NCwd4CIiFRZZnUNbd4cSkI0axZqBPXvD9OmKQmIiNRA5rQI3n8/FImbMgU6dAgd7/vtF3dUIiIZLzMSwaefhtLQ224bZgV16hRGYkVEpMYyIxGsWAEdO8JDD8HPfx53NCIidYq5e9wxVMrMioBlcccB7AZ8FncQNZDp8UPmvwfFH79Mfw9ViT/H3Sutq5MRiSBdmFmhu6fPFptVlOnxQ+a/B8Ufv0x/D8mIP7NmDYmISK1TIhARyXJKBFVTEHcANZTp8UPmvwfFH79Mfw+1Hr/GCEREspxaBCIiWU6JIAFmtouZPWVmb5vZEjM7Ou6YqsrMepjZIjNbaGajzWy7uGOqiJkNM7PVZrawxLFdzWySmb0XXTeKM8bKlPMe+kX/j+ab2TNmtkucMVakrPhL/OxqM3Mz2y2O2BJV3nsws8uif4dFZtY3rvgqU87/oVZmNt3M5ppZoZkdWdPzKBEk5j7gRXc/GGgJLIk5nioxs72Ay4E8d28B1APOiTeqSg0HTi517HpgsrsfAEyO7qez4fz0PUwCWrj7YcC7wA2pDqoKhvPT+DGzfYCTgEzYd284pd6DmbUHOgIt3b05cHcMcSVqOD/9N+gL3OburYCbo/s1okRQCTPbGTgOGArg7uvd/at4o6qW+sD2ZlYfaACsjDmeCrn7VOCLUoc7Ao9Gtx8FzkhpUFVU1ntw95fcfWN0dzqwd8oDS1A5/wYA/YG/AWk/wFjOe7gY6OPu66LHrE55YAkqJ34Hdopu70wt/C0rEVRuX6AIeMTM5pjZEDNrGHdQVeHuHxO+9SwHPgHWuPtL8UZVLXu4+yfR7VXAHnEGUws6AS/EHURVmFlH4GN3nxd3LDVwIHCsmc0ws1fNrE3cAVXRlUA/M/uI8Hdd41alEkHl6gNHAAPd/XDgG9K/S+JHor70joSk9nOgoZmdF29UNeNhulvafyMtj5n1BDYCo+KOJVFm1gC4kdAdkcnqA7sCbYFrgTFmGVXF8mKgh7vvA/Qg6q2oCSWCyq0AVrj7jOj+U4TEkElOBD509yJ33wCMBY6JOabq+NTM9gSIrtO2SV8RM/szcCqQ75k1f3t/wpeJeWa2lNCtNdvMfhZrVFW3AhjrwVvAZkL9nkxxIeFvGOBJQIPFyebuq4CPzOyg6FAHYHGMIVXHcqCtmTWIvvl0IMMGvCPjCX8ERNfjYoylWszsZEL/+unuXhx3PFXh7gvcfXd3z3X3XMIH6hHR30gm+TfQHsDMDgS2IbOK0K0E2kW3TwDeq/ErursulVyAVkAhMJ/wn6hR3DFV4z3cBrwNLAQeA7aNO6ZK4h1NGM/YQPjA6Qw0JswWeg94Gdg17jir8R7eBz4C5kaXQXHHWZX4S/18KbBb3HFW499gG2Bk9LcwGzgh7jirGP//AbOAecAMoHVNz6OVxSIiWU5dQyIiWU6JQEQkyykRiIhkOSUCEZEsp0QgIpLllAgkZcyscVQxca6ZrTKzj6PbX5lZStdmRBUcf1vi/ulmVq0V42a2tKwqnGa2s5mNMLP3zey/ZjYqGRVTK3ovZnarmV1T2+eUukWJQFLG3T9391YeqiYOAvpHt1sRVnfWqqjAXnlaAd9/eLr7eHfvU8shDAU+cPdfuPv+hDUEw2v5HJCa9yJ1mBKBpIt6ZjY4qg//kpltD2Bm+5vZi2Y2y8xeM7ODo+O5ZvZKVNd/spk1jY4PN7NBZjYD6GtmDaOa7m9FRQM7mtk2wN+Bs6MWydlm9mczezB6jT2ivQLmRZdjouP/juJYZGZdK3ozZvYLoDVwe4nDfwdamtlBZna8mU0o8fgHo9ITmNnNZjbTwt4RBVvq4JjZFDO7K3ov75rZsZW9l1Ixlfe7/EN0rnlmNrXq/3SS6ZQIJF0cAAzwUB/+K+D30fEC4DJ3bw1cAzwUHX8AeNRDXf9RwP0lXmtv4Bh3vwroCbzi7kcSygr0A7YmFE77V9RC+VepWO4HXnX3loS6Uoui452iOPKAy82scQXvpxkw1903bTkQ3Z4DHFLJ7+JBd2/jYe+I7Ql1ibaoH72XK4Fb3H19Je+lpPJ+lzcDv47e7+mVxCZ1UEVNZ5FU+tDd50a3ZwG5ZrYDoTjekyWKQ24bXR8N/C66/Rg/3pzjyRIfwCcBp5foJ98OaFpJLCcAF8D3H95rouOXm9mZ0e19CMnr88TeXpW0N7O/EfaN2JWQiJ6Nfral2NgsIDfRF6zkdzkNGG5mY0q8vmQRJQJJF+tK3N5E+Ca8FfBVNI5QFd+UuG3A7939nZIPMLOjqvKCZnY8oYrr0e5ebGZTCEmlPIuBVma2lbtvjl5jK8IOd7MJyahki3y76DHbEb6p57n7R2Z2a6nzbPk9baJqf7/l/i7dvVv0+zgFmGVmrd09GQlO0pS6hiRtufvXwIdm9gcAC1pGP36DH7bbzAdeK+dlJgKXlehnPzw6vhbYsZznTCbUfMfM6lnYpW5n4MsoCRxMqGVfUezvE7qBepU43Iuw1eZyYBnQzMy2tbBvcYfoMVs+9D+LvsWfVdF5EngvW+Ip93dpZvu7+wx3v5mwCdM+CZxT6hAlAkl3+UBnM5tH6CLpGB2/DPiLmc0HzgeuKOf5txPGBOab2SJ+GLz9D+GDeK6ZnV3qOVcQumcWELpgmgEvAvXNbAnQh7DNZGU6AQdEU0eLCMmjG4C7fwSMIVTAHENIGnjYBnVwdHwiMDOB81T0Xkoq73fZz8wWWNgg/Q1CVUvJIqo+KpICFvazeA643N2fjzsekZKUCEREspy6hkREspwSgYhIllMiEBHJckoEIiJZTolARCTLKRGIiGQ5JQIRkSz3/xMoFiei1Y8YAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "probplot = sm.ProbPlot(res, stats.gamma, fit=True)\n", "fig = probplot.qqplot(line='45')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEKCAYAAAASByJ7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VOX1x/HPARc2sYq4VCRx14C4kOJWReuuVev2qzVaFRARF9CqqNRarWgVF9xA2dyIaFUsiFgEFFkUJIDsihSJ4lJRccEoW87vj2eiQ5gkkzCTO5N8369XXsncuXPvyYg582znMXdHRESkQdQBiIhIZlBCEBERQAlBRERilBBERARQQhARkRglBBERAZQQREQkRglBREQAJQQREYnZLOoAqmO77bbz3NzcqMMQEcl8a9bAsmXw/ffMhC/dvWVVL8mqhJCbm0tRUVHUYYiIZK716+GRR+Cmm8AMHn4Yu+KK4mReqi4jEZG6YtEiOPJI6NEDjjgC5s+Hyy9P+uWRJQQz28XM3jCzhWa2wMx6RBWLiEhWW7sW+vSBAw6A996Dp56CMWMgJ6dal4myy2gd8Bd3n2VmWwEzzWycuy+MMCYRkewyaxZ06gRz5sA558BDD8EOO9ToUpG1ENz9M3efFfv5e2ARsHNU8YiIZJUff4QbboAOHeB//4MRI+Bf/6pxMoAMGVQ2s1zgQGB6tJGIiGSBSZOgSxf44APo3Bn69oVtttnky0Y+qGxmzYAXgZ7u/l2C57uaWZGZFa1YsaL2AxQRyRTffRcGiTt2DOMG48bB4MEpSQYQcUIws80JyaDQ3UckOsfdB7p7vrvnt2xZ5TRaEZG66dVXoW1bGDAAevYMM4iOPRaAwkLIzYUGDcL3wsKa3SKyLiMzM2AIsMjd74sqDhGRjPbVV3D11fD007DvvjB1Khx66M9PFxZC165QUhIeFxeHxwAFBdW7VZQthMOBC4Dfmdm7sa+TI4xHRCRzuIdB4n33heHD4eabYfbsDZIBQO/evySDMiUl4Xh1RdZCcPcpgEV1fxGRjPXpp9C9O4wcCe3bh7GC/fdPeOpHHyW+REXHKxP5oLKIiMS4w5AhkJcHY8fC3XfDtGkVJgOA1q2rd7wySggiIplg6dIwSNylS0gAc+fCddfBZpV35PTpA02abHisSZNwvLqUEEREorR+PfTrB/vtx9q3ZtB72wE0nPQGucftmdRsoYICGDgwVKkwC98HDqz+gDIoIYiIpEVSU0EXLIDDD4err+aTPY+iDQu44+tulNLg59lCySaFZcugtDR8r0kyACUEEZGUK5sKWlwchgU2+uO+Zg384x9w4IGwZAkMG8bhK0fzwU+7bHCdms4Wqilz99q72ybKz8937YcgIpkuNzckgfJycmDZ8zNCuYl58+Dcc+GBB2D77WnQICSP8szCJ/9NYWYz3T2/qvPUQhARSbFEUz4bU8IVxdfBIYeExWYjR4b1BdtvD6R2tlBNKSGIiKRY+T/iHZnIHPbnWu4JrYOFC+G00zY4J5WzhWpKCUFEJMXK/rg351sG0I2JHE1DK2X8jRPCFKCtt97oNamcLVRTGkMQEUmDN659hX36Xcr26z9j6FY92arfPzi3U5OqX5gGyY4hZMR+CCIidcaKFdCzJ0c/8wy0aQNDXuSSgw+OOqqkqMtIRCQV3MMgcV4ePP88/P3vYXvLLEkGoBaCiMimW74cLrsMRo8OW1oOGRL2LsgyaiGIiNRUaWkY+W3TBiZMgHvvhbfeyspkAGohiIjUzJIlcMklMHEiHH00DBoEu+8edVSbRC0EEZHqWL8+tATatQtjBAMHhtZBlicDUAtBRCR58+dDp04wYwacemrY33jnnaOOKmXUQhARqcrq1WHW0EEHwYcfhtlEI0fWqWQAaiGIiFRu+vRQbmLBgrBsuF8/2G67qKNKC7UQREQS+eEHuOaasKn9t9+GKaXDhtXZZABqIYiIbOz118MMoqVLoVs3uOsuaN486qjSTi0EEZEy33wTEsExx4StziZODAPH9SAZgBKCiEgwalRYYDZ0aNjcfs4c6Ngx6qhqlRKCiNRvX3wRdi47/XRo0SIMIt9998abE9QDSggiUj+5h0HiffeFl14KexwXFUF+lVWi6ywNKotI/fPxx2GweMyYsKXlkCGhSmk9pxaCiNQfpaVhkLhNmzBg3K8fTJmiZBCjhCAi9cMHH4QidN27hxLV8+ZBjx7QsGFKb1NYCLm5YZJSbm54nC2UEESkblu3LgwSt2sXZg4NGQLjxsFuu6X8VoWF0LUrFBeHIYri4vA4W5KCEoKI1F1z5oQxgl694MQTYeHCUJzOLC23690bSko2PFZSEo5nAyUEEal7Vq+Gm28OM4Y+/hj+9S8YMQJ+/eu03vajj6p3PNMoIYhI3fL223DggXD77XDeeaFVcM45aWsVxGvdunrHM40SgojUDatWQc+ecPjh4ecxY+DJJ8Nis1rSp8/G69maNAnHs0GkCcHMhprZF2Y2P8o4RCTLjRsH++0HDzwA3bvz3N8WkHvZSbU+06egIGyglpMTGiQ5OeFxQUHt3H9TRd1CeAI4MeIYRCRbrVwZ9io4/njYYguYNInCQx+mU4+tIpvpU1AAy5aFJQ/LlmVPMoCIE4K7TwK+jjIGEclSL70UFpQ9+STccEOYUXTEEVk/0ydKKl0hItnl88/hyivhhRfggAPglVfC1pYx2T7TJ0pRdxlVycy6mlmRmRWtWLEi6nBEJCru8NRToVUwalQYqX3nnQ2SAWT/TJ8oZXxCcPeB7p7v7vktW7aMOhwRiUJxMZx0Elx4YahOOmcO3HQTbL75BqcVFoYJRuVl00yfKFWZEMxsBzMbYmavxh7nmVnn9IcmIvVeaSk88gi0bRuK0D34IEyeDPvss9GpZWUjvvpqw+MtWmTXTJ8oJdNCeAIYC5Qt8VsM9EzFzc1sOPA2sLeZLVeiEZGfvf9+2LHsiivgsMNg/vwwdtCgAd27w2abhamdm20W6tUlGkwGaNZMySBZyQwqb+fu/zKzGwHcfZ2ZrU/Fzd39T6m4jojUHc88uZaPe95Lj2/+zk8NmrCw6xMc9uiff15p3L17qGBdZv36DR+Xp8Hk5CXTQvjBzFoADmBmhwDfpjUqEamXxvSZTd7FB9PrmxsZze/Zp3Qhxw27kMJnfik7MXBg9a6pweTkJZMQrgFGAbub2VTgKeDKtEYlIvXLTz/BTTdx/F9/w47+KWfxAufwAv9jx43WEKyvpH8im8tGZIIqE4K7zwI6AocBlwJt3H1uugMTkXpi6tSwnuDOO3maC8hjISM4a4NT4rt9KtrPpmHD7C4bkQkqHEMwszMreGovM8PdR6QpJhGpD77/PkwdfeSR0K8zdiy3dj2elcUbnxrf7dO1a+Ixg65dwx9/JYCaq2xQ+dRKnnNACUFEambs2PAX/OOPw8yhPn2gWTP69AmH42cLle/26d8/fB84MHQfNWwYXlN2XDaBu2fNV/v27V1EMtOwYe45Oe5m4fuwYQlO+uor9wsvdAf3ffZxnzKlZteRagGKPIm/sVVOO43NMLoF+C2hZTAFuM3dv6r0hSJSb5QtCiv7ZF9WYRTiunBeeAEuvzysHOvdG/76V2jUaKNrqdsnOsnMMnoWWAGcBZwd+/m5dAYlItml0gqjn30GZ50Vdi3beWcoKgq7mSVIBhKtZBam7eTu/4h7fLuZ/TFdAYlI9km8+Mv5XfETkHcN/Pgj/POf8Je/hKXFkpGSaSG8ZmbnmlmD2Nf/EUpZiIgAsO22Gz7OYRljOYGhdAo7mc2ZA716KRlkuGQSwiXAM8Ca2NezwKVm9r2ZfZfO4EQkuzRgPVfyIPNpy6G8zXVNH4GJE2HvvaMOTZJQZbp2961qIxARyV5ffw37sIghdOYw3uZVTuRSHmN5SWv6ZnyRfSmT1H8qM2tnZqeZ2ZllX+kOTESiU1gYNqdPapP6tWu5a+s+vMsB7M37XMBTnMwYPqa16ghlmWSmnQ4F2gELgNLYYS1ME6mjkppCWmbmTOjUieu+mcsLDf+Py9c/yBfsAKiOUDZKZoTnEHfPS3skIpIRKptC+nNC+PFHuPVWuOceaNkSXnqJ1T/8gca9wT4KpSb69NF6gmyTTEJ428zy3H1h2qMRkchVuUn9pEnQpQt88AF07gx9+8I221CAEkC2S2YM4SlCUnjfzOaa2TwzU7VTkSyUzNhARf3+ea2+CyuNO3aEdetg/HgYPBi22SadIUstSqaFMAS4AJjHL2MIIpJlkh0bSFRg7owtx/B0STcYsBx69gwrjZs2rb3gpVYk00JY4e6j3P1Ddy8u+0p7ZCKSUpWWl4hTUPDLvgLb8SUvNr2AEatPoekOW8Fbb8H99ysZ1FHJtBBmm9kzwMvA6rKDrv0QRLJKlWMDcQrOcwo2fz5scL9yJdx8c8gcW26Z3iAlUskkhMaERHB83DFNOxXJMq1bh26iRMc38OmnYSf7kSOhffswVtCuXa3EKNFKZqXyxbURiIikT2EhrFq18fEN1gq4w9ChoQDd6tVh9lDPnqo/VI8kszCtEdAZaAP8XK/W3TulMS4RSZHyg8llWrSABx6IDSgvXQqXXAKvvx5mEQ0eDHvsEUm8Ep1kBpWfBnYETgDeBFoB36czKBFJnUSDyQDNmkHBuevDIPF++8GMGfDooyEpKBnUS8m0Bfdw93PM7HR3fzI2wDw53YGJSGokGjcAaFa8AA7vDNOnwymnhGTQqlXtBicZJZkWwtrY92/MrC2wNbB9+kISkVQpLASzDY9tzhpu5jZmcSAsWRJOevllJQNJKiEMNLNtgJuBUcBC4K60RiUi1VLRCuTevcNYcZl8ZlBEPrdxC58eejYsWgTnnbdx1pB6KZlZRoNjP74J7JbecESkuipbgVy2xqAxJdzKLVzDfXzGTpzKKF5+69RoApaMVWELwcxONbOcuMd/M7M5ZjbKzHatnfBEpCo9elS8Arl1a+jIROawP9dxD4PpQhsWMC9HyUA2VlmXUR9gBYCZ/R44H+hE6DZ6NP2hiUi8RN1ChYXw1VeJz/+m+FvG7d6NiRyN4RzN63TjMdY12Vr7FEhClXUZubuXfe44Exji7jOBmWbWPf2hiUiZirqFGjdOfP4pjGZQw27sNPEzFp78F86cdxuLlzchR/sUSCUqSwhmZs2AEuAYoH/cc40Sv0RE0qGiwnTlj23HCh6gB+cxnG92agsvjiCvQwfeq71QJYtV1mXUD3gXKAIWuXsRgJkdCHxWC7GJSExFhel+4ZzLcBaSx9m8wF2N/86v/jsTOnSojfCkjqgwIbj7UKAjoWzFyXFPfQ6ovpFIGlQ0fbSiTWtatIA9Gi1nFKcxnPNYym4c3mgWrQbdAltsUVthSx1R6ToEd//E3We7e2ncsc/cvcrPK8kwsxNjO7EtMbMbUnFNkWxVNk5QXBzWDpSNExQWhn7/Jk02PL9p41JeOX0gC6wNx9oEruE+/tT6LXoObqsxAqkR8/hVK7V5Y7OGwGLgOGA5MAP4U2V7N+fn53tRUVEtRShSu3JzE5eZyMmBZctCYujdO3QfHbHTEp7b+hJ2XDQRjj4aBg2C3Xev5YglW5jZTHfPr+q8ZFYqp0sHYIm7L3X3NcCzwOkRxiMSmcLCimsOlY0fFBTAsiXrKL37Ht78ej92/GRWSAQTJigZSEpUmRDM7OlkjtXAzsDHcY+Xx46Vv1dXMysys6IVK1ak4LYimaWsq6giP48fzJsHhx0G110Hxx0HCxdCly4qOyEpk0wLoU38g1hXT/v0hLMxdx/o7vnunt+yZcvauq1I2pQfOO7WLXF5agjjBnf+fTXccgscdFDoO3r22bCb2c4bfX4S2SQVrkMwsxuBm4DGZvZd2WFgDTAwBff+BNgl7nGr2DGROivRArPKvHj9dE68pzMsWADnnx/2Lthuu/QHKvVSZdNO73T3rYC+7t489rWVu7dw9xtTcO8ZwJ5mtquZbQGcSyiLIVJnVbRZTXlN+IFBW13DibceCt9+C6NHw9NPKxlIWiXTZTTazJoCmNn5ZnZffNG7mnL3dcAVwFhgEfAvd1+wqdcVyWRVLzCDo3mdubSjy/f3h/6kBQvCBjYiaZZMQhgAlJjZ/sBfgP8CT6Xi5u4+xt33cvfd3V3ltqTOq2iBGcDWfMNALuF1jsEaNICJE6F/f2jevNbik/otmYSwzsNihdOBh939EWCr9IYlUjclWmAGcBojWUgenRjKfZtdzzuD54bN7kVqUTIJ4fvYAPP5wCtm1gDYPL1hiWS/RGUoCgpg4MCw2MwMDmr1BTP2OJeR/IEVtOQPO05nhyfu4tyLKyhjKpJGVe6YBvwROA/o7O6fm1lroG96wxLJXoWFYdOa+H0K4ncxKyiAgvP8lxO/WAX/+Af79+rFy5vrs5ZEJ7LSFTWh0hWS6cpPKy0vJweWTf44DBaPGQOHHAJDhkBeXu0GKvVKWktXmNm8mrxOpK6rbFqpUcrJxQOgTZswYNyvH0yZomQgGaOyhWlnVvQUsGN6whHJbhVNK92TxQymC0cyGQ4+Ngwk7KqtySWzVDaG8BxQCCTqU9KOaSIJtG694erjhqzjGu7jVm7hJxrx9iVDOfSxi1R/SDJSZQlhLnCPu88v/4SZHZu+kESyV58+v4whtGMOQ+lEe2YxeoszWH3vI5x1xU5RhyhSocoSQk/guwqeOyMNsYhkvYICmD5pNTsMup3r/Z98zbYMOv55LvnPWWoVSMarMCG4++RKntNUH5EExv79bS4b1Jl9fRFP8meu4T5+mtKCJs+gXcwk40W5QY5I3bFqFfTsyXG3Hk4T/4ETeZWLeJKvaUFJSZh9JJLpklmYJiKVGTcuDBwsW8YALucG7mRVueouyRS1E4lahS0EM+sR+3547YUjkkVWroROneD442GLLWDSJPrmPLxRMoDKi9qJZIrKuowujn1/qDYCEckqL70UFpQ99RTceCPMmQNHHJGweF2TJmH2kUimqywhLDKzD4C9zWxu3Nc8M5tbWwGKpFpZ0Tkz2Gyz8L2s+FyVPv8czjkHzjwTdtwR3nkH7riDwhcbkZsLF1wAjRtDixbhujk5YQ2aBpQlG1Q2y+hPZrYjYQOb02ovJJH0KV9raP368L188bmNuIcdy3r2DC++4w649lrYfPONrvnVV6FV8PTTSgSSXZIqbhfb4nKv2MP33X1tWqOqgIrbyabKza18H+OcnLCP/QaKi+HSS2HsWDjssFCMbp99qrxmwmuJRCBlxe3MrCPwAfAI0B9YbGZHbnqIIrWnrJuoqk3tN5gNVFoKDz8citFNmQIPPQSTJ2+QDDZ6TRLHRTJVMtNO7wOOd/f3AcxsL2A40D6dgYmkSlUlqeP9PBvo/fehc2eYOhVOOAEeeyx85K/gNYkSjWYWSbZJZmHa5mXJAMDdF6Md0ySLVFaSOl6TJnDHrWvhzjth//1h4UJ44gl49dUKkwEk3hZTM4skGyXTQigys8HAsNjjAkAd+ZI1Kuu6adgwDCzn5ED/S2Zz8gOdYfZsOPvs0EW0Y9WV3ssGjnv3Dvdq3TokAw0oS7apclDZzLYELgd+Gzs0Gejv7qvTHNtGNKgsNVHloO9PP8Ftt8Hdd8N220H//mFaqUgdkeygcpUthNgf/vtiXyJZJ74kdZmfu3SmTIEuXcKYwcUXw733wjbbRBarSJRU3E7qvIKCsDgsJ+eXxWJDH/iegrevgCOOgNWrw5TSoUOVDKReU3E7qRcKCuL69MeODU2Gjz+Gq64KTYVmzSKNTyQTVKuFYGYNzKx5uoIRSYeyNQgt7GtebHYhnHhi6DOaMgUeeEDJQCQmmYVpz5hZczNrCswHFprZdekPTWTTla1ByC9+gYXsy2k/PMNdm/Vm+PWzw6pjEflZMi2EPHf/DvgD8CqwK3BBWqMSSZF+vT7jqZKzeIFzWE4rfsMMblh3Ozfe2ijq0EQyTlIL08xsc0JCGBWrY1R1ASSRNCnrAmrQoJIqpe7w+OO89kkep/AKvfgnBzOdORwAqKyESCLJDCo/BiwD5gCTzCwH+C6dQYlUpHwZioRVSj/8MBwcP57FWx7BBasH88HPtRkDlZUQ2ViVLQR3f9Ddd3b3kz0oBo6uhdhENpKoDMXPexavXw8PPght28K0adC/P0sGTeSTJhsmA5WVEEksmUHlHcxsiJm9GnucB1yY9shEEqioq6dJ8aKwpqBHD+jYERYsgMsuo+CCBhutQdCGNSKJJTOG8ARhk5xfxx4vBnqmKyCRihQWhnGDeJuxlpvow2wOCKuNn34aXnllgz6hgoJQoqK0NHxXMhBJLJmEsJ27/wsoBXD3dcD6TbmpmZ1jZgvMrNTMqqyvIVI2drA+7l/eQcykiHz68Fc+P/gPsGgRnH9+aAqISLUlkxB+MLMWxGYWmdkhwLebeN/5wJnApE28jmSBpGYFVSF+7KARP/JPejGdg2nJCt7s+RI5056D7bdPZdgi9U4ys4yuAUYBu5vZVKAlcPam3NTdFwGYPsnVeUnNCkpC2djBEUxiMF3Yiw8YRBeupy8r7/9VaoMWqaeSmWU0C+gIHAZcCrRx97npDqyMmXU1syIzK1qxYkVt3VZSpNJZQdWQ1+o7HqE7k+jIZqzjGMbTlUFsnaNkIJIqFSYEMzuz7As4Ddgb2As4NXasUmY23szmJ/g6vToBuvtAd8939/yWLVtW56VSC6rqDkrJfsNjxjC9pC3deJT7uJr9mMfrHKPpoyIpVlmX0amVPOfAiMou7O7H1igiyRrJdAdt0n7DX34JV18Nw4bRNC+P/1z5Fg8+fgg/fgQ52pVMJOWq3DEtrTc3mwhc6+5JbYOmHdMyS5U7kZF4g/smTapYC+AOzz8PV1wBK1fCTTeFry23TPFvIFI/bPKOaWZ2vrsPM7NrEj3v7jXeQc3MzgAeIgxQv2Jm77r7CTW9nkQjme6gau83/Omn0L07jBwJ+fkwfjy0a5fSuEUkscq6jJrGvm+V4LlNala4+0vAS5tyDYlest1BG2xOUxF3GDIErr027GB2zz1h1fFm2sNJpLZU+H+buz8W+3G8u0+Nf87MDk9rVJIVKt2ruDqWLoVLLoHXXw9lJwYPhj32SGmsIlK1ZBamPZTkMalnEu1VXK06QevXw/33h2J0M2bAY4+FpKBkIBKJysYQDiWsPWhZbhyhOdAw3YFJdkiqOyiRBQugc2eYPh1OOQUefRRatUp5fCKSvMpaCFsAzQhJY6u4r+/YxJXKUo+tWQO33QYHHgj//S888wy8/LKSgUgGqGwM4U3gTTN7IrYHgsimmTEDOnWC+fPhvPOgXz/QYkORjJHMFI4tzWwgkBt/vrv/Ll1BSR1TUgJ/+1sYL9hpJxg1Ck6tbN2jiEQhmYTwPPAoMJhNLHst9dDEidClS+geuvRSuOsu2HrrqKMSkQSSSQjr3H1A2iORuuXbb+H668O0o913D7OHjtbOqyKZLJlppy+bWXcz28nMti37Sntkkr1Gj4Y2bcJ6gmuvhblzlQxEskAyLYSy/ZOvizvmwG6pD0ey2ooVYXXx8OFhbcGIEdChQ9RRiUiSqkwI7r5rbQQiWcwdnn0WrroqdBXdeivccANssUXUkYlINSRVKMbM2gJ5QKOyY+7+VLqCkiyyfDlcdlnoJjr44FCPqE2bqKMSkRqoMiGY2S3AUYSEMAY4CZgCKCHUZ6WlMGgQXHcdrFsH990XWggNtYhdJFslM6h8NnAM8Lm7XwzsD2jeYH22ZAkccwx06wa/+U1YaHb11UoGIlkumYTwo7uXAuvMrDnwBbBLesOSTFNYCLvnrOM6u4ef9tqPNdNnhRbC+PGwm+YXiNQFySSEIjP7FTAImAnMAt5Oa1SSUQoL4YEu8xj+0WH05TrG+vHk+UIKG3cJZU5FpE6o1haaZpYLNHf3uekKqDLaQjMCq1fzwA530P3bO1jJNlzBwzzPOYBtsFWmiGSuTd5CM+5CRyY65u6TahqcZIlp06BzZ3p8u5CnOZ+e9ONrWvz8dEVbaIpIdkpm2mn8grRGQAdC15GK29VVP/wAN98cqpHuvDMXbf8KT35x8kanld8qU0SyW5VjCO5+atzXcUBbYGX6Q5NITJgA++0XKpN26wYLFnDcfSfTpMmGp9Voq0wRyWjJDCqXtxzYN9WBSMS++Sbsa3zssWFj+zffhP79oXnzTd8qU0SyQjJjCA8RahdBSCAHEGYaSV0xcmRYbfzFF9CrF9xyCzRuvMEpNd4qU0SyRjJjCPHTetYBw919aprikdr0xRdhdfFzz8H++4etLNu3jzoqEYlIshvk7BH7+X13X53GeKQ2uIfFBT16wKpVcPvtYe+CzTePOjIRiVCFYwhmtrmZ9QM+Bh4HngCWmtkNsecPqJUIJbU++ghOOQUuuAD23hvefRd691YyEJFKB5XvBZoBue7e3t0PIgwm72ZmA4CXaiNASZHSUhgwIFQiffNNeOABmDwZ9tX8ABEJKusyOhnY0+OWMrv7d2Z2GfAloeqpZIPFi8O+xpMnh1lEAwfCrtrmQkQ2VFkLodQT1LVw9/XACneflr6wJCXWrYO77w4DxvPmwdCh8NprSgYiklBlCWGhmf25/EEzOx9YlL6QJCXmzAkb1vTqBSedBAsXwsUXqxidiFSosi6jy4ERZtaJUKoCIB9oDJyR7sCkhn76KcwauusuaNECXngBzjor6qhEJAtUmBDc/RPgYDP7HVC2J+IYd59QK5FJ9b31FnTuDO+9BxdeGHYx23bbqKMSkSxR5ToEd38deL0WYpGaWrUqTB196CHYZRf4z3/ghBOijkpEskxNahlJJhk3LhSje/BBuPzysJ2lkoGI1EAkCcHM+prZe2Y218xeiu3IJtWxciV06gTHHw9bbhmmlD70EGy1VdSRiUiWiqqFMA5o6+7tgMXAjRHFkZ1GjIC8PHjqKbjxxrDa+Le/TfltCgshNxcaNAjfCwtTfgsRySDJ1DJKOXd/Le7hNODsKOLIOp9/DldcAS++CAccAGPGwIEHpuVWhYXQtSuUlITHxcXhMajqqUhdlQljCJ2AVyt60sy6mlmRmRWtWLGiFsPKIO7w5JOhVTB6NNxxB7zzTtoS4CAtAAANLElEQVSSAYQx6rJkUKakJBwXkbrJEixGTs2FzcYDOyZ4qre7j4yd05uwtuHMRKuiy8vPz/eioqKqTqtbiovh0kth7Fg4/HAYPBj22Sftt23QIOSh8sxCWSQRyR5mNtPd86s6L21dRu5+bGXPm9lFwO+BY5JJBvVOaWnYseyGG8Jf4YcfDpvYNKidRl3r1iEXJTouInVTVLOMTgSuB05z95Kqzq933n8fjjwSrrwyDBbPnx+mlNZSMoCwX7L2URapX6IaQ3gY2AoYZ2bvmtmjEcWRWdauhTvvDMXoFi4M4wavvho2Ma5l2kdZpP6JapbRHlWfVc/Mnh3WFbz7Lpx9dugi2mGHSEPSPsoi9UsmzDKq3376Kawl+M1vwrTSF1+E55+PPBmISP0TSQtBYqZMCcXoFi8OpanvvRe22SbqqESknlILIQrffx8WmB1xBKxZEzatGTpUyUBEIqWEUNvGjoW2bcOU0h49wk5mxx0XdVQiIkoItearr8IeBSeeCE2bwtSp0K8fNGsWdWQiIoASQvq5h13L8vLgmWfgr38NM4oOPTTqyERENqBB5XT67LOwoOyll6B9+zBWsP/+UUclIpKQWgjp4A6PPx5aBa++GvY3njZNyUBEMppaCKn24YehTvT48aH8xKBBsNdeUUclIlIltRBSZf36sI1l27YwfToMGABvvKFkICJZQy2EVFi4ELp0gbffhpNOgsceC5vdi4hkEbUQNsXatXD77WGjmsWLYdgweOUVJQMRyUpqIdTUzJmhGN3cufDHP4buou23jzoqEZEaUwuhun78EXr1gg4dYMUK+Pe/4dlnlQxEJOuphVAdkyaFsYIPPoBLLoG774Zf/SrqqEREUkIthGR89x107w4dO4bZRBMmhN1ilAxEpA5RQqjKmDHQpk2YOXTNNWHM4He/izoqEZGUU0KoyJdfwvnnwymnQPPm8NZbYb+Cpk2jjkxEJC2UEMpzh+eeC2UnnnsObrkFZs2Cgw+OOjIRkbTSoHK8Tz+Fyy6DUaPClpYTJsB++0UdlYhIrVALAUKrYPDg0CoYNw7uuSesOlYyEJF6RC2EpUvDFNLXX4ejjgrF6PbYI+qoRERqXf1tIaxfD/ffH4rRFRWFWUQTJigZiEi9VT9bCPPnQ+fO8M478Pvfh8qkrVpFHZWISKTqVwthzRq49VY46KDQVfTMM2EAOcOSQWEh5OZCgwbhe2Fh1BGJSH1Qf1oIM2aEYnTz58N554UN7lu2jDqqjRQWhv11SkrC4+Li8BigoCC6uESk7qv7LYSSErj2WjjkEFi5El5+OfzVzcBkANC79y/JoExJSTguIpJOdbuFMHFiKEb33//CpZeGvY233jrqqCr10UfVOy4ikip1s4Xw7bchARx9dHj8xhvw6KMZnwwAWreu3nERkVSpewnh5ZfDArPBg0NX0dy5YX1BlujTB5o02fBYkybhuIhIOtWdhLBiRRgsPu00aNECpk2Dvn03/uua4QoKQmXtnBwwC98HDtSAsoikX/aPIbjD8OFw1VVh34Lbbgs7mm2xRdSR1VhBgRKAiNS+7E4Iy5eHYnSjR4dqpEOGhL0LRESk2iLpMjKzf5jZXDN718xeM7NfV+sCpaWh1EReXqhBdP/9MHWqkoGIyCaIagyhr7u3c/cDgNHA35J+5ZIlcMwx0K1b2Oh+3jzo2RMaNkxbsCIi9UEkXUbu/l3cw6aAJ/XC//0vlKTecsswi6hTpzDyKiIimyyyMQQz6wP8GfgWODqpFy1fDqefDv37w6+r18skIiKVM/fkPpxX+8Jm44EdEzzV291Hxp13I9DI3W+p4DpdgVg1H9oC81Mda5bZDvgy6iAipvdA7wHoPYDk34Mcd6+yXk/aEkKyzKw1MMbd2yZxbpG759dCWBlL74HeA9B7AHoPIPXvQVSzjPaMe3g68F4UcYiIyC+iGkP4p5ntDZQCxUC3iOIQEZGYqGYZnVXDlw5MaSDZSe+B3gPQewB6DyDF70HkYwgiIpIZ6k5xOxER2SRZlxA2uexFHWBmfc3svdj78JKZ/SrqmGqbmZ1jZgvMrNTM6s1MEzM70czeN7MlZnZD1PFEwcyGmtkXZlYvp6Cb2S5m9oaZLYz9P9AjVdfOuoTAppS9qDvGAW3dvR2wGLgx4niiMB84E5gUdSC1xcwaAo8AJwF5wJ/MLC/aqCLxBHBi1EFEaB3wF3fPAw4BLk/Vv4OsSwg1LntRh7j7a+6+LvZwGtAqynii4O6L3P39qOOoZR2AJe6+1N3XAM8Spm3XK+4+Cfg66jii4u6fufus2M/fA4uAnVNx7awsf12jshd1VyfguaiDkFqxM/Bx3OPlwMERxSIZwMxygQOB6am4XkYmhKrKXrh7b6B3rOzFFUDCshfZLJnSH2bWm9B8LKzN2GpLsuVPROojM2sGvAj0LNdzUmMZmRDc/dgkTy0ExlAHE0JV74GZXQT8HjjG6+jc4Wr8O6gvPgF2iXvcKnZM6hkz25yQDArdfUSqrpt1YwgqexFmmgDXA6e5e0nU8UitmQHsaWa7mtkWwLnAqIhjklpmZgYMARa5+30pvXa2fbg0sxeBDcpeuHu9+pRkZkuALYGvYoemuXu9Kv9hZmcADwEtgW+Ad939hGijSj8zOxnoBzQEhrp7n4hDqnVmNhw4ilDp83/ALe4+JNKgapGZ/RaYDMwj/B0EuMndx2zytbMtIYiISHpkXZeRiIikhxKCiIgASggiIhKjhCAiIoASgoiIxCghSFqYWSszG2lmH5jZUjN72My2TMN9jjKzw+IedzOzP8d+fsLMzk71PePuNTxWcfbqcseviVWinGtmE8wsp5rXPcrMRic43sTMCs1snpnNN7MpsdWqaWNmq9J5fcksSgiScrGFMyOAf7v7nsCeQGPg7jTc7ijg54Tg7o+6+1NpuM8GzGxH4Dexyrv3l3t6NpAfq0b7Aqn7vXsA/3P3/dy9LdAZWJuia4soIUha/A74yd0fB3D39cDVwJ/NrJmZXWRmD5edbGajzeyo2M8DzKwoVuf91rhzlpnZrWY2K/YJeZ9YYa9uwNWx/TGOMLO/m9m15QMys/Zm9qaZzTSzsWa2U+z4VXGf5p9N8LpGZvZ47J6zzaysmOJrwM5l941/jbu/EbeC/OdqtLFP/hPN7IXYfhaFseRZts/Be2Y2i1DWO5GdiCtV4e7vu/vq2Ov/HfvdFphZ17j4V1nYP2OBmY03sw6xGJaa2Wmxcy6KteYmxlp0CUvBmNl1ZjYj9l7dmugcyW5KCJIObYCZ8QdixbeWAXtU8dre7p4PtAM6mlm7uOe+dPeDgAHAte6+DHgUuN/dD3D3yYkuGKv78hBwtru3B4YCZSt8bwAOjH2aT7Ta+/IQvu8H/Al40swaAacB/63svjGdgVfjHh8I9CTsZ7AbcHjseoOAU4H2JC7oRyzuXmb2tpndXq6MS6fY75YPXGVmLWLHmwKvu3sb4HvgduA44AzgtrjXdwDOIrzv51i5TYfM7HhCS68DcADQ3syOrOT3liykhCCZ5v9in5JnExJL/MYfZUW8ZgK51bjm3kBbYJyZvQv8lV/2kJgLFJrZ+YTKseX9FhgG4O7vEcql7JXMTWPXzAf6xh1+x92Xu3sp8G7s99gH+NDdP4gVKhyW6Hru/i4hifQFtgVmmNm+saevMrM5hBbJLoQ/3gBrgP/Efp4HvOnua2M/58Zdfpy7f+XuPxLe59+Wu/3xsa/ZwKxYzHsidUpGVjuVrLcQ2GAw18yaEz75vk/44xz/YaRR7JxdgWsJffMrzeyJsudiVse+r6d6/3YNWODuhyZ47hTgSMKn895mtl/c5kM1ZmbHAr2BjmXdOjHxP1f398DdVxH+YI8ws1LgZDPbATgWONTdS8xsIr+8b2vjquGWlt3f3UvNLP7e5WvYlH9swJ3u/lh14pXsohaCpMMEoEncbJ+GwL3Aw7FPoMuAA8ysgZntQuiGAGgO/AB8G/sjd1IS9/oe2KqKc94HWprZobF4NjezNmbWANjF3d8AegFbA+Vn7UwGCmKv2wtoHbtehczsQOAxQjXaL5L4Hd4Dcs1s99jjP1Vw3cPNbJvYz1sQWk/FsbhXxpLBPoRtFavrODPb1swaA38AppZ7fizQqWxWk5ntbGbb1+A+ksHUQpCUc3e3UI30ETO7mVCR9Lm4ypxTgQ8JLYlFhC4I3H2Omc0m/IH8mI3/KCXyMvCCmZ0OXFlBPGssTD990My2Jvy770fYj3pY7JgBD7r7N+Ve3h8YYGbzCF1KF7n76thYcEX6EhLL87HzPnL30yo62d1/ig0Ev2JmJYQklCjJ7R6LxQgf5l4h1MTfAuhmZosIyWpaZcFV4J3YtVoBw9y9qFyMr8W6p96O/U6rgPOBZBKeZAlVO5W0s7BOYDhwRtlesJI5LGy2lO/uV0Qdi0RLLQRJO3d/C6jW4iwRqX1qIYiICKBBZRERiVFCEBERQAlBRERilBBERARQQhARkRglBBERAeD/ATdvp7xXEse2AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "x = np.random.normal(loc=8.25, scale=2.75, size=37)\n", "y = np.random.normal(loc=8.75, scale=3.25, size=37)\n", "pp_x = sm.ProbPlot(x, fit=True)\n", "pp_y = sm.ProbPlot(y, fit=True)\n", "fig = pp_x.qqplot(line='45', other=pp_y)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VdW9//H3NwMJQ5gJMpogYYggCmFQqWNVEJVer71OdUSp/akd7nWsWtF6q1XvbW2rUi7lWpwt0hYoFUfUq6gEgUASAiEMCVNCgITMw1m/P3K0KSoc4IR9zj6f1/OcJ2evvZLzXZF8nu3ae69tzjlERMRf4rwuQEREwk/hLiLiQwp3EREfUriLiPiQwl1ExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHwowasP7tmzp0tLS/Pq40VEotKKFSt2O+d6HaqfZ+GelpZGdna2Vx8vIhKVzGxLKP00LSMi4kMKdxERH1K4i4j4kMJdRMSHFO4iIj50yHA3szlmVmpma79hv5nZb8ys0MxyzGx0+MsUEZHDEcqR+3PApIPsnwxkBF/TgWePviwRETkah7zO3Tn3gZmlHaTLVGCua3le3ydm1tXM+jjndoSpRhGRsAoEHDWNzdQ2NFPX+MUrQF3TP9oamgM0B9yXr4BzNAegORBsdy0/pym474tHljoHLvgVwOG+bPuiQ1Zad84Ycsj7kI5KOG5i6gcUt9ouCbZ9JdzNbDotR/cMHDgwDB8tItKioSnAjopadlbUsbOyjl2VdeysqKd0fx0VtY3/9KqsbSTg4eOjbznzhKgI95A552YBswCysrL0ZG4ROWy1Dc3k7aggf8d+isqq2bS7iqLd1RTvqflKYHdoF09qShJdO7SjW4d2pPXoSNcOiXRpn0jn5ESS28WTnBBHcmI8yYnxtE+MJzmxZbtdQhxxZiTEGfFxRlxcy/s4a9mONyM+vuVrXBwYhlnL5xpgZgQ3MWvZPpbCEe7bgAGttvsH20REjkog4Fi3cz+fbSonZ1sFa7dVUFha9WWIt0+MJ71nR0b068Ilo/oyoHsH+nRJ5rjOyfTukkxKUsIxD9VIEY5wXwDcZmavAOOBCs23i8iRKiqrYmlBGZ8UlfPppj1U1DYC0CsliZH9ujBpRB9G9O3Mif260KdzMnFxsRneh3LIcDezl4GzgJ5mVgI8CCQCOOdmAouBC4FCoAa4oa2KFRH/aQ44sjfv4Z11pbydv4uismoABnbvwAUn9mbCoB6MH9SDfl3be1xpdAnlapkrD7HfAbeGrSIR8T3nHDklFfxl1TYW5eygbH89ifHGhEE9uO7UNM4ZlsqA7h28LjOqebbkr4jEntLKOl5dXsz8ldvYtLuadvFxnD2sFxeP6suZQ3qRkpzodYm+oXAXkTblnOPjjeW8+OkW3szdRVPAceqgHtxy5iAmjehDl/YK9LagcBeRNlHb0MyfVhTz3EebKdpdTbcOidw4MZ0rxw0kvWdHr8vzPYW7iIRVeVU9c5dtYe6yzeytaeSUgV351eWjmDyiD8mJ8V6XFzMU7iISFtv21fL79zfyWnYxdY0BzsvszffPGERWWnevS4tJCncROSo7Kmp55r2NvLJ8KwCXntKfm88YxODUTh5XFtsU7iJyRHZV1vHs0o289OlWHI5/yxrArWcPpq+uR48ICncROSwVNY08vbSQP368maaA47tj+nPr2YN1XXqEUbiLSEjqm5p5ftkWfvtuIZV1jfzLKf340bkZHN9DV75EIoW7iByUc45FOTt4fMk6ivfU8q2Mntw7eTiZfTt7XZochMJdRL7Rp0Xl/GJxPqtLKhh2XApzbxzX5uuQS3go3EXkK4r31PDI3/JYkruL4zon88RlJ3Hp6P7EawXGqKFwF5Ev1TY08+z7G/n9+xuJM+OO84cwbeIg2rfTzUfRRuEuIjjneGPtTh75Wz7b9tVy8ai+/PTCYfTpossao5XCXSTGbdi1nxkLc/mosJxhx6XwyvQJTBjUw+uy5Cgp3EViVGVdI0+9vYE/fryZDu3ieeiSE7l6/EAS4uO8Lk3CQOEuEmMCAcfrn5fwyzcKKK+u54qxA7jj/KH06JTkdWkSRgp3kRiSU7KPBxfksnLrPk4Z2JU512dxUv+uXpclbUDhLhIDyqvqeWJJAa9mF9OjYxJPfncUl57STw+X9jGFu4iPNTUHeOGTLfz3W+upaWhm2unp/PDbGXTW4+x8T+Eu4lOfFJUzY0Eu63buZ+Lgnsy4JJPBqSlelyXHiMJdxGd2VNTyn3/LZ1HODvp1bc/M743mghOPw0xTMLFE4S7iE/VNzcz+cBO/e7eQgHP86NwMbjnzBN1dGqMU7iI+8O66XTy0MI8t5TVccGJv7p+SqfXVY5zCXSSKFe+p4aGFubydX8oJvTry/LRxfCtDqzaKwl0kKtU1NjPz/Y08u3Qj8XHGvZOHccPp6bRL0N2l0kLhLhJl3l23ixkL8ti6p4aLTurDfVOGa4Ev+QqFu0iUaD0FMzi1Ey/eNJ7TB/f0uiyJUAp3kQinKRg5Egp3kQj2Tn7LVTCagpHDFVK4m9kk4CkgHpjtnHvsgP1dgBeAgcGf+aRz7n/DXKtIzNhaXsPDi/4xBfPSTeM5TVMwchgOGe5mFg88DZwHlADLzWyBcy6vVbdbgTzn3MVm1gsoMLMXnXMNbVK1iE99MQXzzNKNJGgKRo5CKEfu44BC51wRgJm9AkwFWoe7A1Ks5f7mTsAeoCnMtYr42jv5u5ixMJfiPbWagpGjFkq49wOKW22XAOMP6PM7YAGwHUgBLnfOBcJSoYjPbd9Xy4wFubyZt0tTMBI24TqhegGwCjgHOAF4y8w+dM5Vtu5kZtOB6QADBw4M00eLRKem5gDPfbyZ/35rPQHnuHvSMKZN1BSMhEco4b4NGNBqu3+wrbUbgMeccw4oNLNNwDDgs9adnHOzgFkAWVlZ7kiLFol2q4r38dP5a8jbUcnZQ3vx8NQRWgtGwiqUcF8OZJhZOi2hfgVw1QF9tgLnAh+aWW9gKFAUzkJF/KCyrpEn3ijghU+3kJqSxDNXj2byCC3HK+F3yHB3zjWZ2W3AElouhZzjnMs1s1uC+2cCPweeM7M1gAF3O+d2t2HdIlHFOcfCnB38fFEe5VX1XHdqGv9x/hBS9EQkaSMhzbk75xYDiw9om9nq/Xbg/PCWJuIPW8qruf8va/lww25G9uvCnOvGMrJ/F6/LEp/THaoibaSpOcDs/9vEr95aT2J8HDMuzuSaU9OI10Op5RhQuIu0gdztFdz9eg5rt1VywYm9eeiSERzXJdnrsiSGKNxFwqiusZnfvruBme8X0a1DO569ejSTR/bxuiyJQQp3kTBZvnkPd7+eQ1FZNZeN6c/9U4bTtUM7r8uSGKVwFzlKVfVNPP7GOuYu20L/bu2Ze+M4zhiiR92JtxTuIkfhvYJS7pu/hh2Vddxwehp3nD+Ujkn6sxLv6V+hyBGoqGnkoUW5zP98GxmpnZh3y2mMOb6b12WJfEnhLnKY3iso5Z7Xc9hd1cAPzxnMrecMJikh3uuyRP6Jwl0kRPvrGnlkUT6vZhczpHcnZl+rm5EkcincRULwUeFu7pqXw46KWn5w1gn8+NsZOlqXiKZwFzmI6vomHv17Pi98spVBvToy7wenMXqg5tYl8incRb7Bp0Xl3Dkvh+K9NUybmM6dFwwlOVFH6xIdFO4iB6htaObxJet47uPNDOjWgVenn8q49O5elyVyWBTuIq2sLt7HT15dRdHuaq6ZcDz3TB6m69YlKulfrQgtKzg+u3QjT72zgV4pSbx403hO13NMJYop3CXmbSmv5ievruLzrfu4ZFRffj51BF066CEaEt0U7hKznHP8KbuEhxbmEhdnPHXFyUw9uZ/XZYmEhcJdYlJ5VT33zl/Dm3m7mDCoO//1byfTr2t7r8sSCRuFu8Sc9wpKufNPOVTWNnLfhcOZNjGdOD0dSXxG4S4xo7ahmV8szuf5T7YwtHcKz08bx/A+nb0uS6RNKNwlJuRur+CHL69kY1k1N01M5w7dkCQ+p3AXX3PO8dzHm3l08Tq6dUzUJY4SMxTu4lt7qhu4a95q3s4v5dxhqTzx3VF076jH3klsULiLLy3bWM6PX13J3upGHrw4k+tPS8NMJ00ldijcxVeamgM89c4GfvdeIek9OvKH68Yyop/WXJfYo3AX3yjZW8OPX1lF9pa9fHdMf2ZccqLWhZGYpX/54gtv5u7kjj+tJuDQnaYiKNwlyjU1B3hiSQG//6CIkf268LurTuH4Hh29LkvEcwp3iVqllXXc9tJKPtu8h+9NGMgDF2Xq0XciQQp3iUrLNpZz+8srqa5v4teXn8x3TtE0jEhrCneJKoGAY+YHG3lySQHpPTvy8s3jyeid4nVZIhEnLpROZjbJzArMrNDM7vmGPmeZ2SozyzWz98NbpghU1jUy/flsHn+jgCkn9WXBbRMV7CLf4JBH7mYWDzwNnAeUAMvNbIFzLq9Vn67AM8Ak59xWM0ttq4IlNhWVVXHz3Gy2lNcw4+JMrtNNSSIHFcq0zDig0DlXBGBmrwBTgbxWfa4C5jvntgI450rDXajErqUFpdz+8koS4+N44abxTBjUw+uSRCJeKNMy/YDiVtslwbbWhgDdzGypma0ws2u/7geZ2XQzyzaz7LKysiOrWGKGc47fv7+RG59bTv9uHVhw2+kKdpEQheuEagIwBjgXaA8sM7NPnHPrW3dyzs0CZgFkZWW5MH22+FBdYzP3vJ7DX1ZtZ8rIPjzx3ZPo0E7n/0VCFcpfyzZgQKvt/sG21kqAcudcNVBtZh8Ao4D1iBym8qp6bp6bzedb93HH+UO49ezBml8XOUyhTMssBzLMLN3M2gFXAAsO6PNXYKKZJZhZB2A8kB/eUiUWFJVVcemzH5O7vZJnrh7NbedkKNhFjsAhj9ydc01mdhuwBIgH5jjncs3sluD+mc65fDN7A8gBAsBs59zatixc/Gf55j3cPDebODNeunkCY47v5nVJIlHLnPNm6jsrK8tlZ2d78tkSeRas3s4dr62mf7f2/O8NY7U+jMg3MLMVzrmsQ/XTGSrx3HMfbWLGwjzGpnVj1jVZdNPTkkSOmsJdPOOc4zfvFPKrt9dzfmZvfnPlKXpotUiYKNzFE4GA45G/5TPno0386+j+/PJfR5IQH9JqGCISAoW7HHNNzQHumb+GeStKuOH0NB6YkklcnK6IEQknhbscU43NAW5/aSVv5O7k388bwu3n6Bp2kbagcJdjpnWwP3BRJtMmpntdkohvKdzlmGhsDvDDl1uC/WcXZXKjgl2kTekMlrS5xuYAP3plJX9f23LErmAXaXsKd2lTzQHHT15dxeI1O7l/ynBNxYgcIwp3aTPOOe7/y1oW5ezg3snDuOlbg7wuSSRmKNylzTz5ZgEvf7aVW88+ge+feYLX5YjEFIW7tInZHxbx9HsbuXLcQO44f6jX5YjEHIW7hN3rK0p45G/5XDjyOB75zghdxy7iAYW7hNUH68u46/UcJg7uya8uP5l43Xkq4gmFu4TN+l37ufXFz8lI7cTMa8aQlKBFwES8onCXsNhdVc+Nzy0nuV08c64fS6ck3R8n4iWFuxy1usZmbp6bze6qemZfm0Xfru29Lkkk5unwSo6Kc4475+Wwcus+Zn5vNKMGdPW6JBFBR+5ylJ5ZupGFq7dz96RhTBrRx+tyRCRI4S5H7P31ZTz5ZgFTT+7LLWfq7lORSKJwlyNSvKeGH768kqG9U3j00pG6ll0kwijc5bDVNjTz/edX4Jzj99eMoUM7nboRiTT6q5TDdt9f1pC/s5I5143l+B4dvS5HRL6GjtzlsLy+ooT5n2/j9nMyOHtYqtfliMg3ULhLyIrKqnjgr2sZn96dH52b4XU5InIQCncJSX1TM7e9tJKkhDh+fYXWjBGJdJpzl5A8ungdeTsqmX1tFn266A5UkUinI3c5pLfydvHcx5u54fQ0vp3Z2+tyRCQECnc5qNLKOu6at5oT+3bmnsnDvC5HREKkcJdv5Jzj3vlrqGlo5jdXnqIlfEWiSEjhbmaTzKzAzArN7J6D9BtrZk1mdln4ShSv/GlFCe+sK+XuScM4oVcnr8sRkcNwyHA3s3jgaWAykAlcaWaZ39Dvl8Cb4S5Sjr2SvTU8vDCPCYO6c/1paV6XIyKHKZQj93FAoXOuyDnXALwCTP2afrcDrwOlYaxPPBAIOO6al4NzjicuG0WcLnsUiTqhhHs/oLjVdkmw7Utm1g/4F+DZ8JUmXnn+ky18vLGcBy7KZED3Dl6XIyJHIFwnVH8N3O2cCxysk5lNN7NsM8suKysL00dLOG0pr+bRv+dz1tBeXD52gNfliMgRCuUmpm1A67/y/sG21rKAV4LLvvYELjSzJufcX1p3cs7NAmYBZGVluSMtWtqGc46f/nkNiXFxPHbpSVrGVySKhRLuy4EMM0unJdSvAK5q3cE5l/7FezN7Dlh0YLBL5Ju3ooSPCst55DsjOK5LstfliMhROGS4O+eazOw2YAkQD8xxzuWa2S3B/TPbuEY5BnZX1fOfi/MZm9aNq8YN9LocETlKIa0t45xbDCw+oO1rQ905d/3RlyXH2sML86ipb+bRS0fq6hgRH9AdqsJ7BaUsWL2d/3f2CQxOTfG6HBEJA4V7jKuub+L+P69lcGonfnDWCV6XIyJhoiV/Y9xv3y1k275a5t1yqtaOEfERHbnHsE27q/nD/xXxr6P7k5XW3etyRCSMFO4x7JFFebSLj+PuSUO9LkVEwkzhHqOWFpTyzrpSbj83g9TOuqZdxG8U7jGooSnAw4vySOvRgRtOT/O6HBFpAwr3GDR32WaKyqp54KJMnUQV8SmFe4wp21/PU29v4MwhvThnWKrX5YhIG1G4x5gnlxRQ29jMAxdlamEwER9TuMeQNSUVvLaimOtPS2Nwqh6bJ+JnCvcY4ZxjxsJcenRsxw+/neF1OSLSxhTuMeKvq7azYste7rxgKJ2TE70uR0TamMI9BlTXN/Ho3/MZ2a8L3x2jpyuJxAKtLRMDnllayK7Kep65erSW8xWJETpy97mt5TX8z4eb+M7JfRlzvNaPEYkVCnefe+RveSTEGfdMHu51KSJyDCncfez/Nuzmzbxd3Hr2YD0TVSTGKNx9qrE5wEMLcxnYvQPTJqYf+htExFcU7j71widb2FBaxX1ThpOcqPVjRGKNwt2Hyqvq+dVb65k4uCfnZ/b2uhwR8YDC3Yf+6631VDc08+DFWj9GJFYp3H0md3sFL3+2lWsmHE9G7xSvyxERjyjcfcQ5x0ML8ujaPpGffHuI1+WIiIcU7j6yKGcHn23ew50XDKNLB60fIxLLFO4+UdvQzKOL8zmxb2cuH6v1Y0RincLdJ559fyPbK+p48OITidf6MSIxT+HuAyV7a/j9+xu5eFRfxqVr/RgRUbj7wi8W52MG904e5nUpIhIhFO5R7uONu1m8Zif/76zB9O3a3utyRCRCKNyjWFNzgIcX5tG/W3umnzHI63JEJIKEFO5mNsnMCsys0Mzu+Zr9V5tZjpmtMbOPzWxU+EuVA7346VbW7dzPfRdq/RgR+WeHDHcziweeBiYDmcCVZpZ5QLdNwJnOuZHAz4FZ4S5U/tnuqnqefLOAb2X0ZNKI47wuR0QiTChH7uOAQudckXOuAXgFmNq6g3PuY+fc3uDmJ0D/8JYpB3r8jXXUNjTz4MUnav0YEfmKUMK9H1Dcarsk2PZNpgF//7odZjbdzLLNLLusrCz0KuWfrNy6l9eyS5g2MZ3BqZ28LkdEIlBYT6ia2dm0hPvdX7ffOTfLOZflnMvq1atXOD86ZjQHHD/7ay6pKUncfm6G1+WISIQKJdy3Aa3vZ+8fbPsnZnYSMBuY6pwrD095cqDXsotZs62C+6YMp1NSgtfliEiECiXclwMZZpZuZu2AK4AFrTuY2UBgPnCNc259+MsUgH01DTz+xjrGpXfnklF9vS5HRCLYIQ/9nHNNZnYbsASIB+Y453LN7Jbg/pnAz4AewDPBk3tNzrmstis7Nj35ZgGVdU08dIlOoorIwYX0//XOucXA4gPaZrZ6fxNwU3hLk9ZWbNnLi59u5frT0hjep7PX5YhIhNMdqlGgoSnAT+evoU/nZP7j/KFelyMiUUBn5KLA/3xYRMGu/cy+NksnUUUkJDpyj3Cbdlfz1DsbmDKyD9/O7O11OSISJRTuEcw5x0/nryEpIY4HLz5wxQcRkW+mcI9gr2UXs6yonHsnDye1c7LX5YhIFFG4R6jiPTU8vDCPCYO6c4WeiSoih0nhHoECAcdd83IAeOKyUcTpmagicpgU7hFo7rLNLCsq54GLMhnQvYPX5YhIFFK4R5iisioee2MdZw3txeWajhGRI6RwjyB1jc3c9tJKkhPjeezSk7TEgIgcMd0RE0F+sTifvB2V/OG6LI7roqtjROTI6cg9Qryxdgdzl23hponpnDtcNyuJyNFRuEeA4j013Dkvh1H9u3DXpGFelyMiPqBw91h1fRM3z83GgN9eOZp2CfpPIiJHT3PuHgoEHP/+2irW79rPczeMY2APXfYoIuGhw0QPPfXOBpbk7uKnFw7njCF6pqyIhI/C3SN/XbWNp97ZwGVj+jNtYrrX5YiIzyjcPfD++jL+47XVjE/vziPfGaHr2UUk7BTux9jKrXv5wQsryOidwv9cl0VyYrzXJYmIDyncj6GVW/dy7ZzP6NkpiT/eOJbOyYlelyQiPqVwP0ZWbNnDNX/4jG4d2vHy9AmkpugOVBFpOwr3Y2BpQSnX/uEzenZqx6vfn0C/ru29LklEfE7h3sZe/mwr0/6YzfE9OvLq90+lTxcFu4i0Pd3E1EYamgI89vd1zPloE2cO6cXTV4+mU5J+3SJybCht2sD2fbXc+tLnrNy6j+tPS+P+KcNJiNf/JInIsaNwDyPnHPNWlPDzRXk0Bxy/u+oULjqpr9dliUgMUriHyabd1cxYkMv768sYm9aNxy8bRXrPjl6XJSIxSuF+lPZWN/DbdwuZu2wzSQlxzLg4k2tPTdNDrUXEUwr3I7SjopbZH27ipU+3Ut/UzOVjB/CT84bo+nURiQgK98PQ2Bzg3XWl/Cm7mPcKygC4ZFRffnDWCQzpneJxdSIi/xBSuJvZJOApIB6Y7Zx77ID9Ftx/IVADXO+c+zzMtXpiX00DH2zYzbv5u1i6vox9NY2kpiQx/YxBXDVuIAO6aw12EYk8hwx3M4sHngbOA0qA5Wa2wDmX16rbZCAj+BoPPBv8GlX21TRQWFrF+l1V5JTsY8WWvWworQKge8d2nDMslSkj+3DmkF66tFFEIlooR+7jgELnXBGAmb0CTAVah/tUYK5zzgGfmFlXM+vjnNsR9opD1BxwNDYHqKpvoqquiar6JvYHv1bWNrJrfx27KurYVVnPzso6SvbWsruq/svv79I+kdEDu3LJqL6cNrgHJw/oRrxOkopIlAgl3PsBxa22S/jqUfnX9ekHhD3clxaU8vNFeTQFHE3NLQHe8j7wj7ZAAOcO/bNSkhM4rnMyvTsnc/bQXmT07kRGagqDUzvRr2t7XfEiIlHrmJ5QNbPpwHSAgQMHHtHPSElOZNhxnUmINxLi4kiMtwPex5EY1/I1Id5ISUqgU3ICnZIS6ZSUQEpyAp2SEkjtnESHdjqfLCL+FEq6bQMGtNruH2w73D4452YBswCysrJCOLb+qjHHd2PM8d2O5FtFRGJGKGcFlwMZZpZuZu2AK4AFB/RZAFxrLSYAFV7Ot4uIxLpDHrk755rM7DZgCS2XQs5xzuWa2S3B/TOBxbRcBllIy6WQN7RdySIicighTTo75xbTEuCt22a2eu+AW8NbmoiIHCldrC0i4kMKdxERH1K4i4j4kMJdRMSHFO4iIj5kLpT79Nvig83KgC1H+O09gd1hLCfSabz+FUtjhdgab1uN9XjnXK9DdfIs3I+GmWU757K8ruNY0Xj9K5bGCrE1Xq/HqmkZEREfUriLiPhQtIb7LK8LOMY0Xv+KpbFCbI3X07FG5Zy7iIgcXLQeuYuIyEFEXbib2SQzKzCzQjO7x+t6wsHM5phZqZmtbdXW3czeMrMNwa/dWu27Nzj+AjO7wJuqj4yZDTCz98wsz8xyzexHwXbfjdfMks3sMzNbHRzrQ8F23421NTOLN7OVZrYouO3L8ZrZZjNbY2arzCw72BY5Y3XORc2LliWHNwKDgHbAaiDT67rCMK4zgNHA2lZtjwP3BN/fA/wy+D4zOO4kID34+4j3egyHMdY+wOjg+xRgfXBMvhsvYECn4PtE4FNggh/HesC4/x14CVgU3PbleIHNQM8D2iJmrNF25P7lw7qdcw3AFw/rjmrOuQ+APQc0TwX+GHz/R+A7rdpfcc7VO+c20bKG/rhjUmgYOOd2OOc+D77fD+TT8rxd343XtagKbiYGXw4fjvULZtYfmALMbtXs2/F+jYgZa7SF+zc9iNuPert/PM1qJ9A7+N43vwMzSwNOoeWI1pfjDU5RrAJKgbecc74da9CvgbuAQKs2v47XAW+b2Yrg86EhgsaqJ0RHAeecMzNfXdZkZp2A14EfO+cqzezLfX4ar3OuGTjZzLoCfzazEQfs981YzewioNQ5t8LMzvq6Pn4aLzDRObfNzFKBt8xsXeudXo812o7cQ3oQt0/sMrM+AMGvpcH2qP8dmFkiLcH+onNufrDZt+MFcM7tA94DJuHfsZ4OXGJmm2mZMj3HzF7Ap+N1zm0Lfi0F/kzLNEvEjDXawj2Uh3X7xQLguuD764C/tmq/wsySzCwdyAA+86C+I2Ith+h/APKdc//dapfvxmtmvYJH7JhZe+A8YB0+HCuAc+5e51x/51waLX+b7zrnvocPx2tmHc0s5Yv3wPnAWiJprF6fcT6CM9QX0nKFxUbgPq/rCdOYXgZ2AI20zMVNA3oA7wAbgLeB7q363xccfwEw2ev6D3OsE2mZq8wBVgVfF/pxvMBJwMrgWNcCPwu2+26sXzP2s/jH1TK+Gy8tV+ytDr5yv8iiSBqr7lAVEfGhaJuWERGRECjcRUR8SOEuIuJDCncRER9SuIuI+JDCXUTEhxTuIiI+pHAXEfGh/w/lz1YuAAAABElEQVSNLOsysRpNwgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "nobs = 300\n", "np.random.seed(1234) # Seed random generator\n", "dens = sm.nonparametric.KDEUnivariate(np.random.beta(0.5,1.0,size=nobs))\n", "dens.fit()\n", "plt.plot(dens.cdf)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAH39JREFUeJzt3Xl8VPW9//HXJ5ONEEKQhLALQgARATWKdcVaBNRWe61LtdeqtdaqrV0t1tb2tvX3u1Zbb+9DKqXWaitqN61oEdxwa5WyKEsggbAHCElYAmSd5Xv/yGinEcwAk5yZM+/n45HHnPM9J5nPN4/Jm8NZvl9zziEiIv6S4XUBIiKSeAp3EREfUriLiPiQwl1ExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHxI4S4i4kOZXr1xUVGRGzZsmFdvLyKSkpYuXVrvnCvubD/Pwn3YsGEsWbLEq7cXEUlJZrY5nv10WkZExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHyo03A3s0fMrNbMVh1iu5nZ/5pZlZmtMLOTE1+miIgcjniO3B8Fpn3E9ulAafTrJuChoy9LRESORqf3uTvn3jCzYR+xyyXA71z7fH3vmFmhmQ1wzu1IUI0iIkfNOUdbOEJzW5im6Ff7coi2cIRQ2BGKOELhCMHoayjsCEair+EIkei0pM6B++DV8f5spc65Q297/xuBsmHHcM6oTp9DOiqJeIhpELA1Zr062vahcDezm2g/umfo0KEJeGsRSVctwTA1DS3U7Gth574Wave1sqepjT1NQRqa29jTGGRvc5CGpjb2t4RoCoYJR7yfM9oMbj53REqEe9ycc7OB2QBlZWXe/5ZFJKk1tYVYu/MAG+sPsLG+iY31jWysP0D1nmb2NgU/tH9mhlGYl01hXhaFPbIYVJjLCQMLyM/JpGdOgLzsTHpkBcjLDtAj+1/r2ZkZZAaMrIzoa8DI/GA5g8wMIzOQQSDDADDaQ9owrL3p39bbt9u/9nt/p26UiHDfBgyJWR8cbRMRiVtbKMKK6r0s27KH8u37WLWtgQ31jR+c1sgwGNSnB8OL8pk4pJABvXtQUpBL/4Jc+vfOoV9BLr1yMj0J0mSUiHCfC9xmZk8Bk4AGnW8Xkc6EI453t+zhrap6Fm3YzbIte2gNRQAY0DuXEwb25uLxAxk7sIARxT0ZckweOZkBj6tOHZ2Gu5k9CUwGisysGvgBkAXgnJsFzAMuBKqAJuD6ripWRFLbgdYQb66t4+U1tSysrGV3YxtmcHz/Aq6eNJRJw/tSNqwPRfk5Xpea8uK5W+aznWx3wK0Jq0hEfCUYjvDmujqeeXc7L62uoSUYoXePLM4bXcz5x5dwTmkxvfOyvC7Tdzwb8ldE/G3dzv3MWbSF55ZvZ1djG33ysrj8lCFceOIATh3Wh8yAHpDvSgp3EUmYYDjCS6t38ru3N/HOht1kBzKYMraES08axLmjisnOVKB3F4W7iBy1A60hnli0mUfe2kTNvhYGFfbgO9PGcEXZYPrq/LknFO4icsR2N7bx6N838tjbm2loDnLGiL785NJxnDem3wf3hIs3FO4ictjqD7Ty0GvreWLRFpqDYaaeUMItk0cyYUih16VJlMJdROK2vyXIr9/cyMNvbqA1FOGSiQO5ZfIIRvbr5XVp0oHCXUQ61RoK8/u3NzNzYRV7moJcNH4A35wyiuOK870uTQ5B4S4iH2lhZS0/em41G+sbObu0iDumjuHEwb29Lks6oXAXkYPaVN/Ij59fzSsVtRxX1JPfXn8q543u53VZEieFu4j8m7ZQhFmvr+fBV6vIChh3Th/D9WcO1z3qKUbhLiIfWFG9lzv+vIKKmv1cPH4Ad188ln4FuV6XJUdA4S4itIUiPPDyWn71+nqK8nP49bVlTBlb4nVZchQU7iJpbkPdAW5/6j1WbmvgirLB3HXRWHr30EBeqU7hLpKmnHP8aWk1P5xbTnZmBrM+dwrTxvX3uixJEIW7SBpqCYb57jMreXrZNk4/7hgeuHIiA3r38LosSSCFu0ia2ba3mZt/v5SV2xr42idK+crHSzUOjA8p3EXSyNvrd3HbE8toC0V4+NoyPqGLpr6lcBdJE39aspU7n17JsX3zmH1tGSM0dICvKdxFfM45x8yFVdz/4lrOHNmXhz53CgW5uhvG7xTuIj4WCke4e245TyzawqdPGsS9l43Xk6ZpQuEu4lNtoQhfffJd5pfX8OXJI7hj6mjMdOE0XSjcRXyoNRTmlseX8UpFLXdfPJYbzhrudUnSzRTuIj7TEgzzpd8v5fW1dfzk0nF87vRjvS5JPKBwF/GRlmCYGx9bwt/X13PvZSdy5alDvS5JPKJwF/GJYDjCLXOW8ff19dz3mQl85pTBXpckHtJlcxEfCEcc3/jjcl6tqOWeS09UsIvCXSTVOef4/rOreG75dmZMH8PVk3QqRhTuIinvvgWVPLFoC7dMHsHN547wuhxJEgp3kRQ2Z9Fmfvnaeq6eNJRvTx3tdTmSRBTuIilqYWUtdz9bznmji/nRp07QA0ryb+IKdzObZmaVZlZlZjMOsr23mT1nZsvNrNzMrk98qSLyvvLtDdw2Zxlj+vfiwatPJjOg4zT5d51+IswsAMwEpgNjgc+a2dgOu90KrHbOTQAmAz8zs+wE1yoiwI6GZm54dDEFPbJ45LpT6ZmjO5rlw+L55/40oMo5t8E51wY8BVzSYR8H9LL2/xfmA7uBUEIrFZEPnj490BLit9efSklBrtclSZKKJ9wHAVtj1qujbbEeBI4HtgMrgdudc5GEVCgiQPstj3c9s4oV1Q08cOVExvQv8LokSWKJOlE3FXgPGAhMBB40sw998szsJjNbYmZL6urqEvTWIunht3/fxF+WVfO1T5RywQmayFo+Wjzhvg0YErM+ONoW63rgadeuCtgIjOn4g5xzs51zZc65suLi4iOtWSTt/KOqnnvmreGCsSV89eOlXpcjKSCecF8MlJrZ8OhF0quAuR322QKcD2BmJcBoYEMiCxVJV1t3N3HrE8s4rqgnP79yIhmazFri0OlldudcyMxuAxYAAeAR51y5md0c3T4L+DHwqJmtBAz4jnOuvgvrFkkLraEwtz6xjFDEMfvaMvJ1Z4zEKa5PinNuHjCvQ9usmOXtwAWJLU1E/v+8ClZUN/Cr/zyF4UU9vS5HUoiefBBJUi+s3MGj/9jEDWcOZ6ouoMphUriLJKHNuxq5488rmDCkkBnTP3RvgkinFO4iSaYl2H6e3QxmXn0S2Zn6M5XDp6szIknmv1+oYNW2ffz62jIG98nzuhxJUTokEEkiCytqPzjPPmVsidflSApTuIskifoDrXz7z8sZ078Xd0zT2OxydHRaRiQJOOf4zp9XsK8lxJwbTyc3K+B1SZLidOQukgTmLNrCKxW1zJg2htH9e3ldjviAwl3EY1W1B/jJ31ZzdmkR150xzOtyxCcU7iIeCoYjfP0P79EjK8D9l0/QuDGSMDrnLuKhX72+npXbGvjlNSdr4g1JKB25i3ikomYfv3hlHReNH8CFJw7wuhzxGYW7iAeC4Qjf+tNyCnKz+NGnTvC6HPEhnZYR8cCvXl/Pqm37eOiak+mbn+N1OeJDOnIX6Waxp2Om63SMdBGFu0g3CoUjfPtPK3Q6RrqcTsuIdKNH/7GJldsamHm1TsdI19KRu0g3qd7TxM9eXMv5Y/px4YmafEO6lsJdpBs457j72XLM4EeXjsNMDytJ11K4i3SDF1bV8GpFLd+YMopBhT28LkfSgMJdpIvtawnyw7nlnDCwQGPHSLfRBVWRLnbf/ErqD7Ty8OfLyAzoeEq6hz5pIl1o6eY9PL5oM58/YxjjBxd6XY6kEYW7SBcJhSPc9cxK+hfk8s0LNLOSdC+Fu0gXmbNoCxU1+7n74rHk5+gMqHQvhbtIF9h1oJWfvVjJmSP7Mm2c7mmX7qdwF+kC979YSVNbmB9+8gTd0y6eULiLJNiK6r08tXgr150xjNISzYcq3lC4iyRQJNL+JGrfnjnc/olSr8uRNKZwF0mgvyyr5r2te5kxfQy9crO8LkfSmMJdJEH2tQS5d34FJw0t5D9OGuR1OZLm4gp3M5tmZpVmVmVmMw6xz2Qze8/Mys3s9cSWKZL8/vfldexqbONHnxpHRoYuooq3Or351swCwExgClANLDazuc651TH7FAK/BKY557aYWb+uKlgkGW2qb+Sxtzdx+SmDOXFwb6/LEYnryP00oMo5t8E51wY8BVzSYZ+rgaedc1sAnHO1iS1TJLndO7+CrEAG39KTqJIk4gn3QcDWmPXqaFusUUAfM3vNzJaa2bUH+0FmdpOZLTGzJXV1dUdWsUiSWbxpNy+squFL54ygX0Gu1+WIAIm7oJoJnAJcBEwFvm9mozru5Jyb7Zwrc86VFRcXJ+itRbwTiTh+8vxqSgpy+OI5w70uR+QD8Qx4sQ0YErM+ONoWqxrY5ZxrBBrN7A1gArA2IVWKJKnnVmxneXUD918+gbxsjR8jySOeI/fFQKmZDTezbOAqYG6HfZ4FzjKzTDPLAyYBaxJbqkhyaQmGufeFCk4YWKBbHyXpdHqo4ZwLmdltwAIgADzinCs3s5uj22c559aY2XxgBRABHnbOrerKwkW89pu3NrK9oYX7r5igWx8l6cT1/0jn3DxgXoe2WR3W7wPuS1xpIsmrbn8rD722nk8cX8IZI4q8LkfkQ/SEqsgReODltbQEw9x54RivSxE5KIW7yGGqqj3AU//cwjWThjKiON/rckQOSuEucph+/lIluVkBvnK+Rn2U5KVwFzkMy7fuZd7KGm48+ziK8nO8LkfkkBTuIofhvgWV9MnL4otn64ElSW4Kd5E4vbWunreq6rn1vJEaq12SnsJdJA7OOX66oIKBvXP53OnHel2OSKcU7iJxmL+qhhXVDXx9yihyswJelyPSKYW7SCdC4Qj3vVhJab98/uPkwV6XIxIXhbtIJ/6yrJoNdY18a+poAhpmQFKEwl3kI7QEwzzw0jomDinkgrElXpcjEjeFu8hH+N3bm6jZ18J3po3BTEftkjoU7iKH0NAcZObC9ZwzqpiPjejrdTkih0XhLnIIv35jAw3NQe6YqnlRJfUo3EUOonZ/C795ayOfnDCQcYN6e12OyGFTuIscxIOvVhEMR/jmlA9NBSySEhTuIh1s2dXEE4u2cOWpQxhW1NPrckSOiMJdpIOfv1RJZsD4qob0lRSmcBeJsXr7Pp5dvp3rzxxOSUGu1+WIHDGFu0iM+1+spFdOJjefM8LrUkSOisJdJOqfG3fzakUtX548kt55GtJXUpvCXYT2IX3vnV9BSUEO150xzOtyRI6awl0EeGVNLUs37+H280fRI1tD+krqU7hL2gtHHPctqGR4UU8uL9OQvuIPCndJe8++t43Knfv55gWjyAroT0L8QZ9kSWutoTA/f2kt4wYVcOG4AV6XI5IwCndJa08s2kL1nmbumDqGDE3EIT6icJe0daA1xIOvVnHGiL6cXVrkdTkiCaVwl7T1mzc3squxjTs0EYf4kMJd0tKuA63MfmM9007oz8QhhV6XI5JwcYW7mU0zs0ozqzKzGR+x36lmFjKzzySuRJHEm7lwPc3BMN/SRBziU52Gu5kFgJnAdGAs8FkzG3uI/e4FXkx0kSKJVL2nicff2czlpwxhZL98r8sR6RLxHLmfBlQ55zY459qAp4BLDrLfV4C/ALUJrE8k4R54aR0Y3P4JDekr/hVPuA8CtsasV0fbPmBmg4BPAw8lrjSRxKus2c/T71Zz3RnDGFjYw+tyRLpMoi6o/g/wHedc5KN2MrObzGyJmS2pq6tL0FuLxO++BZXkZ2fy5XM1pK/4W2Yc+2wDhsSsD462xSoDnoreTlYEXGhmIefcX2N3cs7NBmYDlJWVuSMtWuRILN28m5fX7ORbF4yiT89sr8sR6VLxhPtioNTMhtMe6lcBV8fu4Jwb/v6ymT0KPN8x2EW85Jzj3hcqKcrP4Yazhnf+DSIprtPTMs65EHAbsABYA/zROVduZjeb2c1dXaBIIrxWWcc/N+3m9vNHkpcdzzGNSGqL61PunJsHzOvQNusQ+1539GWJJE4k0j4Rx9Bj8rjy1KFelyPSLfSEqvjeX9/bRkVN+5C+2Zn6yEt60CddfK0lGOb+BZWcOKg3nxw/0OtyRLqNwl187bd/38T2hha+e+HxGtJX0orCXXxrd2Mbv1xYxflj+vGxEX29LkekWyncxbf+95V1NLaFmDF9jNeliHQ7hbv40sb6Rh5/ZzNXnTaU0pJeXpcj0u0U7uJLP51fQXZmBl/T4GCSphTu4jtLN+/mhVU1fOmcEfTrlet1OSKeULiLrzjnuOdva+jXK4cvnqNhBiR9KdzFV+avqmHZlr18Y8ooDTMgaU3hLr7RFopw7/wKRpf04vKyIZ1/g4iPKdzFN3739iY27WpixoVjCOiBJUlzCnfxhfoDrfzilXVMHl3MeaP7eV2OiOcU7uILP3txLc1tYb530YfmbhdJSwp3SXmrt+/jD4u3cO3HhjGyX77X5YgkBYW7pDTnHD96vpzCvGxuP18PLIm8T+EuKW3+qhre2bCbb0wZRe+8LK/LEUkaCndJWS3BMPfMW8OY/r246lTd+igSS+EuKes3b22kek8zd188lsyAPsoisfQXISlp+95mZi6sYuoJJZwxssjrckSSjsJdUtKPn19NxDm+f7FufRQ5GIW7pJzXKmt5YVUNX/l4KYP75HldjkhSUrhLSmkJhvnB3HKOK+7JF88+zutyRJKWhs2TlDLr9fVs3tXEnBsnkZ2pYxORQ9Ffh6SMzbsa+eVr6/nkhIGcqYuoIh9J4S4pwTnHD+aWkx3I4HsXHe91OSJJT+EuKeG5FTt4rbKOr08ZRUmBps4T6YzCXZLe7sY2/mtuOROGFHLdGcO8LkckJeiCqiS9Hz+/mobmIHMuO1GTcIjESUfuktQWVtbyzLvbuOW8kYzpX+B1OSIpQ+EuSetAa4i7nl7JyH753HreCK/LEUkpcYW7mU0zs0ozqzKzGQfZfo2ZrTCzlWb2DzObkPhSJd389wtr2LGvhXsvG09OZsDrckRSSqfhbmYBYCYwHRgLfNbMOg7osRE41zl3IvBjYHaiC5X0srCylsff2cINZw7nlGP7eF2OSMqJ58j9NKDKObfBOdcGPAVcEruDc+4fzrk90dV3gMGJLVPSyZ7GNu748wpGleTz7amjvS5HJCXFE+6DgK0x69XRtkP5AvDCwTaY2U1mtsTMltTV1cVfpaQN5xzffWYle5vaeODKieRm6XSMyJFI6AVVMzuP9nD/zsG2O+dmO+fKnHNlxcXFiXxr8Yln3t3GC6tq+MaU0ZwwsLfX5YikrHjuc98GxM5hNjja9m/MbDzwMDDdObcrMeVJOtm6u4kfPFvOqcP6cNM5GvFR5GjEc+S+GCg1s+Fmlg1cBcyN3cHMhgJPA//pnFub+DLF79pCEW57YhkY/PyKiXpYSeQodXrk7pwLmdltwAIgADzinCs3s5uj22cBdwN9gV+aGUDIOVfWdWWL3/y/eWtYXt3ArM+dwpBjNAGHyNGKa/gB59w8YF6HtlkxyzcCNya2NEkXL6zcwaP/2MT1Zw5j2rj+Xpcj4gt6QlU8tbG+kTv+vIIJQwq5c7qG8hVJFIW7eGZfS5AvPLaYzIDx4GdP0sxKIgmkUSHFE+GI46tPvsuWXU08fuMknWcXSTCFu3jip/MreK2yjns+PY7Tj+vrdTkivqP/B0u3e2LRFn71xgb+8/RjuWbSsV6XI+JLCnfpVgvKa/jeX1dy3uhi7v5kx/HnRCRRFO7SbRZv2s1Xn3yX8YMLmXnNyWQF9PET6Sr665JuUb69gS88uphBfXrwyHWnkpetyz0iXUnhLl2ufHsD1zy8iPycTB67/jSO6ZntdUkivqdwly61evs+rnl4EXlZAZ666WO65VGkmyjcpcusqN7LNQ+/Q4+sAE/edDpD+yrYRbqLwl26xBtr67hq9jv0zMnkyS+ezrF9e3pdkkha0VUtSbhn39vGN/+4nNKSXjx2/an0K8j1uiSRtKNwl4SJRBy/eGUdv3hlHacfdwyzry2jIDfL67JE0pLCXRLiQGuIb/7xPRaU7+Sykwdzz6fHaf5TEQ8p3OWoVdXu55Y5y6iqPcD3Lx7LDWcOIzppi4h4ROEuR8w5xx8Wb+WHz5XTMzuT390wibNKi7wuS0RQuMsR2t3YxvefXcXfVuzgzJF9eeCKibpwKpJEFO5yWJxz/G3lDn7wbDkNzUG+PXU0Xz53BBma0FokqSjcJW4797Xw/b+u4sXVOxk/uDdzvjiJMf0LvC5LRA5C4S6dagmG+c1bG5m5sIpwxHHn9DF84azhZGpUR5GkpXCXQ3LOMX9VDffMW0P1nmYuGFvCXRcdr6dNRVKAwl0+xDnH62vreOCltSyvbmB0SS/m3DiJM0fqThiRVKFwlw8453irqp7/eXkdSzfvYVBhD+697EQuO3mwTsGIpBiFu9AaCvPc8h08/OYGKmr2078gl59cOo4ryoaQnalQF0lFCvc0VtPQwp+WbOX372ymdn8ro0t6cd9nxvOpiQPJydTQASKpTOGeZtpCEV6t2MkfFm/l9bV1RBycXVrE/ZdP4OzSIg0bIOITCvc0EApHeGfDbuat2sGCVTXsamyjpCCHL08ewRVlQ3T3i4gPKdx9qqE5yNvr61lYUceLq2vY0xQkLzvAeWP6cdnJgzintFgXSUV8TOHuE/tagqzY2sA/N+3mzXV1LN+6l4iD/JxMzj++H9PHDWDy6GINwyuSJuIKdzObBvwCCAAPO+f+u8N2i26/EGgCrnPOLUtwrUL77Yp1+1tZV3uAdTv3U759H+9u3cv6ugM4BxkGE4YUctt5Izl7VDEThxSSpSN0kbTTabibWQCYCUwBqoHFZjbXObc6ZrfpQGn0axLwUPRVDlMoHGFvc5C6/a1s39vM9oYWtu9tZsfeZrbuaWbdzv3sawl9sH+fvCxOGtqHT00YyElDC5kwpFCzH4lIXEfupwFVzrkNAGb2FHAJEBvulwC/c8454B0zKzSzAc65HQmv2EPhiCMYjhAMRwiF25fbYpaD0ddQJEJLMEJTW5imthDNbWGa2sI0B8MxyyEaW8PsbQ6yt6mNPU1t7G0Ksj8muN+XmWH0753LoMIefHLCQEr75VNa0ovSfvkU98rRHS4i8iHxhPsgYGvMejUfPio/2D6DgISH+2uVtfz4+dU4AAeO9lMV7a/gcO2v7hDt0XX+bT1mv4MshyKOUDhCxB19/RkGedmZ5GYF6JkToLBHFn3ysjmuqCeFedkU5rWvF+XnMLAwl4GFPSjKzyGgIXVF5DB06wVVM7sJuAlg6NChR/QzeuVmtQ8za2DtPzP62nHdYtrAiK4bwCG2Rb8PICPm+wMBIzuQQWZGBlmZRlZGBlkBIzOQQXagvS0zI4OswL/aczMzyMvOpEd2gB7ZAfKy2l9zMjN0pC0iXS6ecN8GDIlZHxxtO9x9cM7NBmYDlJWVHdFx8CnH9uGUY/scybeKiKSNeG6jWAyUmtlwM8sGrgLmdthnLnCttTsdaPDb+XYRkVTS6ZG7cy5kZrcBC2i/FfIR51y5md0c3T4LmEf7bZBVtN8KeX3XlSwiIp2J65y7c24e7QEe2zYrZtkBtya2NBEROVJ6ukVExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHzI2m908eCNzeqAzUf47UVAfQLLSXbqr3+lU18hvfrbVX091jlX3NlOnoX70TCzJc65Mq/r6C7qr3+lU18hvfrrdV91WkZExIcU7iIiPpSq4T7b6wK6mfrrX+nUV0iv/nra15Q85y4iIh8tVY/cRUTkI6RcuJvZNDOrNLMqM5vhdT2JYGaPmFmtma2KaTvGzF4ys3XR1z4x2+6M9r/SzKZ6U/WRMbMhZrbQzFabWbmZ3R5t911/zSzXzP5pZsujff2vaLvv+hrLzAJm9q6ZPR9d92V/zWyTma00s/fMbEm0LXn66pxLmS/ahxxeDxwHZAPLgbFe15WAfp0DnAysimn7KTAjujwDuDe6PDba7xxgePT3EfC6D4fR1wHAydHlXsDaaJ9811/aJ/LKjy5nAYuA0/3Y1w79/gbwBPB8dN2X/QU2AUUd2pKmr6l25P7BZN3OuTbg/cm6U5pz7g1gd4fmS4DHosuPAZfGtD/lnGt1zm2kfQz907ql0ARwzu1wzi2LLu8H1tA+367v+uvaHYiuZkW/HD7s6/vMbDBwEfBwTLNv+3sQSdPXVAv3Q03E7Ucl7l+zWdUAJdFl3/wOzGwYcBLtR7S+7G/0FMV7QC3wknPOt32N+h/gDiAS0+bX/jrgZTNbGp0fGpKor906QbYcGeecMzNf3dZkZvnAX4CvOef2xU4a7qf+OufCwEQzKwSeMbNxHbb7pq9mdjFQ65xbamaTD7aPn/oLnOWc22Zm/YCXzKwidqPXfU21I/e4JuL2iZ1mNgAg+lobbU/534GZZdEe7HOcc09Hm33bXwDn3F5gITAN//b1TOBTZraJ9lOmHzezx/Fpf51z26KvtcAztJ9mSZq+plq4xzNZt1/MBT4fXf488GxM+1VmlmNmw4FS4J8e1HdErP0Q/TfAGufcz2M2+a6/ZlYcPWLHzHoAU4AKfNhXAOfcnc65wc65YbT/bb7qnPscPuyvmfU0s17vLwMXAKtIpr56fcX5CK5QX0j7HRbrgbu8ridBfXoS2AEEaT8X9wWgL/AKsA54GTgmZv+7ov2vBKZ7Xf9h9vUs2s9VrgDei35d6Mf+AuOBd6N9XQXcHW33XV8P0vfJ/OtuGd/1l/Y79pZHv8rfz6Jk6queUBUR8aFUOy0jIiJxULiLiPiQwl1ExIcU7iIiPqRwFxHxIYW7iIgPKdxFRHxI4S4i4kP/Bxe/rHqwrxmaAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VOXd9/HPL5ONhCwkhJAdAkEIW4CIiCxi1QIuuINLtbUWaWttre1T+3i391192tq7vVu7aBWprba1iFaUKkpdWWQNQtghIYQsLNkTkpBkMnM9f2TwHjGQASY5s/zer1deM3POdWZ+FzBfTq5zznXEGINSSqngEWJ1AUoppfqWBr9SSgUZDX6llAoyGvxKKRVkNPiVUirIaPArpVSQ0eBXSqkgo8GvlFJBRoNfKaWCTKjVBXRn4MCBZsiQIVaXoZRSfmPr1q01xpgkT9r6ZPAPGTKEgoICq8tQSim/ISKHPW2rQz1KKRVkPAp+EZktIvtFpFhEHulm/TwR2SEi20WkQESmua0rFZGdp9Z5s3illFLnrsehHhGxAU8BVwEVwBYRWWGM2ePW7H1ghTHGiMg4YBkw0m39LGNMjRfrVkopdZ482eOfDBQbY0qMMR3AUmCeewNjTLP53/mdowGd61kppXyUJ8GfBpS7va5wLfsMEblRRPYBbwH3uq0ywHsislVEFp7pQ0RkoWuYqKC6utqz6pVSSp0zrx3cNcYsN8aMBG4AHndbNc0YkwfMAb4pIjPOsP1iY0y+MSY/KcmjM5KUUkqdB0+CvxLIcHud7lrWLWPMGiBbRAa6Xle6HquA5XQNHSmllLKIJ+fxbwFyRGQoXYG/ALjDvYGIDAcOug7uTgQigFoRiQZCjDEnXM+vBh7zag+U6kVlta0UHK7jSMNJRISMhCguGZpAcmyk1aUpdd56DH5jTKeIPACsAmzA88aY3SKyyLX+GeBm4G4RsQMngfmu/wSSgeUicuqzXjLGvNNLfVHKKxxOw8qdR1my7hCF5Q2fWy8CM0ck8fBVFzE2Pc6CCpW6MOKLN1vPz883euWussLOikb+7/Kd7KxsZOjAaO68JJMZI5LITIgCoLiqmX/vOc5fN5TScNLOopnDePiqEYTa9FpIZS0R2WqMyfekrU9O2aBUXzPG8OyaEv77nX0k9o/gtwvyuG5cKiEh8pl2Y9LiGJMWx9emD+Wnb+3ljx8dZN/RJv5wx0SiI/TrpPyD7qaooNdmd/DNlz7hibf3MWdMCu99dybz8tI+F/ruYiLDeOLmcfzsxrGsKarh/r9upb3T0YdVK3X+NPhVUDvRZuee5zfz9q5j/N+5I/nDHROI6xfm8fZ3XJLJL24ex7riGr77ciFOp+8NnSp1Ov3dVAWt5vZOvvSnzeyqbOTJ+XnMy/vcdYkeuWVSOnUt7fxs5T5GpcTwwBU5Xq5UKe/SPX4VlNrsDha+WMDOykaeunPieYf+KV+bns28vFR+/e4BNhys9VKVSvUODX4VdIwxfO+VQtYfrOWXt4zji6MHX/B7igg/vXEsQxKj+fbSbTSetHuhUqV6hwa/CjpPf3SQN3cc5QezR3LTxHSvvW//iFB+u2ACNc3t/HzlXq+9r1LepsGvgsqH+6r45ar93JCXyqKZ2V5//7HpcXxtejZLt5Sz/qDORK58kwa/ChpVTW08/EohIwfH8MTN43BdUe5137lyBFmJUfzH67uwO5y98hlKXQgNfhUUnE7Dw68U0trRyR/umEBkmK3XPqtfuI0fXZNLSXULf9/o8W1QleozGvwqKLy0uYy1RTX8+NrRDB8U0+uf94VRg7hseCJPvl9EY6se6FW+RYNfBbxjjW088fY+LhueyO2TM3rewAtEhEfn5tJ40s5v3y/qk89UylMa/Crg/eiNXXQ6nfzsxrG9Nq7fndzUWObnZ/DihlIO17b02ecq1RMNfhXQ1hyo5t09x/n2F0aQlRjd55//0FUjsIUIv3u/uM8/W6kz0eBXAavT4eT/vbWHrMQo7p02xJIakmMjuWtKFsu3VVBS3WxJDUqdToNfBaylW8o5cLyZH84ZSURo753F05NFM4cREWrjdzrWr3yEBr8KSG12B799v4iLhwzwypQMFyIpJoK7p2bxRuERio6fsLQWpUCDXwWolzaVUX2ine9edVGfHtA9k/tnDKNfmI3ff6Bj/cp6Gvwq4LTZHfxx9UGmZCdw6bBEq8sBICE6nDsvyeTNHUf0DB9lOQ1+FXD+sblrb/87V46wupTPuG96NqEhITy7psTqUlSQ8yj4RWS2iOwXkWIReaSb9fNEZIeIbBeRAhGZ5um2SnlTp8PJn9Yd4uIhA5iS7Rt7+6ckx0Zy86R0Xi2ooKqpzepyVBDrMfhFxAY8BcwBcoHbRST3tGbvA+ONMXnAvcCSc9hWKa9Ztfs4FfUnuW+692fe9IZFM7PpdHb956SUVTzZ458MFBtjSowxHcBSYJ57A2NMszHm1M1GowHj6bZKedOSdSVkJUZx5ahkq0vpVlZiNNeMS+VvGw/rHD7KMp4EfxpQ7va6wrXsM0TkRhHZB7xF116/x9u6tl/oGiYqqK6u9qR2pT7jk7J6tpU1cO9lQ7GFWH8mz5l8feYwWjocvLih1OpSVJDy2sFdY8xyY8xI4Abg8fPYfrExJt8Yk5+UlOStslQQeWlTGdHhNm6Z5L27avWG3NRYrhg5iOc/PkRrR6fV5agg5EnwVwLuUxqmu5Z1yxizBsgWkYHnuq1S56upzc6bO45wfV4q0RGhVpfTo69fPoz6Vjuvbq2wuhQVhDwJ/i1AjogMFZFwYAGwwr2BiAwX11UyIjIRiABqPdlWKW94Y/sR2uxOFlycaXUpHsnPGsDEzHieW1tCp96lS/WxHoPfGNMJPACsAvYCy4wxu0VkkYgscjW7GdglItvpOotnvunS7ba90REV3F7eUsaolFjGpcdZXYpHRISFM4ZRXneSt3cds7ocFWQ8+p3YGLMSWHnasmfcnv8C+IWn2yrlTbsqG9lV2cRj80b7xPQMnroqN5nsgdEsXlPCteNS/Kp25d/0yl3l95ZuKSMyLIR5ed2eMOazbCHCfdOz2VnZyIaDtVaXo4KIBr/ya3aHkzd3HOXq3MHE9QuzupxzdtPENAb2D9dpHFSf0uBXfm1dUQ0NrXauH59qdSnnJTLMxpenDmH1gWr2Hm2yuhwVJDT4lV9bUXiEuH5hzBjhv9d+3DUli6hwG8/pXr/qIxr8ym+d7HDw793HmDNmMOGh/vtPOT4qnPkXZ7Ci8AhHGk5aXY4KAv77bVFB74N9VbR0OPx2mMfdV6cNxQDP6+Rtqg9o8Cu/taKwkkExEVziY9Mvn4/0AVFcOy6Ff2wuo/GkTt6mepcGv/JLJ9rsfLi/mmvGpfj0hGznYuGMbFo6HPx902GrS1EBToNf+aWP9lfT0elk7tgUq0vxmtGpcUzPGcifPy6lvdNhdTkqgGnwK7/07z3HGdg/nImZA6wuxavunzGM6hPtvL5N5zJUvUeDX/md9k4HH+6r4spRyQEzzHPKZcMTyU2J5dk1JTidpucNlDoPGvzK72wsqaO5vZOrR/vmXbYuhIhw/8xsSqpbeG/vcavLUQFKg1/5nX/vPkZUuI2pwwZaXUqvuGZsCmnx/VisF3SpXqLBr/yK02l4d89xLr8oicgwm9Xl9IpQWwj3TR9KweF6th6us7ocFYA0+JVf2VHZSNWJdq7OHWx1Kb1q/sUZxEeF8exq3etX3qfBr/zK+3uPYwsRZl00yOpSelVUeChfmpLFu3uPc7C62epyVIDR4Fd+ZU1RDXkZ8cRF+d8UzOfqnqlDCLOFsGSt7vUr79LgV36jobWDHRUNTM8JzIO6pxvYP4JbJqXzz62VVJ1os7ocFUA0+JXfWH+wFmNgeo7/TsF8rr42PRu708kL60utLkUFEI+CX0Rmi8h+ESkWkUe6WX+niOwQkZ0isl5ExrutK3Ut3y4iBd4sXgWXtUXVxESGMt5PbqjuDUMHRvPF3MH8dcNhmts7rS5HBYgeg19EbMBTwBwgF7hdRHJPa3YImGmMGQs8Diw+bf0sY0yeMSbfCzWrIGSMYc2BGqYOSyTUFly/qN4/M5umtk5e3lJudSkqQHjyDZoMFBtjSowxHcBSYJ57A2PMemNMvevlRiDdu2WqYFda20plw8mgGuY5ZULmACYPSeBPa0uwO5xWl6MCgCfBnwa472pUuJadyVeBt91eG+A9EdkqIgvPvUSluoZ5AGYEYfBD117/kcY23txxxOpSVADw6u/MIjKLruD/gdviacaYPLqGir4pIjPOsO1CESkQkYLq6mpvlqUCwJoDNWQmRJGZGGV1KZaYddEghg/qz7OrSzBGJ29TF8aT4K8EMtxep7uWfYaIjAOWAPOMMbWnlhtjKl2PVcByuoaOPscYs9gYk2+MyU9KCs69OtU9u8PJxpLaoDmNszshIcLCGdnsO3aCNUU1Vpej/Jwnwb8FyBGRoSISDiwAVrg3EJFM4DXgS8aYA27Lo0Uk5tRz4Gpgl7eKV8Fhe3kDze2dQTm+725eXiqDYiJYvOag1aUoP9dj8BtjOoEHgFXAXmCZMWa3iCwSkUWuZj8GEoGnTzttMxlYJyKFwGbgLWPMO17vhQpoaw9UYwsRLh3m//fWvRARoTbunTaUj4tr2VXZaHU5yo+JL44X5ufnm4ICPeVfdbnhqY8JEXjtG5dZXYrlmtrsTP35B8waOYjf3z7B6nKUDxGRrZ6eMh9cJ0Qrv9PYandN0xDcwzynxEaGccclmby14wjlda1Wl6P8lAa/8mnrD9bgNDBjRPAe2D3dVy4bQogIf1p3yOpSlJ/S4Fc+bU1RDTERoYxPj7e6FJ+REtePeXlpvLylnPqWDqvLUX5Ig1/5rK5pGqq5NAinaejJwhnZnLQ7+OvGw1aXovyQfpuUz/p0moYROr5/uosGxzDroiT+sr6UNrvD6nKUn9HgVz5r3afTNOj4fnfunzmMupYOXtlaYXUpys9o8CuftaaohoyEfmQlRltdik+6ZGgC49PjWLK2BIfT907LVr5Lg1/5JLvDyYaDtXoa51mICPfPHMbh2lZW7T5mdTnKj2jwK59U6JqmQYd5zu6LoweTlRjFs6sP6uRtymMa/MonrSmqIUTg0mEa/GdjCxHum55NYUUjmw7VWV2O8hMa/MonrS2qZnxGPHH9wqwuxefdOimdhOhwFq8psboU5Sc0+JXPaWy1U1iu0zR4KjLMxj2XDuGDfVUcOH7C6nKUH9DgVz5nQ4lrmgYd3/fYly7NIjIsRPf6lUc0+JXPWVNUQ/+IUMZn6DQNnkqIDmd+fgZvbK/kWGOb1eUoH6fBr3yK+zQNYTpNwzm5b3o2DqfhL+tLrS5F+Tj9Zimfcri2lYr6kzrMcx4yEqK4clQyywrKae/UaRzUmWnwK5+ytrjrfrLT9MDueblrShZ1LR28vVMv6FJnpsGvfMraA9WkD+jHkMQoq0vxS9OGDyQrMYq/6ayd6iw0+JXPcJ+mQUSsLscvhYQId12SRcHhevYebbK6HOWjNPiVzygsb+CETtNwwW6ZlE54aIju9asz8ij4RWS2iOwXkWIReaSb9XeKyA4R2Ski60VkvKfbKnXKWtc0DVN1moYLMiA6nGvHpfD6tkqa2zutLkf5oB6DX0RswFPAHCAXuF1Eck9rdgiYaYwZCzwOLD6HbZUCuqZpGJceT1yUTtNwob40JYuWDgfLt1VaXYryQZ7s8U8Gio0xJcaYDmApMM+9gTFmvTGm3vVyI5Du6bZKATSetLO9vEGHebwkLyOe0amx/H3jYZ21U32OJ8GfBpS7va5wLTuTrwJvn+e2KkhtOFiL0+hpnN4iItw1JYt9x07wSVmD1eUoH+PVg7siMouu4P/BeWy7UEQKRKSgurram2UpP7C2qJrocBsTMnWaBm+5bnwqUeE2lm0p77mxCiqeBH8lkOH2Ot217DNEZBywBJhnjKk9l20BjDGLjTH5xpj8pCTd6ws2a4tquHTYQJ2mwYv6R4Ry7bgU/rXjiB7kVZ/hybdsC5AjIkNFJBxYAKxwbyAimcBrwJeMMQfOZVulDte2UFbXyowROr7vbfMvzqC1w8FbO45YXYryIT0GvzGmE3gAWAXsBZYZY3aLyCIRWeRq9mMgEXhaRLaLSMHZtu2Ffig/trbINU3DcA1+b5uYOYBhSdG8rMM9yk2oJ42MMSuBlacte8bt+X3AfZ5uq5S7tUXVpMX3Y+jAaKtLCTgiwoKLM/npyr0UHT9BTnKM1SUpH6ADqspSdoeT9cW1zBgxUKdp6CU3TkwjNER0r199SoNfWWpbWdc0DTNH6AH93jKwfwRXjkrmtW2VdHQ6rS5H+QANfmWp1QeqsIUIU3V8v1fNvziDupYO3t973OpSlA/Q4FeWWn2gmkmZA4iN1GkaetOMEUkMjo3k5QId7lEa/MpCNc3t7Kps0tM4+4AtRLg1P53VB6o50nDS6nKUxTT4lWXWFnVdoT1zxCCLKwkOt+VnYAz8c2uF1aUoi2nwK8us3l9NYnQ4o1NjrS4lKGQkRDF1WCKvbK3A6dSJ24KZBr+yhNNpWFtUw4wRSYSE6GmcfeW2/AzK6lrZeKi258YqYGnwK0vsPtJEbUuHju/3sdljBhMTGcorBTrcE8w0+JUlVh+oAmC6TsPcpyLDbMzLS2XlzqM0tdmtLkdZRINfWWLNgRrGpsUxsH+E1aUEndvyM2jvdLJiu07cFqw0+FWfa2qzs7WsXq/WtcjYtDhGDo7hFT2nP2hp8Ks+93FRDQ6nYYYGvyVEhNvyMyisaGTfsSary1EW0OBXfe69vVXE9Qtjot5tyzI3TEgjzCYs26IHeYORBr/qUw6n4cP9Vcy6KIlQvduWZRKiw7kqN5nl2yp04rYgpN881ae2ldVT19LBlbnJVpcS9G7Lz6C+1a4TtwUhDX7Vp97de5wwm+j4vg+YnqMTtwUrDX7Vp97bc5xLhibqbJw+wBYi3DIpnTUHqjnaqBO3BRMNftVnDtW0cLC6hStH6aRsvuLW/HScBl77pNLqUlQf0uBXfebUWPIXRun4vq/ISoxmSnYCywrKdeK2IOJR8IvIbBHZLyLFIvJIN+tHisgGEWkXke+dtq5URHaKyHYRKfBW4cr/vLvnOCMHx5CREGV1KcrNbfkZHK5tZXNpndWlqD7SY/CLiA14CpgD5AK3i0juac3qgAeBX53hbWYZY/KMMfkXUqzyX/UtHRQcrudK3dv3OXPGpBATEcoyPcgbNDzZ458MFBtjSowxHcBSYJ57A2NMlTFmC6CzPqlufbCvCofT8AUd3/c5/cJtXKcTtwUVT4I/DXDfFahwLfOUAd4Tka0isvBcilOB4+1dR0mJi2R8ul6t64tuy8+gze7kzcKjVpei+kBfHNydZozJo2uo6JsiMqO7RiKyUEQKRKSgurq6D8pSfaWpzc6aAzXMGZOiN13xUePT47goOUbP6Q8SngR/JZDh9jrdtcwjxphK12MVsJyuoaPu2i02xuQbY/KTkvTinkDy/t7jdDicXDNusNWlqDMQEW67OIPC8gZ2VTZaXY7qZZ4E/xYgR0SGikg4sABY4cmbi0i0iMSceg5cDew632KVf3prxzEGx0YyIWOA1aWos7hlYjr9wmz8dcNhq0tRvazH4DfGdAIPAKuAvcAyY8xuEVkkIosARGSwiFQA3wX+Q0QqRCQWSAbWiUghsBl4yxjzTm91RvmeE2121hRVM2fsYB3m8XFxUWHcMCGNNworaWjtsLoc1YtCPWlkjFkJrDxt2TNuz4/RNQR0uiZg/IUUqPzbB/uq6Oh0MndsitWlKA/cfWkW/9hcxisFFXxtRrbV5aheolfuql715o6jJMdGMClTh3n8waiUWCYPSeCvGw/rlbwBTINf9Zr6lg4+2l/FdeNSdZjHj9w9NYuyulZWH9Cz6wKVBr/qNW/uPIrdYbhhwrlc9qGs9sXRgxkUE8ELG0qtLkX1Eg1+1WuWf1LBiOT+jE6NtboUdQ7CbCHceUkWH+2vprjqhNXlqF6gwa96xeHaFj4pa+DGCemI6DCPv7lrSiaRYSE8t+aQ1aWoXqDBr3rF8m2ViMANE1KtLkWdh8T+Edw6KYPl2yo53tRmdTnKyzT4ldcZY3h9WyWXZieSEtfP6nLUebpv+lA6nU7+/HGp1aUoL9PgV163pbSe0tpWbtSDun4tKzGaOWNT+PvGw5zQWTsDiga/8rp/bC4jJjKUa8fpMI+/u39GNifaO3lRp3EIKBr8yqvqWzp4a+dRbpyQRr9wm9XlqAs0Lj2eWRcl8dzaEt3rDyAa/MqrXttWSUenk9snZ1pdivKSh64aQUOrnb/oWH/A0OBXXmOM4R+by5iQGc+oFD13P1CMS4/nylHJPLe2hMaTutcfCDT4lddsKa2nuKpZ9/YD0ENX5dDU1smf1ul5/YFAg195zQsbSl0HdXUmzkAzOjWOuWMHs2RtiZ7XHwA0+JVXVNS38vbOo9wxOZOocI9m+1Z+5gezR9LpMPxq1X6rS1EXSINfecUL60sREe6ZOsTqUlQvyUqM5suXDeHVTyr09ox+ToNfXbDm9k6Wbi5n7tgUUuP1St1A9s1ZwxkQFc7jb+7BGJ2v319p8KsLtmxLOSfaO/nqtKFWl6J6WVy/MB66agSbDtWxovCI1eWo86TBry6I3eHk+Y8PkZ81gLyMeKvLUX3gjsmZjE+P4/E39+i9ef2UBr+6IMu3VVJRf5JvzBpmdSmqj9hChJ/dNJb6VjtPvL3P6nLUefAo+EVktojsF5FiEXmkm/UjRWSDiLSLyPfOZVvlvzodTp7+sJgxabHMumiQ1eWoPjQ6NY77pg1l6ZZyNpXUWl2OOkc9Br+I2ICngDlALnC7iOSe1qwOeBD41Xlsq/zUv3YcobS2lW9dkaM3WwlC374yh8yEKL7/6g6a2zutLkedA0/2+CcDxcaYEmNMB7AUmOfewBhTZYzZApx+PXeP2yr/5HAafv9BMSMHx3DVqGSry1EWiAoP5de3jaeivpXH/rXb6nLUOfAk+NOAcrfXFa5lnvB4WxFZKCIFIlJQXV3t4dsrq6worKSkuoUHrhhOSIju7Qer/CEJfP3yYSwrqGDV7mNWl6M85DMHd40xi40x+caY/KSkJKvLUWfR3ungV6sOMDo1lrljdHqGYPftL4xgTFosP3xtJ1UndDoHf+BJ8FcCGW6v013LPHEh2yof9beNZVQ2nOSROSN1b18RHhrCk/PzaGnv5LsvF+Jw6oVdvs6T4N8C5IjIUBEJBxYAKzx8/wvZVvmgpjY7f/igiOk5A5meo7+ZqS7DB8Xwk+tHs664hj98UGx1OaoHPc6mZYzpFJEHgFWADXjeGLNbRBa51j8jIoOBAiAWcIrId4BcY0xTd9v2VmdU73vmo4PUt9r5weyRVpeifMz8izPYfKiOJ98/QP6QAVw2fKDVJakzEF+cbyM/P98UFBRYXYY6TWlNC1f/Zg3XjEvhN/PzrC5H+aDWjk6u/8PHNLR2sPLB6QyKjbS6pKAhIluNMfmetPWZg7vK9z3+5h7CbMIP5+jevupeVHgof7xzIi3tDr71j210OpxWl6S6ocGvPPLBvuO8v6+KB7+Qo3tx6qxykmP46Y1j2HSojiffK7K6HNUNDX7Vo/ZOB4/9aw/ZSdF85TKdgVP17KaJ6czPz+Cpj4pZfUCvy/E1GvyqR0vWHqK0tpX/um404aH6T0Z55ifzRnNRcgwPvbydo40nrS5HudFvsTqrQzUt/Pb9IuaMGcyMEXr6pvJcZJiNp+6cSJvdwYM63u9TNPjVGTmdhkf+uYOI0BB+cv1oq8tRfmhYUn9+ftNYtpTW86t/H7C6HOWiwa/O6OWCcjYdquPRuaP0gK46b/Py0rh9cibPrD7IB/uOW12OQoNfncHxpjZ+tnIvU7ITmH9xRs8bKHUW/3ldLqNSYvnuskIqG3S832oa/Kpb//nGbjo6nfz8pnE61766YJFhNp6+cyL2TiffeukT7DrebykNfvU5b+44wju7j/HtK3MYOjDa6nJUgBg6MJonbh7HJ2UN/GrVfqvLCWoa/Oozqpra+I/XdzE+I56F07OtLkcFmOvGp3L75EwWry1h/cEaq8sJWhr86lPGGB55bScnOxz8z63jCbXpPw/lfT+6dhRDE6N5eFkhja2n37RP9QX9ZqtPLSso54N9Vfxg9kiGD+pvdTkqQEWFh/LkgjyqT7Tz6Os78cWJIgOdBr8CoLyulcf+tYcp2Ql8eeoQq8tRAW5cejzfuTKHN3cc5fXtem+mvqbBr3A6Dd9/tRAR4Ze3jNe7aqk+8fXLh5OfNYAfv76b8rpWq8sJKhr8ir+sL2VjSR0/unYUGQlRVpejgoQtRPjN/DwM8PAyvWVjX9LgD3LFVc384p19XDFyELfl64Vaqm9lJETxk+tHs7m0jmdWH7S6nKChwR/EOh1OHn6lkH7hNp64aaxeqKUscdPENOaOHcyT7x1gV2Wj1eUEBQ3+IPb7D4opLG/g8XljdC4eZRkR4ac3jGVAVDgPvbydNrvD6pICngZ/kNp6uI7ff1DETRPSuG58qtXlqCA3IDqcX946niLX0KPqXR4Fv4jMFpH9IlIsIo90s15E5Heu9TtEZKLbulIR2Ski20VE76DuA0602fn20u2kDejHT+bpdMvKN8wckcQ9l2bx549LWVekV/X2ph6DX0RswFPAHCAXuF1Eck9rNgfIcf0sBP542vpZxpg8T+8Ar3rXf67YzZGGkzw5P4+YyDCry1HqU4/MGcWwpGi+90ohDa0dVpcTsDzZ458MFBtjSowxHcBSYN5pbeYBL5ouG4F4EUnxcq3KC1YUHuG1Typ54IocJmUlWF2OUp/RL9zGk/MnUNPczo/e2G11OQHLk+BPA8rdXle4lnnaxgDvichWEVl4pg8RkYUiUiAiBdXVenPm3lDZcJJHl+9kQmY8D14x3OpylOrW2PQ4vnNlDv8qPMIbelVvr+iLg7vTjDF5dA0HfVNEZnTXyBiz2BiTb4zJT0rSe7t6m8NpeOjl7Tidhifn5+kEbMqnLZo5jElZA/iP13dxRG/c4nWefPtindG8AAAMEklEQVQrAfcre9JdyzxqY4w59VgFLKdr6Ej1sT9+VMzmQ3X8ZN4YshJ1jn3l20JtIfz6tvE4nYaHlxXi1Kt6vcqT4N8C5IjIUBEJBxYAK05rswK423V2zxSg0RhzVESiRSQGQESigauBXV6sX3lgU0ktv373ANeNT+XmiaeP0inlm7ISo/nxdblsKKllyboSq8sJKKE9NTDGdIrIA8AqwAY8b4zZLSKLXOufAVYCc4FioBX4imvzZGC564rQUOAlY8w7Xu+FOqOa5nYeXLqNrMRofq5X5yo/c1t+Bh/uq+YX7+xnYuYA8ofoCQneIL44F3Z+fr4pKNBT/i+U02m458+b2XSojte/cRm5qbFWl6TUOWtqs3Pd79fRbnfy1oPTSOwfYXVJPklEtnp6yrwe4Qtgf1x9kLVFNfzXdaM19JXfio0M4+k7J1LX2sF3Xt6us3h6gQZ/gFpXVMP//Hu/6x6nOuum8m+jU+N47PrRrHX9u1YXpscxfuV/Dte28M2XPiFnUIzOuqkCxvyLMyisaODpjw6SndSfWyalW12S39LgDzDN7Z187cUCROC5u/OJjtC/YhUYRITH5o2hrK6VH762g4wB/bgkO9HqsvySDvUEkK5znrdTXNXMH26fSGai3k1LBZYwWwhP3zGJzIQo7v/bVkqqm60uyS9p8AcIYwyPv7WHVbuP8+g1uUzLGWh1SUr1irioMP785cnYRLhrySa9X+950OAPEEvWHuLPH5fylcuGcO9lQ6wuR6lelZkYxYtfnUxzeyd3LtnEscY2q0vyKxr8AeDVrRX8dOVerhmXwo+uydWDuSoojE6N48WvXkJdSwcLFm/QPf9zoMHv5/65tYLvv1rItOED+Z9bxxMSoqGvgkdeRjwv3DuZupYObv7jevYebbK6JL+gwe/H/rm1gu+9WsjUYYk8d3c+kWE2q0tSqs9NyhrAq1+fSogItz27gdUHdFr3nmjw+yFjDE9/VMzDr3SF/pK7L6ZfuIa+Cl4jkmP45zemkhbfjy//eTNPfVisM3qehQa/n7E7nDz6+i7++539XD8+lee/rKGvFEBafD9e+8ZUrh+fyi9X7efeF7ZwvEkP+nZHg9+PlNe1cuszG3hpUxlfv3wYT87PIyJUQ1+pU6LCQ3lyfh6PzRvNxpJarv7NGpZvq8AXJ6O0kga/HzDG8Mb2Sq753VoOVjXz1B0T+cHskXogV6luiAh3XzqElQ9OZ1hSNA+9XMhtz25gZ0Wj1aX5DJ2W2ceV1rTwozd2sbaohvEZ8fxuQZ7eQUspDzmchmUF5fxq1X7qWju4aUI6D35heEB+h85lWmYNfh9VUd/KUx8e5JWCciLDbHz/ixdx15QsbLqXr9Q5a2qz8/v3i3hhw2E6HU7m5aXxjcuHkZMcY3VpXqPB76ecTsPGklr+vrmMVbuOESLCgskZPDBrOINiI60uTym/V9XUxnNrS/jbxjJO2h1cmp3InVMyuTp3MOGh/j3yrcHvRzo6nWwprePfu4/x3t4qKhtOEtcvjJsnpnPf9KGkxvezukSlAk5dSwdLt5Tx0qYyKupPMrB/BDdOSOW68amMTYvzy6vfNfh9lDGGY01t7D3axNbD9WwprWdHRQNtdieRYSFMz0li7tjBzBmTohdjKdUHHE7DmqJqXtpUxkf7q7A7DFmJUVwzNoUrRg4iLyOeUJt//Cbg9eAXkdnAb+m62foSY8wTp60X1/q5dN1s/cvGmE882bY7/h787Z0OyutOUl7XSmltC0VVzRw4doL9x09woq0TgNAQYXRqLJOyEpiSncD0nCQ9H18pCzW22lm1+xj/2nGE9QdrcTgNsZGhTM9JYuaIJC4emsCQxCif/W3Aq8EvIjbgAHAVUAFsAW43xuxxazMX+BZdwX8J8FtjzCWebNsdXw/+jk4nVSfaON7UxpGGNsrqWjlc20JZXStlta0cbWrD/Y81NjKUkYNjGTG4Pxclx3DR4FjGpsVp0Cvloxpb7awrrmH1gSpWH6jmeFM7AAOiwsjLiGdC5gBGDo5h+KD+ZCZE+cRvBecS/J7cnmkyUGyMKXG9+VJgHuAe3vOAF03X/yIbRSReRFKAIR5sawljDHaHoa3TQVuHg5P2rp8TbZ3Ut3TQcNJOQ2sHDa126lvtVJ9o41hTG8ca26lpbv/c+yXFRJCZEMWU7EQyE6PITIgiKzGKjIQokvpH+OxeglLq8+KiwrhmXArXjEvBGENRVTOfHK7nk7J6tpU18OH+/50PKNwWwpCBXd/55NhIBsdGkhwXSVL/CGIiQ4mJDCO2X9djdLjNJ7LAk+BPA8rdXlfQtVffU5s0D7f1mmt+t5aTdgfGdI3dOY3BGHAa43rd9byj08lJuwOHB3N5hIYI8VFhDOwfQUpcJGPT4j79yx0cF0lKXD8yEvoRFa63OFQqEIkII5JjGJEcw4LJmUDX6aHFVc0UVzVzsLqZg1XNVNSfZOvheupb7Wd5LwgLCSHMJoSFhhAaEkK467ktREiMDueVRVN7vU8+k1YishBYCJCZmXle7zEiOQa7w0mICCFC12PI559HhNqIDAuhX5iNSNfPqecxkaEMiAonPiqM+Kgw+keE+sT/0Eop3xEbGcbEzAFMzBzwuXVtdgfHm9qobengRFsnJ9rsnz42tzuwO5zYO51dj06DvdNJh8OJw2mIieybSPbkUyqBDLfX6a5lnrQJ82BbAIwxi4HF0DXG70Fdn/Ob+Xnns5lSSnlNZJiNrMRon7462JMjEluAHBEZKiLhwAJgxWltVgB3S5cpQKMx5qiH2yqllOpDPe7xG2M6ReQBYBVdp2Q+b4zZLSKLXOufAVbSdUZPMV2nc37lbNv2Sk+UUkp5RC/gUkqpAHAup3Naf/KpUkqpPqXBr5RSQUaDXymlgowGv1JKBRkNfqWUCjI+eVaPiFQDh89z84FAjRfL8WXB1FcIrv4GU19B++sNWcaYJE8a+mTwXwgRKfD0lCZ/F0x9heDqbzD1FbS/fU2HepRSKsho8CulVJAJxOBfbHUBfSiY+grB1d9g6itof/tUwI3xK6WUOrtA3ONXSil1FgET/CIyW0T2i0ixiDxidT3eICLPi0iViOxyW5YgIu+KSJHrcYDbuh+6+r9fRL5oTdXnR0QyRORDEdkjIrtF5Nuu5YHa30gR2Swiha7+/sS1PCD7C1337xaRbSLyput1IPe1VER2ish2ESlwLfOd/hpj/P6HrimfDwLZQDhQCORaXZcX+jUDmAjsclv238AjruePAL9wPc919TsCGOr687BZ3Ydz6GsKMNH1PAY44OpToPZXgP6u52HAJmBKoPbX1YfvAi8Bb7peB3JfS4GBpy3zmf4Gyh7/pzeEN8Z0AKdu6u7XjDFrgLrTFs8DXnA9fwG4wW35UmNMuzHmEF33RpjcJ4V6gTHmqDHmE9fzE8Beuu7ZHKj9NcaYZtfLMNePIUD7KyLpwDXAErfFAdnXs/CZ/gZK8J/pZu+BKNl03d0M4BiQ7HoeMH8GIjIEmEDXXnDA9tc19LEdqALeNcYEcn+fBP4P4HRbFqh9ha7/xN8Tka2u+4mDD/XXZ262rs6dMcaISECdliUi/YF/At8xxjS53+g+0PprjHEAeSISDywXkTGnrQ+I/orItUCVMWariFzeXZtA6aubacaYShEZBLwrIvvcV1rd30DZ4/fkhvCB4riIpAC4Hqtcy/3+z0BEwugK/b8bY15zLQ7Y/p5ijGkAPgRmE5j9vQy4XkRK6RqGvUJE/kZg9hUAY0yl67EKWE7X0I3P9DdQgj+Ybuq+ArjH9fwe4A235QtEJEJEhgI5wGYL6jsv0rVr/ydgrzHm126rArW/Sa49fUSkH3AVsI8A7K8x5ofGmHRjzBC6vpsfGGPuIgD7CiAi0SISc+o5cDWwC1/qr9VHv714FH0uXWeCHAQetboeL/XpH8BRwE7XuN9XgUTgfaAIeA9IcGv/qKv/+4E5Vtd/jn2dRte46A5gu+tnbgD3dxywzdXfXcCPXcsDsr9ufbic/z2rJyD7StfZhYWun92n8siX+qtX7iqlVJAJlKEepZRSHtLgV0qpIKPBr5RSQUaDXymlgowGv1JKBRkNfqWUCjIa/EopFWQ0+JVSKsj8f6BHGTz87h8mAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x=np.random.normal(size=nobs)+2*np.random.uniform(size=nobs)\n", "dens = sm.nonparametric.KDEUnivariate(x)\n", "dens.fit()\n", "plt.plot(dens.cdf)\n", "plt.show()\n", "plt.plot(dens.density)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.5758280021977709" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dens.entropy" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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.6.2" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false } }, "nbformat": 4, "nbformat_minor": 2 }