{ "metadata": { "name": "Bootstrap" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": "import seaborn\nseaborn.set()\ncolors = seaborn.color_palette()\nimport moss\nimport random\nimport numpy as np\nimport scipy, scipy.stats\nfrom numpy.matlib import repmat", "language": "python", "metadata": {}, "outputs": [], "prompt_number": 8 }, { "cell_type": "heading", "level": 4, "metadata": {}, "source": "\u041e\u0431\u0449\u0438\u0435 \u043f\u0440\u0438\u043c\u0435\u0440\u044b bootstrap." }, { "cell_type": "markdown", "metadata": {}, "source": "\u0421\u0433\u0435\u043d\u0435\u0440\u0438\u0440\u0443\u0435\u043c \u0434\u0430\u043d\u043d\u044b\u0435." }, { "cell_type": "code", "collapsed": false, "input": "n_obs = 100\nX = column_stack((randn(n_obs, 4), ones(n_obs)))\nw = rand(5)\nnoise = randn(100) * 5\ny = dot(X, w) + noise", "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "markdown", "metadata": {}, "source": "\u041e\u0446\u0435\u043d\u043a\u0430 \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432 \u043c\u0435\u0442\u043e\u0434\u043e\u043c \u043d\u0430\u0438\u043c\u0435\u043d\u044c\u0448\u0438\u0445 \u043a\u0432\u0430\u0434\u0440\u0430\u0442\u043e\u0432." }, { "cell_type": "code", "collapsed": false, "input": "ols_fit = lambda X, y: dot(dot(inv(dot(X.T, X)), X.T), y)\nw_ols = ols_fit(X, y)", "language": "python", "metadata": {}, "outputs": [], "prompt_number": 9 }, { "cell_type": "markdown", "metadata": {}, "source": "\u0421\u0440\u0430\u0432\u043d\u0435\u043d\u0438\u0435 \u0440\u0435\u0430\u043b\u044c\u043d\u044b\u0445 \u0432\u0435\u0441\u043e\u0432 \u0441 \u043e\u0446\u0435\u043d\u043a\u043e\u0439." }, { "cell_type": "code", "collapsed": false, "input": "bar(arange(5) + .1, w, .4, label=\"actual weights\")\nbar(arange(5) + .5, w_ols, .4, color=colors[1], label=\"estimated weights\");", "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAewAAAFXCAYAAABgJ33WAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGQhJREFUeJzt3X1sk9ehx/GfidHKwIEGnCdZQ4rqvqwoG6sWaR0rs4pH\nylWatJC0slCZVkXapmmwkrZsoQpCpClaW2WbtEooNxMTFMF66caLEm3cmSlRtVFE1cpaGVd3FFJl\nxCbUSUNSNhrH94/eRctCEjcnjn3s7+cvHB/7nJy6+dqP4yeuRCKREAAAyGjz0r0AAAAwPYINAIAF\nCDYAABYg2AAAWIBgAwBgAYINAIAFjILd0NCg1atXq6qqaspx4XBYK1eu1MmTJ02mAwAgZxkFu6am\nRm1tbVOOicfjevnll7VmzRrxkW8AAGbGKNjl5eXKz8+fcsyBAwf00EMPqaCgwGQqAAByWkrfw45G\nowqFQtq0aZMkyeVypXI6AACyVkqD3dzcrGeeeUYul0uJRIJD4gAAzJA7lXf+7rvvatu2bZKk/v5+\ndXV1ye12KxAITHqbRCLBK3EAAP5NSoMdCoXG/t3Q0KAHH3xwylhLnxw27+u7lsplQZLX62GfU4w9\nTj32eG6wz6nn9XqmHWMU7Pr6ep05c0YDAwPy+/3asmWLRkZGJEnBYNDkrgEAwL9wZeKf1+SZXOrx\njDn12OPUY4/nBvucesm8wuZMZwAAWIBgAwBgAYINAIAFCDYAABYg2AAAWIBgAwBgAYINAIAFCDYA\nABYg2AAAWIBgAwBgAYINAIAFCDYAABZI6Z/XzAXxeFyXLr2X7mWMs2LFHcrLy0v3MgAAs4hgG7p0\n6T1tP75TC5P4SytzYbjvml6s3i2f7650LwUAMIsI9ixY6PXI87kl6V4GACCL8R42AAAWINgAAFiA\nYAMAYAGCDQCABQg2AAAWINgAAFiAYAMAYAGCDQCABQg2AAAWINgAAFiAYAMAYAGCDQCABQg2AAAW\nINgAAFiAYAMAYAHjYDc0NGj16tWqqqq66fXHjx9XdXW1qqqqFAwGdf78edMpAQDIOcbBrqmpUVtb\n26TXL1++XAcPHtSJEyf0ve99Tzt37jSdEgCAnGMc7PLycuXn5096/X333SePxyNJWrVqlSKRiOmU\nAADknDl9D/vIkSPy+/1zOSUAAFnBPVcTnT59Wq+//roOHTo07Viv1zMHK5od/f2L0r2ECQoKFiW1\nhzbts63Y49Rjj+cG+5x+cxLs8+fPq7GxUW1tbVq8ePG04/v6rs3BqmZHLDaU7iVMEIsNTbuHXq/H\nqn22EXuceuzx3GCfUy+ZJ0QpPyR++fJlbdmyRS+99JJuv/32VE8HAEBWMn6FXV9frzNnzmhgYEB+\nv19btmzRyMiIJCkYDOqVV17R4OCgdu3a9cmEbreOHDliOi0AADnFONgtLS1TXt/c3Kzm5mbTaQAA\nyGmc6QwAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECw\nAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQ\nbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsYBbuhoUGrV69WVVXVpGOe\nf/55VVRUqLq6WufOnTOZDgCAnGUU7JqaGrW1tU16fWdnp7q7u3Xy5Ek1NTVp165dJtMBAJCzjIJd\nXl6u/Pz8Sa8PhULasGGDJGnVqlUaHBzU1atXTaYEACAnpfQ97CtXrqioqGjsclFRkSKRSCqnBAAg\nK6X8l84SicS4yy6XK9VTAgCQddypvPPCwsJxr6gjkYgcx5n2dl6vJ5XLmlX9/YvSvYQJCgoWJbWH\nNu2zrdjj1GOP5wb7nH4pDXYgENCrr76qyspKvfPOO8rPz9eyZcumvV1f37VULmtWxWJD6V7CBLHY\n0LR76PV6rNpnG7HHqccezw32OfWSeUJkFOz6+nqdOXNGAwMD8vv92rJli0ZGRiRJwWBQfr9fnZ2d\nWrdunRYsWKA9e/aYTAcAQM4yCnZLS8u0Y3bu3GkyBQAAEGc6AwDACgQbAAALEGwAACxAsAEAsADB\nBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxA\nsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAAL\nEGwAACxAsAEAsADBBgDAAsbB7urq0vr161VRUaHW1tYJ18diMdXV1emRRx7Rww8/rF//+temUwIA\nkHOMgh2Px9XU1KS2tja1t7ervb1dFy5cGDfm4MGDWrlypY4dO6b9+/frxz/+sUZGRowWDQBArjEK\ndjgcVmlpqUpKSjR//nxVVlYqFAqNG+P1ejU0NCRJGh4e1pIlS+R2u02mBQAg5xgFOxqNqri4eOyy\n4ziKRqPjxjz++OP661//qgceeEDV1dXasWOHyZQAAOQko2C7XK5px+zdu1ef//zn9cYbb+jYsWPa\nvXv32CtuAACQHKNj047jqLe3d+xyJBKR4zjjxrz99tv67ne/K0ljh88vXryoL3zhC5Per9frMVnW\nnOrvX5TuJUxQULAoqT20aZ9txR6nHns8N9jn9DMKdllZmbq7u9XT06PCwkJ1dHSopaVl3Jg77rhD\nf/rTn/TlL39ZV69e1cWLF7V8+fIp77ev75rJsuZULJZ5RwtisaFp99Dr9Vi1zzZij1OPPZ4b7HPq\nJfOEyCjYbrdbjY2Nqqur0+joqGpra+Xz+XT48GFJUjAY1He+8x3t2LFD1dXVSiQSevbZZ7VkyRKT\naQEAyDnGv67t9/vl9/vHfS0YDI79u6CgQHv37jWdBgCAnMaZzgAAsAAfiMZNxeNxXbr0XrqXMWbF\nijuUl5eX7mUAQNoQbNzUpUvv6QcvHddnFxemeyn66MMr+tmz1fL57kr3UgAgbQg2JvXZxYVadOtt\n6V4GAEAEG0AWy7S3diTe3sHMEWwAWSuT3tqReHsHZgg2gKzGWzvIFnysCwAACxBsAAAsQLABALAA\nwQYAwAIEGwAACxBsAAAsQLABALAAwQYAwAIEGwAACxBsAAAsQLABALAAwQYAwAIEGwAACxBsAAAs\nQLABALAAwQYAwAIEGwAACxBsAAAsQLABALAAwQYAwAIEGwAACxBsAAAsQLABALCAcbC7urq0fv16\nVVRUqLW19aZj3nzzTT366KN6+OGHtXnzZtMpAQDIOW6TG8fjcTU1NWnfvn1yHEe1tbUKBALy+Xxj\nYwYHB7V792794he/UFFRkWKxmPGiAQDINUavsMPhsEpLS1VSUqL58+ersrJSoVBo3JgTJ06ooqJC\nRUVFkqSCggKTKQEAyElGwY5GoyouLh677DiOotHouDHd3d368MMPtXnzZm3cuFFHjx41mRIAgJxk\ndEjc5XJNO2ZkZETnzp3TL3/5S12/fl3BYFBf+tKXtGLFiklv4/V6TJY1p/r7F6V7CRMUFCxKag+n\nGpNp31ey31OmsXHNtrHpcSzxWMbMGQXbcRz19vaOXY5EInIcZ9yYoqIi3Xrrrbrlllt0yy23qLy8\nXOfPn58y2H1910yWNadisaF0L2GCWGxo2j30ej1Tjsm07yuZ7ynTTLfHMGfb41jisYybS+YJkdEh\n8bKyMnV3d6unp0c3btxQR0eHAoHAuDGBQEBvvfWW4vG4rl+/rnA4rDvvvNNkWgAAco7RK2y3263G\nxkbV1dVpdHRUtbW18vl8Onz4sCQpGAzK5/NpzZo1qq6u1rx58/TYY48RbAAAPiWjYEuS3++X3+8f\n97VgMDjucl1dnerq6kynAgAgZxkHe67F43FduvReupcx5v33u9O9BABADrAu2JcuvacfvHRcn11c\nmO6lSJI+6PmLbvuPdK8CAJDtrAu2JH12caEW3XpbupchSfrow6ikq+leBgAgy/HHPwAAsADBBgDA\nAgQbAAALWPkeNgAgM2TaJ3ckacWKO5SXl5fuZcw6gg0AmLFM++TORx9e0c+erZbPd1e6lzLrCDYA\nwEgmfXInmxFsZLzE6GjGnaAmWw+5AchcBBsZ7/q1Pv38nd9q4d8y48/7Dfdd04vVu7PykBuAzEWw\nYYWFXo88n1uS7mUAQNrwsS4AACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAAL\nEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxgHOyuri6tX79e\nFRUVam1tnXRcOBzWypUrdfLkSdMpAQDIOUbBjsfjampqUltbm9rb29Xe3q4LFy7cdNzLL7+sNWvW\nKJFImEwJAEBOMgp2OBxWaWmpSkpKNH/+fFVWVioUCk0Yd+DAAT300EMqKCgwmQ4AgJxlFOxoNKri\n4uKxy47jKBqNThgTCoW0adMmSZLL5TKZEgCAnGQU7GTi29zcrGeeeUYul0uJRIJD4gAAzIDb5MaO\n46i3t3fsciQSkeM448a8++672rZtmySpv79fXV1dcrvdCgQCk96v1+uZ9Lr+/kUmS84JBQWLptzD\nf2KfZ2429hizw7bHcbKPnUwz2ZrZ47ljFOyysjJ1d3erp6dHhYWF6ujoUEtLy7gx//qedkNDgx58\n8MEpYy1JfX3XJr0uFhsyWXJOiMWGptxD6ZP/+djnmZuNPYY5Gx/HyTx2Ms1U+8wez45knmAYBdvt\ndquxsVF1dXUaHR1VbW2tfD6fDh8+LEkKBoMmdw8AAP6fUbAlye/3y+/3j/vaZKHes2eP6XQAAOQk\nznQGAIAFCDYAABYg2AAAWIBgAwBgAYINAIAFCDYAABYg2AAAWIBgAwBgAeMTpwAAkpMYHdX773en\nexnjrFhxh/Ly8tK9DCSBYAPAHLl+rU8/f+e3Wvi3zPjDFMN91/Ri9W75fHeleylIAsEGgDm00OuR\n53NL0r0MWIj3sAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADB\nBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAgQbAAALEGwAACxAsAEAsADBBgDAAsbB7urq0vr161VR\nUaHW1tYJ1x8/flzV1dWqqqpSMBjU+fPnTacEACDnuE1uHI/H1dTUpH379slxHNXW1ioQCMjn842N\nWb58uQ4ePCiPx6Ouri7t3LlTr732mvHCAQDIJUavsMPhsEpLS1VSUqL58+ersrJSoVBo3Jj77rtP\nHo9HkrRq1SpFIhGTKQEAyElGwY5GoyouLh677DiOotHopOOPHDkiv99vMiUAADnJ6JC4y+VKeuzp\n06f1+uuv69ChQyZTAgCQk4yC7TiOent7xy5HIhE5jjNh3Pnz59XY2Ki2tjYtXrx42vv1ej2TXtff\nv2hmi80hBQWLptzDf2KfZ2429hizg8exGdPHcibucbLfk22Mgl1WVqbu7m719PSosLBQHR0damlp\nGTfm8uXL2rJli1566SXdfvvtSd1vX9+1Sa+LxYZMlpwTYrGhKfdQ+uR/PvZ55mZjj2GOx7E508dy\nJu5xMt9TpknmCYZRsN1utxobG1VXV6fR0VHV1tbK5/Pp8OHDkqRgMKhXXnlFg4OD2rVr19htjhw5\nYjItAAA5xyjYkuT3+yf8IlkwGBz7d3Nzs5qbm02nAQAgp3GmMwAALECwAQCwAMEGAMACBBsAAAsQ\nbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMAC\nBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCwAMEGAMACBBsAAAsQbAAALECwAQCw\nAMEGAMACBBsAAAsQbAAALECwAQCwgHGwu7q6tH79elVUVKi1tfWmY55//nlVVFSourpa586dM50S\nAICcYxTseDyupqYmtbW1qb29Xe3t7bpw4cK4MZ2dneru7tbJkyfV1NSkXbt2mUwJAEBOcpvcOBwO\nq7S0VCUlJZKkyspKhUIh+Xy+sTGhUEgbNmyQJK1atUqDg4O6evWqli1bZjI1AAATJEZH9f773ele\nxjgrVtyhvLw84/sxCnY0GlVxcfHYZcdxFA6Hx425cuWKioqKxi4XFRUpEokYBfujD6/M+Laz7fq1\nmIb7rqV7GWNmcy2Zss/ZvMcXLvzvrN2XqXg8Pis/VGbLbP2Qy5THsZS9j+VM2uPY5f/Rj//7v7Tg\n1oXpXook6Xr/sF6p+4l8vruM78so2C6XK6lxiUTiU93O6/VMel1BwSod2LMoqXnnQjz+oCRl1A86\nn8+X1Hps2eds3WNJ6u/PjD2WpIsXL2rHkRcy4gfd9f5hHXr2P3X33XdPO9aWx7GUnY9l9nh6ye7x\ndIyC7TiOent7xy5HIhE5jjNuTGFhoSKRyJRj/l3fNM/6br21eMrrc10s9tG0Y7xeD/tsIBv3OD+/\nUK8s/km6lzEmP79w2v2zbY8z0Ww8ltnjqSW7x9MxCnZZWZm6u7vV09OjwsJCdXR0qKWlZdyYQCCg\nV199VZWVlXrnnXeUn5/P+9dABsrLy5uVw3YAUsMo2G63W42Njaqrq9Po6Khqa2vl8/l0+PBhSVIw\nGJTf71dnZ6fWrVunBQsWaM+ePbOycAAAcokr8e9vMGeA6Q5xwVwyhxJhhj1OPfZ4brDPqZfMIXHO\ndAYAgAUINgAAFiDYAABYgGADAGABgg0AgAUINgAAFiDYAABYgGADAGABgg0AgAUINgAAFiDYAABY\ngGADAGABgg0AgAUINgAAFiDYAABYgGADAGABgg0AgAUINgAAFiDYAABYgGADAGABgg0AgAUINgAA\nFiDYAABYgGADAGABgg0AgAUINgAAFiDYAABYgGADAGABgg0AgAXcM73hwMCAtm3bpsuXL+u2227T\nT3/6U+Xn548b09vbq+3btysWi8nlcunxxx/XN7/5TeNFAwCQa2b8Cru1tVWrV6/W7373O91///1q\nbW2dMMbtdmvHjh1qb2/Xr371Kx08eFAXLlwwWjAAALloxsE+deqUNmzYIEnasGGDfv/7308Y4/V6\nde+990qSFi5cKJ/PpytXrsx0SgAActaMg/3BBx9o2bJlkqRly5bpgw8+mHJ8T0+P/vKXv+iLX/zi\nTKcEACBnTfke9pNPPqmrV69O+PpTTz017rLL5ZLL5Zr0foaHh7V161Y999xzWrhw4QyXCgBA7poy\n2Pv27Zv0uqVLl6qvr09er1dXrlxRQUHBTcd9/PHH2rp1q6qrq/WNb3wjqUV5vZ6kxsEM+5x67HHq\nscdzg31OvxkfEl+7dq1+85vfSJKOHj160xgnEgk999xz8vl8+ta3vjXjRQIAkOtciUQiMZMbDgwM\n6KmnnlJvb++4j3VFo1E1NjaqtbVVZ8+e1RNPPKF77rln7JB5fX29vv71r8/qNwEAQLabcbABAMDc\n4UxnAABYgGADAGABgg0AgAVmfC7xVOjq6tILL7yg0dFR1dbW6tvf/na6l5R1Ghoa1NnZqaVLl+rE\niRPpXk5W4hz6qfePf/xDTzzxhG7cuKGPP/5YgUBATz/9dLqXlZXi8bhqampUVFSkvXv3pns5WWnt\n2rVauHCh8vLy5Ha7deTIkZuOy5hgx+NxNTU1ad++fXIcR7W1tQoEAvL5fOleWlapqanR5s2b9cMf\n/jDdS8la/zyH/r333qvh4WFt3LhRX/va13gsz6LPfOYz2r9/vxYsWKCRkRFt2rRJZ8+eVXl5ebqX\nlnX2798vn8+n4eHhdC8lqx04cEBLliyZckzGHBIPh8MqLS1VSUmJ5s+fr8rKSoVCoXQvK+uUl5dP\n+KtqmF2cQ39uLFiwQNInJ2eKx+PT/rDDpxeJRNTZ2anHHnss3UvJesl8YCtjgh2NRlVcXDx22XEc\nRaPRNK4IMMc59FNndHRUjzzyiFavXq2vfOUruvPOO9O9pKzzwgsvaPv27Zo3L2NSkZVcLpeefPJJ\nbdy4Ua+99tqk4zLmv8JU5yIHbMQ59FNr3rx5OnbsmLq6unT27Fm9+eab6V5SVvnDH/6gpUuXauXK\nlUm9+sPMHTp0SEePHlVbW5sOHjyos2fP3nRcxgTbcRz19vaOXY5EInIcJ40rAmZuJufQx8x4PB75\n/X79+c9/TvdSssrbb7+tU6dOae3atXr66ad1+vRpbd++Pd3LykqFhYWSpIKCAq1bt07hcPim4zIm\n2GVlZeru7lZPT49u3Lihjo4OBQKBdC8L+NQ4h37qxWIxDQ4OSpL+/ve/649//KNWrlyZ5lVll/r6\nenV2durUqVNqaWnR/fffrxdffDHdy8o6169f19DQkCTpo48+0htvvKG77777pmMz5rfE3W63Ghsb\nVVdXN/axLn6rdvbV19frzJkzGhgYkN/v19atW1VTU5PuZWWVt956S8ePH9c999yjRx99VBLn0J9t\nfX19+tGPfqTR0dGx97K/+tWvpntZwKd29epVff/735f0yaelqqqq9MADD9x0LOcSBwDAAhlzSBwA\nAEyOYAMAYAGCDQCABQg2AAAWINgAAFiAYAMAYAGCDQCABQg2AAAW+D/ceM+BduJ7CgAAAABJRU5E\nrkJggg==\n", "text": "" } ], "prompt_number": 10 }, { "cell_type": "markdown", "metadata": {}, "source": "\u0421\u0440\u0430\u0432\u043d\u0438\u043c \u0431\u0443\u0442\u0441\u0442\u0440\u0430\u043f\u043e\u0432\u0441\u043a\u0438\u0435 \u043e\u0446\u0435\u043d\u043a\u0438 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432 \u043c\u043e\u0434\u0435\u043b\u0435\u0439 \u0441 \u0440\u0430\u0437\u043d\u044b\u043c \u0443\u0440\u043e\u0432\u043d\u0435\u043c \u0448\u0443\u043c\u0430. " }, { "cell_type": "code", "collapsed": false, "input": "X = column_stack((randn(1000, 4), ones(1000)))\nw = rand(5)\n\n# Large sample high noise; N = 1000, noise std = 5\ny_large_n = dot(X, w) + randn(1000) * 5\n\n# Small sample low noise; high signal to noise \ny_lownoise = dot(X[:100], w) + randn(100)\n\n# Small sample, high noise\ny_noisy = dot(X[:100], w) + randn(100) * 5\n\n# Bootstrap \nn_boot = 1000\nw_boot1 = moss.bootstrap(X, y_large_n, n_boot=n_boot, func=ols_fit)\nw_boot2 = moss.bootstrap(X[:100], y_lownoise, n_boot=n_boot, func=ols_fit)\nw_boot3 = moss.bootstrap(X[:100], y_noisy, n_boot=n_boot, func=ols_fit)\n\nw_model1 = mean(w_boot1, axis=0)\nw_model2 = mean(w_boot2, axis=0)\nw_model3 = mean(w_boot3, axis=0)", "language": "python", "metadata": {}, "outputs": [], "prompt_number": 29 }, { "cell_type": "markdown", "metadata": {}, "source": "\u0427\u0435\u043c \u043c\u0435\u043d\u044c\u0448\u0435 \u0443\u0440\u043e\u0432\u0435\u043d\u044c \u0448\u0443\u043c\u0430 \u0438 \u0431\u043e\u043b\u044c\u0448\u0435 \u0434\u0430\u043d\u043d\u044b\u0445 \u0432 \u0432\u044b\u0431\u043e\u0440\u043a\u0435, \u0442\u0435\u043c \u043e\u0446\u0435\u043d\u043a\u0430 \u043b\u0443\u0447\u0448\u0435." }, { "cell_type": "code", "collapsed": false, "input": "barx = linspace(0, 1, 6)[:-1]\nmodels = [w, w_model1, w_model2, w_model3]\ncis = [None] + map(seaborn.ci_to_errsize, [ci1, ci2, ci3], models[1:])\nfor i, model in enumerate(models):\n bar(barx + i, model, 0.2, yerr=cis[i], color=colors[i], ecolor=\"gray\")\nxticks([.5, 1.5, 2.5, 3.5], [\"model\", \"large N\", \"low noise\", \"noisy\"]);", "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAFXCAYAAABtOQ2RAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X101OWd9/HPZEK7eQ4Jk4BG1hIRERZBgtul9ESCRI8h\nIRFrcw6FLc+noFi0WsHDg0bFasW2siwnS8u24V6syy2ChRothbhsEZUC2Rt0dwGBBkMeSAJ5oAiZ\n6/4jh7QRSOaRuSbzfv2Vmbnmur7zuzLzmd/D/H4OY4wRAACwUlSoCwAAANdGUAMAYDGCGgAAixHU\nAABYjKAGAMBiBDUAABaL9ufJ1dXVevLJJ9XQ0CCHw6GHHnpI06dP79Jm7969mj9/vm666SZJUm5u\nrubPn+/PsAAARAy/gjo6OlpLlizR0KFD1draqgceeEDf+MY3lJmZ2aXdmDFjtHbtWr8KBQAgEvm1\n6dvlcmno0KGSpLi4OGVmZqq2tjYghQEAgADuo66qqtInn3yiESNGdLnf4XBo//79Kigo0Jw5c3Tk\nyJFADQkAQK/nCMQpRFtbWzVt2jTNnz9f99xzT5fHWlpa5HQ6FRMTo4qKCr3wwgsqLy/3d0gAACKC\n32vUFy9e1MKFC1VQUHBFSEtSfHy8YmJiJEnZ2dm6ePGimpqauu3z0qV2f8sCAKBX8OtgMmOMnn76\naWVmZuq73/3uVdvU19crNTVVDodDlZWVkqTk5ORu+21sbPOnLOu5XAmqq2sOdRnwAXMX3pi/8NXb\n587lSrjmY34F9b59+7R161YNGTJEhYWFkqRFixapurpaklRcXKzy8nJt3Lixc/P3qlWr/BkSAICI\nEpB91IHWm781Sb3/m2FvxtyFN+YvfPX2uetujZozkwEAYDGCGgAAixHUAABYjKAGAMBiBDUAABYj\nqAEAsBhBDQCAxQhqAAAsRlADAGAxghoAAIsR1AAAWIygBgDAYgQ1AAAWI6gBALAYQQ0AgMUIagAA\nLEZQAwBgMYIaAACLEdQAAFiMoAYAwGIENQAgIMrK1qmsbF2oy+h1CGoAACxGUAMAYDGCGgAAixHU\nAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxQhqAAAs5ldQV1dXa9q0acrLy9OkSZP0q1/96qrtnnvu\nOeXm5qqgoECHDx/2Z0gAACJKtF9Pjo7WkiVLNHToULW2tuqBBx7QN77xDWVmZna2qaio0IkTJ/Tu\nu+/q4MGDWrFihd544w2/CwcAIBL4tUbtcrk0dOhQSVJcXJwyMzNVW1vbpc2OHTtUVFQkSbrjjjt0\n7tw51dfX+zMsAAARI2D7qKuqqvTJJ59oxIgRXe6vra1V//79O2/3799fp0+fDtSwAAD0an5t+r6s\ntbVVCxcu1NNPP624uLgrHjfGdLntcDi67a9v31hFRzsDUZq1XK6EUJcAHzF34Y35C56oqI7P9mAt\n40idO7+D+uLFi1q4cKEKCgp0zz33XPF4WlpalzXo06dPKz09vds+Gxvb/C3Lai5XgurqmkNdBnzA\n3IU35i+43O6OlbJgLOPePnfdfQnxa9O3MUZPP/20MjMz9d3vfveqbSZMmKC33npLknTgwAElJiaq\nX79+/gwLAEDE8GuNet++fdq6dauGDBmiwsJCSdKiRYtUXV0tSSouLlZ2drYqKio0ceJExcTEaOXK\nlf5XDQBAhPArqLOysvTpp5/22G7ZsmX+DAMAQMTizGQAAFiMoAYAwGIENQAAFiOoAQCwGEENAIDF\nCGoAACxGUAMAYDGCGgAAixHUAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxQhqAIgwZWXrVFa2LtRl\nwEMENQAAFiOoAQCwGEENAIDFCGoAQMQKh/31BDUAABYjqAEAsBhBDQCAxQhqAAAsRlADAGAxghoA\nAIsR1AAAWIygBgDAYgQ1AAAWI6gBALAYQQ0AgMX8DurFixdr7Nixys/Pv+rje/fu1ejRo1VYWKjC\nwkKtWbPG3yEBAIgY0f52MGXKFE2bNk0//OEPr9lmzJgxWrt2rb9DAQAQcfxeo87KylJiYmIgagEA\nAF8S9H3UDodD+/fvV0FBgebMmaMjR44Ee0gAAHoNvzd99+T222/Xrl27FBMTo4qKCi1YsEDl5eXB\nHhYAgF4h6EEdHx/f+Xd2draeeeYZNTU1KTk5+ZrP6ds3VtHRzmCXFlIuV0KoS4CPmLvwxvxJUVEO\nSYFfFsHq97Jg9BvsmgMh6EFdX1+v1NRUORwOVVZWSlK3IS1JjY1twS4rpFyuBNXVNYe6DPiAuQtv\nzF8Ht9tIUsCXRbD6lYI3d8Gs2RvdfVHwO6gfe+wxffjhh2pqalJ2drYeeeQRXbp0SZJUXFys8vJy\nbdy4UU6nUzExMVq1apW/QwIAEDH8Duqegnfq1KmaOnWqv8MAABCRODMZAAAWI6gBALAYQQ0AgMUI\nagAALEZQAwAQYGVl61RWti4gfRHUAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxYJ+UQ4AAL6svb1d\nx48f87h9Y2O8GhpaPGp7882D5HT2niswEtQAgOvu+PFjWvPy20pOSg9ov01nazT/iXxlZg4OaL+h\nRFADAEIiOSld/fre2GM7E///JEmOluHBLslK7KMGAMBiBDUAABYjqAEAsBhBDQCAxQhqAAAsRlAD\nAGAxghoAAIsR1AAAWIygBgDAYgQ1AAAWI6gBALAYQQ0AgMUIagAALEZQAwBgMYIaAACLEdQAAFjM\n76BevHixxo4dq/z8/Gu2ee6555Sbm6uCggIdPnzY3yEBAIgYfgf1lClTtG7dums+XlFRoRMnTujd\nd99VSUmJVqxY4e+QAGCNsrJ1Kiu79mcg4C+/gzorK0uJiYnXfHzHjh0qKiqSJN1xxx06d+6c6uvr\n/R0WAICIEPR91LW1terfv3/n7f79++v06dPBHhYAgF7huhxMZozpctvhcFyPYQEACHvRwR4gLS2t\nyxr06dOnlZ6e3u1z+vaNVXS0M9ilhZTLlRDqEuAj5i68BXr+oqIcQek3mIJVszf9NjbGB3Tsv5aS\nEu/xa7NhWfQk6EE9YcIEbdiwQXl5eTpw4IASExPVr1+/bp/T2NgW7LICpr29XcePH/PqOSkp8Wpo\naPGo7c03D5LT2fu+tFw++GbatNkhrsQ7LleC6uqaQ10GfBSM+XO7O7YYhtP/RbBq9qZfTz8DfdHQ\n0OLxa7NhWUjdB7rfQf3YY4/pww8/VFNTk7Kzs/XII4/o0qVLkqTi4mJlZ2eroqJCEydOVExMjFau\nXOnvkFY5fvyYHn15q2KT0gLed9vZWv30iQJlZg4OeN8AgPDgd1CvWrWqxzbLli3zdxirxSalKb7v\njaEuA4g43m7RamxkaxbCT9A3fQNAsBw/fkx/WLRQA2JjPWr/mYf9Vre1Sa/+jK1ZsAJBDSCsDYiN\n1cD48DmQC/AW5/oGQowzWwHoDkENAIDFCGoAACzGPmpENI4aBmA7ghoR7fjxY3py6zLFBfisRK11\nzXqp4FmOGgbgN4IaES/OlaCEG5JDXQYAy3mzBe7ixYuSpKNH/9ej9i7Xndd8jKAGAMADx48f05qX\n31ZyUvfXq5Ck5BsvSJL+rfTDHts2na3R1/8vQQ0AgN+Sk9LVz4MzUZqoRknyqG1POOobAACLEdQA\nAFiMoAZ6Mc56BoQ/9lFHIF+uoe0pfjsMf4TrdcqBYCKoI1CwrqFty/Wz+bAH0JsQ1BGKa2gD1+Y2\nRidPnvCorbe/l2WrE7xFUAPAl9Scb5NefUXtHlzn+tJdoyVJnz39VI9tuc41fEFQA8BVeHqd65OO\njmNyPWnrzZq6L1hb750IagC4TrxZU/cWa+sd3G63V1+GvNl1EcwvWd0hqAHgOvJ0TR2+Odtcp9/8\nuk7JSTUetffmVJ8nqw5pYMYwv+rzBUFtMePlN0NPhepbIQBcD56e5lPy7lSfjWc9C/9AI6gtdr65\nTq/8ul6xSdUB7fdM1SdKzRga0D4BAMFBUFsuGD+jagvRt0IAgPc4hSgAABYjqAEAsBibvoEg8OZA\nQG/PbCXxe9lIsHvYbZKkcYc+DXElCDWCGgiCtvoWrT7zc8Wd6vlnOKMuZEqSntnzskd9t9Y166WC\nZyP+97JApCCogSCJcyUo4YbkHts5jnTsgfKkLYDIwz5qAAAsRlADAGAxv4P6/fff13333afc3FyV\nlpZe8fjevXs1evRoFRYWqrCwUGvWrPF3SAAAIoZf+6jb29tVUlKi9evXKz09XQ8++KAmTJigzMzM\nLu3GjBmjtWvX+lUoAACRyK816srKSg0cOFAZGRnq06eP8vLytGPHjkDVBsBSZWXrVFa2LtRleGX3\nsNs6f/IEhBO/grqmpkYDBgzovJ2enq6amq6np3Q4HNq/f78KCgo0Z84cHTlyxJ8hAcAq4w59ym+d\nEVR+bfp2OBw9trn99tu1a9cuxcTEqKKiQgsWLFB5eXm3z+nbN1bR0eFxMofGxvhQl2CVlJR4uVw9\n/3Y4Kqrjf8eTtt7ypu9wnD/jduvs2TqPane72yVJjY2eXdglMzPToxOpBGv+vO23sTFenwW0gvAW\n6vdfb3/vhYpfQZ2enq7q6r98AJw+fVrp6eld2sTH/2UysrOz9cwzz6ipqUnJydf+zWhjY5s/ZV1X\nDQ0toS7BKg0NLaqra+6xndttJMmjtt7ypu9wnL+2+hY9X/Ga4g57cDKVto7jRR7dtrzHtt6cSCVY\n8+dtv+E4f8EU6vdfb3/vhYpfQT18+HCdOHFCVVVVSktL0/bt27Vq1aouberr65WamiqHw6HKykpJ\n6jakAfSMk6kAkcOvoI6OjtbSpUs1a9Ysud1uPfjgg8rMzNTrr78uSSouLlZ5ebk2btwop9OpmJiY\nK4IcAABcm9+nEM3OzlZ2dnaX+4qLizv/njp1qqZOnervMAAARCTOTAYAgMUIagAALEZQAwBgMYIa\nAACLWXk96qNH/zco/d588yCPTuYAAIAtrAzqR1/eqtiktID22Xa2Vj99osCjkzkg+KZNmx3qEgAg\nLFgZ1LFJaYrve2OoywAiinG7dfLkCY/aXrx4UZLnW7/YmgX4zsqgBmw0+IhLkvS/t9SFuJLgaKtv\nUdWmV9QeG9tj20t3jZYkffb0Uz22rW5rk179GVuzAB8R1ECI2RT8A2JjNTC+53OIn3R0HIfqSVsA\n/iGoAViDYxeAK/HzLAAALMYaNYCgchvj8UFq3gpWv4BNCGoAQVVzvk161bOD1LxVeeaMRqSmBrxf\nwCYENYCg8/QgNW9Vt7UGvE/ANuyjBgDAYgQ1AAAWI6gBALAY+6gBwELjDn0a6hJgCdaoAQCwGGvU\nCBhvLurgCy7sACASEdQImPPNdVp94B3FnQr8z3Baas7pkTtna+DAv+2xrTdXduKEGQBsR1AjoOJc\nCUq4ITng/bbWntPqAz/36EvAqAuZkqRn9rzcY9u6T6vlum2A3/UBQLAQ1Agbnn4JcBzpOPTCk7at\ntef8rstmNl2ZC4BvOJgMAACLsUaNXoe1SFwNP3dCuIqYoA7WEckcjASgN2tvb9fx48c8asuBnMER\nMUF9vrlOr/y6XrFJ1QHt90zVJ0rNGBrQPgHAFsePH9Oal99WclJ6j22Tb7wgSfq30g97bHuy6pAG\nZgzzu75IEDFBLUmxSWmK73tjQPtsO1sT0P4AwDbJSenq58Fnp4lqlCSP2jby2ekxDiYDAMBifgf1\n+++/r/vuu0+5ubkqLS29apvnnntOubm5Kigo0OHDh/0dEgCAiOFXULe3t6ukpETr1q3Ttm3btG3b\nNh09erRLm4qKCp04cULvvvuuSkpKtGLFCn+GBAAgovgV1JWVlRo4cKAyMjLUp08f5eXlaceOHV3a\n7NixQ0VFRZKkO+64Q+fOnVN9fb0/wwIAEDH8CuqamhoNGPCX0y+mp6erpqbrAQK1tbXq379/5+3+\n/fvr9OnT/gwLAEDE8Ouob4fD4VE7Y4xXz2s7W+tzTddyvrlBkmf12tBvMPsOZr+tdc0B71eS2hpa\nJQ//34LVr2lPlSQ1f94U0H69FcxlUd12yaO2A/d+JEk66UHbuvN/VrDeI8HqO9z6laRTra1yevjb\nZG9/79zk4RHaybEd/z9Njad6bNvcXO/xkghWv8Hs25t+e1q+fgV1enq6qqv/8rvk06dPKz2962/t\n0tLSuqxBX63Nl60v+baiowN7OcP29vGSFPDLJAar32D2HW79SlL7fe1B6dubfrdu3SpJ+mneMwHt\n11s2LAtvjGwP3rIIVt/h1q8kpXz2mf5P6QdKTuo5VL35vXNDU7WmzRurr33taz22vfwe+c60nB7b\ntrdnS/LuvRfofoPZtzf99sSvoB4+fLhOnDihqqoqpaWlafv27Vq1alWXNhMmTNCGDRuUl5enAwcO\nKDExUf369eu2X5crw5+yrOdyJaguSGueCJ6oKKeiohzq25erbYWr3vzea2hoCcrvnSUpKcnl0f/9\nP/7jPI/681Yw33tRUR2hG+i+A9mvX0EdHR2tpUuXatasWXK73XrwwQeVmZmp119/XZJUXFys7Oxs\nVVRUaOLEiYqJidHKlSv9LhoAgEjh95nJsrOzlZ2d3eW+4uLiLreXLVvm7zAAAEQkzkwGAIDFCGoA\nACxGUAMAYDGCGgAAixHUAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxQhqAAAsRlADAGAxghoAAIsR\n1AAAWIygBgDAYgQ1AAAW8/t61ECkmDZttlyuBNXVNYe6FAARhDVqAAAsRlADAGAxghoAAIsR1AAA\nWIygBgDAYgQ1AAAWI6gBALAYQQ0AgMUIagAALEZQAwBgMYIaAACLEdQAAFiMoAYAwGI+Xz2rqalJ\nixYt0ueff64bb7xRP/nJT5SYmHhFu5ycHMXFxcnpdCo6OlqbNm3yq2AAwNU1na3xqF1y7KWO9o2n\nAtYngsfnoC4tLdXYsWM1Z84clZaWqrS0VD/4wQ+u2rasrEzJyck+FwkA6N7NNw/S/CfyPWpbUfE7\nSdL9c+/yuG+Ejs9B/fvf/14bNmyQJBUVFWnatGnXDGpjjK/DAAA84HQ6lZk52KO2f/hDhSR53B6h\n5fM+6jNnzqhfv36SpH79+unMmTNXbedwODRjxgw98MADeuONN3wdDgCAiNTtGvWMGTNUX19/xf3f\n//73u9x2OBxyOBxX7WPjxo1KS0tTQ0ODZsyYoUGDBikrK8uPkgEAiBzdBvX69euv+Vhqaqrq6urk\ncrlUW1urlJSUq7ZLS0uTJKWkpGjixImqrKzsMaj79o1VdLSzp9rDmsuVEOoS4CPmLrwxf1JUVMeK\nVbgsi2DWG6y+A9mvz/uoc3JytHnzZs2dO1dvvfWW7rnnnivanD9/Xu3t7YqPj1dbW5t2796thx9+\nuMe+GxvbfC0rLLhcCaqraw51GfABcxfemL8ObnfHcUPhsizcbqOoKEdQ6g3WsvC23+4C3ed91HPn\nztUf/vAH3Xvvvfrggw80d+5cSVJNTU3n3/X19Zo6daomT56shx56SOPHj9e4ceN8HRIAgIjj8xp1\ncnKy/vVf//WK+9PT01VaWipJuummm7RlyxafiwMABN60abNDXQK8wJnJAACwGEENAIDFCGoAACxG\nUAMAYDGCGgAAixHUAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxQhqAAAs5vO5vgEAwNUF8nzqrFED\nAGAxghoAAIsR1AAAWIx91ACAiBUO1+ZmjRoAAIsR1AAAWIygBgDAYgQ1AAAWI6gBALAYQQ0AgMUI\nagAALEZQAwBgMYIaAACLEdQAAFiMoAYAwGIENQAAFiOoAQCwGEENAIDFfA7q3/72t8rLy9PQoUN1\n6NCha7Z7//33dd999yk3N1elpaW+DgcAQETyOahvvfVWrV69WllZWdds097erpKSEq1bt07btm3T\ntm3bdPToUV+HBAAg4kT7+sTMzMwe21RWVmrgwIHKyMiQJOXl5WnHjh0ePRcAAAR5H3VNTY0GDBjQ\neTs9PV01NTXBHBIAgF6l2zXqGTNmqL6+/or7Fy1apJycnB47dzgcPhXVt2+soqOdPj03XLhcCaEu\nAT5i7sIb8xd+oqI6siRS567boF6/fr1fnaenp6u6urrz9unTp5Went7j8xob2/wa13YuV4Lq6ppD\nXQZ8wNyFN+YvPLndRlFRjl49d919CQnIpm9jzFXvHz58uE6cOKGqqip98cUX2r59uyZMmBCIIQEA\niAg+B/V7772n7OxsHTx4UPPmzdPs2bMldeyXnjt3riQpOjpaS5cu1axZs5SXl6f777+fA8kAAPCC\nw1xrdTiEevPmDYnNb+GMuQtvzF94Kitbp6goh6ZOnRXqUoIm6Ju+AQBAcBDUAABYjKAGAMBiPp+Z\nDACA62HatNkRfXwBa9QAAFiMoAYAwGIENQAAFiOoAQCwGEENAIDFCGoAACxGUAMAYDGCGgAAixHU\nAABYjKAGAMBiBDUAABYjqAEAsBhBDQCAxQhqAAAsRlADAGAxghoAAIsR1AAAWIygBgDAYgQ1AAAW\nI6gBALAYQQ0AgMUIagAALEZQAwBgMYIaAACLEdQAAFgs2tcn/va3v9Xq1at17Ngxbdq0ScOGDbtq\nu5ycHMXFxcnpdCo6OlqbNm3yuVgAACKNz0F96623avXq1Vq2bFmPbcvKypScnOzrUAAARCyfgzoz\nM9PjtsYYX4cBACCiBX0ftcPh0IwZM/TAAw/ojTfeCPZwAAD0Kt2uUc+YMUP19fVX3L9o0SLl5OR4\nNMDGjRuVlpamhoYGzZgxQ4MGDVJWVpZv1QIAEGG6Der169f7PUBaWpokKSUlRRMnTlRlZWWPQe1y\nJfg9ru0i4TX2VsxdeGP+wlekzl1ANn1fax/0+fPn1dLSIklqa2vT7t27deuttwZiSAAAIoLPQf3e\ne+8pOztbBw8e1Lx58zR79mxJUk1NjebOnStJqq+v19SpUzV58mQ99NBDGj9+vMaNGxeYygEAiAAO\nwyHZAABYizOTAQBgMYIaAACLEdQAAFiMoLZATk6Ompqa/G6Dno0aNeq6jvfmm29q6NCh+u///u/O\n+yZNmqTPP//8utYRbq73PH3Z66+/rrfeeiukNaCrn/3sZ9qzZ0+oywgJn08hCvR2l4+zdDgcfvXT\nv39/rV27Vq+++mpA+kPwFRcXh7oEfMnChQtDXULIENQ+qqqq0uzZszVq1Cj98Y9/1PDhw1VUVKTV\nq1ersbFRL7/8sgYOHKglS5aoqqpKMTExevbZZzVkyBA1Njbq8ccfV21trUaOHNml3y1btmjDhg26\nePGiRowYoRUrVigqig0fgdba2qoFCxbo3Llzunjxor7//e9rwoQJqqqq0qxZszRy5EgdOnRIpaWl\n2rx5s95++22lpKRowIABGjZsmGbOnKmTJ0/q2WefVUNDg2JiYlRSUqJBgwZ1GcfhcOjuu+/Wxx9/\nrM8++0xf+9rXQvSKw5MxRi+99JJ2794tSfre976n+++/X88884y++c1vKicnRwsWLFBSUpJeeOEF\nbdq0SX/605+0aNGiLv2MGjVK06dP165du/Q3f/M3WrNmjVJTU1VVVaUlS5aoqalJKSkpWrlypQYM\nGKDXXntNcXFxmjlzpn71q1/p17/+tZxOp2655RatWrVKbW1tKikp0ZEjR3Tp0iU9/PDDmjBhQigW\nUdiqqqrSnDlzlJWVpf379ys9PV1r1qzRsWPHtHz5cv35z3/WwIED9cILLygxMVFPPfWUxo8fr3vv\nvVc//vGPtXPnTjmdTo0bN04LFixQQUGBysvLFR0drZaWFk2ePFnvvvuunE5nqF+q/wx88qc//cnc\nfvvt5n/+53+M2+02RUVFZvHixcYYY373u9+Z+fPnm5KSErN69WpjjDF79uwxkydPNsYYU1JSYv7p\nn/7JGGPMrl27zJAhQ0xjY6M5cuSImTdvnrl06ZIxxpjly5ebzZs3G2OMGT9+vGlsbLzeL7PXGTly\npDHGmEuXLpnm5mZjjDFnzpwxEydONMZ0zOttt91mDh48aIwx5uDBg2by5MnmwoULpqWlxeTm5ppf\n/OIXxhhjpk+fbo4fP26MMebAgQNm+vTpV4z35ptvmmeffdZs3rzZ/PCHPzTGGDNp0iRz6tSp4L7Q\nMHd5nt555x0zY8YM43a7TX19vbn77rtNbW2t2bZtm/nRj35kjDFmypQp5tvf/rYxxpinnnrK7N69\n+4r+hgwZYnbu3GmMMeall14ya9asMcYYM2/evM732KZNm8z8+fONMca89tprnfM8btw488UXXxhj\nTOf/zCuvvGK2bNlijDHm7NmzJjc317S1tQV8OfRmlz9DP/nkE2OMMY8++qjZsmWLyc/PNx999JEx\nxpif/vSn5vnnnzfGdMxteXm5aWhoMPfee29nP5fn5KmnnjLvvfeeMcaY119/3bz44ovX8+UEFatq\nfsjIyNDgwYPlcDg0ePBgjR07VlLHJUCrqqq0b98+TZ48WZL09a9/XU1NTWppadHHH3+sgoICSVJ2\ndraSkpJkjNGePXt06NAhTZkyRYWFhfrggw9UVVUVstfXm7ndbq1atUoFBQWaOXOmamtrdebMGUnS\nDTfcoBEjRkiS/vjHP+qee+7RV77yFcXFxWn8+PGSOs60t3//fj366KMqLCzU8uXLr3pe/Mvy8/N1\n4MAB5tNL+/bt06RJk+RwOJSamqoxY8bov/7rvzR69Gjt27dPR48e1eDBg5Wamqq6ujodOHDgqvu3\n+/Tpo7vvvluSNGzYMJ06dUqSdODAAeXn50uSCgoKtG/fviueO2TIED3++OPaunVr59at3bt361/+\n5V9UWFio6dOn6+LFi6qurg7SUui9MjIydNttt0nqmJeTJ0+qubm58zTTRUVF+vjjj7s8JzExUV/9\n6le1ZMkSvffee/rqV78qSfrWt76lN998U5K0efNmTZky5Tq+kuBi07cfvvKVr3T+HRUVpT59+kjq\n2NzpdrvldDqveXrVa91fVFSkxx57LPDFoou3335bjY2N2rx5s5xOp3JycnThwgVJUmxsbGc7h8PR\nZa4u/+3xOG3PAAADlUlEQVR2u5WYmOjxAUdOp1MzZ85UaWlpAF9F73e15e9wOJSenq5z587pP/7j\nP5SVlaWzZ89q+/btio2N7TJ/l0VH/+WjLioqSu3t7V36vJrL95eWluqjjz7Szp07tXbtWr399tuS\npNdee00333xzIF5mxPrrz1Cn06nm5uYuj395bowxcjqd+vd//3ft2bNH5eXl2rBhg375y1/qzjvv\n1KlTp7R37161t7frlltuuS6v4XpgjTqIRo8e3fmm3rt3r1JSUhQfH6+srCz95je/kSRVVFTo7Nmz\ncjgc+od/+AeVl5eroaFBktTU1MTRwUHS0tKi1NRUOZ1OffDBB9dcznfeead27typL774Qq2trdq1\na5ckKT4+XhkZGXrnnXckdXyAfPrpp1c8/68/aIqKirRnz57O+UXPsrKytH37drndbjU0NGjfvn2d\nWztGjhypX/7yl7rrrruUlZWlX/ziFxozZoxX/Y8aNUrbtm2T1PHl7fLzL8+bMUaff/65/v7v/16P\nP/64mpub1dbWpm9+85sqKyvr7Ofw4cOBeLkRLyEhQUlJSZ1r0Vu2bNFdd93VpU1bW5uam5uVnZ2t\nxYsXd3nfFRYW6gc/+EGvWpuWWKMOqC8fzfvwww9ryZIlKigoUGxsrF588cXO+x9//HFNmjRJo0aN\n0g033CBJyszM1KOPPqqZM2fK7XYrOjpaK1as6Hwc/rs8R/n5+fre976n/Px8DR8+XJmZmVdt/3d/\n93fKyclRfn6++vXrpyFDhig+Pl6S9OMf/1grVqzQP//zP+vSpUvKy8vr3Iz31+NdHrNPnz6aPn26\nXnjhhSC+wt7h8jKbOHGi9u/f37kL6YknnlBqaqqkji/C//mf/6mbbrpJ/fv317lz5655Zb6/fm/+\n9d9Lly7V4sWL9fOf/1ypqalauXJlZxuHw6H29nY9+eSTnWt606dPV0JCgubPn6/nn39e+fn5MsYo\nIyNDa9euDfyCiEAvvviili9frvPnz2vgwIGdcyJ1zEtra6vmz5+vCxcuyBijxYsXdz4+adIk/eQn\nP9GkSZNCUXrQcK5voAdtbW2KjY3V+fPn9Z3vfEfPPfechg4dGuqyAHzJO++8o507d+pHP/pRqEsJ\nKNaogR4sXbpUR48e1YULF1RUVERIAxYqKSnR7t27e+VxIKxRAwBgMQ4mAwDAYgQ1AAAWI6gBALAY\nQQ0AgMUIagAALEZQAwBgsf8P5r4ZoZ+Ll9UAAAAASUVORK5CYII=\n", "text": "" } ], "prompt_number": 30 }, { "cell_type": "heading", "level": 4, "metadata": {}, "source": "\u041e\u0446\u0435\u043d\u043a\u0430 \u043a\u043e\u0432\u0430\u0440\u0438\u0430\u0446\u0438\u043e\u043d\u043d\u043e\u0439 \u043c\u0430\u0442\u0440\u0438\u0446\u044b \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432." }, { "cell_type": "code", "collapsed": false, "input": "mas = []\nfor n_boots in xrange(100):\n w_boot = moss.bootstrap(X, y_large_n, n_boot=n_boot, func=ols_fit) \n w_est_boot = []\n for indices in w_boot:\n boot_X = X[indices, :]\n boot_y = y[indices]\n w_est = w_ols\n w_est_boot.append(w_est.tolist())\n w_est_boot = np.array(w_est_boot)\n w_end = w_est_boot - repmat(w_est_boot.mean(axis=0), w_est_boot.shape[0], 1)\n cov_estimation = np.dot(np.transpose(w_end), w_end) / (n_boots - 1)\n diff = np.linalg.norm(cov_estimation - cov)\n mas.append(diff)", "language": "python", "metadata": {}, "outputs": [], "prompt_number": 12 }, { "cell_type": "code", "collapsed": false, "input": "plt.plot(range(2, 100), difference_bootstrap, 'g-')\nplt.legend(('n_boots', 'Bootstrap'))", "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 13, "text": "" }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEACAYAAABcXmojAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXl4VNXdx793sq+ThWQgCyDZIBATFomCSxBBBUFFsYLF\nVKogvEK1Koq1r6hVqNa6oRa0orQWqu2rQg0ISAJSDIskgIBAIIEkhEBWkkyS2e77x3BO7p25syUz\nSSb8Ps/Dwyx37j0zmTnf81uPIIqiCIIgCIIAoOrpARAEQRC9BxIFgiAIgkOiQBAEQXBIFAiCIAgO\niQJBEATBIVEgCIIgOHZFoa2tDdnZ2cjKykJ6ejqWLl2qeNzixYuRkpKCzMxMFBUV8ceXL1+O4cOH\nIyMjA7Nnz0Z7e7vsdW+88QZUKhXq6urc8FYIgiCIrmJXFAIDA5Gfn4/i4mIcOnQI+fn52LVrl+yY\nvLw8lJSU4OTJk1i9ejUWLFgAACgrK8OHH36IAwcO4PDhwzAajVi/fj1/XXl5ObZu3YpBgwZ54G0R\nBEEQncGh+yg4OBgAoNPpYDQaERUVJXt+w4YNyM3NBQBkZ2ejoaEB1dXVCA8Ph5+fH7RaLQwGA7Ra\nLeLj4/nrfvvb3+K1115z53shCIIguohDUTCZTMjKyoJGo8GECROQnp4ue76yshKJiYn8fkJCAior\nKxEVFYUnn3wSAwcORFxcHCIiInDLLbcAAL7++mskJCTg6quvdvPbIQiCILqCQ1FQqVQoLi5GRUUF\ndu7ciYKCAqtjlDplnDp1Cm+99RbKyspw7tw5NDc347PPPoNWq8Wrr76KF1980e7rCYIgiO7H19kD\n1Wo1pk6div379yMnJ4c/Hh8fj/Lycn6/oqIC8fHxKCgowLhx4xAdHQ0AmDFjBnbv3o3MzEyUlZUh\nMzOTHz969Gjs3bsXsbGxsmsmJyfj1KlTXXl/BEEQVxRJSUkoKSnp9OvtWgo1NTVoaGgAALS2tmLr\n1q0YOXKk7Jjp06dj7dq1AIDCwkJERERAo9EgLS0NhYWFaG1thSiK2LZtG9LT0zFixAhUV1ejtLQU\npaWlSEhIwIEDB6wEATBbG6Io0j9RxAsvvNDjY+gt/+izoM+BPgvb/7q6kLZrKVRVVSE3Nxcmkwkm\nkwlz5szBxIkTsWrVKgDA/PnzMWXKFOTl5SE5ORkhISFYs2YNACArKwsPPvggxowZA5VKhVGjRmHe\nvHlW1xAEoUtvgCAIgnAfdkUhIyMDBw4csHp8/vz5svsrV65UfP2SJUuwZMkSuwM4ffq0ozESBEEQ\n3QRVNHsJ0jjOlQ59Fmboc+iAPgv3IYii2GtTfwRBQC8eHkEQRK+jq/Om09lHBEEQrhAVFYX6+vqe\nHkafJTIy0iMtgshSIAjCI9Dv17PY+ny7+rlTTIEgCILgkCgQBEEQHBIFgiAIgkOiQBAEQXBIFAiC\nIOxQUFAg6wTd1yFRIAiC6AHKysqgUqlgMpl6eigySBQIgiB6kN6WtkvFawRBdDvCi+5rhCm+4Nqk\nOnjwYCxatAhr167FmTNncNttt+HTTz9FQECA3dctX74cf/7znxEaGopXXnkFs2fPBgA0NjZi0aJF\n2Lx5M4KDg/HII4/gueee4/UCr7zyCj766CO0trbitttuw7vvvovw8HDceOONAICIiAgAwLZt2xAd\nHY1f//rXOHjwIPz8/DBx4kTZNsbdAVkKBEFcUQiCgC+++ALffvstSktLcejQIXzyySd2X3P+/HnU\n1tbi3Llz+PTTTzFv3jycOHECALBo0SI0NTWhtLQUO3bswNq1a3m36DVr1uDTTz9FQUEBTp8+jebm\nZjz22GMAgO+//x6AWVSampqQnZ2N3//+97jtttvQ0NCAyspKLF682HMfhC3EXkwvHx5BEHborb/f\nwYMHi5999hm/v2TJEvHRRx+1eXx+fr7o6+srarVa/th9990nvvzyy6LBYBD9/f3FY8eO8edWrVol\n5uTkiKIoijfffLP4wQcf8OeOHz8u+vn5iUajUSwtLRUFQRCNRiN//sEHHxTnzZsnVlRUOHwftj7f\nrn7uZCkQBHHF0b9/f347KCgIzc3Ndo+PjIxEUFAQvz9o0CBUVVWhtrYWer0egwYN4s8NHDgQlZWV\nAMx70lg+ZzAYUF1drXid1157DaIoYuzYsRgxYgS3OLoTEgWCIK5onNnoq76+Hlqtlt8/c+YM4uLi\n0K9fP/j5+aGsrIw/d/bsWSQkJAAA4uLirJ7z9fWFRqNRvK5Go8Hq1atRWVmJVatWYeHChd2+5wyJ\nAkEQVzSik9k/L7zwAvR6Pb7//nt88803mDlzJlQqFe677z787ne/Q3NzM86cOYM333wTv/zlLwEA\ns2bNwptvvomysjI0Nzfjueeew/333w+VSoWYmBioVCrZ9plffPEFKioqAJgD0IIgQKXq3mmaso8I\ngriiEQTBrrUgCAIGDBiAyMhIxMXFISQkBKtWrUJqaioA4N1338WiRYswZMgQBAYGYt68eXjooYcA\nAHPnzsW5c+dw4403oq2tjWcfAUBwcDB+97vfYfz48TAYDNi0aRP279+PJ554Ao2NjdBoNHjnnXcw\nePBgj38GsvcrOiuTPQC13iUI74V+v56FWmcTBEEQHodEgSCIK55XX30VYWFhVv+mTp3a00Prdsh9\nRBCER6Dfr2ch9xFBEAThcUgUCIIgCA6lpBIE4REiIyOdKgwjOkdkZKRHzksxBYIgiD4ExRQIgiAI\nt+F1oqA36tGsa4beqO/poRAEQfQ5HIpCW1sbsrOzkZWVhfT0dCxdulTxuMWLFyMlJQWZmZkoKiri\njy9fvhzDhw9HRkYGZs+ejfb2dgDA008/jWHDhiEzMxMzZsxAY2OjUwOe+o+pCFsehvyyfKeOJwiC\nIJzHoSgEBgYiPz8fxcXFOHToEPLz87Fr1y7ZMXl5eSgpKcHJkyexevVqLFiwAIB5D9IPP/wQBw4c\nwOHDh2E0GvkuQpMnT8aRI0dw8OBBpKamYvny5U4N2N/HHwCgM+pceqMEQRCEY5xyHwUHBwMAdDod\njEYjoqKiZM9v2LABubm5AIDs7Gw0NDSguroa4eHh8PPzg1arhcFggFarRXx8PABg0qRJvPtfdnY2\n7wzoCBIFgiAIz+GUKJhMJmRlZUGj0WDChAlIT0+XPV9ZWYnExER+PyEhAZWVlYiKisKTTz6JgQMH\nIi4uDhEREbjllluszv/xxx9jypQpTg2YiQLFFAiCINyPU6KgUqlQXFyMiooK7Ny5EwUFBVbHKKVA\nnTp1Cm+99RbKyspw7tw5NDc347PPPpMd88orr8Df359vgu0IshQIgiA8h0vFa2q1GlOnTsX+/fuR\nk5PDH4+Pj0d5eTm/X1FRgfj4eBQUFGDcuHGIjo4GAMyYMQO7d+/GAw88AAD45JNPkJeXh++++87m\nNZctW8Zv5+TkwE/lB4BEgSAIAgAKCgoUF+qdxtEmzhcvXhTr6+tFURRFrVYr3nDDDeK2bdtkx3zz\nzTfi7bffLoqiKP7www9idna2KIqiWFRUJA4fPlzUarWiyWQSH3zwQXHlypWiKIripk2bxPT0dPHi\nxYs2r600vEc3PipiGcT3977vaOgEQRBXHE5M63ZxaClUVVUhNzcXJpMJJpMJc+bMwcSJE7Fq1SoA\nwPz58zFlyhTk5eUhOTkZISEhfLPprKwsPPjggxgzZgxUKhVGjRqFefPmAQAWLVoEnU6HSZMmAQCu\nu+46vP/++w5FjNxHBEEQnsPr2lw8veVp/OmHP+H1Sa/jqXFP9dDICIIgeidXXJsLshQIgiA8B4kC\nQRAEwfE6UfDzoewjgiAIT+F1okCWAkEQhOcgUSAIgiA4XisK1OaCIAjC/XitKJClQBAE4X68ThR4\nmwsTiQJBEIS78TpRIEuBIAjCc5AoEARBEByvFQUKNBMEQbgfrxUFshQIgiDcD4kCQRAEwfE6UaA2\nFwRBEJ7D60SBLAWCIAjPQaJAEARBcLxWFPQmyj4iCIJwN14rCmQpEARBuB+vEwXe5oJEgSAIwu14\nnSiQpUAQBOE5SBQIgiAIDokCQRAEwfFaUaDeRwRBEO7H60TBV+ULwJySKopiD4+GIAiib+F1oiAI\nAs9AoloFgiAI9+J1ogBQXIEgCMJTkCgQBEEQHK8WBQo2EwRBuBevFgWyFAiCINyLXVFoa2tDdnY2\nsrKykJ6ejqVLlyoet3jxYqSkpCAzMxNFRUX88eXLl2P48OHIyMjA7Nmz0d7eDgCoq6vDpEmTkJqa\nismTJ6OhocGlQdOeCgRBEJ7BrigEBgYiPz8fxcXFOHToEPLz87Fr1y7ZMXl5eSgpKcHJkyexevVq\nLFiwAABQVlaGDz/8EAcOHMDhw4dhNBqxfv16AMCKFSswadIknDhxAhMnTsSKFStcGjRZCgRBEJ7B\nofsoODgYAKDT6WA0GhEVFSV7fsOGDcjNzQUAZGdno6GhAdXV1QgPD4efnx+0Wi0MBgO0Wi3i4+Ot\nXpObm4uvvvrKpUGTKBAEQXgGh6JgMpmQlZUFjUaDCRMmID09XfZ8ZWUlEhMT+f2EhARUVlYiKioK\nTz75JAYOHIi4uDio1WrccsstAIDq6mpoNBoAgEajQXV1tUuDJlEgCILwDL6ODlCpVCguLkZjYyNu\nvfVWFBQUICcnR3aMUmXxqVOn8NZbb6GsrAxqtRozZ87EZ599hgceeEB2nCAIEATB5vWXLVvGb+fk\n5CAnJ4c22iEIgrhMQUEBCgoK3HY+h6LAUKvVmDp1Kvbv3y8Thfj4eJSXl/P7FRUViI+PR0FBAcaN\nG4fo6GgAwIwZM7B792488MAD0Gg0OH/+PPr374+qqirExsbavK5UFBhkKRAEQZhhi2XGiy++2KXz\n2XUf1dTU8Myg1tZWbN26FSNHjpQdM336dKxduxYAUFhYiIiICGg0GqSlpaGwsBCtra0QRRHbtm3j\nrqfp06fj008/BQB8+umnuOuuu1watDMb7YiiiBZdi0vnJQiCuNKxKwpVVVW4+eabkZWVhezsbEyb\nNg0TJ07EqlWrsGrVKgDAlClTMGTIECQnJ2P+/Pl4//33AQBZWVl48MEHMWbMGFx99dUAgHnz5gEA\nnn32WWzduhWpqanYvn07nn32WZcG7chSEEURczfMhXqFGoeqD7l07r5OXWsdHtnwCPZV7uvpoRAE\n0QsRxF7calQQBMV4xYx/zsCXP3+Jf9/3b8wYNsPq+bcL38bj3z4OAFh9x2o8MvoRj4/VW1h7cC1y\nv8rF/SPux7p71vX0cAiCcDO25k1n6XMVzd+f+R5PbX2K369ucS2zSconxZ9gxS7Xaih6O7XaWgDA\n+ebzPTwSgiB6I14tCpa9j6qaqnDfv+6DwWTAQPVAAF2b/J749gks/W4p6lrrOj/YXkZDmzlGdKHl\nQg+PhOjNlNaXokZb09PDIHoArxQFW4HmOV/Owfnm88gZnINXb34VgHOiYBJNVo8ZTAY+gbL/+wLs\nvVQ3d96CIvo2zbpmZHyQgUl/m9TTQyF6AK8UBSX3UZuhDd+VfgcfwQfr71mP+HBz9bQj99Gq/asQ\nsSICxeeLZY9LhaCpvcldQ+9xGtrN76u2tRYGk6GHR0P0Rsoby9Gib8Gxi8dod8MrkD4jCq36VgBA\nqH8oNKEa9A/tD8CxpfD92e/RpGvC7vLdssfrW+v57WZds1vG3RuQvi9yDxBK1Laa407txna06Cmt\n+0qj74iCwSwKgb6BAOC0KLBzWMYN6ts6Js8mnbWl0NDWoOh26u1ILSCKKxBKSBcLtHC48vBqUZC2\nuWgztAEAgvyCAADqADX8ffzRrGu2W8TWbrzczttSFCQrakv3UXljOTR/0uDhDQ934V30DCQKhCNI\nFK5svFoUlNxHQb5mURAEgVsL9uIK7BxSy8DyvqWlcKzmGHRGHX6s+rGzb6HHcEUUjCajp4dD9EKk\nQnCx5WIPjoToCbxSFJQ22WHuI2YpAM65kGy6j+xYCpfaLwHwzlWUVBTsZSAdrzmOiD9GYPn3y7tj\nWDLONJzhnzHR/bBaFsC177jRZOxT8bcrFa8UBXuWAospAIAmxNye297k50xMwfKLLhUFb8rOMJgM\nMqvHnqWQX5aPZl0ztpVu646hcU7Xn0bayjTc8/k93XpdooOa1s65j5ZsXYLo16JxuPqwJ4ZFdBN9\nRhR4TMHXA5aCTtlS0Bl1XpWdYbn6ticKp+tPA+j+eoZNJzeh3diO705/RzGPHqIzMQWjyYhPDn4C\nnVGHz4987qmhEd1AnxGFzrqP2g02As1tjt1HgHe5kKRCBwAXtLYn3VP1pwB0rU1IZ8gvywcAiBDx\nbcm3Tr/uy2NfYvGmxRQHcQOdEYW9lXv5b2jr6a0eGRfRPXi1KEjbXFgGmgGJ+8iJQHNda53MFWQv\n0OytomBZme2MpVCr7b4iN5NoQkFZAb+/qWST069dsm0J3t37Lnad3eX4YMIu0pjCRa1zgea8k3n8\n9r5z+6wWIIT34JWiwNtcmKzdR9KYgivuI51RB61eyx93xn0EeKcosL5QtkRBFEUuCiLEbstAOXLh\nCGpbaxHqHwoA2Fyy2amVv1avxak6s2XDxk10ns5YCnklZlEI9Q+FSTRxi89Vmtqb0NjW2KnXEu7B\nK0XBrvtIIabgjKUAyF1IzrqPpKsqT1DWUIYZ/5zhln0hmCikRacBsC0Kta21svfYXS4kNpHcNfQu\nJEUmob6tHnsq9zh83c81P0OE2cojUegaBpNB9t13RhSqmqpwoOoAgnyDsHjsYgDA1lOuu5BEUcT4\nj8cj8y+ZtKtiD9J3REFvHVPQhJrdR85YCoCFKNhpc9GdlsLr/30dX/78Jd7f936Xz8VEIT48HkG+\nQdDqtYophJYTa3cFm5koTBg8AVNSpgAwB54dceTCEX77dAOJQldgvwGVYJ4anPl+by7ZDACYOGQi\npqVNA4BOZa1VNVfh8IXDONN4hlt+RPfTd0TBjqVwvvm8zdRRVtEM2LEUesh9JIoiNp7YCAA403im\ny+djohAZGInYEPO+2ErWgpUodIOlYBJN2FG2A4BZFG5Pvh1Ah1vCHkcuSkSBLIUuwSzfwRGDzfdb\nax22c/nm5DcAgCnJUzAmbgzUAWqU1JWgrKHMpWtLxf3nmp9lz5lEE94qfIvSXbuBPiMKSjGFUP9Q\nhPiFoM3Qpti/yPIcTBSMJqNs4u+p7KPDFw6j/FI5ALj8A1OCiUJEYIRdUbBcpXWHpXDw/EHUt9Vj\nkHoQroq8CjmDcxDoG4gDVQcc9q+SikJpfamnh9ptiKKINUVrZJOlp2Hf5/6h/RERGAGTaLIbNNYb\n9dhyagsA4PaU2+Gr8sWEqyYAcN2FJP07WorCllNb8MS3T+Cuf95FGWYexitFgVU0K2YfSdxHgGMX\nklQUmHXAJk9mQtuzFFhHSU+w8fhGfvtMw5kuF8qx9+dIFNhq+6qIqwB0zy5t3HV0eUIJ8gvChMHm\n28w9YYufLvzEb1e3VNvtdeVNHKg6gLkb5mLBNwtcet3z259Hzic5/DfhCkwU+gX3Q7/gfrLHlPhv\n+X/RpGtCekw6ty4mDTHvw+Bqaqr07/hzrVwUDlQdAGD+bv7fsf9z6byEa3ilKDjrPgLsZyCZRJMs\n3ZJZCmzyjAuLA2C2FKQTslQkPGkp/Ofkf/jtVkNrl68ltRRYuq6ipXC5RmFc4jgA3eM+ksYTGDyu\nYCc1tVnXjLKGMvip/LiIlTb0DWuh4lIFAOtVsz2a2pvw+u7XsePMDuw7t8/la3JRCHJOFFgq6pTk\nKfwxJgrflX7n0qrenqVw+EKH2+iP//2jV3US8Db6jijYsBR4BpKCC8RyO08uCpfN5diQWPj7+MMo\nGnnsQRTFbnEfXWi5gD0VexDgE4DkqGQAXXchOes+YpbCdQnXAfC8KBhMBuw8sxOAXBRYXGHLqS02\nayWOXTwGAEiNTkVaP3NWVV+JKzAr9KL2otM9hTaXbOa/i5O1Jzt9zX7B/RATHAPASVFI6RCF5Khk\nDFIPQl1rHYrOFzl1XVEUZW6y4zXHZRM/y77zVfnix6ofO53ySjimz4hCm9E6pgB0FLApWQrSIDNg\nbSlEBkYizD8MQEdcoc3QJpugPCUK35z4BiJE3HzVzUiPSQfQ9WCzUqDZUizbDe2ouFQBH8EH18Rf\no3iMuymqKsKl9ktIikxCojqRP54UlYSUqBQ0tDVgX6XyqpetLkfEjsCQiCEA5KIgiiJe2fkK3i58\n2+t2mpOmOzsbK/nq+Ff89sk610WBfZ+jg6MdWgrnm8/jyMUjCPELwfiB4/njgiBwa0Fa1GaP8kvl\naNI1oV9wP0QFRaGxvZEvRtoN7ThecxwCBCwZtwSA2Vq4UjlRewJPbH4Ceyv3euT8fUYUlCqaAfvu\nI8tcaEtLITIoEmEBl0XhssuIWQnhAeEAzCsrT5iyzHU0LXUaBqkHATDHFbqCoqVg0eqirKEMIkQM\nVA9EQngCAGVLwZ0bDLEqZqmVwLj5qpsBADvO7FB8LfNDD48ZjqsiL7uPJBPozzU/4/n85/H4t49j\n/MfjuWXhDUjjVc64xPRGPb458Q2/X1JX4vI1XYkpFFYUAgCyE7L5b5LBUlNf+f4Vp9qVMCtheMxw\nXkfDXEjHao7BKBqREp2CJ8c9iRC/EGw5tcVqC92+TpuhDcsKliHjgwy8tectzNs4zyPX8WpRkG6y\no9T7CLBfwGZTFOxYCkwUYkNiEeIXAp1R53K7YFEU8cG+DzDi/RGKbRnaDe08o2Nq6lQewLN0Hx25\ncAQfHfjIab+tM+4jtsoeEjlE5j6QXuODfR8gYkUEfjznnv0kDlYfBABcm3Ct1XM5g3MAQNb+Qgqz\nFIbHDseQyMuWgqRWgU1cgLk/z8hVI/F24dvuGLbHkVoKli6xHWU7kPmXTNn3Z8eZHWhsb+QLo65Y\nClJRsNXqgn2218Zb/92mpU7DorGLoDPqcNc/70J+qX13j9TiG9pvKIAOUWBpqBmxGYgKisIjox4B\nALz239dcem/ezLGLx3D1B1fjxR0vQmfUwU/lh4PVB3H04lG3X8srRYG3uXDQOhuw7z5yaCkE2rcU\nnAnEWdKia8GcL+dgYd5CHLl4BP858R+rY3ac2YFmXTMyNZkYqB7YYSlYuI8e/eZRPLLxEbyz5x2n\nru2MKLAgc1JkEvx8/BAdFA2TaJK9x6+Pf40mXZNDv67eqHdqZc6+2MxNJuWmQTcBAHad3WUVAwI6\nVpgjYkd0iEK9tSg8f8PzmJs1F+3Gdjz+7ePYf26/w3H1NDJLwcJ99OnBT3Go+hB+veHX/Hv89c9f\nAwAeynoIgNlScNWKlcYUHH2/WbW5kpgLgoC3b3sb80bNQ5uhDXesu8NuXyou7jHDrUSBxROu1lwN\nAHjiuifgq/LF50c+d0uqtjewbMcynKw7iWH9hmHHr3YgNzMXALDu8Dq3X8srRcHZ1tmAc+4jdj4r\nSyEokvfhsbQUwgPCER0cDcB5UThRewLZH2Xjs8Of8ccsayCAjlTUO1LvAAAMirAWBaPJyNP0fp//\ne5xtPGv32nqjHi36FqgEFUL9Q3mqrj1LAehI6ZVaWsdrj5vH48Cd9afdf0L6++l2q7FNoon/+IfF\nDLN6fkDYAKRFp6FF32K1092l9ksov1SOAJ8AJEUmdWQf1ZfyybCw0iwKtybfir/e+Vf8YvgvAJjr\nIno7UlGwrNRmk+iJ2hN4Z887EEURXx83i0JuVi76BfeDVq/FuaZzLl2TxxSC7McUDCYDj/NkJ2Qr\nnksQBHxwxwfIzcyFVq/FzC9m2nQ7cjdgrLUosMyjjNgMAObeXfePuB9G0Yg3f3jTpffnKqIo4tPi\nT3u80eKJ2hMAgE/u+gQ3DroRszJmAQD+8dM/3O6+7jOi0Bn3EWubzY5RtBQuu4+Yi0jJUnC2VuGh\nrx/CkYtHkBadhmfGPwPAugYCAHZX7AYA3Jp0KwAouo9K6kp4A78WfQsWbVpk99pSK0EQBNkPXuoa\nYqKQFJUEwHqjojZDGxeDssaO8SjxQ8UPAID/zf9fm03OzjScQauhFXFhcYgIjFA8hlkLli4kZiUM\n7TcUPiofqAPViAqKQquhFdUt1Whqb8JPF36Cr8oXowaM4scCHRaRJ7nYcrFLNR62As2iKMrcBi/t\neAmbSjah/FI54sLiMCZuDFKiUgC47kKSuo/sZR8duXAELfoWXBVxFbc6lVAJKvx1+l8RGRiJ883n\nFbPdTKKJvx9nLAUAeOq6pwAAHxV9ZNX23p3858R/8Kuvf4UHv3zQY9dwhLRBJVv43DToJgwIHYDT\n9afdHnDuO6JgI9DMV7rN1VaKyl4fHRQNH8EHLfoWtBvaZZaCu9xHJtHEV/a75u5Cdny27HxS2AQ+\nIGwAH1+wXzAutV/izzE//DVx1yDMPwwbjm/AVz9/ZXUuy3OyiddX5ctdQ1JRY5OlLUuhpK6EN59z\nZCmwc9W21trMFmGTwbB+1lYCw1ZcQRpPYEhdSPvP7YdJNCFTk4lgv2AAZreYdGyeQhRFjPt4HEat\nGsWtWFexDDSz72/5pXI065oRGxKLO1LvQJOuCff/634AwJ1pd0IlqJASfVkUXEhLNZgMaGhrgEpQ\nISIwwu73257ryBIflQ9f2Ch9Z840nIFWr0X/0P6IDo7GVRFXwU/lh7ONZ1HeWI6q5iqE+IXwRAIA\nyOyficlJk6HVa/HBvg+cfo+uYDQZ8dz25wCYP39pF+XupL6tHpfaLyHEL4T/TXxUPrh/hPlvvu4n\n97qQ7IpCW1sbsrOzkZWVhfT0dCxdulTxuMWLFyMlJQWZmZkoKjLnJR8/fhwjR47k/9RqNd55x+z7\n3rt3L8aOHYuRI0fimmuuwb59rhXZONvmgt1XB6ihN+ll/Yykrw/wDUBUUBQA8x/AmUBzuH84ooOc\ndx9VNVWhzdCGmOAY9Avux7OXlETBMsNJEASrDCTm/picNBnLJ5r3UX4s7zFFdxQgT0dlWMYVpCsS\nLgoWlsLxmuP89WUNZTZNV5Nokvn23yp8C5WXKq2OO1ZjjjkoxRMYNw02Wwr/Lf+vLK7A4wkxI/hj\nUlFg8QT/2TBiAAAgAElEQVRWbwF0WECebrhW2VSJkroSVDVXdSrjSRRFbimE+odCq9fyvxN73+kx\n6fjz5D/DT+XHFy13pt0JADYthUPVhxT/DkCHpRwVFAUflY/dQDMPMjshCkCHtauUVi2NJwDmjgXJ\nUckQIeLfx/4NwBwzYh0GGMxaeHfvu50WXnus+2mdrMq6M9lc7oBZiUMih0AQBP74rBFmF9L6n9a7\ntfWHXVEIDAxEfn4+iouLcejQIeTn52PXLrlvLS8vDyUlJTh58iRWr16NBQvMJflpaWkoKipCUVER\nfvzxRwQHB+Puu+8GACxZsgQvv/wyioqK8NJLL2HJkiUuDdpX5QvAvLJhk5It9xFgO64gjSlwUWit\nl6ek+ju2FJxpn80DuJcnJUsLRIqlKADWLiRmKWRqMvHomEcxNn4sKpsq8VbhW4rXt7QUAGtRuNBy\nAVq9FlFBUfw4y42KWDyBjd1y4x7GuaZzXATvTb8XrYZWvFDwgtVx9oLMjLiwOKRGp6JZ18ytLUDZ\nUmDm9en60zyeIJ24ustSkLp3OtP2vFnXDL1Jj2C/YJ6iyURWOommRKfgt9f9FoD5+8LahLCCR6ko\nlDWUYczqMcj5NEdxEpHGEwBAHaiGj+CDS+2XrJIyXBUFe2nVbOIdEdsh7qwQ8YujXwDoiCdIuWXI\nLcjqn4Xqlmr8/dDfnRqHs+iMOvw+//cAwOOKnSkGdAcsHVlqKQHAmLgxSI5KRnVLtVuL+Ry6j4KD\nzWa3TqeD0WhEVFSU7PkNGzYgN9ccCc/OzkZDQwOqq+X++23btiEpKQmJiebCpAEDBqCx0exjbmho\nQHx8vEuDFgSBZyCxtFRb7iPAdv8jJVGoa62TWQr2As2uuI/YKoNNSrYshXZDO085C/AJ4I9bZiCx\nHO3M/pnwUfngyeueBACrYCxDSRQsg82WriPpMUwUWMCLYSv7g63Ek6KS8OrNr8JX5Ys1xWusUuic\ncR8BQM6gHAByF5LlClM6dqmlIJ24WCpxQ1uDR33R0upcJuCuwFxH0UHRHfUXlycHy/f9uxt+h/uG\n34flE5dzK5pbCpKJ7NuSb6E36VFSV6KY9SaNJwDmeABLppAufBraGnCs5hgCfAKQ1T/LqfejlCzB\nUPo7Do02xxV2l5vja9J4AkMQBG4t/Gn3n1yunVHKZmOs/nE1yhrKkB6Tjl+P/DUA6+9+d8EsBbbg\nYQiCgNkjZgNwbxaSQ1EwmUzIysqCRqPBhAkTkJ4uX9FVVlbyyR4AEhISUFFRITtm/fr1mD17Nr+/\nYsUKPPnkkxg4cCCefvppLF++3OWBS11Ioig6ZSlYVe9ermi2EgWF4jV7geaaVseiwCdJC1Gw1YE1\nPCBcZiryH1XDGdRqa1HZVIkQvxB+Pva/rSInRUshWG4p8CDz5XMB1p8dsxTUAWoAdkRBktqaEp2C\neaPmwSSasKxgGT9GFEWn3EeAJK5wpgCA2X12rukcgnyDZCsoJgr5Zfm40HIB/YL7yUROEASnXEha\nvVbmKnOVrloKbBKODo62qtTmhV6XLaSwgDD8895/YuE1C/nrWUzhVP0pPllK9zh4d++7Vte0FAUA\nisFmlnU0csBIq6I1W9hKq1Z6P0BHQgAjQ2NtKQDAfcPvQ2J4Io7XHnfJWqi4VIHBbw/G+I/HW80L\nzbpmvLzzZQDAKze/wsdiL2ivN+qxKG8RVu5d6fQYnMXSpSuFZSH9+9i/3eZCcygKKpUKxcXFqKio\nwM6dO1FQUGB1jKVfWTqZ6XQ6bNy4ETNnzuSP/frXv8Y777yDs2fP4s0338TcuXNtXn/ZsmX8n/Ta\nUlEwmAwwiSb4CD7ctSSlf4h991GAT0dM4aL2IhrbzVaMOkBt7T7SSVJSXYgpWLmPLp/X0lJQch0B\nEvdRYxlfeWZoMuCj8gHQYVqerj+t6Od3xn3EJkmZpSBxH4miyCfKiUMmArDdeoOdi7kxlt5gjkd9\nc/Ib/uU913QOl9ovITooGjEhMYrnYbC4wq6zu/Cvo//C9WuuB2Au7pP6mtnYWYrutQnXyr6P0jHZ\ncyH9ZtNvMPS9oTz331WO1nSIwsHqg52uF5BZCvWlVpk6tggPCEdsSCzaDG2ovFQJk2jC9tLtAMx1\nPt+VfmcV62BCJBUFJWvYXtGaLaSLGilGk5EvDGSWgqUoKLiPAHP84fc3mt088zbOczp19L297+Fc\n0znsLt+NcR+P4xbV4erDuGXtLbjQcgHZ8dm4M+1OpEanArBvKfxh5x+wct9KPLPtGbdW+wMS95GF\npQCYP6fkxmSE7A7BM797BsuWLevy9ZzOPlKr1Zg6dSr275cX/cTHx6O8vJzfr6iokLmDNm3ahNGj\nRyMmpuNHv3fvXh5fuPfee7F3r+2UKqko5OTk8MelomDPSgBccx8xU00doIaPysdm9lFYQFjnYgqX\nV+Eh/iEQIKBF3yLz7zJBshQFqU+Wu440mfz5yMBIhAeEo1nXrJgiK22bzZD2PzKJJnx/9nsANtxH\nzdWo0dagvq0eYf5hfEKwZSmU1MvdZQnhCcjqnwWtXsub3zkTT2DEhcUhJSoFzbpmzPxiJpp1zZid\nMRtr71orOy4xPFEmEkoTFxuTvcDh3nPm7+RTW59yeWtIacqov48/arQ1LqemSnsQSSu1zzaeRYu+\nBZoQDXft2EIabC4+X4y61joMUg/C3JHmRZjlqtYypgBAMdisFKtxhC1L4XT9abQZ2pAQngB1oJo/\nzmIKgPlvb++9PjzqYSwcsxDtxnZMXzfdoYXXqm/Fhwc+BGBeIJyuP41xH4/DorxFGLV6FPZU7kFc\nWBw+mPoBBEFwmN5bWFGIV75/BYDZwmTdbd2FrZgCo3hFMSo3VOLtP77teVGoqalBQ4N5hdna2oqt\nW7di5MiRsmOmT5+OtWvNP8zCwkJERERAo9Hw59etW4dZs2bJXpOcnIwdO8y9bLZv347U1FSXB85b\nXRj1duMJgO1aBSVRYJN3ZJA5S8dm9pGLMQWpjx0ALyID5Nt9svNLfyCA3CfLLAWpP1cQBFnxliV2\nLQXtBby842V8V/odwgPCcVvybdbHtFzgeeNp/dLsZpMovV+go70ya5LmrOuIwVxIAgSsmLgCf7/7\n71YLAT8fPwxUD+T3lSYuR8FmURS52JXUlbi8FWpVcxUa2hoQGRjJU49djStw91FQhyiU1pfKMo8c\nIU1L3Xba7Dq6ZcgteGzsYwDMVdHS+hEl95Hld1wUReypcD4dlREVFIUQvxBZWjWgHE8AzN9T9rtV\niidIEQQBb9/+NqalTkN9Wz1u/+x2u00c1/+0HrWttRg9YDSK5hfhtuTbUKOtwcp9K2E0GbFwzEIc\nXXgUIweY5zq2fe2FlgtW9TbNumbM+XIOjKIRAswWqaNW59+f+d5pl6JJNPHvIvvNWRLiH+LUuZzF\nrihUVVXh5ptvRlZWFrKzszFt2jRMnDgRq1atwqpVqwAAU6ZMwZAhQ5CcnIz58+fj/fc7fkAtLS3Y\ntm0bZsyYITvv6tWrsWTJEmRlZeH555/H6tWrXR4422jHGUtBmm4qxa4oXE7d5IFmhewjaUWzPfdA\nfas5zTXEL4S7Y9g5pOe0PL+U/qH9+aqTBd+klgIAq4CkFHuB5h1lO7BsxzIIELD+nvW8ER77bCID\nI2EUjdxtkBqdykXKmZgCg7VXZqLgbJCZ8djYxzA5aTI2zNqAZ65/xsotxGCTqACBd3qV4iimwPLC\nGS/teMmloDR378QO538jV+MKUvfRQPVACBBQfqmcW4n2XEcM6QpXKgojYkcgZ3AOWvQt+PTgp/x4\nFhuzJwqn6k+htrUWmhCNTHwdIQiCogtJqTCNwVxItlxHUnxVvlh3zzqMiRuD0oZSPLzxYcXjRFHE\nO3vNqfGLxi5CqH8oNty/AYvHLsYNA2/Arrm78N7U92SLMpWgUszmAoAnv30SJXUlyIjNwIOZ5gI3\nW5bKxZaL+MW/foEbP7kR1350rVMprueazkFn1CE2JJbPRZ7G2gEvISMjAwcOHLB6fP78+bL7K1cq\nB1dCQkJQU2O9ih4zZgz27NnjyjitkLqPBJN5crCsUWCwx1kFM4Pdl4lCnYWlEGDbUgj0DUSIXwha\n9C1o0jVZTeQMaTxBOpGFB4SjsqnSKVFQCSoMVA9ESV0JSupKIECwCr45Yyko1Skwl9Vrk17D7Sm3\nW71WE6pBfVs9dp41u33SotPsFiPVtdahoa0Bof6hsmrX7IRsRARG4GTdSZysPemS+wgwTxzf/tJx\nx80hEUOwHdsxPHa44t/EkaXAhG5E7AhoQjT4rvQ7vLzjZbx5m3MtFfj76peOzP72RUEURTz41YO4\n0HIBebPzeIxIGmj29/FHojoRZxvPYvMp8y500qCsLdhE9tOFn7ivnXWdXTR2EQrKCvDevvfw2NjH\noBJUTsUUmOtPKVbjiEHqQTh68SjONJ7hn4tlCwspOYNyUFBWgFuG3OLU+UP8Q7Dh/g0Y9NYg5J3M\nw/nm89zaYPy3/L8oPl+MmOAY/GKEueWJn48f3r7dfpPE1OhUHL5wGCdqT2BM3BgAwPbS7Vh9YDX8\nffzx2YzPuPAqWQpfHPkCC/MW8s+x1dCKhzc8jO25263qL6TYyjzyJF5Z0QxYxBQcuI+YKFhG55Us\nBbZCY5OnvTYXAJyKK1hmHjGUahWkxXGWSM3H5Khkq5UDdzM4aSlIJ+w5V8/haa2WMOuGTSxp0WmI\nCY5BkG+Q1aoakKffSicOX5Uvb92xqWSTy+4jZ2FuE2nRmpREdSJ8Vb4413ROcctK6Q/xjclvQICA\n9/a953SeulTs2ArYlvvoq5+/wt8P/R1bTm2RFftJLQU2FqAjRdMVS2Hb6W1oNbTias3V/G8+PW06\nEsMTcaL2BA9AO5N99EnxJwCAqSlTHV7fEqVaBXuWwtIbluLEYycwOWmy09cYEDYAt6fcDpNowudH\nPrd6nmVdzRs9z+YiUgmlFF+WBrpk3BJkaDJ4HMRyK9Etp7bgvn/dhxptDW6+6mbse2QfYoJjsOPM\nDqz+0b6XxFE8wRP0DVFw4D5yJArS7CMGFwUH+yk4E1dQcqVIz+GMpQB0/KgA8JWWFHvbUSqJQph/\nGGaNmIW7ht6F1dNW21z5MTcTO0dqdKpNdwCgHE9gMBfS2oNrUaOtQZh/GN/21F3MHz0fz45/Fi/c\nZF0sB5jFiQms0i5tUh9uZv9MPJT1EPQmPd4sdM5SYH7y9Jh0jIgdAQECfq752cpS1Rv1ePa7Z62u\nC0hE4bKLkgk+y2xxxVIwiuZEhluu6lhx+6p88fAos4vlowMfAZAHtxnS7/exi8fw/dnvEeofylss\nuIJlrYJWr8XJ2pPwEXysso0A82+cCbwrsNx9aeNJwJyG+u+j/4aP4INHxzzq0jl5BlKdOQNJFEVs\nOW1ub3/3MHPSDHsPlu4j5i5dOGYhts3ZhjFxY/DelPcAAEu2LrHbzNKy51F30CdEwVaHVIYzloLU\nrQIoB5rZtXwEH34tZzql8pWzxSRpGcQGXBAFjYIoSNJSLVESBUEQ8I97/oEvf/Gl3VWTNA4CdPxA\n2Hgs4wq2RBAAD2KzIrv0mHSX3RCOUAeqsfyW5YgPt10Uac+FZJkCyH70zrRplm4rOTx2OIL9gpES\nnQKDycAtI8Zfi/4qS3OUirk00CwdC2COL1kuYpQICwiTuU8s3TAPZT0EAQK+/PlL1Ghr7AaaL2ov\n8oyd2SNm88WSK1hmIB29eBQiRAztNxQBvgH2XuoS09KmIdQ/FHsr98r89m/sfgNG0YgZw2bI4mbO\nYNlL6kTtCZxtPIt+wf14wscg9SAE+ASgsqlS9ptme2XfOfRO/l2/N/1e3D30bjTpmvDIxkdwpuGM\nYlySfSeUahQ8hdeKgrSi2dZeCgxX3EcMJhJ+PubKYqNo5Pn80sKy7rQUpO4jpUpSqZ/fso2BUkqq\ns0hFISE8gWc72MpAsqzelhIbEst9soD7XUfOwmsVFILNltkejjabkXKh5QLq2+qhDlBjQKi5oaFS\nsLlZ18wL+djfUslSYNeWTgqufGbM7eGn8sMNg26QPZeoTsRtybdBZ9RhTdEaNLY38mZ4DHb9c03n\neFB63ujO7fhlaVnacx11hWC/YNw91CzkzMXz04Wf8O7ed6ESVFh6vXIPN3tIaxVEUcTW01sBAJOG\nTOIxAR+VDxcPJvZ6o563ZpF+7wVBwHtT3kNkYCS2nNqCwW8PRviKcFz70bX47vR3/DiKKbiAO9xH\n0opmy8mSWQpARwYSayQmnbD7BTlun23LneKy+yjCvqUQ7BcMTYgGepNe1ke/zdCGNkMbfFW+vFuo\nKzD3EQDehwdwbCmwidcSlpoKOJ955G7sWQrs/TDLi/nVL7Y4FgVpPIEtHHhcQbKHwxu730B1SzXG\nxo/F49mPA7BhKVy2RKU+ZWfiCQwmCtclXqeYvcJcSG/tMffMig6KlgU+mShcaLmAutY6jBowCqPj\nRjt9fSmWloJ0RzV3Mzujw4UkiiL+J+9/YBSNeHT0ozzV1BVigmMQHhCOxvZG1Ghr+M6IlvEOy7bf\nRy4eQZuhDclRyVYLzwFhA/DV/V/h5qtuRmxILJp1zdhTuQdPbX2KH8PdRxRTcIy7A80+Kh+ZMEjd\nScxUZhOtTBQcWAqt+lZUNlXCV+VrlcJnWS0NdGQCWdYpAOZJVoCA2JBYm+avUloqy61meym4itRS\nYCsmwLalYC+mAHTEFYCesxTY2CzTAkVR5J8de3+s2toZS0Ep756JwqEL5pXxidoTeH336wCA1255\njVsBTIx0Rh2adE3wEXx4OxGppeCKKFw/0Fz5fc+wexSfvyP1DsSGxPLvtmWRWLBfsMwCnzeq8/sC\nDwgbAD+VHy60XECrvpV/Hu62FACzqywmOAbHa4/j6a1PY+eZnegX3A9/uPkPnTqfIAj8u3/k4hHe\ngG7SkEmy4yz3l2Z7HVwTZ50aDQA3DroR3z34Haqfqsb5J88jKigKxeeLUXy+GO2GdpxrOgcfwQeJ\n4YmKr/cEfUIUuhpTYP5MqZJLLQU2eVc2WVsKjmIKbIIZpB5k1YLDVUshLiwOn8/8HP+a+S+bk7tS\nWqpSOqor2LIUlDb/adG1oKq5Cn4qP5tf5DFxY9A/tD9UgsojE4Iz2LIUarQ10Oq1iAiM4IuEMP8w\n+Pv4Q6vXOuypr5RmK3Uf7a3ci/Efj0eLvgV3Db0LNw2+yepzlLawZn9nTYiGf7+dCTIzcrNycfDR\ng1g0VnkTJn8ff761IyCPJwDmyZBZSiF+IbzXTmdQCSokqs3fibONZ7n7yFZfo67gq/Llu+y98cMb\nAMwCLP1duwqzutYeXItmXTPSY9Kt4lY82Hy5RxjrEzU2fqzD82tCNTxI/knxJzjTeAYiRCSqE3ld\nVnfQJ0SBuY+6ElMALERBwVJQdB852H3NsgeQFHZeZ0UBMAeoLH3DUpQykJSCzK4gtRSk7QeUso+Y\nuTs4YjDPubfER+WDb2Z/g42zNvJJoruRrs4NJgN/3NJKAGC1U509lERhoHog1AFqXGi5gJs+uQk1\n2hrcnnw7/nb33wCYxd5P5YfzzefRqm+1ch2xMdw19C6kRKVgZH/n3R9MeO1ZiKwLKGAtCtLHZo2Y\nZfN76SzMhbSncg9qtDVQB6g9tgpmLiTAXFeRm5Vr52jHMEuBbWozeYh1qqyl+4i1S7FlKVjyq6xf\nATC7vVgWU3fGE4C+Igp6+zEF1oK63dgui/DbFQUnLQVHk4W9TBzeKVWpTqGTPz6lzeu7LAqhyu4j\nVmV9UXsRLboWAI7jCYxRA0bJ3EjdTZBfEOLC4mAwGVDe2NG7i8cTLH6IzsYVlPZ4EASBW0RthjY8\nlPUQvr7/a+7j91H5cNfimcYzVjUKjH/c8w+cWHTC7W0N0vql4YaB5oUGi5FJGRs/FoG+gbw9Rldg\nC4mNJ8z7kGdoMtyefca4NuFaDOs3DL4qX7w35T27RWLOwCwFtrhUqp9gv4+TdSfR1N6EIxeOwEfw\ncTqOMWrAKGTEZqBGW8NrKkgUnIRnHxn1HYFmG+4jQRBkwsCQBpoB25YCDzQruY8cdEq151931X3k\nDEoxha6KQqBvIK5NuBbpMemytFhWZQ10dCW1VajXG1FyIdnqM+NMXOFiy0VeexEfJncrMN/z8zc8\nj79O/6uVO4Bdr7S+VNFS8DTPjH8GvipfRSt05ZSVqHiiQrE2xlXY9+fbEnNl+tWxnnMfCoKA7bnb\n8dOCn/ge3V1BuiDy9/HHjYNutDomPCAccWFxaDO04aufv4JRNGJE7AinEzwEQeDWAstw6s50VMCL\nRUEppmAv117JhWRlKQR2iIKsyMsJ91FnLAVX6xScwV5MobOiAAC7HtqFg48etHIJWfrDLVuE92aU\neiCxz81KFOxsYs9g+ehKtRfP3/g8Ljx1AS/f/LLiypj93coaymxaCp5kaupUtDzXwvv3SPFV+bpN\noJgoMOvY0zGl/qH9ZS7PriAtpBufON6mxcbibn87ZHYPOus6Yvzy6l/K4o/dmXkE9BFRcOQ+AuyL\nArMimKUQHhAum/zY5K2UfSTdmUqp+MRW4Zr0PEwIlIrjXCVRnQgfwQfnms7xCtqu1CgwfFTKe1VY\nphnaq1HobSRHKmxZ2VgGwNpk57UKNtxHWr0WT3z7BABlt4IgCHb3jOCWQkOpVeFad+HshjldQZpW\nDXgmyOwpIgIj+OLAXusNFlf4rtRcb6DUlNEesSGxMtcquY+cRLFOwc5E6pSlcFkULLN0LFNHpaIQ\n6BuIUP9Q6E16nk7KMJqMfAWtZAJaigKzGNSB6k77WX1VvkhUJ0KEyCdqd1gKtmCT2a6zu/Dyjpex\np9Lc6NBRTKE3wPy8zEwHHFsKttxHS7YuwYnaE0iPScdzNzzn8ljYalBmKXSj+6i7kLofAfm+zN7A\n+IHj4afyw/S06TaPYaLAWpI4k3lkyUNZD/Hb5D5yEleK1wAXRcEibc2ypN/StcN8jUu2LpFZC8Xn\ni6E36TEgdICiT9FRX6XOYulC8qQosB/5Z4c/w/8W/C8utV9CUmSSV7iPJl41EeEB4ThUfYhXqjIh\ntRlTULAUvi35Fu/tew9+Kj/8/e6/u9RojSG1FJQ2u+krJKoT+b4DgyMGd/m73t18cucnOLLwiN36\nGmnadqBvoEt1JYwpKVN47yxp48ruwGtFgQXq9CZ9p2MK0tbZgG1LwbIS1PKL/Jepf0GgbyA+PPAh\n383qpws/Yeo/zJ0kbbX+tbQUbO265iqWaaldrVOwx/iB43lDwYeyHsJ/Zv0HRxYe6RZXRFcJ8A3A\nnWl3AgD+dfRfON98Hm2GNkQHRVstBHhMwWI/7rrWOjz0tXlV92LOi52qlgVsxBT6oKXg7+OPAWHm\n9h89VaPSFdSBaodN+qTN/Ub2H9mpGgN/H38cmHcAxfOLPZadZQu7+yn0ZlypaAacsxRyBudgWuo0\nHv1nMPcRw3LSvib+Gnw8/WPM/r/ZePzbx2EUjfjDzj+gtrUWtwy5BX+54y+KYwryDYJKUKHN0Aa9\nUe8+S0HSGG/LqS34+rh5n2FP1AQMiRyChmcb4CP4dGuBjbuYmT4Tfzv0N3xx9AtMGDwBgHJgz1ZM\n4eUdL6OquQrjEsdhyfglnR6HJlSDAJ8A1GhreN1HX7QUALN1ea7pnEfaW/QGEtWJCPINQquhtVOu\nI4Y7mwS6gtdaCu50H7EPPywgDBtmbcCMYfKd4hy5jwBgVsYs/O6G38EkmvDEt0+gtrUWU1KmYOOs\njTbT0QRBkNUquEsUmA/y6+NfY9q6aWgztOHhkQ9jfOL4Lp3XFoG+gV4pCAAwKWkSwvzDUHy+mG+S\norTtoa2U1APnzc3OXrjpBZvFes6gElQ8CMtqHfqipQCAF9+xFhx9DZWg4i5lVzOPegN9QhQctbkA\nnLMUbOHIUmC8NOEl3p3xzrQ78X/3/Z9D/7LUheTumMKJ2hPQGXVYOGYhVk1b1e1mqDcQ6BvIg4b2\nioVsFa+xRAJ3ZFux67IK675qKbw68VXs/NVO2V7gfY2nxz2NaanTMC1tWk8PxWX6lPuoq3UKtnDG\nUgDMK4TPZ36OA1UHMHrAaKdWjtJaBXu7rrmCNFvh8ezH8edb/0yCYIeZ6TPx2eHPuBWgZClEBUVB\ngID6tnoYTAb4qnxhMBlQeakSAgS3uOYsr9tXLQV1oNpuq5a+wANXP4AHrn6gp4fRKfqGKHTSfWRZ\n0WwLR4FmKb4qX5f8iEqWglKHVFfQhGrw+qTXEegbiP+55n9IEBwwOWkyQv1D+ZarSpaCj8oHUUFR\nqG2tRa22FppQDSouVcAoGhEfFu+WwLr0uqH+oV4RrCf6Hl7rPlLaZKe73EdKfek7i7QpnrvcRwDw\n1Lin8NjYx0gQnCDILwjTUjvMfCVLAbCOK7CAsGVBVmeRXrevuo6I3o/XioJiTKGLFc22kLqPwvzD\nutxYS4onAs2E69ybfi+/bVMULOIKtvokdRZp1lNfdR0RvZ8+5T7yWExBYim4e8Jm8QN3WwqEa9ye\nfDuSo5LRP7S/zcUFsxRYcRkTBcsq3c5ClgLRG+gbotAJ95EoilwUHKVTsn2a243tbp+wpe4jdxWv\nEa4T5BeEIwuPKPZ3YrC20tx9ZKP6ubPEBMcg2C8YWr2WLAWix/B691GrvhV6kx4CBLsrfktR0Jv0\nAMyBYWfcQSyO4HZLgbmP2sl91NP4+/jb/S5Ytrpwt/tIEAR+LrIUiJ7Ca0WBre7ZRBrkF2Q3qGop\nCs66jhhsRe92S8HfM4Fmwv1YNsVjloK73EcASBSIHsdrRYFN5mwidVQkZksUHAWZGWzy9pilQIHm\nXo80pmA0GfnGQmyjIXeQ3s/caM1dGU0E4SpeH1NgfnhH+w/0VktBsU4hoGt1CoRn4P2PtBdxrukc\nDCYDNCEau1lvrrL0hqXI6p+Fe9Lvcds5CcIV7FoKbW1tyM7ORlZWFtLT07F06VLF4xYvXoyUlBRk\nZmxvCmUAABEfSURBVGaiqKgIAHD8+HGMHDmS/1Or1XjnnXf4a959910MGzYMI0aMwDPPPOPywC0t\nBUc/zC6LgocsBU/VKRDuR5qS6u4gMyMqKAoPXP1Ap9pvE4Q7sGspBAYGIj8/H8HBwTAYDLj++uux\na9cuXH99RyOrvLw8lJSU4OTJk9izZw8WLFiAwsJCpKWlcYEwmUyIj4/H3Xeb+wLl5+djw4YNOHTo\nEPz8/HDxov3N0JVgk7kzbbOlz7PjLdtmO8LTgeb6tnpo9VqoBJXT+7kS3YvUfcTTUcnNQ/QxHMYU\ngoPNE5ROp4PRaERUVJTs+Q0bNiA3NxcAkJ2djYaGBlRXV8uO2bZtG5KSkpCYaO4P88EHH2Dp0qXw\n8zMHi2NibG9TaAvLydxb3UfMApHu/0xVyL0T6X7cPPNIPbjnBkQQHsChKJhMJmRlZUGj0WDChAlI\nT5fvOFRZWcknewBISEhARUWF7Jj169dj9uzZ/P7Jkyexc+dOXHvttcjJycH+/ftdHjhrc8HwtPvo\nmrhroBJUvO2vu2Aio7T/M9G7kG69erD6IACyFIi+h8NAs0qlQnFxMRobG3HrrbeioKAAOTk5smMs\nN6yXrnR1Oh02btyIP/7xj/wxg8GA+vp6FBYWYt++fbjvvvtw+vRpxesvW7aM387JyeHXdpel4OxG\nFguvWYg5V8+x6pjaVZgIGEWj7D7RO4kJjkGzrhn7z5kXMu6OKRCEqxQUFKCgoMBt53M6+0itVmPq\n1KnYv3+/TBTi4+NRXl7O71dUVCA+Pp7f37RpE0aPHi1zESUkJGDGDPNGNtdccw1UKhVqa2sRHW2d\nmy0VBSmWotDZlFRXOlG6WxCUzkmi0LuJCYlBaUOp2wvXCKKzSBfLAPDiiy926Xx23Uc1NTVoaDDv\n79va2oqtW7di5Ei5+2T69OlYu3YtAKCwsBARERHQaDT8+XXr1mHWrFmy19x1113Yvn07AODEiRPQ\n6XSKgmAPK0vBRfeRs22zPU2AT4DMFUbpqL0bloHEcGfhGkH0BuxaClVVVcjNzYXJZILJZMKcOXMw\nceJErFq1CgAwf/58TJkyBXl5eUhOTkZISAjWrFnDX9/S0oJt27bhww8/lJ137ty5mDt3LjIyMuDv\n789FxRW6O9DsKQRBQFhAGOpa6wCQpdDbYcFmdjvEP6QHR0MQ7seuKGRkZODAgQNWj8+fP192f+XK\nlYqvDwkJQU1NjdXjfn5++Nvf/ubKOK2wbFzmraIAmIWARME7kFoKZCUQfRGvbXMhCILM7eLpNhee\nRCoEJAq9G1arAFA8geibeK0oAPJVvqdTUj2JJ/drINwLWQpEX6fviIKL7iNXK5o9CVkK3oM0pkCW\nAtEX6Tui4MBSYG6idmO7bIOd3iAK0rRUEoXeDbmPiL6OV4uCdMc0RzEFQRBkwtCbRIFtyQlQSmpv\nR+Y+ompmog/i1aLgivsIkLuQKNBMdAappUAxBaIv4rX7KQCuuY8Asyg0tjfKRKE3WArkPvIewvzD\nMDVlKvx8/KAOJKuO6Hv0HVFw0VLoLRXNAFkK3oQgCPjP7P/09DAIwmP0GfeRM5uSKLmPeoMoUEoq\nQRC9hT4jCs66j4DeJwpkKRAE0VvwalGQVjR3OtDsZOtsT8KEQIBAvXQIguhRvFoU+oqlwALNYQFh\nUAle/SchCMLL8eoZqCsxhd4YaKYaBYIgepo+IwqddR/1BlFIikxCYngibhp8U08PhSCIK5y+k5Lq\nxe6jEP8QnP7NafgIPj09FIIgrnC8WhRcaXMhPaa3VTQD1vtDEARB9ATe7T5S9Y06BYIgiN6Cd4vC\n5Qk9wCfAqawdWaC5F7XOJgiC6C30CVFwJp4AkKVAEAThiD4hCs64jqTHkSgQBEEo0ydEwZl0VKD3\nVjQTBEH0FrxaFFj2EbmPCIIg3INXi0JXLIXeVNFMEATRW+gTokAxBYIgCPfQJ0SB3EcEQRDuoW+I\nQhcCzSQKBEEQHXi1KLD9FLpiKfSWNhcEQRC9Aa8WhYTwBABAYniiU8czUWg1tFJFM0EQhAJ2RaGt\nrQ3Z2dnIyspCeno6li5dqnjc4sWLkZKSgszMTBQVFQEAjh8/jpEjR/J/arUa77zzjux1b7zxBlQq\nFerq6jo1+ElJk7B77m68NOElp45notCia4EIESpBBR8VdSYlCIJg2G3NGRgYiPz8fAQHB8NgMOD6\n66/Hrl27cP311/Nj8vLyUFJSgpMnT2LPnj1YsGABCgsLkZaWxgXCZDIhPj4ed999N39deXk5tm7d\nikGDBnV68CpBhesSr3P6eCYKl9ovASArgSAIwhKH7qPg4GAAgE6ng9FoRFRUlOz5DRs2IDc3FwCQ\nnZ2NhoYGVFdXy47Ztm0bkpKSkJjY4eb57W9/i9dee63Lb8AVSBQIgiDs41AUTCYTsrKyoNFoMGHC\nBKSnp8uer6yslE32CQkJqKiokB2zfv16zJ49m9//+uuvkZCQgKuvvrqr43cJS1GgIDNBEIQch6Kg\nUqlQXFyMiooK7Ny5EwUFBVbHiKIouy8IAr+t0+mwceNGzJw5EwCg1Wrx6quv4sUXX7T5ek8hDTQD\nZCkQBEFY4vR2X2q1GlOnTsX+/fuRk5PDH4+Pj0d5eTm/X1FRgfj4eH5/06ZNGD16NGJiYgAAp06d\nQllZGTIzM/nxo0ePxt69exEbG2t13WXLlvHbOTk5smu7imXlM4kCQRDeTkFBgeJivbMIop1lek1N\nDXx9fREREYHW1lbceuuteOGFFzBx4kR+TF5eHlauXIm8vDwUFhbi8ccfR2FhIX/+/vvvx+23387j\nDpZcddVV+PHHH61iFYDZ4nCnFdFmaEPQKx01DWnRafj5sZ/ddn6CIIiepqvzpl1LoaqqCrm5uTCZ\nTDCZTJgzZw4mTpyIVatWAQDmz5+PKVOmIC8vD8nJyQgJCcGaNWv461taWrBt2zZ8+OGHdt9Ad2EZ\nQyBLgSAIQo5dS6GncbelAACBfwjkHVLHxI3Bvkf2ufX8BEEQPUlX502vrmjuDNK4AlkKBEEQckgU\nCIIgCA6JAkEQBMEhUSAIgiA4V7QoUEUzQRCEnCtaFMhSIAiCkEOiQBAEQXBIFAiCIAgOiQJBEATB\nuaJFgQLNBEEQcq5oUSBLgSAIQg6JAkEQBMEhUSAIgiA4JAoEQRAE54oWhQBfCjQTBEFIuaJFgSwF\ngiAIOSQKBEEQBIdEgSAIguCQKBAEQRAcEgWCIAiCc0WLArW5IAiCkHNFiwJZCgRBEHJIFAiCIAgO\niQJBEATBIVEgCIIgOFe0KFCbC4IgCDlXtCiQpUAQBCGHRIEgCILgOBSFtrY2ZGdnIysrC+np6Vi6\ndKnicYsXL0ZKSgoyMzNRVFQEADh+/DhGjhzJ/6nVarzzzjsAgKeffhrDhg1DZmYmZsyYgcbGRje+\nLduQKBAEQdjGoSgEBgYiPz8fxcXFOHToEPLz87Fr1y7ZMXl5eSgpKcHJkyexevVqLFiwAACQlpaG\noqIiFBUV4ccff0RwcDDuvvtuAMDkyZNx5MgRHDx4EKmpqVi+fLkH3p7C+/FSUSgoKOjpIfQa6LMw\nQ59DB/RZuA+n3EfBwcEAAJ1OB6PRiKioKNnzGzZsQG5uLgAgOzsbDQ0NqK6ulh2zbds2JCUlITEx\nEQAwadIkqFQq/pqKioquvRMnkVYxe1NFM33pO6DPwgx9Dh3QZ+E+nBIFk8mErKwsaDQaTJgwAenp\n6bLnKysr+WQPAAkJCVaT/Pr16zF79mzF83/88ceYMmWKq2PvFIIgICkyCZoQjcxqIAiCIJwUBZVK\nheLiYlRUVGDnzp2KqiyKouy+IAj8tk6nw8aNGzFz5kyr173yyivw9/e3KRieYO8je3FowSH4qHy6\n7ZoEQRBegegiL730kvj666/LHps/f764bt06fj8tLU08f/48v//VV1+Jt956q9W51qxZI44bN05s\nbW1VvFZSUpIIgP7RP/pH/+ifk/+SkpJcndZl+MIBNTU18PX1RUREBFpbW7F161a88MILsmOmT5+O\nlStX4v7770dhYSEiIiKg0Wj48+vWrcOsWbNkr9m8eTNef/117NixA4GBym6ckpISR8MjCIIg3Igg\nihZ+HwsOHz6M3NxcmEwmmEwmzJkzB08//TRWrVoFAJg/fz4A4LHHHsPmzZsREhKCNWvWYNSoUQCA\nlpYWDBo0CKWlpQgLC+PnTUlJgU6n40Hr6667Du+//75H3iRBEAThHA5FgSAIgrhy6JUVzZs3b8bQ\noUORkpKCP/7xjz09nG6lvLwcEyZMwPDhwzFixAhe7FdXV4dJkyYhNTUVkydPRkNDQw+PtPswGo0Y\nOXIkpk2bBuDK/SwaGhpw7733YtiwYUhPT8eePXuu2M9i+fLlGD58ODIyMjB79my0t7dfMZ/F3Llz\nodFokJGRwR+z996XL1+OlJQUDB06FFu2bHF4/l4nCkajkbuijh49inXr1uHYsWM9Paxuw8/PD2++\n+SaOHDmCwsJCvPfeezh27BhWrFiBSZMm4cSJE5g4cSJWrFjR00PtNt5++22kp6fzjLYr9bP4zW9+\ngylTpuDYsWM4dOgQhg4dekV+FmVlZfjwww9x4MABHD58GEajEevXr79iPouHHnoImzdvlj1m670f\nPXoU//znP3H06FFs3rwZCxcuhMlksn+BLoWpPcDu3btlmUrLly8Xly9f3oMj6lnuvPNOcevWrbKM\nrqqqKjEtLa2HR9Y9lJeXixMnThS3b98u3nHHHaIoilfkZ9HQ0CBeddVVVo9fiZ9FbW2tmJqaKtbV\n1Yl6vV684447xC1btlxRn0Vpaak4YsQIft/We3/11VfFFStW8ONuvfVW8YcffrB77l5nKSgVwlVW\nVvbgiHqOsrIyFBUVITs7G9XV1TyjS6PRWFWM91WeeOKJ/2/njllah+IogB9BJxdRLNVGoRSrRNEq\n1cHR0rFBrYMOFURcnAQ/gwVx0S8gFB3qKqIOIoKgDhKcLZJCSnURHSpCFP5vkF6ez6e+Qdry7vlt\nKeFy74HkQNIbrK2tqd3vALTMwnEctLa2Ym5uDkNDQ1hYWMDT05OWWTQ3N2N5eRmdnZ1ob29HU1MT\n4vG4llmUfbb2YrEIwzDUef9yP625Uvh905vOSqUSkskkNjY23v1rC3jLSIec9vb24PP5MDg4+GFz\nZJkuWby+vsK2bSwuLsK2bTQ2Nn54PKJLFjc3N1hfX0c+n0exWESpVML29va7c3TJ4m++W/t3udRc\nKQQCAbiuq45d133XdDp4eXlBMplEKpXC+Pg4gLf2v7u7AwDc3t7C5/NVc4oVcXZ2ht3dXQSDQczM\nzOD4+BipVErLLAzDgGEYGB4eBgBMTU3Btm34/X7tsri8vMTo6ChaWlpQX1+PyclJnJ+fa5lF2WfX\nxJ/300KhgEAg8OVYNVcK0WgUuVwO+XwenudhZ2cHlmVVe1oVIyKYn5+HaZpYWlpSv1uWhUwmAwDI\nZDKqLP5n6XQaruvCcRxks1mMjY1ha2tLyyz8fj86OjpwfX0N4O0Dk729vUgkEtpl0dPTg4uLCzw/\nP0NEcHR0BNM0tcyi7LNrwrIsZLNZeJ4Hx3GQy+UwMjLy9WA//QLkJ+zv70s4HJZQKCTpdLra06mo\n09NTqaurk4GBAYlEIhKJROTg4EDu7+8lFotJV1eXxONxeXh4qPZUK+rk5EQSiYSIiLZZXF1dSTQa\nlf7+fpmYmJDHx0dts1hdXRXTNKWvr09mZ2fF8zxtspienpa2tjZpaGgQwzBkc3Pzy7WvrKxIKBSS\n7u5uOTw8/HZ8bl4jIiKl5h4fERFR9bAUiIhIYSkQEZHCUiAiIoWlQERECkuBiIgUlgIRESksBSIi\nUn4B/V0UDOPUEk4AAAAASUVORK5CYII=\n", "text": "" } ], "prompt_number": 13 } ], "metadata": {} } ] }