{ "metadata": { "name": "bootstrap" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "import pandas as pd\n", "import numpy as np\n", "import pandas.rpy.common as com\n", "from scipy import stats\n", "from statsmodels.formula.api import ols" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 39 }, { "cell_type": "code", "collapsed": false, "input": [ "# emulates abline function, only possible to basic line styles and width\n", "def abline(a, intercept, gradient, *args, **kwargs):\n", " #a = gca()\n", " xlim = a.get_xlim()\n", " ylim = a.get_ylim()\n", " \n", " if args:\n", " sty = args[0]\n", " else:\n", " sty = 'r'\n", " \n", " if kwargs:\n", " lw = kwargs['linewidth']\n", " else:\n", " lw = 5\n", "\n", " a.plot(xlim, [intercept + gradient * x for x in xlim], sty, linewidth=lw)\n", " a.set_xlim(xlim)\n", " a.set_ylim(ylim)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 40 }, { "cell_type": "code", "collapsed": false, "input": [ "np.random.seed(333)\n", "x = np.random.normal(size=30)\n", "\n", "bootMean = [np.random.choice(x, size=30, replace=True).mean() for i in range(1000)]\n", "sampledMean = [np.random.normal(size=30).mean() for i in range(1000)]\n", "\n", "bootDensity = stats.kde.gaussian_kde(bootMean)\n", "sampledDensity = stats.kde.gaussian_kde(sampledMean)\n", "\n", "x = np.arange(-0.8, 0.8, .01)\n", "plot(x, bootDensity(x), 'k')\n", "plot(x, sampledDensity(x), 'r');" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD9CAYAAABHnDf0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8jXf/x/FXZCARe8eOkNWMlqZWJPYIUsEttNEKVS29\n9We0dye97yqlVW4dVFG7RpFWQmucEoRbxSgtMUIEsWs1sr6/Py5S2YfknCs5+Twfjzwq57rOOW+a\nfHLlc32HlVJKIYQQwqKU0TuAEEKIoifFXQghLJAUdyGEsEBS3IUQwgJJcRdCCAskxV0IISxQvsU9\nISGBwMBAPDw88PT0ZNasWTnOMRgMVKpUCV9fX3x9ffnPf/5jsrBCCCGMY5PfQVtbW2bMmIGPjw+3\nb9/mqaeeonPnzri5uWU5r3379kRERJg0qBBCCOPle+Veu3ZtfHx8AKhQoQJubm6cP38+x3kyD0oI\nIYqXfK/cHxYfH09sbCx+fn5ZHreysmLXrl14e3vj5OTE9OnTcXd3z3GOEEKIR/fYF8/KCLdu3VJP\nPfWUWrt2bY5jN2/eVHfu3FFKKRUZGalcXFxynGPk2+ju/fff1zuCUUpCzpKQUSnJWdQkZ9EqTO0s\ncLRMamoqISEhPPfccwQHB+c47ujoiL29PQDdu3cnNTWVa9euPd5PGiGEEEUi3+KulCI8PBx3d3fG\njBmT6zlJSUmZvzbs3bsXpRRVq1Yt+qRCCCGMlm/PfefOnSxZsgQvLy98fX0BmDx5MmfPngVgxIgR\nrF69mi+//BIbGxvs7e1ZsWKF6VObSEBAgN4RjFIScpaEjCA5i5rkLD6slDL9UBcrKysZUSOEEI+o\nMLVTZqgKIYQFkuIuhBAWSIq7EEJYICnuQghhgaS4CyGEBZLiLoQQFkiKuxBCWCAp7kIIYYGkuAsh\nhAWS4i6EEBZIirsQQlggKe5CCGGBpLgLIYQFkuIuhBAWSIq7EEJYICnuQghhgaS4i1IjJSWF06dP\nc+nSJb2jCGFyUtyFxdu3bx/t27enQoUKtG/fnubNm+Pq6sq4ceM4d+6c3vGEMAkp7sJiKaUYN24c\nvXr1IiwsjLt373L27FmuXLnCihUrUErh5eXFyJEjuXz5st5xhShSUtyFxXr//fcxGAz8/vvvhIeH\nY2Oj7QdvbW2Nj48Pn3zyCcePH6ds2bK4u7sza9YsUlNTdU4tRNGQDbKFRfrqq6/49NNPiY6OpmbN\nmgWef+TIEcaMGcOFCxf47LPP6NSpkxlSCpG/wtROKe7C4pw6dYqnn36aPXv24OzsbPTzlFKsX7+e\ncePG0aBBAyZNmkTbtm2xsrIyYVoh8ibFXYj7lFIEBQXh7+/PG2+88VivkZqaypIlS/jwww9xdHRk\nxIgR9O/fn2rVqhVxWiHyJ8VdiPvWrVvHW2+9xYEDB7CzsyvUa2VkZLB582a+/vprfvrpJ/z8/OjX\nrx/BwcFGtXqEKCwp7kIAaWlpNG/enLlz59KxY8cife07d+4QFRXFmjVriIqKwt/fn9GjR9OpUydp\n2wiTKUztlNEywmJ8//331KlTJ//Cvm8fvPwyNGsGjo7g7Q3h4XDoUL6v7eDgQL9+/Vi+fDmJiYkE\nBQXx+uuvExgYSFxcXBH/TYQoPCnuwiIopZg+fTrjx4/P/YTUVJg4EXr2BGdnWLUKEhJgwQJo2hS6\ndoU+feD8+QLfy8HBgZdeeomDBw8SHBxMq1atmD9/ftH+hYQoJGnLCIuwfft2hg0bxh9//EGZMtmu\nWdLSIDgYUlJg4UKoWzfnCyQnw5Qp8NVX8OWX8OyzRr/38ePH6dy5MxMmTODVV18t3F9EiIcUpnba\nFHEWIXQxffp0xo4dm7OwKwWvvALp6bBhA9ja5v4C5cppV/Zdu8LAgfDHH/Dmm2BEP71Zs2YYDAY6\ndOiAtbU1L7/8cuH/QkIUkly5ixIvISEBb29vEhMTKV++fNaD06bB8uXwyy9aj90Y589Djx7QujXM\nng3Zf2Dk4eTJkzzzzDP89NNP+Pr6PuLfQoic5IaqKNWWLl1K//79cxb2336Djz+G9euNL+ygtW22\nb4ejR+GllyAjw6inOTs7M3PmTEJDQ7lz584j/A2EKHpy5S5KNKUUnp6ezJ07lzZt2vx9ID1du/IO\nD9cK9OO4fVu7gndz03rxRg55HDJkCBUqVODzzz9/vPcV4j65chel1v79+0lOTqZ169ZZD8yaBfb2\nMGzY4794hQpan/7gQXj3XaOfNmPGDFauXMkff/zx+O8tRCFJcRcl2qJFiwgLC8s6kejKFfjwQ5gz\nx+h+eZ4cHeGHH+C772DuXKOeUrVqVSZMmMC//vWvwr23EIUgbRlRYqWlpVG3bl12796ddYGwsWPh\nr7/giy+K7s3i4qBdO1i8GDp3LvD05ORkmjdvzrJly7K2i4R4BCZryyQkJBAYGIiHhweenp7MmjUr\n1/Nee+01XFxc8Pb2JjY29rGCCPGoduzYQcOGDbMW9rNntbHs771XtG/m4gLLlsHzz2vvUYBy5crx\nwQcf8M477xRtDiGMlG9xt7W1ZcaMGRw5coSYmBg+//xzfv/99yznREZGcuLECeLi4pg7dy4jR440\naWAhHli3bh3BwcFZH5w4EUaOhNq1i/4NO3TQfisICdEmPRVg0KBBnDhxggMHDhR9FiEKkG9xr127\nNj4+PgBUqFABNzc3zmebnh0REcGQIUMA8PPz48aNGyQlJZkorhAapVTO4n76tDbscdw4073xuHFQ\nrx4YcUVua2vLq6++ysyZM02XR4g8GD1DNT4+ntjYWPz8/LI8npiYSP369TM/r1evHufOnaNWrVpZ\nzps4cWLmnwMCAggICHi8xEIABw8exM7ODnd3978fnD5dG/ZYubLp3tjKSrux6u0NvXpB+/b5nj58\n+HCaNm3KlClTcnxPCJGdwWDAYDAUyWsZVdxv375Nv379mDlzJhUqVMhxPHvDP7clUB8u7kIU1oOr\n9syvtaQkbSZqtrahSdSooY3EefFFbZhkPhOkqlWrxoABA5gzZw7vFfV9AGFxsl/4Tpo06bFfq8Bx\nYqmpqYSEhPDcc8/l7G8CTk5OJCQkZH5+7tw5nJycHjuQEMZYt24dffr0+fuBmTMhNBTMdXXcqxcE\nBBjVnhk9ejRfffUVaWlpps8lxH35FnelFOHh4bi7uzNmzJhcz+nduzeLFi0CICYmhsqVK8uvn8Kk\nzpw5w/nz52nVqpX2wK1b2pW0KXvtufn4Y+23hcOH8z3N09OTevXqsWXLFjMFE6KAtszOnTtZsmQJ\nXl5emQshTZ48mbP3h4KNGDGCHj16EBkZSdOmTXFwcGDBggWmTy1KtU2bNtGlSxesra21BxYu1Eay\nNG5s3iDVq8OkSTBqFBgM+S5PEBYWxqJFi+jatav58olSTSYxiRInJCSE4OBgnn/+eW1Rr2bN4Ntv\nQY/JQunp0KIFvPGGtlRwHq5cuULTpk05e/YsFStWNGNAUZLJ2jKi1EhLS2Pr1q10fjBLdMMGqFpV\nWyRMD9bW8Mkn2qSpfHrq1atXJzAwkDVr1pgxnCjNpLiLEmXv3r00bNiQ2g8mKc2YAWPGGL1io0l0\n6ABOTtrSBPkICwvj22+/NVMoUdpJcRclyqZNm/7uW//2m7ZjUr9++oYC+Pe/4YMPtK388tCzZ08O\nHTqUYyKgEKYgxV2UKA9upgLawmAvvQR2dvqGAmjbVuv95zOgwM7OjqCgINauXWvGYKK0khuqosS4\ndu0ajRo14vLly5S9dw8aNdKGIRaXeRXR0drEpmPH8lxqOCIighkzZrBt2zYzhxMlkdxQFaWCwWCg\nTZs2lC1bFpYsgY4di09hB220TuXK8OOPeZ7SpUsXYmNjuXTpkhmDidJIirsoMbZt20ZgYCAoBZ9/\nDq+8onekrKys4P/+T7vJm4dy5crRvXt31q1bZ8ZgojSS4i5KDIPBoK27sX27VuCL4+Jz/frBiROQ\nz74GISEhrF692oyhRGkkPXdRIly+fBkXFxeuXLmCzaBB4O+vzQwtjqZOhaNHtYlVubhz5w5169Yl\nPj6eKlWqmDmcKEmk5y4snsFgoG3btthcugQ//6ztiFRcDR2qrSv/55+5HnZwcMDf35+NGzeaOZgo\nTaS4ixLBYDBo/fZ587Rp/pUq6R0pbzVqaPusLl+e5ym9evXihx9+MGMoUdpIcRclwrZt2whs21Zb\n/bEkbOUYHq79IMpDUFAQGzduJDU11YyhRGkixV0UexcvXuTixYt4nzunrfzo5aV3pIJ17gyXLmmb\neeSibt26NGnShJ07d5o5mCgtpLiLYm/79u20bdsW63nz4OWX9Y5jHGtrbULTN9/keUpQUBA/5jMm\nXojCkOIuir3o6GiCPD3hf/+DkBC94xhvyBBYsSLP1SKl7y5MSYq7KPaio6PpcfEiPPcclC+vdxzj\nNWmifeSxA9OTTz7J7du3iYuLM3MwURpIcRfF2s2bNzl57BhOmzbB8OF6x3l0AwdqV++5sLKyokuX\nLmzatMnMoURpIMVdFGu7d+9mVOPGWDVuDB4eesd5dAMGaGPek5NzPdy1a1cp7sIkpLiLYi06Opqw\ne/e0pX1Lorp1wdsb8piw1LlzZ3755RdS8lkHXojHIcVdFGvHN2+mcVIS9O+vd5THFxqa54SmatWq\n4erqKkMiRZGT4i6KrZSUFHx+/RU1aFDJupGaXd++2pX7X3/lerhbt27SmhFFToq7KLb2793LUKDs\nq6/qHaVwqleHJ5+EzZtzPSx9d2EKUtxFsXV+/nzuVqsGTzyhd5TCCw6GPLbX8/PzIz4+nqSkJDOH\nEpZMirsotupGRXG5Tx+9YxSN4GD44YdcJzTZ2NjQoUMHfvrpJx2CCUslxV0USyopCbekJOq9/rre\nUYpGw4ZQvz7kceNUWjOiqElxF8XSpc8+Y4uDA3WbN9c7StF59lnIY3u9rl278vPPP5ORkWHmUMJS\nSXEXxY9S2CxezLFWrfROUrSCg7XinsvOOg0bNqRKlSocOHBAh2DCEklxF8VPbCwZf/5Jtb599U5S\ntDw9ISMDfv8918PSmhFFSYq7KH4WLGCprS1t/f31TlK0rKyge/c8Z6tKcRdFSTbIFsXLvXuk16lD\nC6XYf+0aVlZWeicqWuvXw+zZ2j6w2dy5c4fatWtz/vx5HB0ddQgnihvZIFtYjogIrjg50cDf3/IK\nO0CHDrBnD9y+neOQg4MDfn5+bNu2TYdgwtJIcRfFy4IFbKpThzZt2uidxDQcHaFlS8ijgHfq1Ikt\neaz/LsSjkOIuio/z5yEmhi+TkmjXrp3eaUyne3eIisr1UMeOHaW4iyIhPXdRfEydSsrRo1Rds4Zr\n165hZ2endyLTOHIEgoLg1CntJutD0tPTqVGjBkePHqV27do6BRTFhfTchWVYupT9np60bNnScgs7\ngLu7tgzBsWM5DllbW9O+fXu2bt2qQzBhSaS4i+LhyBG4fp0fr1+37JYM/D0kUlozwoTyLe5Dhw6l\nVq1aPJHHqnwGg4FKlSrh6+uLr68v//nPf0wSUpQCy5fDP/7B9uho2rZtq3ca0zOiuEsrUxRGvj33\nHTt2UKFCBcLCwjh8+HCO4waDgU8//ZSIiIj830R67iI/SoGLCymLF1O1c2cuXLhg+eO8b94EJye4\neBEcHLIcUkrh5OTEjh07cHZ21imgKA5M1nNv164dVapUyfcFpGiLQtu3D8qUYW9aGq6urpZf2AEq\nVoQWLXIdEmllZSWtGVFoNoV5spWVFbt27cLb2xsnJyemT5+Ou7t7rudOnDgx888BAQEEBAQU5q2F\nJVm+HEJDid65s3S0ZB540JoJCspxqGPHjkRFRfFSSd0YXDwWg8GAwWAoktcqcChkfHw8vXr1yrUt\nc+vWLaytrbG3tycqKop//vOfHD9+POebSFtG5CU9HRo0gC1b6Dl2LEOHDiUkJETvVOZx+DD06QMn\nT+YYEpmQkMCTTz5JUlISZcrIuIfSSrehkI6Ojtjb2wPQvXt3UlNTuXbtWmFeUpQ2O3ZAzZqku7iw\na9eu0nXl7ukJqakQF5fjUP369alSpUquF1VCGKNQxT0pKSnzp8revXtRSlG1atUiCSZKifstmd9+\n+42aNWtSq1YtvROZj5UVdOkCeWyvJ313URj5FvfQ0FBat27NsWPHqF+/PvPnz2fOnDnMmTMHgNWr\nV/PEE0/g4+PDmDFjWLFihVlCCwuRkgJr1sA//kF0aRkCmZ0Ud2EisvyA0M+GDTB5MuzcycCBA+nW\nrRsvvPCC3qnM68oVcHaGy5ch26zcq1ev0qRJE65cuYKtra1OAYWeZPkBUTLdb8kopdixY4flz0zN\nTfXq0KwZxMTkOFStWjWcnZ3Zu3evDsFESSfFXejj7l348Ufo35/Tp0+TkZFBkyZN9E6ljy5dII8d\nmKQ1Ix6XFHehjw0b4OmnoVYtoqOjadeunWVuzmGMrl3z7Lt36NBBNu8Qj0WKu9DH/ZYMUHpbMg88\n8wwcP67137Np27Yt+/bt46+//tIhmCjJpLgL8/vzT9iyBZ59FoCtW7eW7hnLdnbQvr32b5KNo6Mj\nXl5e7N69W4dgoiST4i7Mb+1abS/RypWJj4/n9u3beHp66p1KX/kMiQwMDJT13cUjk+IuzO+hlszW\nrVvp0KFD6e23P/CguOcy7K1Dhw5S3MUjk+IuzOvSJdizJ3OxrC1bttChQwedQxUDLi5gbQ2//57j\nUKtWrTh8+DC3bt3SIZgoqaS4C/NatQp69gR7e5RSbN26lY4dO+qdSn/5LEVQvnx5WrRoQXR0tA7B\nREklxV2Y10MtmaNHj1KuXLnSO749uwKGREprRjwKKe7CfM6ehT/+0K5QQa7as+vQAaKj4d69XA5J\ncRePRoq7MJ/vvoO+fTPXUJF+ezZVqoCHB+zcmeNQy5YtiYuL4/r16zoEEyWRFHdhPg+1ZO7du8e2\nbdvo3LmzzqGKmTz67nZ2drRu3ZpffvlFh1CiJJLiLszj2DFISgJ/fwCio6Nxc3OjRo0aOgcrZvIZ\n7y6tGfEopLgL81i+HAYM0Ib7ARs2bKBHjx46hyqGnn4aTp3SfhBmI5OZxKOQ4i5MT6ksLRmAyMhI\nevbsqWOoYsrWFgIDYfPmHId8fX1JTEwkKZfCL0R2UtyF6cXGQloatGwJwMmTJ7lx4wa+vr46Byum\n8hgSaWNjg7+/PwaDwfyZRIkjxV2Y3vLlMHCgNlEH7aq9R48elCkjX365ymcpAmnNCGPJd5cwrYwM\nWLEiR0tG+u35aNIEHBzgt99yHJL13YWxpLgL09q1CypXhvurPt64cYNdu3bR5f5EJpGHrl1h48Yc\nD3t6enL9+nUSEhJ0CCVKEinuwrSy3UiNiIggMDCQihUr6hiqBOjRQ9utKpsyZcoQEBAgV++iQFLc\nhemkpcHq1Vq//b6VK1cyYMAAHUOVEIGBsH8/5DIjVca7C2NIcRems2ULNG6s9ZCB69evs2PHDnr1\n6qVzsBLA3h7atct11MyDvrvK5YarEA9IcRemk60ls379ejp27Iijo6OOoUqQoKBcWzPNmjUjNTWV\nU6dO6RBKlBRS3IVpJCfD+vXarNT7Vq5cSf/+/XUMVcL07KndVE1Pz/KwlZWVtGZEgaS4C9OIjARf\nX6hTB4CkpCR2795N0P0dmIQRGjSAWrXgf//LcUiGRIqCSHEXppGtJbNw4UL69u0rLZlHFRQEERE5\nHn4wmUn67iIvUtxF0bt5U7sRGBICgFKKefPmMXz4cJ2DlUB9+sC6dTkebty4MeXLl+f3XPZcFQKk\nuAtTiIjQlvatWhUAg8FAuXLl8PPz0zlYCfT003DjhrZkcjbSdxf5keIuil62lszcuXN56aWXsLq/\ntox4BGXKQHAwrF2b41BgYKD03UWerJQZmnZWVlbSGywtrl4FZ2c4dw4qVCApKQlXV1dOnTpFlSpV\n9E5XMv38M7zzDuzZk+XhxMREvLy8uHz5sizCZqEKUzvlK0IUrVWroFs3qFABgFmzZhEaGiqFvTAC\nAiAuDhITszzs5OREjRo1OHjwoD65RLEmxV0UraVLYfBgAG7dusWcOXMYN26czqFKOFtbbcx7LjdW\npTUj8iLFXRSd+Hj44w9tRUO0XnunTp1ocn/5AVEI/fppvxVlIzdVRV6k5y6KzuTJkJAAX37JvXv3\ncHZ25ocffpAdl4rCvXvahLDDh8HJKfPhy5cv07RpU65evYqNjY2OAYUpSM9d6E+pLC2ZuXPn4uXl\nJYW9qJQtq415X7kyy8M1atSgYcOG/PrrrzoFE8VVvsV96NCh1KpViyeeeCLPc1577TVcXFzw9vYm\nNja2yAOKEuLgQbhzB1q35tatW3z44YdMnjxZ71SWJTRU29UqG2nNiNzkW9xffPFFNuayG8wDkZGR\nnDhxgri4OObOncvIkSOLPKAoIZYuhUGDoEwZPvvsMzp27IiPj4/eqSxLhw7afY1sq0FKcRe5ybe4\nt2vXLt8hbBEREQwZMgQAPz8/bty4QVJSUtEmFMVfero2cWnwYK5cucLMmTP54IMP9E5leWxstBur\n2a7e/f39iYmJ4d69ezoFE8VRoe7AJCYmUr9+/czP69Wrx7lz56hVq1aOcydOnJj554CAAAICAgrz\n1qI42b4dqlcHDw/ee+UVBg0ahLOzs96pLNPgwRAeDv/6F9yf8Vu5cmVcXV3Zs2cP/v7+OgcUhWEw\nGDAYDEXyWoW+vZ79Tm5eU8wfLu7Cwty/kXr48GFWr17NH3/8oXciy9WqFWRkQEyM9uf7HrRmpLiX\nbNkvfCdNmvTYr1Wo0TJOTk5ZdmE/d+4cTg8N0xKlQHIyfP89auBAXn/9dd577z2q3l8wTJiAlRUM\nHQrz52d5WPruIrtCFffevXuzaNEiAGJiYqhcuXKuLRlhwSIjwdubiP37uXDhAi+//LLeiSxfWBis\nWaONTrqvbdu27N+/n9u3b+sYTBQn+bZlQkND+eWXX7hy5Qr169dn0qRJpKamAjBixAh69OhBZGQk\nTZs2xcHBgQULFpgltChGFiwgbfBgxo8fz3//+1+ZSGMOdepAmzZagQ8LA8DBwYGWLVtiMBhktysB\nyAxVURjnz4OHB1+89RY/bN1KVFSU3olKj7VrYcYM7Wb2fR999BEXLlxg1qxZOgYTRakwtVOKu3h8\nU6aQfPQoDTZuZNu2bXh4eOidqPRITYVGjSAqCry8ANi/fz+DBg2SG9oWRJYfEOanFMyfz5yUFPr2\n7SuF3dxsbeHll+HzzzMf8vHx4dq1a5w5c0bHYKK4kCt38Xiio7n3wgs4Xb/OkaNH5Ua6Hi5eBDc3\nbcbq/cmGgwcPJiAgQPartRBy5S7M75tvWF6+POPGj5fCrpfataFHD1i4MPOhLl26sGnTJv0yiWJD\nrtzFo7t1i7S6dXmmcmWi4+IoV66c3olKr9274bnn4PhxsLbmwoULuLu7c/nyZRm5ZAHkyl2YlVqx\ngp12doz56CMp7Hpr1Uq7gv/+ewDq1KlDw4YNiYmJ0TmY0JsUd/HI/pwxg+X29oSGhuodRQBMmABT\np2o3uYGePXuyYcMGnUMJvUlxF49EHT1K6okTBE6ZgrW1td5xBECvXtps1fsLTj2YXChKN+m5i0dy\nKiSErb/8wotJSVLci5P587VdmjZuJD09nVq1ahEbG5tl1VZR8kjPXZiFunOHqhER1Hr3XSnsxc3g\nwXD0KOzZg7W1NV27dpWr91JOirswWtwHH7DPzo4eo0bpHUVkV7YsvP02vPceoPXdpbiXbtKWEcZR\nilOVKvH7iy/Sc+ZMvdOI3KSkQPPmsHgxV93caNy4MZcuXZIRTSWYtGWEyZ369lvS796l40cf6R1F\n5MXOTrtyf/ddqlWrhre3t6zxXopJcRdGufr++5zu3p1y9vZ6RxH5ef55SEyErVsJDg5m3bp1eicS\nOpG2jCjQlZgYaN2aMvHxVG3QQO84oiBLl8IXX3Dy229p07YtiYmJcgO8hJK2jDCpM6+9Roy3txT2\nkmLgQLhxA+cTJ6hZs6bMVi2lpLiLfKUmJNBk3z4af/qp3lGEsaytYeJEePddgvv0kdZMKSXFXeTr\nxD//ybbatfEIDNQ7ingUISGQlsYLFSuydu1aaYuWQlLcRd6uXaPODz9Q7q239E4iHlWZMjBtGo3n\nzIGUFH777Te9Ewkzk+Iu8nRpwgQ22NnRecQIvaOIx9GpE1bNmzPd2ZlVq1bpnUaYmRR3kbtLl7Bf\nsoSrI0dia2urdxrxuKZNo+fBg0QtXSqtmVJGirvI1V+TJrFcKf4xfrzeUURhuLtjM2AAI65c4cCB\nA3qnEWYkxV3kFB+PWriQQ0FBsoWeBbCaNInQ1FR+/uorvaMIM5LiLnJQ48fzVdmyPDdhgt5RRFGo\nVYs/w8PxWLJEWjOliBR3kdWOHSRv387KevV4+umn9U4jikidKVPwTUnh6Ny5ekcRZiLFXfwtPR3G\njGGuszPPjxiBlZWV3olEEbFycGD3s89S8Z13IC1N7zjCDKS4i7/Nnk1quXK8f+QIgwYN0juNKGIt\nPv6Y+Bs3SJ09W+8owgxk4TChOXMGnnqKBcOHsy0xkUWLFumdSJjAUD8/vvzjD8rGxUHNmnrHEQWQ\nhcNE4SgFr76KGjOGaevXM3z4cL0TCRMJHDWKyGrV4M039Y4iTEyKu4CFCyEhgZh27cjIyKBt27Z6\nJxIm0rdvX/557RrpUVGwe7fecYQJSXEv7U6dggkTYMkS5i5cyLBhw+RGqgVzcHCgc0gIGwIC4NVX\ntZvowiJJz700S0uD9u0hJIQ/w8Np1KgRx44do6b0Yi3avn37GNC/PycbNsSqXz+QDc+LLem5i8fz\n7rtgbw9jxrB8+XI6deokhb0UaNGiBdWqV2d7aKi27vvZs3pHEiYgxb20ioiAJUtg2TIoU4Z58+bJ\njdRS5JVXXmHaDz/AP/8JI0dqN9WFRZHiXhqdOgXDhsHKlVCjBrGxsVy5coVOnTrpnUyYyT/+8Q9i\nYmI4PWCAduW+bJnekUQRK7C4b9y4EVdXV1xcXJg6dWqO4waDgUqVKuHr64uvry//+c9/TBJUFJHk\nZOjXD961gkaCAAAWjUlEQVR5B1q1AmDevHmEh4dTpoz8rC8t7O3tGTJkCF/MmwcLFsDrr0Niot6x\nRBHK94Zqeno6zZs3Z/PmzTg5OdGyZUuWL1+Om5tb5jkGg4FPP/2UiIiIvN9EbqgWH8OHw82bsGIF\nWFlx9+5d6tevz8GDB6lXr57e6YQZnTlzhieffJLTp09TccYM2LULNm4EGS1VbJjshurevXtp2rQp\njRo1wtbWloEDB7J+/foc50nhLiHmz4foaJg3L/MbeNWqVbRq1UoKeynUsGFDunbtytdffw1vvQXX\nr8MXX+gdSxQRm/wOJiYmUr9+/czP69Wrx549e7KcY2Vlxa5du/D29sbJyYnp06fj7u6e47UmTpyY\n+eeAgAACAgIKl1w8mn374I03YPt2cHTMfHjevHmMHTtWx2BCT+PGjaNPnz689tpr2C5ZAm3aQNu2\n4O2td7RSyWAwYDAYiubFVD5Wr16thg0blvn54sWL1ahRo7Kcc/PmTXXnzh2llFKRkZHKxcUlx+sU\n8DbC1C5fVqphQ6VWr87y8NGjR1Xt2rVVSkqKPrlEsdChQwe1aNEi7ZMlS5Rq1kypmzf1DSWUUoWr\nnfm2ZZycnEhISMj8PCEhIcev746Ojtjb2wPQvXt3UlNTuXbtWtH85BGFl5YGAwdqHyEhWQ7NnTuX\noUOHyh6ppdwbb7zBRx99RHp6OgweDO3awYgRMjyyhMu3uLdo0YK4uDji4+NJSUnhu+++o3fv3lnO\nSUpKyuy57927F6UUVatWNV1i8WjeeUf7b7ZRTH/99ReLFy+Wse2Czp07U7lyZVatWqU9MGsW/PEH\nzJihbzBRKPn23G1sbJg9ezZdu3YlPT2d8PBw3NzcmDNnDgAjRoxg9erVfPnll9jY2GBvb8+KFSvM\nElwYYc0abVTM//4HNln/V69atYqWLVvSqFEjfbKJYsPKyoqJEycyZswY+vfvj7W9PaxdC888A088\nAZ076x1RPAZZW8ZS/f47+PtDVBS0aJHjcJs2bRg/fjzBwcE6hBPFjVKKNm3a8NprrzFw4EDtwV9+\ngQEDYOtW8PDQN2ApVZjaKcXdEt28CX5+MH48DB2a4/Bvv/1Gt27diI+Px8Ym31/eRCny008/MXr0\naI4cOfL318WSJdowyZ074aGRc8I8ZOEw8TeltIlK/v65FnaAOXPmEB4eLoVdZNG5c2fq1avH/Pnz\n/37wuefgtdegWzdtHLwoMeTK3dIsXgxTp2p99vLlcxy+c+cODRo04MCBA1nmMAgB2nLAffr04fjx\n4zg4OPx9YOxY2LsXfvop168rYRpy5S408fHwf/8HS5fm+Q343Xff0bp1aynsIlctWrSgXbt2fPbZ\nZ1kPTJumtWUGDdKG14piT67cLUV6OgQEQJ8+MG5cnqf5+fnx7rvvEhQUZL5sokQ5efIkfn5+HDp0\niLp16/59ICUFevWCGjXg22/B2lq/kKWEXLkL+PhjsLXVrtzzcODAAS5cuED37t3NGEyUNM7Ozrz0\n0ktMmDAh6wE7O22I5Pnz8NJLkJGhT0BhFLlytwS//grdu2v/zafdMnToUJydnXn77bfNGE6URHfu\n3MHNzY1ly5bl3DD99m3tBquXF3z+uawiaUIyFLI0u3sXnnwS3n8fQkPzPC0pKQlXV1fi4uKoXr26\nGQOKkuq7777jo48+4tdff8U6ewvm5k3o0kWb6DRjhhR4E5G2TGk2fjw89VS+hR3gyy+/ZMCAAVLY\nhdEGDBhAlSpVMmekZ1Gxorb2+44d8Oabsg5NMSRX7iVZZKS2/+XBg1C5cp6nJScn06hRI7Zt25Zl\noxUhCnL48GE6duzI0aNHc78wuHoVOnTQbuR/8IH5A1o4uXIvjS5d0vZB/fbbfAs7wJIlS/D19ZXC\nLh7ZE088QWhoaN73aapVg59/1tYx+vBD84YT+ZIr95JIKQgOBldXbcJSPtLS0nB1dWX+/Pn4+/ub\nKaCwJDdu3MDd3Z3Vq1fTunXr3E+6cAHat9dG0eQzFFc8msLUTpl/XhLNm6ftWL9yZYGnLl++HCcn\nJyns4rFVrlyZWbNmER4eTmxsLOXKlct5Up062gJj7dtD2bIwerT5g4os5Mq9pImLg1attO3yctnO\n8GHp6em4u7vzxRdf0LFjRzMFFJZIKUVISAgeHh78+9//zvvE+HhtMt1bb2lX8aJQ5Mq9tEhN1RZy\nev/9Ags7wMqVK6levTodOnQwQzhhyaysrPj888/x9vamX79+eOe1x2qjRrBli1bg7ezghRfMmFI8\nTK7cS5L334c9e7RRMmXyvxeenJyMu7s78+bNk+Iuisz8+fP54osviImJyX9V0WPHtFE006Zp69GI\nxyKjZUqD3bvhq69g/vwCCzvAzJkz8fLyksIuitSLL75I5cqVmVHQFnzNm2srSI4dC6tXmyecyEKu\n3EuCW7fAx0e7Curbt8DTk5KS8PDwYPfu3bi4uJghoChNTp06xdNPP010dDSurq75n3zwoDaT9euv\nIdv+y6JgsvyApQsP1/77zTdGnf7iiy9SpUoVPv30UxOGEqXZnDlzmDNnDrt376Zs2bL5n7xvH/To\nAYsWaWvSCKNJcbdkK1dqIw9iY8HRscDTf/zxR0aPHs2hQ4dwNOJ8IR6HUoq+ffvi7OzM9OnTC37C\n7t3aLNbly0FGbhlNirulOnFCG/a4caO2fkwBrl69ipeXF8uWLaN9+/ZmCChKs6tXr+Lt7c38+fPp\n0qVLwU/Yvh369dN68DLvwihS3C1RcjK0bq21ZF59tcDTMzIy6N+/Pw0bNpR2jDCbLVu2EBYWRmxs\nLDVr1jTmCdoid+vXaxcuIl9S3C3Rq6/C5cvw3XdGLac6ceJEfvrpJ7Zu3Zr7DEIhTOTNN9/kyJEj\nREREYGXM0r8bN0JYGGzYAC1bmj5gCSZDIS3NypWwaZM2wsCIb5YVK1awcOFC1q5dK4VdmN0HH3zA\nxYsXCx4e+UC3btrggKAgOHDAtOFKMblyL27i4qBNG+3q5sknCzx9zZo1vPLKK/z88894eXmZIaAQ\nOZ05cwY/Pz+WLl1q/FIXa9bAqFHaqpKenqYNWELJlbuluHZNu5r58EOjCvvSpUsZNWoUmzZtksIu\ndNWwYUOWL1/O4MGDOX36tHFPCgmBTz/VxsHHxpo2YCkkV+7FRUoKdO2qjYopYGhZeno67733HosX\nLyYqKgoPDw8zhRQif7Nnz+aLL75g586dVKlSxbgnrVmjbTqzbBl06mTagCWM3FAt6dLT4fnntf1Q\n16yB7PtVPuTixYu88MILJCcns3LlSuNGKAhhRmPHjmXfvn1s2rTJ+HtA27dD//7w2WcFbhlZmkhb\npiTLyIARIyApSZvgkUdhV0qxcuVKvL29adGiBZs3b5bCLoqladOmUbt2bUJDQ0lNTTXuSf7+sHkz\nTJigbbgtCk2u3PWUnq79OnrkiDY6pkKFXE87deoUo0eP5tSpUyxcuBA/Pz8zBxXi0aSkpNCvXz9s\nbW1ZsWIFtra2xj3x7Fno3h3atoVZs7SNP0oxuXIvif76S5utd/o0REXlWtiTk5OZNGkSTz/9NP7+\n/hw8eFAKuygR7OzsWLVqFffu3aNfv37cvXvXuCc2aKAtVXDpkrYmfEKCSXNaMinuejh9WtuOzN5e\nm8hRsWKWw0op1q9fj6enJ4cOHWL//v288cYb2NnZ6RRYiEdXtmxZvv/+eypWrEiHDh24fPmycU+s\nWFG799SnjzbAYNUq0wa1VMoMzPQ2xV9GhlIrVypVo4ZSn36qfZ7N3r17lb+/v/L09FQbN27UIaQQ\nRSsjI0O9/fbbqlGjRmrPnj2P9uS9e5VycVFqwAClEhNNE7AYK0ztlOJuLsePK9W9u1Jubkrt3p3j\n8IkTJ1RoaKiqU6eO+vrrr1VaWpoOIYUwnTVr1qgaNWqo6dOnP9rX9507Sr31llLVqin18cdK3b1r\nupDFTGFqp7RlTO3QIW2Y4zPPQGCgNt36mWcyDx8+fJjBgwfj5+dHs2bNOH78OMOGDcM6n+GQQpRE\nffv2JSYmhoiICJ555hlijZ24ZG+vTeyLjtb68S4uMHMm3Lxp2sAlnBT3hxgMhsK/iFJw+LC2a5Kv\nr3bn38MDTp6E8ePBzg6lFDt27KB379506dIFLy8vTp48ycSJE6mQx4iZIs9pYiUhI0jOolZQziZN\nmrBt2zZGjhxJt27deO6554iLizPuxV1d4fvvYe1a2LlT24x75Ej45Rdt5FkR5rQEBRb3jRs34urq\niouLC1OnTs31nNdeew0XFxe8vb2N/2lcDD3y//A//9R2mVm+XNu8ulcvqF0bgoO1Yv7JJ9rQrjff\nhMqVOXbsGJMnT8bV1ZWXXnqJbt26cerUKd544w0qVapkupw6KAkZQXIWNWNylilThqFDh3LixAlc\nXV1p1aoVwcHBbNq0iXRjinTLltriegcPaqNr/vlPqFMHBg7U9hnetavAq/qS8u9ZGPlsX65Ncx81\nahSbN2/GycmJli1b0rt3b9zc3DLPiYyM5MSJE8TFxbFnzx5GjhxJTEyMyYObnFJw/TqcOaMV6LNn\ns/751CltRmnTptqvic2bw4svwuefQ4MGZGRkEB8fT+y6dWzbto3Nmzdz69YtgoODWbBgAa1atTJu\neVQhLJSjoyPvvPMOr7/+OsuWLePtt98mLCyMZ599lq5du9K+fXuqVq2a9wvUrw//+pf2ceaMtlb8\njh3aipNHj0L16uDuDvXqacX/4Y8//9SGI5cvb76/sJnlW9z37t1L06ZNadSoEQADBw5k/fr1WYp7\nREQEQ4YMAcDPz48bN26QlJRErVq1TJf6cSmlbYJx65Y2jjb7R1QU7NnzdxG3sYGGDbWrgwcfLVpA\ngwbcq1OH246OXLt+ncTERM6dO0diXBzxP//MoUOHOHz4MJUqVcLHxwd/f3+WLl2Kr68vZcpIJ0yI\nhzk4ODB8+HCGDx/OyZMn+f7775kzZw5DhgyhTp06+Pj44OLiQr169TI/qlWrRtmyZSlXrhxly5bF\n2smJlIEDuffss6SkpHDv7l1Sjh1DHT2KSkzE6sQJyuzejd2VK5S9do2UpCTS/vtflJUVd8qV46/y\n5fmrfHmS7e25V6ECaRUrklGpEqpqVahWjTLVq2NTsya2tWphV6sW9tWqYe/ggL29PeXKlSueF2r5\n3W1dtWqVGjZsWObnixcvVqNGjcpyTlBQkNq5c2fm5x07dlT79u3LccdXPuRDPuRDPh7943Hle+Vu\n7E8jlW16bPbnZT8uhBDCtPLtETg5OZHw0PTfhIQE6tWrl+85586dw8nJqYhjCiGEeBT5FvcWLVoQ\nFxdHfHw8KSkpfPfdd/Tu3TvLOb1792bRokUAxMTEULly5eLZbxdCiFIk37aMjY0Ns2fPpmvXrqSn\npxMeHo6bmxtz5swBYMSIEfTo0YPIyEiaNm2Kg4MDCxYsMEtwIYQQ+Xjsbn0+rl69qjp16qRcXFxU\n586d1fXr13M9b/Lkycrd3V15enqq0NBQlZycbIo4hc55/fp1FRISolxdXZWbm5vancvyAcUhp1JK\npaWlKR8fHxUUFGTGhBpjcp49e1YFBAQod3d35eHhoWbOnGmWbFFRUap58+aqadOmasqUKbmeM3r0\naNW0aVPl5eWl9u/fb5Zc2RWUc8mSJcrLy0s98cQTqnXr1urgwYM6pDTu31Mpba0ka2trtWbNGjOm\n+5sxObdt26Z8fHyUh4eHat++vXkD3ldQzsuXL6uuXbsqb29v5eHhoRYsWFDga5qkuI8fP15NnTpV\nKaXUlClT1BtvvJHjnNOnT6vGjRtnFvQBAwaohQsXmiJOoXIqpVRYWJj65ptvlFJKpaamqhs3bpgt\no1LG51RKqU8++UQNGjRI9erVy1zxMhmT88KFCyo2NlYppdStW7dUs2bN1NGjR02aKy0tTTk7O6vT\np0+rlJQU5e3tneM9N2zYoLp3766UUiomJkb5+fmZNNPj5ty1a1fm119UVFSxzfngvMDAQNWzZ0+1\nevXqYpnz+vXryt3dXSUkJCiltCJaHHO+//776s0338zMWLVqVZWamprv65pk0PXDY9+HDBnCunXr\ncpxTsWJFbG1tuXv3Lmlpady9e9fsN2KNyfnnn3+yY8cOhg4dCmitqkeZTWqunKDdzI6MjGTYsGG6\njFAyJmft2rXx8fEBoEKFCri5uXH+/HmT5np4voatrW3mfI28sj88X8OcjMnZqlWrzK8/Pz8/zp07\nZ9aMxuYE+O9//0u/fv2oUaOG2TOCcTmXLVtGSEhI5kCR6tWrF8ucderU4eb9Wbc3b96kWrVq2Njk\n21U3zdoyD09iqlWrVq7fJFWrVmXs2LE0aNCAunXrUrlyZTqZeXNcY3KePn2aGjVq8OKLL/Lkk08y\nfPhw4zceMGNOgNdff51p06bpNlHK2JwPxMfHExsba/INSBITE6lfv37m5/Xq1SMxMbHAc8xdOI3J\n+bBvvvmGHj16mCNaFsb+e65fv56RI0cCxg+rLkrG5IyLi+PatWsEBgbSokULFi9ebO6YRuUcPnw4\nR44coW7dunh7ezNz5swCXzf/0p+Pzp07c/HixRyPf/jhh1k+t7KyyvV/7MmTJ/nss8+Ij4+nUqVK\n9O/fn6VLlzJ48ODHjWSSnGlpaezfv5/Zs2fTsmVLxowZw5QpU/jggw+KVc4ff/yRmjVr4uvra9J1\nMwqb84Hbt2/Tr18/Zs6cadRiaYVRVPM1TO1R3m/btm3Mnz+fnTt3mjBR7ozJ+eD75ME2cXr8JmlM\nztTUVPbv38+WLVu4e/curVq14plnnsHFxcUMCTXG5Jw8eTI+Pj4YDAZOnjxJ586dOXjwII6Ojnk+\n57GL+88//5znsVq1anHx4kVq167NhQsXct3Ied++fbRu3Zpq1aoB2nKgu3btKvLiXticD6Y7t2zZ\nEoB+/foxZcqUIs1YFDl37dpFREQEkZGRJCcnc/PmTcLCwjKHqRaXnKB9Q4WEhPDcc88RHBxcpPly\nU1LmaxiTE+DQoUMMHz6cjRs3UqVKFXNGBIzL+euvvzJw4EAArly5QlRUFLa2tjmGUuuds379+lSv\nXp3y5ctTvnz5zO0szVncjcm5a9cu3n77bQCcnZ1p3Lgxx44do0WLFnm/sCluEIwfPz7zju9HH32U\n6421AwcOKA8PD3X37l2VkZGhwsLC1OzZs00Rp1A5lVKqXbt26tixY0op7cbGhAkTzJZRKeNzPmAw\nGHQZLWNMzoyMDPX888+rMWPGmC1XamqqatKkiTp9+rS6d+9egTdUd+/ercuNSmNynjlzRjk7O5t9\nxNbDjMn5sBdeeEGX0TLG5Pz9999Vx44dVVpamrpz547y9PRUR44cKXY5X3/9dTVx4kSllFIXL15U\nTk5O6urVq/m+rsmGQnbs2DHHkLjExETVo0ePzPOmTp2aORQyLCxMpaSkmCJOoXMeOHBAtWjRQnl5\nealnn33W7KNljM35gMFg0GW0jDE5d+zYoaysrJS3t7fy8fFRPj4+KioqyuTZIiMjVbNmzZSzs7Oa\nPHmyUkqpr776Sn311VeZ57z66qvK2dlZeXl5qV9//dXkmR4nZ3h4uKpatWrmv13Lli2LZc6H6VXc\nlTIu57Rp0zLrkLmG5j5qzsuXL6ugoCDl5eWlPD091dKlSwt8TSulZOEXIYSwNLL+rBBCWCAp7kII\nYYGkuAshhAWS4i6EEBZIirsQQlggKe5CCGGB/h8/qYVpZhfRvAAAAABJRU5ErkJggg==\n" } ], "prompt_number": 41 }, { "cell_type": "code", "collapsed": false, "input": [ "# this is a bit of a shame\n", "import imp\n", "\n", "bootstrap = imp.load_source('scikits.bootstrap', '/Users/erriza/scikits-bootstrap/src/package/scikits/bootstrap/bootstrap.py')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 42 }, { "cell_type": "code", "collapsed": false, "input": [ "np.random.seed(333)\n", "x = np.random.normal(size=30)\n", "sampledMean = [np.random.normal(size=30).mean() for i in range(1000)]\n", "\n", "# this computes the confidence interval\n", "bootstrap.ci(x, statfunction=np.mean, n_samples=100)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "/Users/erriza/scikits-bootstrap/src/package/scikits/bootstrap/bootstrap.py:114: InstabilityWarning: Some values used top 10 low/high samples; results may be unstable.\n", " warnings.warn(\"Some values used top 10 low/high samples; results may be unstable.\", InstabilityWarning)\n" ] }, { "output_type": "pyout", "prompt_number": 43, "text": [ "array([-0.27856777, 0.29389323])" ] } ], "prompt_number": 43 }, { "cell_type": "code", "collapsed": false, "input": [ "nuclear = com.load_data('nuclear', package='boot')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 44 }, { "cell_type": "code", "collapsed": false, "input": [ "nuke_lm = ols('np.log(cost) ~ date', nuclear).fit()\n", "\n", "plot(nuclear['date'], np.log(nuclear['cost']), 'ok')\n", "abline(gca(), nuke_lm.params['Intercept'], nuke_lm.params['date'], 'r', linewidth=3)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD9CAYAAACyYrxEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1UlGX+P/A3AoL6VdMkwYZCR0VAAR/IXFOmUsavloVN\nm0VlWEilgJ2203f3fDvLaXfTOn1/8lC2tKs9bGWW6dYyPi3VpKkISmLbgxrSKigkJoqwgIPX748r\ngfEeYGBm7nl6v87hHLive+b+dDd+uLnu6/58/IQQAkRE5LX6uToAIiJyLiZ6IiIvx0RPROTlmOiJ\niLwcEz0RkZdjoici8nI9Jvr6+noYDAZERUUhOjoaxcXFFuN1dXWYN28e4uPjMXHiRLz55pvOipWI\niPrAr6d19EuWLEFiYiKWLl0Ks9mMxsZGDB06tH08OzsbLS0tWLVqFerq6hAZGYna2loEBAQ4PXgi\nIupZt1f058+fx+7du7F06VIAQEBAgEWSB4CwsDBcuHABAHDhwgVce+21TPJERG6k24xcWVmJkJAQ\npKamory8HFOnTkVubi4GDhzYvk9aWhpuu+02jBo1Cg0NDfjggw8U7+Pn5+f4yImIfIAjihd0e0Vv\nNptRVlaGJ598EmVlZRg0aBBWr15tsc8LL7yA+Ph4nDp1CocOHcLy5cvR0NBgNVh3+vr973/v8hgY\nk3PjSkpKsvq51uv1PnOuGJPnxiSE46rTdJvoNRoNNBoNEhISAAAGgwFlZWUW++zduxf33nsvAECr\n1WL06NE4cuSIwwIk6quWlhar25ubm1WOhMi1uk30oaGhCA8Px9GjRwEARUVFiImJsdhnwoQJKCoq\nAgDU1tbiyJEjGDNmjJPCJbJdUFCQ1e3BwcEqR0LkWj3eNc3Pz0dKSgpaW1uh1Wqxfv16FBQUAADS\n09Pxu9/9DqmpqYiLi8Ply5fx0ksvYfjw4U4P3F46nc7VISgwJtvZEldmZiYqKipQUVHRvk2r1SIj\nI8NlMamNMdnGHWNypB6XVzrkIH5+Dp1vIrKV0WhEfn4+mpubERwcjIyMDCxYsMDVYRHZxFG5k4me\niMidCAHs2QPk5MDvo48ckju54J2IyB20tgIffADk5AAHDzr0rZnoiYhcqa4OKCgAXn0VOH3aKYdg\nUTMiIlf45htg2TIgPBz43/+1TPJBQcBjjznsUJyjJyJSy+XLwI4dcnpm507leGgosHw5kJ4OhIQ4\nLHdy6oaIyNmamoC33wZyc4Hvv1eOT5kCPPUU8OtfA/37O/zwTPRERM5SVSXn3gsKgHPnLMf8/IC7\n75YJ/pZb5M9OwkRPRORoJSVyeubDDwGz2XJs8GA5/75iBaBSFQEmeiIiRzCbgS1bZILfu1c5Pno0\nkJUFpKYCQ4aoGhoTPRGRPerrgb/+FcjPB06cUI4nJgIrVwJ33gn4+6sfH5joiYj65tgxIC8PeOMN\noLHRciwwELj/fnkFP2WKa+LrhImeiMhWQgCffw6sWQMYjfLnzkaMAJ54Qn6FhbkmRiuY6ImIetLc\nDGzYIOffDx9Wjk+cKKdnHngAGDBA/fh6wERPRNSV2lrgtdfk108/KccXLJAJ/vbbnbo80l5M9ERE\nVzt0SF69b9ggi411NnAg8MgjQGYmEBnpkvB6i4meiAgA2trkvPuaNYDJpBzXaICMDCAtDRg2TPXw\n7MFET0S+raFBrpzJywM6dSNrd/PNcnpm0SK5msYDMdETkW/68Ue59v2vfwUuXLAc8/cHDAaZ4G++\n2SXhORITPRH5jk7dm7Bli6wm2dk118jSwStWyPLBXqLHevT19fUwGAyIiopCdHQ0iouLFfuYTCZM\nnjwZEydO9Pomu0TkgVpbgXffBW66CZg1C/joI8skP368LD5WVQW8+KJXJXnAhnr0S5YsQWJiIpYu\nXQqz2YzGxkYMHTq0fby+vh4zZ87Ejh07oNFoUFdXhxEjRlgehPXoicgVzp7t6N506pRyfM4cWT1y\n3jygn/v1YVKlOfj58+cxefJkHD9+vMs3WLt2LWpqavD88893fRAmeiJS07ffytrvb78tH3bqLCgI\nePBBOf8+caJr4rORKo1HKisrERISgtTUVJSXl2Pq1KnIzc3FwIED2/c5duwYLl26hFtvvRUNDQ3I\nysrCQw89pHiv7Ozs9u91Oh2neIjIsYTo6N60Y4dy/KruTe7IZDLBZG1pp526vaI/cOAAZsyYgb17\n9yIhIQErV67EkCFDLK7eV6xYgbKyMnz66adoamrCjBkzYDQaMW7cuI6D8IqevJjRaEReXh5aWloQ\nFBSEzMxMLFiwwNVh+Y6mJuBvf5NX8N99pxyfPLmje1NQkPrx2UGVK3qNRgONRoOEhAQAgMFgwOrV\nqy32CQ8Px4gRIzBgwAAMGDAAs2fPRnl5uUWiJ/JWRqMRWVlZqOi0/vrK90z2TlZd3dG96eefLceu\ndG9auVLefHXj8gRq6PbuQ2hoKMLDw3H06FEAQFFREWJiYiz2ueuuu/Dll1+ira0NTU1N2L9/P6Kj\no50XMZEbycvLs0jygEz0+fn5LorIB5SWAikpQEQEsGqVZZIfPFgm9x9+ADZvBmbP9vkkD9iwjj4/\nPx8pKSlobW2FVqvF+vXrUVBQAABIT0/HhAkTMG/ePMTGxqJfv35IS0tjoief0dLSYnV789U3AMk+\nZjPw97/L8gRddW/KzASWLnVI9yavm44TKlDpMNSNwsJCkZSUJBITE0VSUpIoLCx0dUheISkpSQBQ\nfOn1eleH5h3OnRPi5ZeFuPFGIeTtVsuv2bOF2LxZCLPZYYcsLCwUWq3W4v+nVqt1yb8ZR+VOJnof\n4E4fXG/Dc+skx44JkZEhxKBByuQeGCjEQw8JceCAUw7tTr+8HZU7WQLBB3Q3j+zRf466gSvnLz8/\nH83NzQgODkZGRgbPa18IIatGrlkDFBZa7970+OOye9OoUU4Lwxun45jofYA3fnDdyYIFC5jY7dHc\nDLz/vlz/Xl6uHI+JkTdYU1JU6d4U1MUSzODgYKcf21mY6H2AN35wyQv01L1p/nyZ4OfMUXXlTGZm\nJioqKiz+CtZqtcjIyFAtBkdjovcB3vjBJQ9WXi6v3t97z3r3piVLgKwsh3Zv6s0qGm+cjuuxqJlD\nDsInY13OaDR61QeXPMzly3LePScH+Pxz5bhGI0sDp6UBw4c79NDWHmrTarXIzc11+38DqhQ1cxQm\neiIfdfFiR/emH35Qjt90kyxPcM89TuvepNfrsXPnTqvbt2/f7pRjOooqJRBIHV73cAbRv//d0b3p\n/HnLMX9/mdifekqV7k1cjMBE73KslUJeQwj51GpOjiw/YK17U1qanKK54QbVwuJiBBs6TJFzsVYK\nebxLl+SN1enTgVtuATZtskzy48bJ4mMnTwIvvaRqkgfkYgStVmuxzdcWI/CK3sWqq6utbq+qqlI5\nEqJeOnsWeP114JVXrHdvuv12OT3z3//t0u5N3riKprd4M9bFRowYgbNnz1rdfubMGRdERNSD777r\n6N70n/9Yjl3p3pSVBUya5Jr4vAhvxnqJUaNGWU30YWFhLoiGqAtCADt3yvl3aytVRo7s6N503XXq\nx0fdYqJ3sbCwMHz99deK7aOcWMuDyGZNTcA778gEb617U3y8nJ657z6P697kS5joXYxPrfoej1hO\nW10NrF0ruzdd/Rennx9w112yPAEbe3gEJnoX440i3+L2y2kPHJBX7xs3ymYfnf3XfwGPPgpkZABX\nrWIh98absUQqcsunNM1m4OOPZXngPXuU4xERHd2bhg5VPTxfxpuxRB7I1U9pdp42GtavH/40Zgyi\ni4rkk6xXmzVLTs/cdZd8mpU8Vo+LW+vr62EwGBAVFYXo6GgUFxdb3a+0tBQBAQHYvHmzw4Mk8hau\nfEqzfdpo504s+uILvP3554het84yyQcEyOWRBw4Au3YBixYxyXuBHq/os7KyMH/+fGzatAlmsxmN\njY2Kfdra2vDss89i3rx5nKIh6obLbr4LAVN2Nv6vogJ3wsoV3rXXyu5NTz7p1O5N5BrdztGfP38e\nkydPxvHjx7t9k5ycHPTv3x+lpaW44447cM8991gehHP0RO1ULRnd0gJs2NBl96ZvAGwdPx7PHDqk\nSvcm6h1V5ugrKysREhKC1NRUlJeXY+rUqcjNzcXAgQPb96mursbHH3+Mzz77DKWlpfDrYqlVdnZ2\n+/c6nQ46nc7u4Ik8kSqtB3/6SXZuWrvWavemrQByAPwTgH70aDzDJO8WTCYTTCaT49+4u87hpaWl\nIiAgQJSUlAghhMjKyhLPPfecxT4Gg0EUFxcLIYRYsmSJ2LRpk+J9ejgMETlKebkQqalC9O8vhHye\ntf3LHBQk3hkyREQCAr98abVaUVhY6OqoqQuOyp3dXtFrNBpoNBokJCQAAAwGA1avXm2xz8GDB7F4\n8WIAQF1dHbZt24bAwEAsXLjQ8b+ViEjp8mXAaJTTM599phy//npgxQr4L1uGa/btQ0R+PkL5zIZP\n6TbRh4aGIjw8HEePHsX48eNRVFSEmJgYi306z9+npqbizjvvZJInUsPFi8Cbb8oCYzZ2b1Jl2ojc\nTo+rbvLz85GSkoLW1lZotVqsX78eBQUFAID09HSnB0hEV/n3v2Vp4L/8Rdm9qV8/y+5NLE9A4JOx\nRJ5BCGDfvo7uTW1tluNDhwLLlskKkjfe6JoYyeH4ZCx5DI8o4uWuLl2SHZtycoCSEuX4uHGy9vuS\nJbIWDZEVTPTkVG5fxMtd/fxzR/cma13IbrtNTs/Mn+/S7k3kGTh1Q07llkW83Nn338ubq2+9Zb17\nU0qKvIKPjXVNfKQqTt2QR3B1ES+PIATwz3/K6Zlt25TjI0fK0gSPP87uTdQnTPTkVK4s4uX2/vOf\nju5N336rHI+Lk9MzixezexPZhYmenIodtKw4dUqWJvjzn613b1q4UJYHTkzk8khyCCZ6cip20Ork\n4EHZ3KOr7k1Ll8ruTWPHuiY+8lq8GUvkTG1tHd2bvvxSOX7jjbJ706OPsnsTKfBmLJE7O38eWLcO\nyM8HfvxROX7LLR3dmwL4z5Cci58wIkeqqADy8oD162Utms4CAoD77pMJfto018RHPomJnsheQsi2\ne2vWAJ98In/ubPhwuTRy+XJ2byKXYKIn6quWFuD99+XyyEOHlONRUfLq/cEHgU7NeojUxkRP1Fs/\n/SSXRq5dC9TWKsfnzZPr3+fO5fJIcgtM9ES2+vprefX+7rvyar6zAQOAhx+W5QmiolwTH1EXmOiJ\nunP5MrB1q0zwn36qHP+lexPS0oBrr1U/PiIbMNETWXPxoiwslpsLHDumHE9IkNMzBkN79yYid8VE\nT9TZiRMd3Zvq6y3H+vUDFi2SCX7GDM6/k8dgoicCOro3ffSR9e5NaWlyiobdm8gDMdGT77p0SSb2\nnBxg/37l+NixHd2bBg9WPz4iB+mxNU19fT0MBgOioqIQHR2N4uJii/F3330XcXFxiI2NxcyZM3H4\n8GGnBUu+w2g0Qq/XQ6fTQa/Xw2g0Ou7Nf/4ZePFFYMwY4P77lUn+ttvkg09HjsireCZ58nSiBw8/\n/LBYt26dEEKIS5cuifr6eovxvXv3tm/btm2bmD59uuI9bDgMUbvCwkIRGhoqALR/hYaGisLCQvve\n+LvvhHjiCSEGDhRCPr/a8dW/vxCpqUIcOuSY/wgiB3BU7uy2euX58+cxefJkHD9+3KZfGufOncOk\nSZNQVVVlsZ3VK6k3pkyZgq+++srq9oMHD/buzYQAiopkeQJr3Zuuu66je9PIkX2MmMg5VKleWVlZ\niZCQEKSmpqK8vBxTp05Fbm4uBnbxOPe6deswf/58q2PZ2dnt3+t0Ouh0uj4HTd7tR2vVHiE/jzb7\nz3/kg005OcA33yjHY2Pl6pn772f3JnIbJpMJJpPJ4e/b7RX9gQMHMGPGDOzduxcJCQlYuXIlhgwZ\ngueff16x7+eff47ly5djz549GDZsmOVBeEVPvTB8+HCcO3dOsX3YsGH4+eefu3/x6dPAq6923b3p\nzjtl/Rmdjssjye05Knd2ezNWo9FAo9EgISEBAGAwGFBWVqbY7/Dhw0hLS8Mnn3yiSPJEvRUREWF1\n++jRo7t+UVkZ8NBDcvnjn/5kmeQHDZKdm44ckU1Abr2VSZ58SreJPjQ0FOHh4Th69CgAoKioCDEx\nMRb7nDhxAosWLcI777yDsWyBRg7whz/8AaGhoRbbQkNDlX9JtrUBmzcDs2cDU6fKRtuXLnWM33gj\n8PLLQFWVrBE/bpwK0ZNTV0xRn/TYSrC8vByPPfYYWltbodVqsX79emzcuBEAkJ6ejsceewxbtmzB\nDTfcAAAIDAxESUmJ5UE4dUO9ZDQau+4ze+GC7N6Ul2e9e9PMmXJ65u672b1JZUajEVlZWYpm8Lm5\nub7ZJ9hOjsqd7BlLnuP48Y7uTQ0NlmMBAcCvfy0T/C9TjaQ+vV6PnTt3Wt2+fft2F0Tk2dgzlnzD\nle5NOTlyft1a96b0dNm96frrXRMjtWu5unzzL5qbm1WOhDpjoqdeMRqNyMvLQ0tLC4KCgpCZmemc\nP8lbWoCNG2WCt7Kmnt2b3FNQF0tVg4ODVY6EOmOi93G9SdzW5l+vfO+wZH/mTEf3ppoa5bhe39G9\nqV+PFTxIZZmZmaioqFDM0WdkZLgwKlKlNoGzDlNYWCiSkpJEYmKiSEpKsv8ReR9TWFgotFqtRakB\nrVbb5XlMSkqy2PfKl16vtz+Yw4eFePRRIYKClOUJgoOFWLZMiG++sf845HSFhYVCr9eLxMREodfr\n+e/SDo7KnR6b6HubpEipt4k7MTHR6v6JiYl9C6CtTYjCQiHmzFEmd0CIUaOEeOEFIerq+v4fSeTB\nHJU7PfZv37y8PIs/DwE5jZCfn++iiDxPb2+cOWz+tbFRTs1ERQF33CFr0XQ2bZosX1BZCfz2t2zR\nR2Qnj52j5919+/U2cds9/3rypOze9PrrXXdvWrkS+NWv+OQqkQN5bKK/cOGC1e0NV6+vpi71NnFf\nueHa5YNMXSkulqtnNm1Sdm8aMqSje1MXpQ+IyD4em+jJfn1J3AsWLLBthc2lS7I8wZo11rs3abWy\ne9Mjj7CxB5GTeWyiHzJkiNXtg5k0esXmxG2rc+dkY+38fFlj5mq33iqnZxYsAPz9HXdcIuqSxyZ6\nPpjhZo4cAXJzgbfeApqaLMf69wceeEBewcfHuyY+Ih/msYmeD2a4ASGATz+V0zNbtyrHQ0I6ujdd\nVY2SiNTjsYm+zzcGyUKfShrY0r1p5UrZvYl/YRG5HKtX+rBel5Q9fVquf//zn4G6OssxPz+5Jn7l\nSjb2IHIQlikmu9lcUrasTF69v/++ZWMPQHZvSk0FMjPZ2IPIwVimmOzW7UNnbW3AJ5/IBL9rl3Kn\nG26Q7fkeewy45honR0pE9mCi92HWVi4NBrC4pkZenVdWKl/0q1/J6pHs3kTkMTy21g3Zb8aMGQj4\nJVmPBvD/AFQBePzIEcskHxAgb6zu3w/s2QMYDEzyRB6E/1p92L69ezHDbMZKAHcBUDy+NGxYR/cm\njUb9AInIIXq8oq+vr4fBYEBUVBSio6NRXFys2CczMxPjxo1DXFwcvrLWDYjcS2sr8Le/IW/fPuwC\nsAiWSf7fAwfKlTVVVcCqVXYneaPRCL1eD51OB71eD6PRaNf7EVHv9HhFn5WVhfnz52PTpk0wm81o\nbGy0GN+6dSt++OEHHDt2DPv378cTTzxh9ZcBuYEzZ4CCAuDVV4GaGkReNbwDwBoAfrfcgm3p6Q45\npCpdqYioW91e0Z8/fx67d+/G0qVLAQABAQEYOnSoxT6ffPIJlixZAgCYPn066uvrUVtb66RwqU/+\n9S9ZIfKGG4DnnrNo0dfs54cCADEA5gH4QavFisxMhx2afQOIXK/bK/rKykqEhIQgNTUV5eXlmDp1\nKnJzczGwUzPm6upqhIeHt/+s0WhQVVWFkSNHWrxXdnZ2+/c6nQ46nc4x/wVk3eXLwPbtcnnkP/+p\nHA8LA1aswK6ICGx5+22ENDdD74Sni/vaN0C1JuREbsRkMsFkMjn+jbtrP1VaWioCAgJESUmJEEKI\nrKws8dxzz1nsc8cdd4gvv/yy/efbb79dHDx40GKfHg5DjnTxohBr1woRGWm9Pd/UqUK8844QLS2q\nhNOXPrNsE0kkOSp3djt1o9FooNFokJCQAAAwGAwoKyuz2Of666/HyZMn23+uqqrC9ddf77BfRGSj\nkyeB//kfIDxcFhI7cqRjrF8/4J57gN27gdJSICVFVpRUQWZmJrRarcW2norPcbqHyLG6nboJDQ1F\neHg4jh49ivHjx6OoqAgxMTEW+yxcuBCvvPIKFi9ejOLiYlxzzTWKaRtyov37ZfXIrro3PfaY7N40\nerRLwutL8Tm2iSRyrB5X3eTn5yMlJQWtra3QarVYv349CgoKAADp6emYP38+tm7dirFjx2LQoEF4\n4403nB60zzObO7o3WVvhNGaMrP2emuoW3Zt629yEvQaIHItFzTzJle5Nr7wip2quptPJ6pF33OHR\n3Zt6XVWTyEuxqJkvOXpUdm96803r3Zvuv19ewU+e7JLwHI29Bogci1f07koI4LPP5PSMtSdJQ0KA\nJ56QX+zeROSVeEXvrZqbO7o3/etfyvFJk+T0zAMPsHsTEdmEid5d1NR0dG86c0Y5fscdsjwwuzcR\nUS8x0bvaV1/Jq/cNG5TdmwYO7OjeNH68a+IjIo/HRO8KbW3AP/4hE/wXXyjHw8M7ujcNG6Z+fETk\nVZjo1XThAvDGG0BeHnD8uHJ8xgw5PZOczMYeROQwzCZqqKwE8vOBdetksu/M3x+49155g3X6dNfE\nR0RejYneWYQAvvxSTs/8/e+ymmRn7N5ERCphone01lbggw9kgj94UDkeGSmv3h96CBg0SP34iMjn\nMNE7Sl1dR/em06eV40lJMsHr9bKaJBGRSpjo7fXNN7I8wd/+Jh926iw4WF65Z2UBV1X99FRsCELk\neZjo++LyZWDHDlmeoKvuTcuXA8uWyVIFXoL9X4k8E2vd9EZjo7xyz80Fvv9eOT5lilwe+etfW23s\n4elXw3q9Hjt37rS6ffv27S6IiMi7sdaNmqqqZGng11+XpYI769cPuPtuOf9+yy1dlifwhqthNgQh\n8ky8K9idkhJZAjgiAnjxRcskP3iwvHr/4Qfgo4+AWbO6rUHjDe3x2BCEyDPxiv5qZjOwZYucf9+3\nTzk+ZoysPZOaKlv12cgbroYzMzNRUVGhaAjSXf9XInI9Jvor6utl96b8fOvdmxIT5fTMnXf2qXuT\nN1wNsyEIkWfizdhjxzq6NzU2Wo4FBsqpm5Ur7e7exPZ4RNRbqt2MjYiIwJAhQ+Dv74/AwECUlJRY\njNfV1eHBBx9ETU0NzGYzfvOb3+CRRx6xOzCnutK9KSdHdm+6+kSOGCE7Nz35pMO6N/FqmIhcpccr\n+tGjR+PgwYMYPny41fHs7Gy0tLRg1apVqKurQ2RkJGpraxHQqfqi21zRNzcD770nE/zXXyvHJ06U\nV+8pKezeREQup+ryyu4OFBYWhsOHDwMALly4gGuvvdYiybuFmhrgtdfkl7XuTQsWyBU0t93G7k1E\n5HV6zMh+fn6YM2cO/P39kZ6ejrS0NIvxtLQ03HbbbRg1ahQaGhrwwQcfWH2f7Ozs9u91Oh10Op1d\ngdvk0KGO7k2trZZjAwcCjzwiyxOwexMRuQGTyQSTyeTw9+1x6ub06dMICwvDmTNnMHfuXOTn52PW\nrFnt43/84x9RV1eHnJwcVFRUYO7cuSgvL8fgwYM7DqLm1E1bG1BYKBO8tRPG7k1E5CEclTt7fGAq\nLCwMABASEoLk5GTFzdi9e/fi3nvvBSBXkYwePRpHjhyxO7Bea2iQnZsiI+WTqlcn+ZtvBjZuBCoq\ngGeeYZInIp/RbaJvampCQ0MDAKCxsRE7d+7EpEmTLPaZMGECioqKAAC1tbU4cuQIxowZ46Rwrfjx\nR+Dpp2Xzjqwsmciv8PcHFi+WDz7t2ydr0AQGqhebBzAajdDr9dDpdNDr9TAaja4OiYgcrNs5+tra\nWiQnJwMAzGYzUlJSkJSUhIKCAgBAeno6fve73yE1NRVxcXG4fPkyXnrppS5X6DiMEMCePXJ6ZssW\n692bli2TFSTDw50biwfzhvo7RNQzz3pgqrUV+PBDmeAPHFCOjx8vl0c+/DC7N9mA1SiJ3JtvVa+s\nq5OVI199FTh1Sjk+d65M8PPmsXtTL3hD/R0i6pl7J/pvv5VX79a6NwUFdXRvmjjRNfF5OG+ov0NE\nPXO/RH/5MrBzp6weaWVaAaGhcu49Pd2ruje5AqtREvkG90n0TU3yyj0nx3r3psmT5dOr991ntXsT\n9R7r7xD5BtffjK2u7uje9PPPV7+wo3tTD409iIi8jeffjC0tldMzH34om310Nngw8Oij8glWNdfk\nExF5IXUT/ZXuTTk5wN69yvHRo2X3pqVLe9W9iYiIuqZeon/5Zdm96cQJ5djs2XJ6ZuHCPnVvIiKi\nrqk3R3/1xsBAWZ5g5UpgyhRnh0BE5HE8d45+xAjg8cdl96ZfCqaR6xiNRuTl5aGlpQVBQUHIzMzk\nqhsiL6Neoo+J6ejeNGCA1V2YdNTFWjdEvkG9RP/1190uj2TSUV9eXp7F+QbkOc/Pz+c5J/Ii6hWG\n6WENfHdJx9N4Sulf1roh8g1u82SstyQdT/rLhLVuiHyD25R69Jak40l/mWRmZkKr1VpsY60bIu/j\nNlf03lJgq7q62ur2qqoqlSPpGWvdEPkG19e66cRoNHp80hkxYgTOnj1rdfuZM2dcEBEReSpHraN3\nq0TvDWJjY/H1118rtk+aNAmHDx92QURE5KkclTvdZo7eW4R18RDYqFGjHHocT1nZQ0Su1+McfURE\nBIYMGQJ/f38EBgaipKREsY/JZMJTTz2FS5cuYcSIETCZTM6I1SNkZmbi8OHDqKmpad8WGhrq0HsN\nnrSyh4hcr8dE7+fnB5PJhOHDh1sdr6+vx/Lly7Fjxw5oNBrU1dU5PEiyxAediKg3bJq66W6O6L33\n3sM999x1N4X5AAAL7ElEQVQDjUYDQN509GV5eXkWV/MAUFNT49Dlld7yzAERqcOmK/o5c+bA398f\n6enpSEtLsxg/duwYLl26hFtvvRUNDQ3IysrCQw89pHif7Ozs9u91Oh10Op3dwbsjNZKwtzxzQESW\nTCaTc6a+RQ9OnTolhBDip59+EnFxcWLXrl0W48uXLxczZswQTU1Noq6uTowbN04cPXrUYh8bDuM1\nkpKSBADFl16vd9gxCgsLhVartXh/rVYrCgsLHXYMInI9R+XOHq/or6wiCQkJQXJyMkpKSjBr1qz2\n8fDwcIwYMQIDBgzAgAEDMHv2bJSXl2PcuHGO/p3kEdR48IsPOhFRb3S7jr6pqQltbW0YPHgwGhsb\nkZSUhN///vdISkpq3+f777/HihUrsGPHDrS0tGD69OnYuHEjoqOjOw7iQ+voAe948IuIXE+VxiO1\ntbVITk4GAJjNZqSkpCApKQkFBQUAgPT0dEyYMAHz5s1DbGws+vXrh7S0NIsk74sWLFjAxE5EboNP\nxhIRuSk+GUtERDZhoici8nJM9EREXo6JnojIyzHRExF5OSZ6IiIvx0RPROTlmOiJiLwcEz0RkZdj\noici8nJM9EREXo6JnojIyzHRExF5OSZ6IiIvx0RPROTlmOiJiLwcEz0RkZdjoici8nJM9EREXq7H\nRB8REYHY2FhMnjwZN910U5f7lZaWIiAgAJs3b3ZogEREZJ+Annbw8/ODyWTC8OHDu9ynra0Nzz77\nLObNm8cm4EREbsamqZueknd+fj4MBgNCQkIcEtTVjEYj9Ho9dDod9Ho9jEajU45DROSNbLqinzNn\nDvz9/ZGeno60tDSL8erqanz88cf47LPPUFpaCj8/P6vvk52d3f69TqeDTqezKUCj0YisrCxUVFS0\nb7vy/YIFC2x6DyIiT2AymWAymRz/xqIHp06dEkII8dNPP4m4uDixa9cui3GDwSCKi4uFEEIsWbJE\nbNq0SfEeNhymS0lJSQKA4kuv1/f5PYmIPIE9ubOzHq/ow8LCAAAhISFITk5GSUkJZs2a1T5+8OBB\nLF68GABQV1eHbdu2ITAwEAsXLnTIL6KWlhar25ubmx3y/kRE3q7bOfqmpiY0NDQAABobG7Fz505M\nmjTJYp/jx4+jsrISlZWVMBgMeO211xyW5AEgKCjI6vbg4GCHHYOIyJt1e0VfW1uL5ORkAIDZbEZK\nSgqSkpJQUFAAAEhPT3d6gJmZmaioqLCYo9dqtcjIyHD6sYmIvIHfL/NAzj2In59dyy6NRiPy8/PR\n3NyM4OBgZGRk8EYsEXk9e3PnFR7zZOyV/1gVfi8REXmVHm/GuhqXVxIR2cftr+jz8vIskjwgE31+\nfr6LIiIi8ixun+i5vJKIyD5un+i5vJKIyD5un+gzMzOh1WottnF5JRGR7dz+ZuyVG65cXklE1Dce\nsY6eiMgX+dw6eiIi6hsmeiIiL8dET0Tk5ZjoiYi8HBM9EZGXY6InIvJyTPRERF6OiZ6IyMsx0RMR\neTkmeiIiL+ezid5kMrk6BAXGZDt3jIsx2YYxqa/HRB8REYHY2FhMnjwZN910k2L83XffRVxcHGJj\nYzFz5kwcPnzYKYE6mjv+j2VMtnPHuBiTbRiT+nqsXunn5weTyYThw4dbHR8zZgx27dqFoUOHYvv2\n7Vi2bBmKi4sdHigREfWNTWWKu6ueNmPGjPbvp0+fjqqqKvujIiIih+mxTPGYMWMwdOhQ+Pv7Iz09\nHWlpaV3u+/LLL+Po0aN4/fXXLQ/i5+eYaImIfIwjyhT3eEW/Z88ehIWF4cyZM5g7dy4mTJiAWbNm\nKfb7/PPPsX79euzZs8cpgRIRUd/0eDM2LCwMABASEoLk5GSUlJQo9jl8+DDS0tLwySefYNiwYY6P\nkoiI+qzbRN/U1ISGhgYAQGNjI3bu3IlJkyZZ7HPixAksWrQI77zzDsaOHeu8SImIqE+6nbqpra1F\ncnIyAMBsNiMlJQVJSUkoKCgAAKSnp+P555/HuXPn8MQTTwAAAgMDrV71ExGRiwg7nTt3Ttxzzz1i\nwoQJIioqSuzbt0/cd999Ij4+XsTHx4uIiAgRHx9v9bXbtm0TkZGRYuzYsWL16tX2huKQmG688UYx\nadIkER8fLxISEpwa0/79+8W0adNEfHy8mDZtmigpKbH6WmedJ3vjUvNcHTp0SNx8881i0qRJ4s47\n7xQXLlyw+lo1P1O2xuSM8/T999+3f57j4+PFkCFDRG5urjh79qyYM2eOGDdunJg7d644d+6c1dc7\n4zzZG5OzPk/W4srJyREffvihiI6OFv369RMHDx7s8vVqnavexNTbc2V3on/44YfFunXrhBBCXLp0\nSdTX11uMP/300+IPf/iD4nVms1lotVpRWVkpWltbRVxcnPj222/tDceumIQQIiIiQpw9e9YhcfQU\nU2Jioti+fbsQQoitW7cKnU6neJ0zz5M9cQmh7rmaNm2a2LVrlxBCiPXr14vnnntO8Tq1P1O2xCSE\n887TFW1tbSI0NFScOHFCPPPMM+LFF18UQgixevVq8eyzzyr2d/Znqi8xCeH883R1XN999504cuSI\n0Ol0XSZVtc+VLTEJ0ftzZVeir6+vF6NHj+5y/PLlyyI8PFz88MMPirG9e/cKvV7f/vOqVavEqlWr\n7AnH7piEkCewrq7O7jhsiWnx4sVi48aNQggh3nvvPZGSkqLYx1nnyd64hFD3XA0dOrT9+xMnTojo\n6GjFPmp/pmyJSQjnnKfOduzYIW655RYhhBCRkZGipqZGCCHE6dOnRWRkpGJ/Z36m+hqTEM4/T1fi\nmjlzpsW27pKqWueqNzEJ0ftzZVetm8rKSoSEhCA1NRVTpkxBWloampqa2sd3796NkSNHQqvVKl5b\nXV2N8PDw9p81Gg2qq6vtCcfumAC55n/OnDmYNm0a/vKXv9gdT3cxrV69Gk8//TRuuOEGPPPMM1i1\napXitc46T/bGBah3rhobGxETE4OPP/4YAPDhhx/i5MmTiteq+ZmyNSbAOeeps/fffx/3338/AHlf\nbeTIkQCAkSNHora2VrG/Mz9TfY0JcP55uhLXAw88YPP+ap2r3sQE9P5c2ZXozWYzysrK8OSTT6Ks\nrAyDBg3C6tWr28c3bNjQ5X+Asx6isicmQD438NVXX2Hbtm149dVXsXv3bqfEtGrVKjz66KPIy8vD\niRMnsGbNGixdulTxWmc+bGZPXIB65+rFF1/E+vXrsXbtWkybNg0XL15E//79Fa9V8zNla0yAc87T\nFa2trfjHP/6Be++9VzHm5+dn9Zw4+wHGvsQEOPc89RRXV1x5rrrT23NlV6LXaDTQaDRISEgAABgM\nBpSVlQGQ/zi2bNmC++67z+prr7/+eosroJMnT0Kj0dgTjt0xAbY9N+ComEpKStpXNRkMBqvHctZ5\nsjcuQN1zFRkZiR07duDAgQNYvHix1b/I1P5M2RIT4JzzdMW2bdswdepUhISEAJBXzDU1NQCA06dP\n47rrrlO8xpmfqb7GBDj3PFmLyxZqnytb9fZc2ZXoQ0NDER4ejqNHjwIAioqKEBMT0/59VFQURo0a\nZfW106ZNw7Fjx/Djjz+itbUVGzduxMKFC+0Jx+6YbHluwJExjRs3Dl988QUA4LPPPsP48eMVr3XW\nebI3LrXP1ZkzZwAAly9fxh//+Mf25bydqf2ZsiUmZ52nKzZs2NA+RQIACxcuxFtvvQUAeOutt3D3\n3XcrXuPMz1RfY3L2ebIWV2eii6f31T5XtsTUp3PVx/sH7Q4dOiSmTZsmYmNjRXJycvsKl0ceeUQU\nFBRY7FtdXS3mz5/f/vPWrVvF+PHjhVarFS+88IK9odgdU0VFhYiLixNxcXEiJibG6TGVlpaKm266\nScTFxYmbb75ZlJWVKWISwnnnyZ641D5Xubm5Yvz48WL8+PHit7/9bfu+rvxM2RKTM8/TxYsXxbXX\nXmuxrPPs2bPi9ttvVyxlVOs89TUmZ56nruLavHmz0Gg0Ijg4WIwcOVLMmzdPEZcQ6p4rW2Lqy7nq\nsagZERF5Np/tMEVE5CuY6ImIvBwTPRGRl2OiJyLyckz0RERejomeiMjL/X+f4lF1vEUhLgAAAABJ\nRU5ErkJggg==\n" } ], "prompt_number": 46 }, { "cell_type": "code", "collapsed": false, "input": [ "# bootsrapping explicitly\n", "f, ax = subplots(ncols=3)\n", "\n", "for axis in ax:\n", " nuclear0 = nuclear.ix[np.random.choice(range(1, len(nuclear)+1), size=len(nuclear), replace=True)]\n", " nuke_lm0 = ols('np.log(cost) ~ date', nuclear).fit()\n", " axis.plot(nuclear0['date'], np.log(nuclear0['cost']), 'ok')\n", " abline(axis, nuke_lm0.params['Intercept'], nuke_lm0.params['date'], 'r', linewidth=3)\n", " \n", "f.set_size_inches(11, 3)\n", "f.tight_layout();" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAxAAAADRCAYAAABcpVRxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xtc1FX++PEXOJO39ZZjakKpo6Z4T1m1tqQ2xRW1r4rZ\narVhEhYI9t3KXGudLIpqLwG6ZW1u7pZmuXYTQ35maGpeCoUyzet+zdx8AAmaJjAyvz9Ow8zAcJ+Z\nz2dm3s/Hw8cDzozzeQ96DvP+nHPeJ8Rms9kQQgghhBBCiAYI1ToAIYQQQgghhP+QBEIIIYQQQgjR\nYJJACCGEEEIIIRpMEgghhBBCCCFEg0kCIYQQQgghhGgwSSCEEEIIIYQQDVZvAlFSUkJsbCwDBgwg\nIiKCXbt21XhObm4uw4cPZ9CgQURFRXkjTiGEDsh4IIRwVt+YUFRUxIQJExg2bBiDBg3i9ddf1yZQ\nIYRHhdR3DsTvfvc7xo4dy5w5c7BarVy4cIEOHTpUPV5SUsKNN97Ipk2bCAsLo6ioCJPJ5PXAhRC+\nJ+OBEMJZfWOCxWKhrKyMZ599lqKiIq677jrOnDmDwWDQMGohRHPV2YNLS0v59NNPWbVqlXqyweAy\nMACsXr2a6dOnExYWBuD2w0JISIin4hVCNIA3zoeU8UAI/6XVmNC9e3cKCgoAOHfuHJ07d66RPMiY\nIIRveWI8qHMJ04kTJ+jSpQtxcXFcf/31xMfHc/HiRZfnHDlyhB9++IFbbrmFkSNH8q9//avWYL31\nZ8mSJV59fV9cQ96DPq7RrNc/cABbz57YwPFnwQJsVqtP34O3+Mt4oMX/Sy2uFajvbcOGDYwfP55r\nr72W8ePHs2HDhoB5bx673v/7f9iGD3cdawCb2Yzt7bexVVa6PF/LMSE+Pp4DBw5w9dVXM3ToUNLT\n0zUdE/zu31qn16rrehs2bMBsNrv8+5rN5mb1Zb28N3+/ls3mufGgzgTCarWSl5fHgw8+SF5eHm3b\ntiUtLc3lORUVFeTl5bFx40Y2bdrEU089xZEjRzwWoBB+4eOP4YYb4D//Ud+HhkJmJvz1r9Cihaah\neYqMB8LbsrKySElJIScnh//7v/8jJyeHlJQUsrKytA5NH/LzYcIEGDcO9u1ztJtMkJEBX38NM2aA\nj+7oN2RMeOaZZxg2bBinT59m//79JCYmcv78eZ/EJ7SRkZHBsWPHXNqOHTtGZmamRhEJb6gzgQgL\nCyMsLIzIyEgAYmNjycvLc3lOeHg448ePp3Xr1nTu3Jmbb76Z/Px870UsdMdisWAymUhLS8NkMmGx\nWLQOybdWrlS/1EtL1fdt28IHH0BSkrZxeZiMB8Lb5INHLU6ehN/9DoYPh02bHO2tW8PixXDsGMyf\nD1dc4dOwGjIm7Ny5kxkzZgDqLnSvXr345ptvfBqn8K2ysjK37ZcuXfJxJMKb6kwgunXrRnh4OIcP\nHwZg8+bNDBw40OU5t99+O9u3b+fy5ctcvHiR3bt3ExER4b2I3fBFpRdvX8Nf34PFYiE1NZXi4mLK\nysooLi4mNTXVa0mErv4dKivhD3+A++4Dq1W1XX01bN8OMTGeuYaO+Mt4UJ0vf96+/rcNtPem1QcP\n3f67nT0Ljz4K/frBP/8J9uUHoaEwdy4cOQJPPw3t23st1ro0ZEzo378/mzdvBuDMmTN888039O7d\n2+ex2un239rPrlXX9Vq2bOm2vVWrVh6/lrcE8r+bx9jqsX//ftvIkSNtQ4YMsU2dOtV29uxZ28sv\nv2x7+eWXq57zwgsv2CIiImyDBg2ypaen13iNBlxG+KnOnTvXWIYL2Ewmk9ahedfFizbbHXfYbOpX\nuvozbJjN9u23Wkfm1f4m44F3bdiwwTZ+/Hjb2LFjbePHj7dt2LBB65B8avz48W7Hk+joaK1D861L\nl2y2P//ZZuvUyXWMAZtt0iSb7auvGvVyWo4JhYWFtkmTJtmGDBliGzRokO3NN9/0aXzC9zZs2GAz\nm83VtueYg2480ytP9bd6y7h6QkhIiEc3bgj96NixI6X2pTtOOnToQElJiQYR+UBhIdx+O3z2maMt\nJgbWrIF27bSL62d67296j08r9vX/zkt4zGYz6enpxNQxoxVIgv5nUFmpxpHFi+H//s/1schIeOEF\nGDu20S+r9z6n9/hE42VlZZGZmcmlS5do1aoV8+fPD44+7Ac81d8kgRDNYjKZKC4udtteWFioQURe\nduiQShaOH3e0JSWpzdKNrGtusVhYtmwZVqsVg8FAUlKSR5Z+6b2/6T0+rURHR5OTk+O2PTs7W4OI\ntBG0Hzw2b1bLlZw3RwOYzfDssxAb2+TN0Xrvc3qPT4hA4qn+Jie5iGZJSkoiNTUVq30PAKoWeGJi\nooZReUluLkydCvaZlZAQePFFSE5u9EvZ9444/9xSU1OrHhPBRzYeKjExMcGRMNjl58PCha6bo0FV\nVvrjHyEhweebo4UQoj51bqIWoj4Wi4XFixdjMpno0KEDJpOJxYsXB96H4FWrYPx4R/LQpg28916T\nkgegaubBmdVqZfny5c2NVPgpb2w8FDqm08pKQgjRELKESYi62GywZAk89ZSjrXt3+PBDGDGiyS/r\nzb0jeu9veo9PK0G//j9YnD2rliRlZIDzrFNoKMyZAxYL9OjR6JfNysoiIyODsrIyWrZsSXJyctX/\nG733Ob3HJ0QgkSVMQnjbpUuqROvq1Y62wYMhKwvCw5v10oZa9ksYjcZmva7wX/YPe0G5/j8YlJXB\n8uWq7OrZs66PTZoEaWlQrQRqQ7lLPu1fy/8fIYQ3yAyEEO4UFan9Dtu3O9omTIC1az1Sc93dHgiD\nweCR5V967296j08Ij7JXVnr8ccdJ9XbNqKzkrL4N+Hrvc3qPT4hA4qn+JnsghKju8GEYM8Y1eZg3\nTy1b8tCBTUGzdyRI2E9j79ixY3Cexi7c27wZRo6Eu+5yTR7MZnUzYvfuZicPIBvwhRC+J0uYhHC2\nbZuaefjhB/V9SAj86U/w0ENNLqFYG4vFIh80A4BU1BI1+LiykmzAF0L4msxACGH3xhtw222O5KF1\na/j3v+F//9fjyYMIHFJRS1TRqLJScnIyZrPZpc1sNjN//nyPXkcIIexkBkIImw2WLlXVT+y6dlVL\nliIjNQtL+Icff/zRbfv58+d9HInQzNmzahN0enrNykpxcfDkk02qrNRQsgFfCOFrsolaBLeyMpg7\nV80+2A0cqCotXXutdnE1g977m97ja6yQOmanAul9Cje8WFnJk/Te5/QenxCBRMq4CtFcP/yg9jts\n2+ZoGzcO3nkHOnTQLi7hV6644grKy8vdtosA5YPKSkIIoWeyB0IEp6NHVaUl5+QhPl7NPEjyIBqh\nXbt2btvbe6hil6gpKyuL6OhooqKiiI6OJisry3cX91FlJSGE0DOZgRDBZ8cOuP12KC52tD3/PDz8\nsGyWFo2WlJTk9kyPxMREDaMKXJodmubjykpCCKFnsgeiHllZWWRkZFBWVkbLli1JTk6WjWn+bM0a\nuPdesC85adUK/vUviI3VNCxP0nt/03t8TWGxWFi+fDkVFRUYjUYSExOlhKuX1HdoWnNZLJaqyloG\ng4HFd9/NQz/8oMYJ5/+3rVurCm2PPuqx82G8Re99Tu/xCRFIZA+ED2h2p0t4ns0Gzzyj1izbXXUV\nfPABjBqlXVwiIMiZHr7jzUPTnM/06AA8Ajzw4ouuT/JRZSUhhNAz2QNRh4yMDJfkAVQCkZmZqVFE\noknKy2HOHNfkYcAA2LVLkodGKikpITY2lgEDBhAREcGuXbvcPm/v3r0YDAbWr1/v4whFoPPmoWnL\nli0j1GrlIeA4sBBwedVJk6CgAP7+d0kehBBBTWYg6uDNO13CR86ehenT4ZNPHG233qoOiOvYUbu4\n/FRKSgoTJ05k3bp1WK1WLly4UOM5ly9fZuHChUyYMEGWJQiPS05O5tixYy43dzxyaFplJVMvXuQP\nQK9qD33RogUjPv5YNkcLIfxbZaXHXkoSiDp4806X8IHjxyEmBg4dcrTNmQMvvSSbHZugtLSUTz/9\nlFWrVgFqo3AHNxWrMjMziY2NZe/evb4OUQQBrxya9vHH8MgjvPrTTy7Nx4BFwCcdO1IoyYMQwp99\n/TXcf7/HXk4SiDp47U6X8L7PPlOVlgoLHW3PPAOPPSaVlproxIkTdOnShbi4OPLz8xkxYgTp6em0\nadOm6jnfffcd77//Plu2bGHv3r21HrLmvF8gKiqKqKgoL0cvAklMTIxn9qHVUlmpEFgKrABsBgOL\nk5Kafy0fys3NJTc3V+swhBB6UFYGzz6rPgNVVHjsZaUKUz2ysrI8e6crSPm0mtXbb8M996hOA9Cy\nJaxaBTNneud6OuOt/vb5558zZswYdu7cSWRkJAsWLKB9+/YsXbq06jkzZszg4YcfZtSoUdx7771M\nnjyZ6dOn+yQ+IRrs5El44gm3lZW2jRzJ7w4c4OzlywFTUUvvfU7v8Qnht7ZvV7MOBw9WNYWAR/qb\nJBDC69xVszKbzaSnp3s2ibDZ4LnnYNEiR5vJBO+/Dzfc4Lnr6Jy3+tv333/PmDFjOHHiBADbt28n\nLS2NDRs2VD2nd+/eVdcuKiqiTZs2vPrqq0yZMsXr8QlRr5ISdScuPd1xgwECvrKSN/tcSUkJc+fO\n5cCBA4SEhLBy5UpGjx7t8pzc3FweeughKioqMJlMNWZHZEwQwsNKS9WKi5dfdm0fNYqQ3buljKvw\nD3VVs/JYAlFRAQ8+qKqj2F13nTpZ2mz2zDWCXLdu3QgPD+fw4cP069ePzZs3M3DgQJfnHD9+vOrr\nuLg4Jk+e7JI8CKGJsjJYvhyefloVVnA2aRKkpUG1/8uiYeorrFBSUkJiYiKbNm0iLCyMoqIijSIV\nIki8+y4kJcHp0462X/xC3Tx54AEweOajvyQQwuu8Xs2qtFQdBLd5s6Nt7FhYvx6uvNIz1xCA2rg6\ne/ZsysvLMZvNrFy5khUrVgCQkJCgcXRCVFNZqQ6PfPxx+M9/XB+LjIQXXpDKSs3QkMIKq1evZvr0\n6YSFhQFgMpncvpbsixKimb77DubPVwmEk9wxY8gdPRqKiuCppzx2OVnCJLzOGyfH2vdUdCotJe2r\nr+jpfNfrnnvg1VeDttKS3vub3uMTAeLnykrs2+fa3ru3uhM3Y0bQFFTwVp/bv38/CQkJRERE1FpY\nwb506cCBA5w/f56UlBTuvvtun8QnRFCorIQVK9SSpXPnHO1du0JmprrB6jTWeaq/yUFywuuSk5Mx\nV1tG1JxqVvY9FWdzcnhx927X5GHpUnj99aBNHoQIegUF8JvfwG23uSYPJhNkZKjNhHfcETTJgzdZ\nrVby8vJ48MEHycvLo23btqSlpbk8p6Kigry8PDZu3MimTZt46qmnOHLkiEYRCxFgvv4abr5ZLeF2\nTh7mzlVjnRdvlMgSJuF1nq7bnpGRwZBjx3gDsN/nKgNeHDyYhU884ZGYhRB+5uRJ+OMf4Z//rFFZ\niYcegkcfBTfnloimCwsLIywsjMjISABiY2NrJBDh4eGYTCZat25N69atufnmm8nPz6dv375ahCxE\nYKitNGu/fvDKKz5ZmlnvDERJSQmxsbEMGDCAiIgIdu3a5fZ5e/fuxWAwsH79eo8HKfxfTEwM2dnZ\n5Obmkp2d3fTN0zYb/3P0KOtwJA/FwG3AR7LfQYjgU1KiznLo10+Va7YnD6GhcN99cOQIpKZK8uAF\nzoUVALeFFW6//Xa2b9/O5cuXuXjxIrt37yYiIkKLcIUIDNu3w/DhqmqcPXkwGNRer/x8n+3rqncG\nor4KCwCXL19m4cKFTJgwQdYxCu+xWiEpiQecKv0cASYCR4FoOSFciOBhr6yUmgo//OD6mFRW8pn6\nCiv079+fCRMmMGTIEEJDQ4mPj5cEQoimqK006+jRat/noEE+DafOTdSlpaUMHz7cpTSjOy+++CJX\nXHEFe/fuZdKkSW4PjlqyZEnV91JhQTTauXNq3bLTibGfAv8D/ICXzpXwE9VPnX3yySd1ncjLhknR\nLFJZqdH03uf0Hp8QmquvNGuLFg1+KU/1tzoTiIZUWPjuu++466672LJlC3PmzGHy5MlMmzbNK8GK\nIHXypLqj+OWXVU3fRUUxz2jkfHm5nBBejd77m97jEzomlZWaRO99Tu/xaclecbCsrIyWLVuSnJws\nv+uCyenTKnGoVpqVyZPVDGx4eKNf0lP9rc4lTPYKC8uWLSMyMpIFCxaQlpbG0qVLq55jb7MHJIOA\n8KgvvlDJw/ffO9qWLKHHkiV8KB8UhAgOBQVqn0P1ss8mk9o4nZAglddEwLFXHHQ+iNX+tSQRAa6y\nUm2GXriwQaVZtVDnDMT333/PmDFjOHHiBADbt28nLS2NDRs2VD2nd+/eVUlDUVERbdq04dVXX3U5\nfVbuLogmef99mDULLl5U3xuN6qTpe+7RNi6d03t/03t8QkekspJH6L3P6T0+rXjjDCXhB77+Gu6/\nH3bscG2fOxeefx46dWrWy/tkBsK5wkK/fv3cVlhw3h8RFxfH5MmTXZIHIRrNZoP0dPjf/3V8aOjU\nSU3hydpmIQJfSYlakpSerjZL24WGQlycqj7So4d28QnhA2XO//edXLp0yceRCJ/QQWnWxqi3ClN9\nFRaE8CirFRYsUGv77Hr3ho0b4brrtItLCOF9UllJiCotW7Z0295KKg4Gnu3b1azDwYOONoNBVV1a\nvBh0+G9e5xImj11EpidFQ5w/D3feqZIFuxtugPfegy5dtIvLz+i9v+k9PqGBykp46y31i1IqK3mc\n3vuc3uPTirs9EMFccTAgebE0q8ViYdmyZVitVgwGA0lJSVgsFt8sYRLCZ06dUncY8/MdbTNnwuuv\n6zLzFkJ4yMcfq70MeXmu7VJZSQQ5e5KQmZnJpUuXpOJgoPFgadbqLBYLqampWK3WqrbU1NTmRFuD\nzED4mYAs6bZvn0oenDvR4sWwdKla8ywaRe/9Te/xCR+Ryko+o/c+p/f4hPAoL5Rmrc5kMlFcXOy2\nvaioSGYggk1AlnTbsEEtW7KfcG4wqM1CcXHaxiWE8I5vv4UnnpDKSiKgBeTNPtE8PizN6jzz4KzC\neXN2M0kC4UcyMjJckgdQCURmZqZ/DkyZmWrDdGWl+r5DB1i/Hm69Vdu4hBCeJ5WVRJAIyJt9onm8\nXJq1OoPB/cd7o9HosWvI+hA/EjAl3S5fhpQUSE52JA89e8Jnn0nyIESgKSuDv/wFzGb1i9J5HJs0\nSS1l+vvfJXkQAaOum30iyJSVgcUCw4a5Jg/9+kFurtoo7eHkASApKalGEmEwGEhMTPTYNSSB8CMB\nUdLtxx9h6lTIyHC0jR4Nu3fDgAHaxSWE8KzKSli9Gvr3h9//3rUsa2QkfPIJfPihlGUVAee///2v\n2/bTzvv8RODbsQOGD1ezq/alQwaD2uOZn+/VynIWi4XFixdjMpno0KEDJpOJxYsXY7FYPHYNSSD8\nSHJyMmaz2aXNbDYzf/58jSJqpNOn4eab1YcGuxkzYMsWuOoq7eISDVZSUkJsbCwDBgwgIiKCXbt2\nuTz+5ptvMnToUIYMGcKNN95IQUGBRpEKTX38sUoSZs92LcvauzesXatuGERFaRWdEF5VW6JQW2Ih\nAkxpqaqi9KtfuZ7rMGqUqjb39NM+qS5psVgoLCykpKSEwsJCjyYPIHsg/Ipfl3TLz1fLFU6dcrQ9\n9pg6MEoqLfmNlJQUJk6cyLp167BarVywb37/We/evdm2bRsdOnQgOzub+++/v0aSIQKYVFYSgm7d\nurmtgNO1a1cNohE+5cXSrHojZVyF9330Edxxh1q+BKoDvfyy2jwkPM5b/a20tJThw4dz/PjxBj3/\n7NmzDB48mFPOSaMX4xMakspKuqb3Pqf3+BorOjqanJwct+3Z1ZNrERh8UJrVU+QgOeEfXnpJdSr7\nZun27WHdOhg3Ttu4RKOdOHGCLl26EBcXR35+PiNGjCA9PZ02bdq4ff5rr73GxIkT3T7mPJUaFRVF\nlCxn8U9SWUmXcnNzyc3N1TqMoJWcnMyxY8dqnCDtN8uN/dCsWbN45513sNlshISEMGPGDFavXu39\nC/uwNKveyAyE8I7Ll9Vdx7/8xdF27bWQlSWbJr3MW/3t888/Z8yYMezcuZPIyEgWLFhA+/btWbp0\naY3nfvLJJyQmJrJjxw46VaswIeNBACgrU3fVUlNdN0eDWqqYlib9XEf03uf0Hl9TZGVl+edyYz80\na9Ys1qxZU6P9t7/9rXeTiIMHVWnW7dtd271UmtVTPNXfJIEQnnfhgto8+f77jrbISPjgA+jWTbu4\ngoS3+tv333/PmDFjOHHiBADbt28nLS2NDRs2uDyvoKCAadOmkZ2dTZ8+fXwWn/CBykp46y1VRcR5\nczSoPv7887I5Wof03uf0Hp/QN6PR6PbgNKPRSHl5uecvWFamZl6fecZRXQlUadZXXvFqdSVP8FR/\nk92rwrP++1/VeZyTh6lTVb1jSR78Wrdu3QgPD+fw4cMAbN68mYHV7jKfPHmSadOm8cYbb7hNHoQf\nk8pKQggdqu3DcKV96XQzWCwWTCYTHTt2xGQy8dqcOZqVZtUbmYEQnvPVVxATAydPOtoefhiee04q\nLfmQN/tbfn4+c+fOpby8HLPZzMqVK1m7di0ACQkJzJ07l3fffZdrrrkGUHeA9uzZ47P4hBfUVVnp\niSdg3jyprKRzeu9zeo9P6Ju3ZiAsFgupqalYrVbaA2nAA9WfNGqUOgxu8OAmX8fXZAmT0JecHHWm\ng30TUYsWsGyZ+nAhfErv/U3v8fkji8XCsmXLsFqtGAwGkpKSml/zWyorBQy99zm9xyf0zVt7IEwm\nE8XFxfwPsAxwKQfhx6VZZQmT0I9XX4WJEx3JQ7t2sGGDJA9C+ID9LllxcTGlpaUUFxeTmpra9ASi\npETNOPTtC6tWOZKH0FC47z44ckRtnpbkQVD/4ZJ2e/fuxWAwsH79eh9HKALd6tWr+e1vf4vRaKRF\nixYYjUaPbKA2lZfzb+BdXJOHjQYDfP21qjDpZ8mDJ8kMhGi6ykpYtEhtnLQLC1OVloYM0S6uIKf3\n/tbc+Lxyt92P2e+SuWsvLCxs+AvVVVkpJkZVVho0qJnRCi14c0z43e9+x9ixY5kzZ07V4ZIdqiWX\nly9fZty4cbRp04a4uDimT5/us/iEaLSfS7Oee/BB2jv9v/wemA/kdu5MYVGRZuE1l5wDIbT1009w\n993w73872q6/Hj78EK6+Wru4REBzXpNql5qaWvVYMHK39hegwrk6SF2kspJootLSUj799FNWrVoF\ngMFgqJE8AGRmZhIbG8vevXtrfS05G0boglNp1vZOza8CjwI/GgwsTkrSKLim8da5MDIDIRrvzBm4\n/XZVdcVuyhRYvRrattUuLgHov781Jz6P3W0PIM36mXz8sdrLkJfn2t67t1rfO2NGwB6CFEy8NSbs\n37+fhIQEIiIiaj1c8rvvvuOuu+5iy5YtzJkzh8mTJzNt2jSfxCdEg9VSmrWoc2fus1rZitqUnZiY\n6Pc3q2QPhNDG11/D6NGuycOCBbB+vSQPwuuafbc9ACUlJWEwuE4mGwwGEhMTa/9LBQXwm9/Abbe5\nJg8mkzpV+uBBuOMOSR5EnaxWK3l5eTz44IPk5eXRtm1b0tLSXJ6zYMEC0tLSqj60SKIgdGfHjlpL\ns5pOneL9khJKSkooLCz0++TBk2QJk2i4jz+G6dOhtFR9HxqqPmz42XSe8F+XL1922+6Jet/+yv4L\nbfny5VRUVNR9l0wqKwkPCgsLIywsjMjISABiY2NrJBBffPEFd955JwBFRUV89NFHGI1GpkyZ4vN4\nhXBRWgqPPQYvv+za7oelWbUgCYRomJUrISEB7HeA27ZVh0fFxGgblwgqV111FT/++KPb9mBmsVjq\nvjNWUqKm59PT1VS9XWgoxMWpO289etT+94Vww/lwyX79+rk9XPL48eNVX8fFxTF58mRJHoT23n1X\n3fw8fdrR9otfqCVMDz4Y1NWVGkoSCFG3ykp4/HH14cPu6qtVpaVhw7SLSwSl8PBwlw8kdmFhYRpE\n4wfKyuBvf4Onn5bKSsIrMjMzmT17tsvhkitWrADU4ZJC6Mrp0ypxePdd1/ZJk9RYGR6uTVx+SDZR\ni9r99BPcey+8/bajbdgwVWlJPrDplt77W3Pii46OJicnx217dvWTkoOZVFYSTgJ5TBBNl5WVRUZG\nBmVlZbRs2ZLk5GRiAnVVwc+lWVm40HFmFUDXrpCZCbGxQbPnS8q4Cu8qLFSVlj77zNEWEwNr1qiD\n4oTQQHJyMseOHePYsWNVbWazmfnz52sYlc5IZSUhRD2ysrJISUlxGUvtXwdcEuFUmtXF3LnqZkqn\nTtrE5edkBkLUdOiQShacl4okJcFf/6oqEwhd03t/a258WVlZZGZmcunSJVq1asX8+fMD7xeekwbf\nJSwoUHfXqs/EdO4Mf/yjOhn+iit8E7TQlUAfE0TjBcVsblmZWqb5zDNQXu5o79tXzUYE6SyszEAI\n78jNhalT1aZLUHcqX3wRkpM1DUsIu5iYmIBOGJw16C6hVFYSQjRSmXMxBSeXLl3ycSResmMHxMer\n2Qc7g0HdZHn8cWjVSrvYAkS950CUlJQQGxvLgAEDiIiIYNeuXS6Pv/nmmwwdOpQhQ4Zw4403UlBQ\n4LVgRcPMmjULo9GIwWDAaDQya9ashv3FVatg/HhH8tCmDbz3niQPQmgkIyPDJXkAlUBkZmaqfrpw\nobqbtmqVI3kIDYX77oMjRyA1VZIHIUQNLVu2dNveyt8/WJeWwgMPwK9+5Zo8jBqllnU+/bQkDx5S\n7wxESkoKEydOZN26dVitVi5cuODyeO/evdm2bRsdOnQgOzub+++/v0aSIXxn1qxZrFmzxqXN/v3q\n1avd/yWbDZYsgaeecrR17642S48Y4a1QhRD1cHeX8Apg0tGjYDZLZSUhRJME5H4yKc3qU3XugSgt\nLWX48OFuyya6c/bsWQYPHsypU6dcLxISwpIlS6q+j4qKIipI1555m9FodHtar9FopNx5DaDdpUvq\nbqVzcjEEd40vAAAgAElEQVR4sCrTKuXM/EJubi65ublV3z/55JO6Xk8s650bzmw2V42/IcCdQCrQ\nq/oTpbKSqIPe+5ze4wtUAbOfTEqzNoqn+ludCcT+/ftJSEggIiKC/Px8RowYQXp6Om3atHH7/D/9\n6U8cPnyYV155xSvBivoZDAa3p/W2aNGiZmJRVKT2OzhXJpgwQR0Q1769lyMV3qL3/qb3+PSkffv2\nnD9/nluB54Ea84FSWUk0gN77nN7jEzpVV2nWjAwZF2vhqf5W5x4Iq9VKXl4eDz74IHl5ebRt27bG\nMfV2n3zyCStXruS5555rdlCi6UJq6SyhodX+qQ8fhjFjXJOHefPUsiVJHoTQhcE2GxuBj3FNHopD\nQtSp0gcPwh13yC9JIURwOXgQxo5V+x2ck4e5c2Vc9JE6E4iwsDDCwsKIjIwEIDY2lrzqtcWBgoIC\n4uPj+eCDD+gk9XQ1NWPGDLftsbGxjm+2bVPJw9Gj6vuQEPjzn9VUn5RpFXWor6gCqLW1ffv2ZejQ\noezbt0+DKAPAt99CXByf/vgjv3Fq/gm1hCmyUydV3EDKsgohgklZGTz5pDrU1vkGaN++8Mkn8Oqr\ncq6Dj9T5abFbt26Eh4dz+PBh+vXrx+bNmxk4cKDLc06ePMm0adN444036NOnj1eDFfWzb5Ret24d\nlZWVhIaGEhsb69hA/cYbMGcOVFSo71u3hjffVEuZgkhQncDpQfUVVdi4cSNHjx7lyJEj7N69mwce\neECKKjRGSYnaBJ2eDpcuVd3huQy8DiwBzhgMLPbnjY5CCNEUUppVV+o9SC4/P5+5c+dSXl6O2Wxm\n5cqVrF27FoCEhATmzp3Lu+++yzXXXAOozbp79uxxvYisb9SezQZLl4LF4mjr2lUtWfp5hilYuKut\nbzabSU9PD4gkwlv9rSFFFebNm8ctt9zCzJkzAejfvz9bt26la9euLvFJUYVqysrUDODTT9eorHS4\nb1/uKyzkS5sNo9FIYmIiFud+LEQ1UlhBBJTSUnjsMXj5Zdf2UaPUjMPgwdrE5ad8sonaU2Rw0FhZ\nmVoX+MYbjraBA1WlpWuv1S4ujQT6CZze6m8NKaowefJkFi1axA033ADAbbfdxnPPPccIp3LAMh44\nqayEt96CxYvhP/9xfWzkSHjhBamsJJpN731O7/EJDUlpVo/zySZqEQB++EEdDuecPIwbp6YCgzB5\ngCA4gdNLGlpUofrAVNvG/qD38cdq9m/2bNfkoXdvlVTs3i3JgxAiOJ0+DdOmqT/OycOkSfD11zB/\nviQPGpMEIpAdPao2S2/b5miLj1czD0F8Om3AnsDpZQ0pqtCjRw++/fbbqu9PnTpFjx49fBqn7hUU\nwG9+A7fdpk5Gtevc2VFZaeZMdaK0EEIEk8pKtVRpwADXcx26dlUl5j/4QM510An5DaUzWVlZREdH\nExUVRXR0NFlZWU17oR07YPRoVa7V7vnnYcUKMBo9E6yfSk5Oxmw2u7T5/QmcPuBcVAFwW1RhypQp\n/POf/wRg165ddOzY0WX/Q1D7ubISw4aB81K5Vq3gD3+AY8ekspIQInhJaVa/IjU7dcTd5l77143a\n3LtmDdx7L9hPnm7VCv71L3Au5RrE7D/LgDiB08cyMzOZPXu2S1GFFStWAKqowsSJE9m4cSN9+vSh\nbdu2/OMf/9A4Yh2oVlmpSmio6qdPPglhYZqFJ4QQmiorU2PkM884PreAKs36yiuylFOnZBO1j9VV\nPrTZm3ttNtUBH3/c0XbVVWrKb9QoT70FoXN67296j89j6qisREyM+oU5aJA2sYmgovc+p/f4hBdJ\naVaf81R/0/UMhMViYdmyZVitVgwGA0lJSX5dvrC+GYZmbe4tL4eEBHj9dUfbgAFqv0OvXs0JWwjR\nGFJZSQgh6ialWf2ebvdAWCwWUlNTKS4uprS0lOLiYlJTU/06gcjIyHBJHkAlEJmZmUAzNveePQsT\nJrgmD7feCjt3SvIghC9t2SKVlYQQoi7vvgsREa7Jwy9+ARkZakZCkge/oNsEwj7z4MxqtbJ8+XKN\nImq++mYYmrS59/hxuOEGdYS7XVwcfPQRdOzY7JiFEA1gr6z0619LZSUhhHBHSrMGFN0uYaqePNhV\nVFT4OBLPqW+GodGbe3ftgilToLDQ0ZaaCosWSaUCIXzh22/hj3+EVavUHiS7Vq3gf/8XHn00qEsm\ni8BXUlLC3LlzOXDgACEhIaxcuZLRo0dXPf7mm2/y/PPPY7PZaNeuHS+99BJDhgzRMGLhc5WVajP0\nwoWu1ZW6dlWzDjNmyGcWP6TbBMJgcB+a0Y9LkCYnJ3Ps2DGXZUzVZxhiYmIaVg3onXfgnnscVV1a\ntlQfYmbO9HTYQojqpLKSEACkpKQwceJE1q1bh9Vq5cKFCy6P9+7dm23bttGhQweys7O5//772bVr\nl0bRCp87eBDuvx+2b3dtnztXlZbv1EmbuESz6TaBGD9+PGvWrKnRPm7cOA2i8QyPlA+12VSne+wx\nR5vJBO+/r5YyCSG8RyorCVGltLSUTz/9lFWrVgHqxl+HajNuY8aMqfp61KhRnDp1yqcxCo1IadaA\np9sEori42G37D9V/afuZBs8wuFNRAQ8+CH//u6PtuutUpaVqeyeEEB4klZWEqOHEiRN06dKFuLg4\n8vPzGTFiBOnp6bRp08bt81977TUmTpzo9jHnAilRUVFESX/yX1KaVVdyc3PJzc31+Ovq9hyIqKgo\ntm7dWqN97NixXvlB6F5pqToIbvNmR9vYsbB+PVx5pXZxCd3Re011vcdXw5Yt8MgjrpujQVU4e/ZZ\ntX5XNkcLHfNWn/v8888ZM2YMO3fuJDIykgULFtC+fXuWLl1a47mffPIJiYmJ7Nixg07Vlq343Zgg\n3CstVXswX3rJtV1Ks+qKp/qbbn/rNbmkaSD6z3/U8iTn5OGeeyAnR5IHIbzlyy9h4sTaKysdOqSb\nykpZWVlER0cTFRVFdHQ0WVlZWockgkBYWBhhYWFERkYCEBsbS171RBsoKCggPj6eDz74oEbyIAKE\nvTSrc/IgpVkDmm6XMDVkw3FQ2LNHVVo6c8bRtnSpmgaUqgVCeJ6fVVaq74BKIbylW7duhIeHc/jw\nYfr168fmzZsZOHCgy3NOnjzJtGnTeOONN+jTp49GkQqvOX0akpJUAuFs0iS1Xyw8XJu4hNfpdgkT\nqF+Mzdpw7O/Wr4e77oKfflLfX3EFrFypDqkSohZ6Xw6g2/j8tLJSdHQ0OTk5btuzs7M1iEjojTf7\nXH5+PnPnzqW8vByz2czKlStZu3YtAAkJCcydO5d3332Xa665BlCVFPfs2eOz+ISXSGlWv+Wp/qbr\nBCJo2Wzw5z+rO532n9uVV8J778FNN2kbm9A9vfc33cXn55WVZL+YqI/u+lw1eo9PVFNbadb77lMF\nJWSZmq4F/B6IQGWxWDCZTHTs2BGTyeRSeQIAqxUeeEBt2rT/A/fpow6Nk+RBCM+prITVq6F/f7U0\nyTl5GDlSne6+YYOukweQ/WJCCB8pK1MzscOGuSYPffuq8fLvf5fkIYhIAuFDFouF1NRUiouLKS0t\npbi4mNTUVEcSce6cWje4YoXjL/3qVyp56NtXk5iFCEhbtkBkpFoO6FyWtVcvVa51926/KcuanJyM\nuVoZ56DcLyaEqDJ27FhCQkKq/owdO7Z5L7hjBwwfDhaL41wHg0GVti4o8JvxUniOLGHyIZPJ5PZ8\nC5PJROEXX6jk4csvHQ/Mng2vvaZOmRaigfTe3zSN78sv1Zrdjz5ybe/cWW2cnjdP7TXyM0G/X0zU\nScaE4DJ27Fi2bdtWo/3mm292u9yxTlKaNeDIHgg/1LFjR0pLS2u03/yLX7D1F7+A7793NC5Zov7I\nJiTRSHrvb5rEV1dlpYceUkmFjiorCeFJMiYEl5A6Pjc06uf87ruqwtLp0462tm3V6dKJidCiRTOi\nFFrxVH/TbRnXQGQw1PxxTwFW//gj/PijajAa1TrCe+7xbXBCNEDPnj1p3749LVq0cFtNpaioiLvu\nuovvv/8eq9XKww8/zL333qtNsOC3lZWEEEIzdZVmXb4cfq6oJYKb7IHwoaSkJFo4ZewpwLtAW3tD\np07w//6fJA9Ct0JCQsjNzWXfvn01kgeAZcuWMXz4cPbv309ubi6///3vsVqtvg+0rAz++lcwm+G5\n51yTh5gYyM9XywMleRBCCKWyEl5+GQYMcE0eunaFtWvhgw8keRBVJIHwocjISNq1a0cLIBN4Ead/\ngN694bPPoLkbnYTwsrqmPrt37865n2uCnzt3js6dO7udefOaykpYs8bvKysJIURT3XzzzY1qB1Rp\n1rFjVRVI53Md7rtPPXbHHbKkWriQJUw+lJGRgbWkhPcB5+2NBzp2ZOCuXdCli1ahCdEgISEh3Hbb\nbbRo0YKEhATi4+NdHo+Pj+fWW2/l6quv5vz587z99ttuX8e5fHFUVBRRnqjgsWWLOjvliy9c23v1\ngmefVQcbhco9ExH4cnNzg+oMEIvFwrJly7BarRgMBpKSkmqWSA8iW7durbGRutYN1GVlapnnM884\nqiuBqvz4yitSXUnUSjZR+9CMMWNYvGsXw5za3gJW3nQTOW4qJgjRFN7sb//973/p3r07hYWFjBs3\njszMTG5yOp/k6aefpqioiBdffJFjx44xbtw48vPzadeunffiC9DKSkJ4it5/BzcnPnt5dOelkgaD\ngcWLFwd1EtEgO3ZAfLyaYbAzGNSNmMcfh9attYtNeI3PDpIrKSkhNjaWAQMGEBERwa5du2o8Jzk5\nmb59+zJ06FD27dvX7KAC0r59vLRvn0vy8DQwCwht06bBLzNo0CCX2s6DZCmG8KHu3bsD0KVLF6ZO\nnVpjH8TOnTuZMWMGoM4i6NWrF9988413gvn2W4iLg6FDXZOHVq1U2cFjxyA52e+Th1mzZmE0GjEY\nDBiNRmbNmqV1SELohn3mwZnVamX58uUaReQHSkvhwQfVOVPOycOoUZCXB6mpkjyIetWbQKSkpDBx\n4kQOHjxIQUEBAwYMcHl848aNHD16lCNHjvDKK6/wwAMPeC1Yv7VhA9x0E6ayMgAqgDjgCaB3Iw58\nGjRoEAcOHHBpO3DggCQRwicuXrzI+fPnAbhw4QI5OTkMrlYDvH///mzevBmAM2fO8M0339C7d2/P\nBlJSAo89Bv36weuvO8qyhobCnDlw5Iiajg+AsqyzZs1izZo1WK1WLl++jNVqZc2aNZJECPGz2oo0\nVFRU+DgSP/HeexAR4XquQ9u2qlLdjh1yroNosDqXMJWWljJ8+HCOHz9e6wvMmzePW265hZkzZwLq\nA8TWrVvp2rWr4yI6nz71qsxMWLBAbe4EKtq25fEBA9jdtm2jD3zyWG1nEdC81d9OnDjB1KlTAfVL\ne/bs2SxatIgVP5+cnpCQQFFREXFxcZw8eZLKykoWLVpU48Nuk+MrK4O//Q2eftp1czSoykppaQG3\nOdpoNLr9gGQ0Gil3Xq8sRB30/ju4OfHVeUBrYWFzQwsctZVmjYlR46pUVwoaPjkH4sSJE3Tp0oW4\nuDjy8/MZMWIE6enptHFacvPdd98RHh5e9X1YWBinTp1ySSDAS5sm9ezyZVUFJiPD0dazJ8aNG3mu\n2iyOEM3hqw2TvXr1Yv/+/TXaExISqr42mUx8+OGHnr1wZaUqIfiHP8B//uP62MiR8PzzcMstnr2m\nTtQ2yFf+fENCiGCXlJTkdg9EYmKihlHpSGWl2gy9cKFrdaWuXdXnkxkzpLqSaJI6Ewir1UpeXh7L\nli0jMjKSBQsWkJaWxtKlS12eV/2XnLs75UG1menHH2HWLHD+IDV6NLz/Plx1lXZxiYBUPSF/8skn\ntQvG04K8slJts46hAfyehWgM+2eL5cuXU1FRgdFoJDExMbg+c9Tm0CG1SXr7dtf2++6DF15QZ08J\n0UR1/hYKCwsjLCyMyMhIAGJjY8nLy3N5To8ePfj222+rvj916hQ9evTwQqh+4vRpuPlm1+Rhxgz1\nQaiZycPAgQMb1S6E3/ryS5g4EX79a9fkoXNntVb30CGYOTOgkwegakN6dbGxsT6ORAj9slgsFBYW\nUlJSQmFhYcAnD/UWVigrgyefVAUmnJOHvn3VWTh//7skD6LZ6vzt261bN8LDwzl8+DAAmzdvrvFh\ndcqUKfzzn/8EYNeuXXTs2LHG8qWgkZ+vqhg4V6J67DF46y2PVDT46quvavz8Bw4cyFdffdXs1xZC\nF4KkslJDrV69mt/+9rcYjUZatGiB0Wjkt7/9LatXr9Y6NCGEBuotrLBjBwwfDhaL41wHg0EtAc3P\nl3MdhMfUew5Efn4+c+fOpby8HLPZzMqVK1m7di3gWPuclJREdnY2bdu25R//+AfXX3+960V0voHL\nIz76SJ3U+OOP6vsWLdSR8HPnahuXCDp6729u4yspgeeegxdfhEuXHO2hoXDvvepuWliYT+MUIlD4\n5Zgg3KqtsEJng4Gi+HjX6koAv/wlvPoqDBniowiF3nmqv8lBcp7w0kswf77aOA3Qvj2sWwfjxmkb\nlwhKeu9vLvGVlan+89RTQVNZSQhf86sxQdTJYDBw2f5Z42e3A8sBl8XjbduqctaJieqGphA/80kV\nJlGPy5fVBs+//MXRds01kJUlH3qEqEsQV1YSQoimci6s0B1YBkyr/iQpzSp8ILB3IHrThQsQG+ua\nPIwcCbt3S/IgRH1++UtVqcw5eejVS+0X2r1bkgchhHDjhhtuIARIAA5SLXm46io1hn74oSQPwutk\nBqIpvv8eJk+Gzz93tE2dCm+8AU5nZAghalG9stIf/wjz5gXN5mghhGiKPpcv8zRwU7X27B49mFBQ\nAFdeqUVYIgjJDERjffWVqrTknDz8/vfwzjuSPAjRGLVUVrJYLJhMJjp27IjJZAr4koz1ycrKIjo6\nmqioKKKjo8nKytI6JCGqlJSUEBsby4ABA4iIiGDXrl01npOcnEzfvn0ZOnQo+5yrFIqG+7k060s7\nd7okD4eBKCCtTx9JHoRPyQxEY+TkqDMd7Kc5tmgBy5apO6dCiIabM8dtZSWLxVLjVNnU1NSqx4JN\nVlYWKSkpHDt2rKrN/nVMTIxWYQlRJSUlhYkTJ7Ju3TqsVisXLlxweXzjxo0cPXqUI0eOsHv3bh54\n4AG3SYaow44d6kC4gwexz9FWAM8DTwOXgOhWrTQLTwQnqcLUUK++Cg884Ki01K4dvP02TJigbVxC\nVKP3/lZXfCaTieLiYrfthYWF3g5Nd6Kjo8nJyXHbnp2drUFEwh95a0woLS1l+PDhHD9+vNbnzJs3\nj1tuuYWZM2cC0L9/f7Zu3epyXpTexyzNlJaqWdpqpVn3t2zJPWVlfPnz92azmfT0dLmpIBpEqjD5\nSmWl6sDPP+9oCwtTlZakrrIQHuWuvjlARUWFjyPRh7KyMrftl5zPyhBCIydOnKBLly7ExcWRn5/P\niBEjSE9Pp43Tct7vvvuO8PDwqu/DwsI4depUjQNnnWcYo6KiiAr2A8/ee0+VYD192tH2c2nW73r2\n5Oq//Y0rL12iVatWzJ8/X5IHUavc3Fxyc3M9/rqSQNTlp5/g7rvh3/92tF1/vapwcPXV2sUlRIAy\nGNwPSUaj0ceR6EPLli3dtreS5QpCB6xWK3l5eSxbtozIyEgWLFhAWloaS5cudXle9budzqVI7YJx\niaJbp0+rc6XWr3dtdyrNGgPETJmiSXjC/1RPyJ988kmPvK5soq7NmTOqlKRz8jBlCmzbJsmDEF6S\nlJRUI4kwGAwkJiZqFJG2kpOTMZvNLm1ms5n58+drFJEQDmFhYYSFhREZGQlAbGwseXl5Ls/p0aMH\n3377bdX3p06dokcPlyPPBKjVDi+/DAMGuCYPUppV6JTMQLjz9dcq23euUb9gAfzpT3KioxBeZL8L\nuXz5cioqKjAajSQmJgbt3Un7soTMzEwuyXIFoTPdunUjPDycw4cP069fPzZv3szAgQNdnjNlyhSW\nLVvGnXfeya5du+jYsWON5UtB79AhtUl6+3bX9jlz4IUXpLqS0CXZRF3dxx/D9Olq8xJAaCikp0NS\nkrZxCdFAeu9veo9PiEDjzT6Xn5/P3LlzKS8vx2w2s3LlStauXQtAQkICoGYWs7Ozadu2Lf/4xz+4\n/vrrfRafrpWXQ1oapKaqr+369oUVK+RATeEVnupvkkA4W7kSEhLAvpGzbVtYu1bNRgjhJ/Te3/Qe\nX3UWi4Vly5ZhtVoxGAwkJSU1eUbEk68lREPpvc/pPT6v2LED7r9frXiwMxjg0Ufh8cehdWvtYhMB\nTaoweVJlpeqwzz7raLv6alVpadgw7eISQmd69uxJ+/btadGiBUajkT179tR4Tm5uLg899BAVFRWY\nTCavVH/wFU+eSyFnXAghaivNyi9/qcrFS3VH4SdkBuKnn+Dee9WZDnbDhqkNS9UOuRLCH3izv/Xq\n1YsvvviCK2tZk1tSUsKNN97Ipk2bCAsLo6ioCJPJ5LP4PM2T51LIGRdCK3rvc3qPz2PqKM1KYqLs\nsRQ+ITMQnlBYCLffDp995miLiYE1a9RBcUKIGuoaeFavXs306dMJ+zn5rp482PlLzXdPnkshZ1wI\nX/FW3XfRRA0ozSqEvwneGYhDh1TndT5BMykJ/vpXtQ5RCD/lzf7Wu3dvOnToQIsWLUhISCA+Pt7l\ncfvSpQMHDnD+/HlSUlK4++67fRafp8kMhAgEeu9zeo+vySor1bKkhQsdhVlAlWbNyIA77gA3Z2II\n4U2e6m/BeQ5Ebi6MGeNIHkJCVKWlzEzNkweLxYLJZKJjx46YTCZZGy10ZceOHezbt4+PPvqI5cuX\n8+mnn7o8XlFRQV5eHhs3bmTTpk089dRTHDlyRKNom8+T51LIGRdCBJFDh2DsWJg3zzV5mDMHDh6E\nmTMleRB+Lfhuta9apeot25cNtGmjlizp4FRH2WQp9K579+4AdOnShalTp7Jnzx5uuummqsfDw8Mx\nmUy0bt2a1q1bc/PNN5Ofn0/fvn21CrlZPHkuhZxxIUQQqK00a58+8MorUppVBIzgWcJks8GSJfDU\nU4627t3VZukRI7SLy4kscRCe4K3+dvHiRS5fvky7du24cOEC48ePZ8mSJYwfP77qOYcOHSIpKYlN\nmzZRVlbGqFGjWLt2LREREV6PTwjhnt77nN7jazApzSr8gGyiboxLl+C++2D1akfb4MGqTGt4uHZx\nVSObLIWenTlzhqlTpwLq/+rs2bMZP348K1asANShUf3792fChAkMGTKE0NBQ4uPjXZIHIYQIOFKa\nVQShwJ+BKCqCqVNdj4ifMEEdENe+vTYx1UJmIIQn6P1unt7jEyLQ6L3P6T2+OklpVuFnZBN1Qxw+\nrDZLOycP8+apZUs6Sx5ANlkKIYQQfuH0aZg+Xd2gdE4eYmLUEqbkZEkeREAL3CVM27apjv3DD+r7\nkBD405/goYd0W/lANlkKIYQQOialWYUAAnUJ0xtvqFJp9r0DrVvDm2+qhEKIAKf35QB6j0+IQKP3\nPqf3+KocOqSqODqvagD1eeOFF+DKK7WJS4hGkCVM7ths8OSTcPfdjuSha1fYulWSByGEEEI0Xnk5\nLF0KQ4e6Jg99+sCWLfDaa7pLHmbNmoXRaMRgMGA0Gpk1a5bWIYkAEzhLmMrKYO5cNftgN3CgqrR0\n7bXaxSWEEEII/7Rzp5p18KPSrLNmzWLNmjUubfbvVztXoxSiGQJjCdMPP6gZhm3bHG3jxsE770CH\nDt67rhA6pPflAHqPT4hAo/c+p8v4SkvhD39QpVmdY3NTmjUrK4uMjAzKyspo2bIlycnJxMTEaBC0\nYjQa3ZaFNxqNlDsfbieCks/OgejZsyft27enRYsWGI1G9uzZ4/J4UVERd911F99//z1Wq5WHH36Y\ne++9t9mBNdjRo6rqweHDjrb4eFi+HIxG38UhhBBCCP/XiNKsWVlZpKSkcOzYsao2+9daJRG1fTis\nrKz0cSQikNU7A9GrVy+++OILrqxlfZ/FYqGsrIxnn32WoqIirrvuOs6cOeNSjtRrdxd27ID/+R91\n1oPdc8/BI49IFQQRtHR5N8+J3uMTItDovc/pJr7Tp2H+fFi/3rU9Jgb+9je45poafyU6OpqcnBy3\n7dnZ2d6KtE4yAyHq4tNN1HVdqHv37pw7dw6Ac+fO0blz5xpnGXjFW2/Br3/tSB5atVJLlh59VJIH\nIYQQQjRMZSWsWAEREa7Jw1VXqc8aH37oNnkAKCsrc9t+6dIlb0TaIDNmzHDbHhsb6+NIRCCr95N+\nSEgIt912Gy1atCAhIYH4+HiXx+Pj47n11lu5+uqrOX/+PG+//bbb13E+yyAqKoqoqKimRWyzqWnE\nxx93tHXpAh98AKNHN+01hfBjubm55Obmah2GEEL4n0OH4P774dNPXdsbWJq1ZcuWbttbtWrlqQgb\nzb5Ret26dVRWVhIaGkpsbKxsoBYeVe8Spv/+9790796dwsJCxo0bR2ZmJjfddFPV408//TRFRUW8\n+OKLHDt2jHHjxpGfn0+7du0cF/HU9GR5OSQkwOuvO9r694eNG6FXr+a/vhABQDfLAWqh9/iECDR6\n73OaxFdeDmlpkJqqvrbr0wdeeQVuuaVBL+NuD4TZbCY9PV3TjdRC1MZnm6i7d+8OQJcuXZg6dSp7\n9uxxSSB27tzJ4sWLAdVpevXqxTfffMPIkSObHZyLs2fVsfGffOJou/VWWLcOOnXy7LWEEEIIEZhq\nK836yCPwxBONKs1qTxIyMzO5dOkSrVq1Yv78+ZI8iIBXZwJx8eJFLl++TLt27bhw4QI5OTksWbLE\n5Tn9+/dn8+bN3HjjjZw5c4ZvvvmG3r17ezbK48fVJqZDhxxtcXHw8stwxRWevZYQQgghGkT3lRqd\nnTsHixY1qDRrY8TExEjCIIJOnUuYTpw4wdSfT3C2Wq3Mnj2bRYsWsWLFCgASEhIoKioiLi6OkydP\nUhkIX64AABJwSURBVFlZyaJFi2qceNis6ZJdu2DKFCgsdLSlpqpBwA83S+utXrQIPLJcQQjhzJt9\nTteVGp299x4kJcF33znaainNKkQg88kSpl69erF///4a7QkJCVVfm0wmPvzww2YH4tY778A994C9\nmkHLlrBqFcyc6Z3reZke60ULIYQQzVFfpcaCggLAx5Ua7ZpQmlUIUT8f9uJGsNng+efhscccbSYT\nvP8+3HCDdnE1U0ZGhkvyACqByMzMlARCCCGE39FdpUa7ykq1LGnhQnWqtN1VV0FGBtxxh1+uYhCi\nsbxVqbHeKkweuUhjpksqKuDBB+Hvf3e0XXcdZGWB2eydAH0kKiqKrVu31mgfO3aslOEUHqP3JUJ6\nj0+IQOPNPqerSo12zSzNKkQg8+lBcj5TWgoTJ7omD2PHqooJfp48gD7rRQvRGD179mTIkCEMHz6c\nX/7yl7U+b+/evRgMBtZXXzYghAgo7io1Otu5c2fVwWbOlRq9orwcli6FoUNdk4c+fWDLFnjtNUke\nhPAQ/SQQ//mPWp60ebOj7Z57ICcnYDp8cnIy5mqJkNlsZv78+RpFJETjhISEkJuby759+2p8ULC7\nfPkyCxcuZMKECR65y2GxWDCZTHTs2BGTyeSy1EEIoZ2LFy9y/vx5gKpKjYMHD3Z5jr1SI+C9So2g\nbjQOHw5LljjOdTAYVMGVgoIGn+sghGgYfeyB2LNHVVo6c8bRtnSpOm06gNYoSr1oEQjqSwoyMzOJ\njY1l7969tT6noeudLRYLqampWK3WqrbU1NQaryGEcPDV6fRnzpypUalx/PjxLpUa//CHPxAXF8fQ\noUOprKzk+eefr7ViU5PUVpo1MlKtZmhiaVYhRN203wOxfj3cdRf89JP6/oorYOVKmD3b22EJEZC8\nud65d+/edOjQodYNk9999x133XUXW7ZsYc6cOUyePJlp06Y1OT6TyURxcbHb9kLn0s5CiFrpfd9R\nk+OT0qxCNJrPTqL2GpsN/vxnePRRx12DK69UA4LTBiwhhH7s2LHDZcNk//79XTZMLliwgLS0tKoB\nqrmDlPPMg7OKiopmva4Qwo9JaVYhNKdNAmG1qrsGP09zAmqT08aN0LevJiEJIernbsOkcwLxxRdf\ncOeddwLqBNqPPvoIo9HIlClTmnS92urFG43GJr2eEMKPSWlWIXTD95uoz52DSZNck4df/UqdOC3J\ngxC61ZANk8ePH+fEiROcOHGC2NhYXnrppSYnDwBJSUk1kgiDwUBiYmKTX1MI4YcOHYKoKJg3zzV5\nmDMHDh5UB8xK8iCEz/h2BuLkSZU8fPmlo232bFVarZYSp0IIfWjIhklPs1gsHD58mHfeeQebzUZI\nSAgzZsyQDdRCBIvyckhLg9RUR3UlUKsWXnlFqisJoRHfbaL+/HOVPHz/veOBJUvUH7lrIITHBNKG\nyaysLFJSUlxOcDebzaSnp0v1MiEayG/HhJ07IT4evv7a0WYwwCOPwBNPQOvWvgtSiADhqfHAdwlE\nmzZw8aJqMBpVebV77vH2pYUIOn77YcGN6OhocnJy3LZnZ2d7OjQhApLfjQlSmlUIr/G/Kkz25KFT\nJ3j3XXXCtBBC1KGsrMxt+6VLl3wciRDCJ95/X5VgrV6aNTVVFV+R0qxC6IJv90D07q0qLV13nU8v\nK4TwTy1r2RvVqlUrH0cihPCq06chORn+/W/X9okTVWnWa6/VJi4hhFu+q8J0ww2q0lKQJw8WiwWT\nyUTHjh0xmUyyGVSIOiQnJ2M2m13azGYz8+fP1ygiIYRXRES4Jg9XXQVvvQUbNkjyIIQO+W4PxE8/\nQbW7hllZWWRkZFBWVkbLli1JTk4O6I2RFouF1NRUl8OxDAYDixcvlkRCeIzfrXeuR1ZWFpmZmVy6\ndIlWrVoxf/78gB4nhPA0vxgTnBvmzIEXXlCHywohPMr/NlFXu0wwVlcxmUwUFxe7bS8sLNQgIhGI\n/OLDgo7jEyLQ6L3PVSUQUppVCK/z1Hjg+4PkfpaRkeGSPAAcO3aMzMxMjSLyPueZB2cVFRU+jkQI\nIYTQkUWLoKCgRvIwa9YsjEYjBoMBo9HIrFmzNApQP7KysoiOjiYqKoro6GiysrK0DkkEId9uonYS\njNVVqp+oa2c0Gn0ciRBCCKEjzzxTo2nWrFmsWbPGpc3+/erVq30Slt64W71h/zpQV28IfdJsBiIY\nq6skJSXVSCIMBgOJiYkaRSSEEELo0zvvvOO2fd26dT6OxEHru//BuHpD6JNmMxDJyckcO3asxh4I\nX1ZX8fUmbvtG6eXLl1NRUYHRaCQxMVE2UAshhBDV1LZOu7Ky0seRKHq4+x+MqzeEPmmWQNg7m1bV\nVbQaCCwWiyQMQgghRD1CQkLctoeGarN4oq67/7767BKMqzeEPmmWQID6oK7Vmj09DARCCCGEcG/G\njBk19kAAxMbGahCNPu7+62H1hhCgcQKhJT0MBEIIIYRwz75Ret26dVRWVhIaGkpsbKxmG6j1cPdf\n69UbQtgFbQKhh4FACCGEELVbvXq1biou6eXuv5arN4SwC9oEQi8DgRBCCCH0T+7+C+GgWRlXT8rN\nzW3034mJiSE9PZ3o6GjGjh1LdHR0nadgN+UajeHt1/fFNeQ96OcawsGXP29f/9vKe/O/a2lxPeHQ\n3J99TEwM2dnZ5Obmkp2dXW/yUNf1PF0SNpD/H8t70596E4iePXsyZMgQhg8fzi9/+Uu3z8nNzWX4\n8OEMGjSIqKgoT8dYr6b+8BszEMgHV+1f3xfXCIT34G31jQlvvvkmQ4cOZciQIdx4440UFBRoEKVD\nIP8ikPfmf9fS4nre5A+fEarHoofr2StB5uTksHXrVnJyckhJSWlWEqGX9+bv1/L19fx1PKh3CVNI\nSAi5ublceeWVbh8vKSkhMTGRTZs2ERYWRlFRkceDFELoR31jQu/evdm2bRsdOnQgOzub+++/n127\ndvk4SiGEL8hnhKaRSpDC3zVoD0Rth7mA2uA0ffp0wsLCADCZTJ6JTAihW3WNCWPGjKn6etSoUZw6\ndcoXIQkhNCKfERpPKkEKfxdiq6vno+4mdujQgRYtWpCQkEB8fLzL4w899BAVFRUcOHCA8+fPk5KS\nwt133+16kVoOgxFCeEc93bpZ6hsTnP3pT3/i8OHDvPLKK1VtMh4I4XveGhPkM4IQ/scT40G9MxA7\nduyge/fuFBYWMm7cOPr3789NN91U9XhFRQV5eXl8/PHHXLx4kTFjxjB69Gj69u3r0UCFEPpQ35hg\n98knn7By5Up27Njh0i7jgRCBQz4jCBGc6t1E3b17dwC6dOnC1KlT2bNnj8vj4eHhjB8/ntatW9O5\nc2duvvlm8vPzvROtEEJz9Y0JAAUFBcTHx/PBBx/QqVMnX4cohPAR+YwgRHCqM4G4ePEi58+fB+DC\nhQvk5OQwePBgl+fcfvvtbN++ncuXL3Px4kV2795NRESE9yIWQmimIWPCyZMnmTZtGm+88QZ9+vTR\nIkwhhA/IZwQhgledS5jOnDnD1KlTAbBarcyePZvx48ezYsUKABISEujfvz8TJkxgyJAhhIaGEh8f\nL4ODEAGqIWPC0qVLOXv2LA888AAARqPR7SyFEMK/yWcEIYKYrZnOnj1rmz59uq1///62AQMG2D77\n7DPbzJkzbcOGDbMNGzbM1rNnT9uwYcPc/t2PPvrIdt1119n69OljS0tL8/jrX3vttbbBgwfbhg0b\nZouMjGzUe9i9e7dt5MiRtmHDhtlGjhxp27Nnj0ffQ0NfvznvYf/+/bbRo0fbBg8ebJs8ebLt3Llz\nHn0PDX39hryHQ4cOVf2bDhs2zNa+fXtbenq6rbi42HbbbbfZ+vbtaxs3bpzt7NmzTXoPzX39pr6H\nF1980fbOO+/YIiIibKGhobYvvvjC7d9t6ntozOs39P9Sc3h7PPDk9Rr78/D2OOHJ63nivXly/PDk\n9Rr73rw9tnjyWp54b54cczzBl2OCL8eD2q7nrTHBl+NBbdfz1pgQqOOBJ67XmPfn6/Gg2QnEPffc\nY3vttddsNpvNVlFRYSspKXF5/Pe//73tqaeeqvH3rFarzWw2206cOGErLy+3DR061Pb111977PVt\nNputZ8+etuLi4ia9h7Fjx9qys7NtNpvNtnHjRltUVJRH30NDXr+572HkyJG2bdu22Ww2m23lypW2\nJ554wqPvoSGv35j3YHf58mVbt27dbCdPnrQ98sgjtueee85ms9lsaWlptoULFzb5PTT19Zv7Hg4e\nPGj75ptvbFFRUbV23ua8h4a8flPeQ1N4ezzw1PVstsb/PLw9Tnjqep56b54cPzx1vaa8N2feHlua\ncy1PvjdvjDlN5csxwZfjQW3X89aY4MvxoLbreWtMCIbxoCnXs9ma/v58MR40K4EoKSmx9erVq9bH\nKysrbeHh4bajR4/WeGznzp226Ojoqu+fffZZ27PPPuux17fZ1A++qKioSe/hzjvvtK1du9Zms9ls\nq1evts2ePduj76Ehr9/c99ChQ4eqr0+ePGmLiIjw6HtoyOv///bOJhTaLwzj138jZSlGDBuZSVOm\nGFLsWMhCFPlY+cjGVhb28rESRUkpK6RIipSU7NCwJhFNSHZYSO7/Qia888x7nnPucxav+7cyZp77\nmuvp8avDOI9qh6/s7OxQbW0tERGFw2G6u7sjIqLb21sKh8NaHUzm63aoqan59r10P7w6HfzMJ/Lf\nwS+2fcCZR+TvfNj2BGceEU83Ln9w5hGZXce23WKSRWTezaZzdHDpBJc+SJdnwwkufZAuz4YTfosP\ndPKI9Pu58MFfd2FKx+XlJXJyctDT04Py8nL09/fj5eUl+fzBwQECgQCKi4v/ODaRSKCwsDD5OBgM\nIpFIsM0HPvaWrq+vRywWw/z8vK8O4+PjGBwcRFFREYaGhjA2NsbaQWW+SYfn52dEIhFsbGwAAFZX\nV3Fzc8PWQXW+aoevLC8vo7OzE8DHZ2wDgQAAIBAI4P7+XquDyXzdDl1dXX99nUkHP/MB/x38YtsH\nnHmAv/Nh2xOceRzdOP3Bmee3209su8UkCzDvZtM5Orh0gksfpMuz4QSXPvDKs+WE3+IDnTxAv58L\nHxgtIN7e3hCPxzEwMIB4PI6srCyMj48nn19aWvIsoHLjGJP5wMf+1CcnJ9je3sbMzAwODg6UMsbG\nxtDX14fp6WlcX19jcnISvb29bB1U55t0mJiYwMLCAmZnZxGLxfD09ISMjAy2DqrzVTt88vr6is3N\nTbS1taV8r6ner5+bEOnM5+zgBVeHdPjpoINtH3DmAf7Oh21PcOZxdOP0B2ee325fse0W0yzATjcv\nXNy8zaUTXPrAK8+WE1z6wCvPlhN+gw908wC9fq58YLSACAaDCAaDqKysBAC0trYiHo8D+Lgo1tfX\n0d7envLYgoKCbyvKm5ub5K3uOeYDavvVe2UcHh4md5dobW1NeaxJB5X5ph3C4TB2dnZwfHyMjo6O\nlL95MemgMl+1wyfb29uoqKhATk4OgI+V+d3dHQDg9vYWubm5Wh1M5pt2UMGkgyp+Ouhg2weceYC/\n82HbE5x5XN24/MGZ57fbV2y7xTSLs5sKut384NIJLn2QLs+GE1z6IF2eDSf8Bh/o5un2c+UDowVE\nXl4eCgsLcXZ2BgDY3d1FJBJJfl1aWor8/PyUx8ZiMZyfn+Pq6gqvr69YWVlBU1MT23yV/anTZZSU\nlGB/fx8AsLe3h1AoxNpBZb5ph4eHBwDA+/s7RkZGkttqcnVQma/a4ZOlpaXkn/gAoKmpCYuLiwCA\nxcVFNDc3a3UwmW/a4SvkccdVkw4q8/120MG2Dzjz/J4P257gzOPqxuUPzjyT69i2W0yzOLt9hcM5\nurh0gksfpMuz4QSXPkiXZ8MJv8EHunm6/Zz5wPd/Zvzg9PSUYrEYlZWVUUtLS3LXg+7ubpqbm/v2\n2kQiQY2NjcnHW1tbFAqFqLi4mEZHR1nnX1xcUDQapWg0SpFIxHO+V8bR0RFVVVVRNBql6upqisfj\nrB1U5pt2mJqaolAoRKFQiIaHh1OeJ5MOKvP9dHh6eqLs7Oxv27U9Pj5SXV3dH1ud6XTQnW/aYW1t\njYLBIGVmZlIgEKCGhgbWDirz/XQwwbYPuPJ0zodtT3DlcXXj9AdXnu51bNstHFmc3TidY4pLJ7j0\ngVeeLSe49IFXni0n/Ms+MMnT6efSB/8ReSxHBEEQBEEQBEEQfmD0ESZBEARBEARBEH4XsoAQBEEQ\nBEEQBEEZWUAIgiAIgiAIgqCMLCAEQRAEQRAEQVBGFhCCIAiCIAiCICgjCwhBEARBEARBEJT5Hw0W\nNXryeXEkAAAAAElFTkSuQmCC\n" } ], "prompt_number": 47 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }