Skip to content

Instantly share code, notes, and snippets.

Created November 11, 2017 20:12
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save anonymous/5caa5f448c0dec3d7066efd447794232 to your computer and use it in GitHub Desktop.
Save anonymous/5caa5f448c0dec3d7066efd447794232 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Imports:"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from scipy.stats import norm"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Functions and simulation:"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"n = 1000\n",
"θ = 0.8\n",
"\n",
"d = np.sqrt(1-theta**2)\n",
"\n",
"def p(x,y):\n",
" \"Stochastic kernel for the TAR model\"\n",
" return norm().pdf((y - θ * np.abs(x))/d)/d\n",
"\n",
"Z = norm().rvs(n)\n",
"X = np.empty(n)\n",
"\n",
"for t in range(n-1):\n",
" X[t+1] = θ * np.abs(X[t]) + d * Z[t+1]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Set up grid, compute kernel"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"grid_size = 200\n",
"ys = np.linspace(-3, 3, grid_size)\n",
"kernel = np.empty(k)"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [],
"source": [
"for (i, y) in enumerate(ys):\n",
" kernel[i] = np.mean(p(X, y))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Plot:"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAAGfCAYAAACdqpz+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl81fWd7/H3NxthFYSwGfY9QIgYgoqIO1itaLW3Wm8r\nnbbWerXOdOrUmfHOfdxp50479TG203qrTvXaVltttaWOUvcqyCIE2fewBAIoAdnXLN/7xyeHExDI\nSXJyvmd5PR+P8zi/nBySD1FO3ue7fL7Oey8AAAC0XlboAgAAANIFwQoAACBOCFYAAABxQrACAACI\nE4IVAABAnBCsAAAA4oRgBQAAECcEKwAAgDghWAEAAMRJTqhv3KNHDz9w4MBQ3x4AACBmixcv3u29\nL2jqecGC1cCBA1VeXh7q2wMAAMTMOVcZy/OYCgQAAIgTghUAAECcEKwAAADiJNgaKwAAEq2mpkZV\nVVU6duxY6FKQpPLz81VYWKjc3NwW/XmCFQAgY1RVValz584aOHCgnHOhy0GS8d5rz549qqqq0qBB\ng1r0NZgKBABkjGPHjql79+6EKpyRc07du3dv1YgmwQoAkFEIVTiX1v7/QbACAACIE4IVAAAJ1KlT\npxb9uXfffVc33nhjwv9sU1r692nsmWee0Y4dO05+/LWvfU2rV69u9dfdsmWLfvOb37T66zQHwQoA\nAAR1erD6xS9+oaKiolZ/XYIVAAAZwnuvBx98UGPGjNHYsWP1wgsvnPPxxhYtWqQLL7xQGzduPOXx\nLVu2aPLkyRo/frzGjx+vefPmnfzcoUOHdNttt2nkyJG688475b2XJC1evFhTpkzRRRddpKlTp2rn\nzp2SpP/8z//UhAkTNG7cON166606cuSIJGnz5s265JJLNHbsWD388MNn/fs9++yzKisrU0lJib7x\njW+orq5OdXV1mjFjxsm/26OPPqoXX3xR5eXluvPOO1VSUqKjR4/qiiuuOHnsXadOnfTggw9q9OjR\nuuaaa7Rw4UJdccUVGjx4sF5++eVz/r0feughzZkzRyUlJXr00UdVV1enBx98UBMmTFBxcbGeeOKJ\nFv23OxfaLQAAMtI3vtE2XzfW39V/+MMftHTpUi1btky7d+/WhAkTdPnll2vevHlnfDxi3rx5uv/+\n+/WnP/1J/fv3P+Vr9uzZU2+++aby8/O1YcMG3XHHHScDypIlS7Rq1Sr17dtXkyZN0ty5czVx4sST\nX6ugoEAvvPCC/vEf/1FPP/20Pve5z+nrX/+6JOnhhx/WU089pfvvv18PPPCAvvnNb+rLX/6yHnvs\nsTP+3dasWaMXXnhBc+fOVW5uru69914999xzGj16tLZv366VK1dKkvbt26euXbvqZz/7mR555BGV\nlpZ+6msdPnxYV111lX70ox/plltu0cMPP6w333xTq1ev1l133aWbbrrprH/vH/zgB3rkkUf0yiuv\nSJKefPJJnXfeeVq0aJGOHz+uSZMm6brrrmtxa4UzIVgBABDA+++/rzvuuEPZ2dnq1auXpkyZokWL\nFp318S5dumjNmjW6++679cYbb6hv376f+po1NTW67777tHTpUmVnZ2v9+vUnP1dWVqbCwkJJUklJ\nibZs2aKuXbtq5cqVuvbaayVJdXV16tOnjyRp5cqVevjhh7Vv3z4dOnRIU6dOlSTNnTtXL730kiTp\nS1/6kr773e9+qo63335bixcv1oQJEyRJR48eVc+ePfXZz35WmzZt0v33368bbrhB1113XZM/p7y8\nPE2bNk2SNHbsWLVr1065ubkaO3astmzZ0uTfu7E33nhDy5cv14svvihJ2r9/vzZs2ECwAgCgtdpg\nFqjN9enTR8eOHdOSJUvOGKweffRR9erVS8uWLVN9fb3y8/NPfq5du3Ynr7Ozs1VbWyvvvUaPHq35\n8+d/6mvNmDFDM2fO1Lhx4/TMM8/o3XffPfm5ploSeO9111136V//9V8/9blly5bp9ddf1+OPP67f\n/e53evrpp8/5tXJzc09+v6ysrJN/j6ysLNXW1jb59z69rp/+9KcnQ2JbYI0VACRQfb1UWSktWCDN\nmSO9+660bJm0f3/oypBokydP1gsvvKC6ujpVV1dr9uzZKisrO+vjktS1a1e9+uqr+vu///tTgk7E\n/v371adPH2VlZenXv/616urqzlnDiBEjVF1dfTJY1dTUaNWqVZKkgwcPqk+fPqqpqdFzzz138s9M\nmjRJzz//vCSd8nhjV199tV588UXt2rVLkvTJJ5+osrJSu3fvVn19vW699VZ9//vf14cffihJ6ty5\nsw4ePNiMn15sf+/Tv+7UqVP185//XDU1NZKk9evX6/Dhwy3+vmfCiBUAJMCmTdLbb0urV0sNa4A/\n5fzzpUmTpMsvl7p0SWx9SLxbbrlF8+fP17hx4+Sc07/927+pd+/eZ3187dq1kqRevXrplVde0fXX\nX6+nn35aEydOPPk17733Xt1666361a9+pWnTpqljx47nrCEvL08vvviivvWtb2n//v2qra3VX//1\nX2v06NH63ve+p4kTJ6qgoEATJ048GVB+8pOf6Itf/KJ++MMfavr06Wf8ukVFRfr+97+v6667TvX1\n9crNzdVjjz2m9u3b6ytf+Yrq6+sl6eSI1owZM3TPPfeoffv2Zxw9a8rZ/t7FxcXKzs7WuHHjNGPG\nDD3wwAPasmWLxo8fL++9CgoKNHPmzGZ/v3NxkV0BiVZaWuojC+oAIF1VVkp//KO0Zk30se7dpUGD\npHbtpOxsadcuacsWKXKKRk6OdNll0i23SGeZ0UALrVmzRqNGjQpdBpLcmf4/cc4t9t5/enX9aRix\nAoA24L30+uvSn/5k03/t2klXXWWBqUePMz9//Xob1Vq+PDpF+KUvSaNHJ7x8AC1EsAKAODt8WHry\nSalh5kZXXSXdeKN0rlkZ56QRI+y2Y4f0y1/aKNZ//Ic0bZp08832HADJjWAFAHF04ID06KMWjjp3\nlmbMkMaMad7X6NtX+u53pTfflGbOlF57TfrkE+muu2yaEK3jvecgZpxVa5dIsSsQAOJk717pkUcs\nVPXpIz38cPNDVURWljR1qnTffTaNuHCh9NOfSg2bmdBC+fn52rNnT6t/eSI9ee+1Z8+es7ZriAWL\n1wEgDg4dkn7wA6m6WurXT3rgARuxiodt22xK8MABadw46Z57LHih+WpqalRVVaVjkZ0CwGny8/NV\nWFio3NzcUx6PdfE6wQoAWqm2Vvrxj6UNGyxUffvbUocO8f0eO3ZIP/qRtWq49FLpy19mzRWQSLEG\nK97zAEAreC8995yFqq5dbeou3qFKsnVX998v5eVJ8+ZJDUefAUgyBCsAaIV337Wgk5sr3Xuvhau2\nMniwHRzsnAWrhnNsASQRghUAtNCOHVLDWa6aMUMaMKDtv+eYMVKk2fVTT0l79rT99wQQO4IVALRA\nba0Fm9paO4amtMmVF/EzbZpUXGzrrZ54wmoAkBwIVgDQAn/8o1RVJRUUSF/4QmK/t3PSV75iHdwr\nK6VZsxL7/QGcHcEKAJpp82Y7eiYrS/rqV63PVKJ16GDTj85Jf/6zBSwA4RGsAKAZ6uttF6D30rXX\n2mHKoQwbZsfl1NdLzzzDlCCQDAhWANAM771nDTvPP1+64YbQ1dgZgj172kL6V18NXQ0AghUAxGj/\nfju7T5Juvz3MFODp8vJsSlCSXn9d+vjjoOUAGY9gBQAxmjlTOnbMduQVF4euJmrIEOvGXlcnvfCC\nTVMCCINgBQAx2L5dmj9fys6W/tt/S77jZD73Oal9e2nVKmn58tDVAJmLYAUAMfjDH2wk6PLLrcVC\nsuncOdo49IUXpJqasPUAmYpgBQBNWL/ejo/Jz0+OBetnM2WKVFho3djfeSd0NUBmIlgBwDl4L730\nkl1fd52NDCWrrCzp1lvt+rXXpMOHw9YDZCKCFQCcw4oV0pYtUpcu0jXXhK6maaNGSSNH2nE3r70W\nuhog8xCsAOAsvI/2hpo6NTnaKzTFOVvILtl04N69YesBMk1Mwco5N805t845V+Gce+gMn7/CObff\nObe04fZP8S8VABJrzRobrercWZo8OXQ1sRswQJowwTqxv/xy6GqAzNJksHLOZUt6TNL1kook3eGc\nKzrDU+d470sabv8c5zoBIKEaj1Zdc01qjFY1Nn26rblasECqrg5dDZA5YhmxKpNU4b3f5L0/Iel5\nSdPbtiwACGvDBqmiwg47vuKK0NU0X0GBdPHFdo7gn/8cuhogc8QSrC6QtK3Rx1UNj53uUufccufc\nn51zo8/0hZxzdzvnyp1z5dW8hQKQxCJh5Oqrrc1CKrr+eltzNX++tHt36GqAzBCvxesfSurvvS+W\n9FNJM8/0JO/9k977Uu99aUEydtgDANmBxqtX2zl8V14ZupqW69lTmjjRRq3YIQgkRizBarukfo0+\nLmx47CTv/QHv/aGG61mScp1zPeJWJQAk0Ftv2f2ll0odO4atpbU+8xkbtZo71xqHAmhbsQSrRZKG\nOecGOefyJN0u6ZR9Js653s7ZyVnOubKGr8s/YQAp58AB6YMPLIxcfXXoalqvVy/bIVhfL739duhq\ngPTXZLDy3tdKuk/S65LWSPqd936Vc+4e59w9DU+7TdJK59wySf8h6XbvOV8dQOp57z1rU1BcbFNp\n6WDqVLufM4du7EBby4nlSQ3Te7NOe+zxRtc/k/Sz+JYGAIlVU2PBSkqNLuuxKiyURo+WVq2yv99n\nPhO6IiB90XkdABqUl0sHD0r9+knDhoWuJr6uu87u33nHAiSAtkGwAoAGkdGqK6+0NVbpZMQIqX9/\nC44LFoSuBkhfBCsAkLR1q7R5szUEnTAhdDXx51x0rdWbb1pneQDxR7ACAEVHqy65xPpXpaPx46Xz\nz5c+/tjWWwGIP4IVgIx39Ki0cKFdX3552FraUlZW9Hied94JWgqQtghWADLeggXSiRPSyJFS796h\nq2lbl10m5ebaiNVHH4WuBkg/BCsAGc17afZsu54yJWwtidCxox3OLEl/+UvYWoB0RLACkNEqK+1s\nwM6drSloJrjqKrufP9+mQQHED8EKQEabO9fuJ06UcmJqmZz6+va1ac/jx6V580JXA6QXghWAjHXi\nRHTR+qRJYWtJtMgi9tmzab0AxBPBCkDGWrJEOnZMGjjQRnEySXGx1LWrLWDfsCF0NUD6IFgByFiR\nabBMG62SpOxs2yEoRXt4AWg9ghWAjLR7t7R2rbUeSMdO67G47DLryL5kiXTgQOhqgPRAsAKQkSLn\n5Y0fL7VvH7aWULp1synBurroIn4ArUOwApBxvJc++MCuJ04MW0tokd5dc+awiB2IB4IVgIxTWSnt\n2iV16SKNGhW6mrCKiqTu3aU9e6R160JXA6Q+ghWAjBOZBiwrs/PzMplzdvC0RE8rIB4y/CUFQKap\nq5MWLbLrTJ8GjIgEqw8/lI4cCVsLkOoIVgAyyqpV0qFDUp8+Ur9+oatJDj16SCNGSDU1Unl56GqA\n1EawApBRGi9ady5sLckk0suL3YFA6xCsAGSM48elZcvsuqwsbC3J5sILpfx8acsWaefO0NUAqYtg\nBSBjrFhh012DB9tOOETl5UUbpbKIHWg5ghWAjLF4sd1fdFHYOpLVpZfa/fz5tsgfQPMRrABkhOPH\nbcRKIlidzaBBUu/e0sGD0sqVoasBUhPBCkBGaDwN2K1b6GqSk3MsYgdai2AFICNEpgFLS8PWkewu\nvtiapq5YYSNXAJqHYAUg7TWeBhw/Pmwtya5LF2nMGKm+PtqaAkDsCFYA0t7y5TYNOGQI04CxiCxi\nnzuXg5mB5iJYAUh77AZsnrFjpU6dpB07pKqq0NUAqYVgBSCtHT8e3eHGNGBscnKia9GYDgSah2AF\nIK0xDdgykc70ixbZeisAsSFYAUhrTAO2zODBdjjzvn3S+vWhqwFSB8EKQNo6diw6DUiwah7noqNW\nCxeGrQVIJQQrAGkr0hR0yBCpa9fQ1aSeSLD68EP7OQJoGsEKQNqiKWjr9Okj9esnHT3KETdArAhW\nANISuwHjIzJqxe5AIDYEKwBpafVqm74aNIhpwNYoK7P1VitWSEeOhK4GSH4EKwBpaelSux83Lmwd\nqa5rV2n4cKm2VlqyJHQ1QPIjWAFIO/X10bMBS0rC1pIOJk60e3YHAk0jWAFIOxUV0uHDUs+eUu/e\noatJfRdeaN3Y162zvlYAzo5gBSDtLFtm9yUltj4IrdOhg50f6L11YgdwdgQrAGnFe9ZXtQWahQKx\nIVgBSCs7d0q7d0udO9uxLIiPsWOl9u2lrVuljz4KXQ2QvAhWANJKZLRq7Fgpi1e4uMnNtbVWUrTx\nKoBP42UHQFppvL4K8RU5b7G8PGwdQDIjWAFIG/v2SVu22OjKqFGhq0k/I0faQvYdO+wG4NMIVgDS\nxvLldl9UJOXlha0lHeXkMB0INIVgBSBtsBuw7UUOtC4vtx2YAE5FsAKQFo4dswaWzknFxaGrSV8j\nRkgdO9rOQKYDgU8jWAFIC6tW2Xl2gwdbqwW0jexsafx4u2YRO/BpBCsAaYHdgInDdCBwdgQrACmv\nri566DLrq9re8OE2Krhrl1RVFboaILkQrACkvI0bpSNHpF697Ia2lZXF7kDgbAhWAFJeZLSKReuJ\nM2GC3TMdCJwqpmDlnJvmnFvnnKtwzj10judNcM7VOudui1+JAHBuK1fa/ZgxYevIJEOHSl26SNXV\ndn4gANNksHLOZUt6TNL1kook3eGcKzrL834o6Y14FwkAZ7Nnj237z8+3X/ZIjKys6O5ApgOBqFhG\nrMokVXjvN3nvT0h6XtL0MzzvfkkvSdoVx/oA4Jwi04BFRdYZHInD7kDg02IJVhdI2tbo46qGx05y\nzl0g6RZJP49faQDQtEiwGjs2bB2ZaOhQ6bzzbNSwsjJ0NUByiNfi9R9L+q73vv5cT3LO3e2cK3fO\nlVdXV8fpWwPIVDU11m1dYn1VCM5FpwM//DBsLUCyiCVYbZfUr9HHhQ2PNVYq6Xnn3BZJt0n6v865\nm0//Qt77J733pd770oKCghaWDABm3ToLVwMG2EJqJF7jYMV0IBBbsFokaZhzbpBzLk/S7ZJebvwE\n7/0g7/1A7/1ASS9Kutd7PzPu1QJAI0wDhjd0qDULra6Wtp/+lhvIQE0GK+99raT7JL0uaY2k33nv\nVznn7nHO3dPWBQLAmXhPsEoGWVnRY4SYDgSkmPbQeO9nSZp12mOPn+W5M1pfFgCc20cf2aLpzp1t\nKhDhjB8vzZljweqmm0JXA4RF53UAKSkyWjV6tC2iRjgjRkgdOkg7d9oNyGQEKwApiWnA5JGdHZ0O\nXLIkbC1AaAQrACnn6FGposLW9xR96hwIhBA5lJl1Vsh0BCsAKWf1aqm+XhoyxKagEF5RkR0rtG2b\ntHt36GqAcAhWAFIO04DJJycn+t+DUStkMoIVgJTivbRypV0TrJILXdgBghWAFLN1q3TwoNS9u9Sn\nT+hq0Njo0VJurrR5s7R3b+hqgDAIVgBSyqpVdk+bheTTrh3TgQDBCkBKiQQrdgMmp8h0IG0XkKkI\nVgBSxtGj0qZN1mZh5MjQ1eBMxo61hewVFdKBA6GrARKPYAUgZaxda20WBg+W2rcPXQ3OJD/fRhO9\nl5YuDV0NkHgEKwApY/Vqux89OmwdODeahSKTEawApATvWV+VKsaNs+nadeukw4dDVwMkFsEKQErY\ntUvas0fq1EkaMCB0NTiXjh1tDVx9vbRsWehqgMQiWAFICZHRqlGjaLOQCiK7AxcvDlsHkGgEKwAp\noXH/KiS/khILwGvW2G5OIFMQrAAkvdpaaf16u2Z9VWro3FkaNkyqq5OWLw9dDZA4BCsASa+iQjpx\nQioslM47L3Q1iBXNQpGJCFYAkh67AVNTpO3CypXS8eNhawEShWAFIOnRvyo1de1qzVxraqLhGEh3\nBCsASW3/fqmqSsrLk4YODV0NmisyasV0IDIFwQpAUouMVo0YYWfQIbVEgtXy5bYJAUh3BCsASY02\nC6mtoMA2HRw7Zmc9AumOYAUgadXXR0esWLieujg7EJmEYAUgaW3damfNde8u9ewZuhq0VKTtwrJl\nFpaBdEawApC0Gu8G5Bib1NWnj9Srl3TokLRhQ+hqgLZFsAKQtOhflR6cY3cgMgfBCkBSOnpU2rRJ\nysqSRo4MXQ1aq3Gw8j5sLUBbIlgBSEpr19p6nMGDpfbtQ1eD1howQOrWTdq3T9q8OXQ1QNshWAFI\nSrRZSC/OcXYgMgPBCkDS8Z5jbNIR04HIBAQrAEln1y5pzx6pUyepf//Q1SBehgyROneWqqul7dtD\nVwO0DYIVgKQTmQYcNYo2C+kkK0sqKbFrpgORrghWAJIO66vSF13Yke4IVgCSSm2ttH69XdO/Kv2M\nGCF16CDt2CF9/HHoaoD4I1gBSCoVFdKJE3Zw73nnha4G8ZaTIxUX2zXTgUhHBCsASYVu6+mPLuxI\nZwQrAEmF9VXpr6hIysuTtmyRPvkkdDVAfBGsACSNfftsG35enjR0aOhq0Fby8qQxY+yaUSukG4IV\ngKQRaQo6YoStxUH6ogs70hXBCkDSoNt65hg71sJzRYV04EDoaoD4IVgBSAr19dFgxcL19Jefbw1g\nvZeWLQtdDRA/BCsASWHrVunwYal7d6lnz9DVIBHYHYh0RLACkBQaTwNyjE1mGDfOjrlZs0Y6ciR0\nNUB8EKwAJIWVK+2eacDM0amTNHy4TQOvWBG6GiA+CFYAgjt6VNq82UYvRo4MXQ0SibMDkW4IVgCC\nW7vWRi2GDJHatw9dDRKppMTuV62Sjh8PWwsQDwQrAMFxjE3m6tpVGjxYqqmJ/n8ApDKCFYCgvKd/\nVaZjOhDphGAFIKhdu6Q9e2whc//+oatBCJFgtWKFVFsbthagtQhWAIKKTP+MGkWbhUxVUCAVFkrH\njtl6OyCVEawABBUJVkwDZrbI2YFMByLVEawABFNTI61bZ9csXM9skenApUtthyiQqghWAIKpqLBw\nVVgonXde6GoQUp8+Uq9edqzRhg2hqwFajmAFIBh2AyLCueh0IGcHIpXFFKycc9Occ+uccxXOuYfO\n8Pnpzrnlzrmlzrly59xl8S8VQLqhfxUaa3wos/dhawFaqslg5ZzLlvSYpOslFUm6wzl3+svg25LG\nee9LJP2VpF/Eu1AA6WXfPmn7dikvTxo6NHQ1SAb9+0vnn2//b2zeHLoaoGViGbEqk1Thvd/kvT8h\n6XlJ0xs/wXt/yPuT7y86SuK9BoBzikwDjhgh5eSErQXJwblTR62AVBRLsLpA0rZGH1c1PHYK59wt\nzrm1kl6VjVoBwFmxvgpn0rgLO9OBSEVxW7zuvf+j936kpJslfe9Mz3HO3d2wBqu8uro6Xt8aQIqp\nr48GK9ZXobEhQ6TOnaXdu22qGEg1sQSr7ZL6Nfq4sOGxM/Lez5Y02DnX4wyfe9J7X+q9Ly0oKGh2\nsQDSw9attq2+e3epZ8/Q1SCZZGUxHYjUFkuwWiRpmHNukHMuT9Ltkl5u/ATn3FDn7DAK59x4Se0k\n7Yl3sQDSQ+Nu6xxjg9NxKDNSWZNLRr33tc65+yS9Lilb0tPe+1XOuXsaPv+4pFslfdk5VyPpqKQv\nNFrMDgCniASrMWPC1oHkNHy41KGDtGOH9PHH1jgUSBUx7cXx3s+SNOu0xx5vdP1DST+Mb2kA0tGR\nI9KmTTblM2JE6GqQjHJypOJiacECmw6cNi10RUDs6LwOIKHWrrXdXkOHSvn5oatBsuJQZqQqghWA\nhFq50u5ps4BzKSqy5rGVldIeVuwihRCsACSM96cuXAfOJjfXpgMldgcitRCsACTMzp12XEmXLlJh\nYehqkOwuusjuy8vD1gE0B8EKQMI0ngakzQKaMmaMTQdu3sx0IFIHwQpAwjANiObIy4tOBy5eHLYW\nIFYEKwAJcfy4VFFhI1UcY4NYlZbaPcEKqYJgBSAh1q2TamulgQOljh1DV4NUMWaM1K6dtGUL04FI\nDQQrAAkROXSZaUA0R+PdgYxaIRUQrAAkBOur0FLsDkQqIVgBaHO7dtmtQwebCgSaIzIdWFkp7d4d\nuhrg3AhWANpcZBqwqMjOCASag+lApBJe4gC0OaYB0VrsDkSqIFgBaFO1tXbwskSbBbTc6NHR6cDq\n6tDVAGdHsALQpioqpBMn7Aibrl1DV4NUlZsrlZTYNaNWSGYEKwBtimlAxEtkdyDBCsmMYAWgTTU+\nHxBojaIiKT9f2rrVdpkCyYhgBaDN7N0r7dhha2OGDAldDVJdbq40bpxdM2qFZEWwAtBmIm0WRo6U\ncnLC1oL0wO5AJDuCFYA2w/oqxFtkOnDbNqYDkZwIVgDaRH29tGaNXROsEC85OewORHIjWAFoE5s3\nS0eOSD17Sj16hK4G6YSzA5HMCFYA2gTTgGgrRUVS+/ZSVZX08cehqwFORbAC0CYiwWrMmLB1IP0w\nHYhkRrACEHcHD9rRIzk50rBhoatBOopMBy5aFLYO4HQEKwBxt2qV5L00fLj1sALibdQoqUMH65O2\nY0foaoAoghWAuFuxwu7Hjg1bB9JXTk501GrhwrC1AI0RrADEVX19tDEo66vQlsrK7H7hQhshBZIB\nwQpAXG3caG0WevWyVgtAWxk2TOraVdqzx9p7AMmAYAUgrpgGRKI4J02YYNdMByJZEKwAxBXBCokU\nmQ4sL7dpaCA0ghWAuPnkE9uhlZ8vDR0auhpkgn79bNr54EFp7drQ1QAEKwBxFBmtGjXKdm0Bbc25\nUxexA6ERrADEDdOACCGyzmrJEqmmJmwtAMEKQFzU1ESnYmizgETq1UsaMEA6diwa7oFQCFYA4mLd\nOgtX/ftL550XuhpkGqYDkSwIVgDigmlAhFRaauutVqyQjh4NXQ0yGcEKQKt5Hw1WTAMihK5d7WzK\n2lpbawWEQrAC0GoffWTdrzt1kgYODF0NMhXNQpEMCFYAWq3xaFUWryoIZPx4KTvbNlHs2xe6GmQq\nXgIBtBrrq5AMOnaUiottappRK4RCsALQKkePShUVNlJVVBS6GmS6Sy6x+/nzLWABiUawAtAqq1bZ\nGW1Dhkjdomo2AAAbk0lEQVQdOoSuBplu9GgbudqxQ6qqCl0NMhHBCkCrLFtm9+PGha0DkOwopUhP\nqwULwtaCzESwAtBidXXSypV2TbBCsrj4Yrv/4AMbTQUSiWAFoMU2bJCOHJH69JF69gxdDWAGDJB6\n95YOHpRWrw5dDTINwQpAiy1fbvfFxWHrABpzLjpqNX9+2FqQeQhWAFrEe2npUrsuKQlbC3C6iRMt\nYC1daqOqQKIQrAC0yI4d1m29c2e6rSP5nH++NGKEHXHz4Yehq0EmIVgBaJHIbsDiYrqtIzlFpgPZ\nHYhE4uUQQIvQZgHJ7sILpbw822Sxe3foapApCFYAmm3fPmnLFik3Vxo1KnQ1wJnl51u4khi1QuIQ\nrAA0W+RswKIiGxEAklXkiJsFCzjiBolBsALQbJHdgEwDItmNGCF17SpVV9uUINDWCFYAmuX4cWnt\nWtvKPnZs6GqAc8vKki691K7nzg1bCzJDTMHKOTfNObfOOVfhnHvoDJ+/0zm33Dm3wjk3zznH+1gg\nTa1ebVvYBw2SunQJXQ3QtEmT7H7xYuno0bC1IP01Gaycc9mSHpN0vaQiSXc454pOe9pmSVO892Ml\nfU/Sk/EuFEByYDcgUk2PHjYlWFMjLVoUuhqku1hGrMokVXjvN3nvT0h6XtL0xk/w3s/z3u9t+HCB\npML4lgkgGdTXR4+xIVghlURGrd5/P2wdSH+xBKsLJG1r9HFVw2Nn81VJf25NUQCS08aN0uHDduBy\n796hqwFiN3681KGDVFkpVVWFrgbpLK6L151zV8qC1XfP8vm7nXPlzrny6urqeH5rAAmwZIndl5TY\n4nUgVeTmSmVlds2oFdpSLMFqu6R+jT4ubHjsFM65Ykm/kDTde7/nTF/Ie/+k977Ue19aUFDQknoB\nBOJ99My1SNNFIJVcdpndf/CBrbcC2kIswWqRpGHOuUHOuTxJt0t6ufETnHP9Jf1B0pe89+vjXyaA\n0Corpb17rSfQoEGhqwGar18/ux05Eu3FBsRbk8HKe18r6T5Jr0taI+l33vtVzrl7nHP3NDztnyR1\nl/R/nXNLnXPlbVYxgCAaj1YxDYhUFRm1oqcV2kpOLE/y3s+SNOu0xx5vdP01SV+Lb2kAkkXjacDx\n48PWArRGWZn04ovSmjXSnj1S9+6hK0K6ofM6gCZt325HgnTuLA0dGroaoOU6dIi+OZg3L2wtSE8E\nKwBNajwNmMWrBlJcpKfV3LnWmw2IJ14iATSJaUCkk+HDpYIC24yxalXoapBuCFYAzmnnTrt16GC/\nkIBU55x0+eV2/e67QUtBGiJYATinxk1Bs7PD1gLEy6WXSjk5NmK1e3foapBOCFYAzommoEhHnTpJ\npaW243X27NDVIJ0QrACc1e7d0rZtUn6+VFQUuhogvqZMsfu5c6Xa2rC1IH0QrACcVWS0qrjYpk2A\ndDJokHViP3RIWrw4dDVIFwQrAGfFbkCkM+eio1bvvRe2FqQPghWAM9q7V9q8WcrLk0aPDl0N0DbK\nymyqe+NGqaoqdDVIBwQrAGcUGa0aM8bCFZCO2rWTLrnErhm1QjwQrACc0cKFdl9aGrYOoK1Felp9\n8IF07FjYWpD6CFYAPmXXLmnLFns3X1wcuhqgbfXta81vjx+XFiwIXQ1SHcEKwKeUl9v9hRdKublh\nawESofEidu/D1oLURrACcArvo9OAEyaErQVIlJISqUsXaccOacOG0NUglRGsAJxi+3Y7G7BjR2nU\nqNDVAImRkxNda/XWW2FrQWojWAE4ReNF65wNiEwyZYoFrOXLbZ0h0BIEKwAneS8tWmTXTAMi03Tp\nYn2tvJfeeSd0NUhVBCsAJ23cKH3yidStmzR0aOhqgMS7+mq7nzdPOnIkbC1ITQQrACc1Hq1yLmwt\nQAiFhdLIkdZ64f33Q1eDVESwAiBJqquLtlkoKwtbCxBSZNTqL3+R6uvD1oLUQ7ACIElau1Y6dEjq\n3dvetQOZauxYqVcvmxaPHO0ExIpgBUDSqb2rmAZEJnNOuuoqu3777bC1IPUQrACopkZassSumQYE\n7GDmDh2kTZvsBsSKYAVAK1bYYt0BA6SePUNXA4TXrp00ebJdM2qF5iBYAdD8+XbPaBUQdeWVUlaW\nrbP65JPQ1SBVEKyADHfggLRypf0CmTgxdDVA8ujWzU4gqK/nmBvEjmAFZLgFC+wXx9ixUufOoasB\nksvUqXY/e7Z08GDYWpAaCFZABvPeOkxL0qRJYWsBklFhoVRcbBs8OOYGsSBYARls82Zp504bqRoz\nJnQ1QHK6/nq7/8tfpKNHw9aC5EewAjJYZLTq4oul7OywtQDJavBgafhwC1XvvRe6GiQ7ghWQoU6c\niJ4NeOmlYWsBkl1k1Oqtt2xaEDgbghWQoZYskY4dkwYOlPr2DV0NkNxGjbI+bwcPSnPnhq4GyYxg\nBWQoFq0DsXMuOmr1+ut2aDlwJgQrIAPt2SOtWyfl5lqfHgBNKymxQ8o/+SR6tiZwOoIVkIHmz7dW\nCyUldh4agKY5J02bZtevvWb/hoDTEayADEPvKqDlysqk88+XPvooenA50BjBCsgw69fbVGC3btKI\nEaGrAVJLdnZ01Orll+3UAqAxghWQYSI7mi65xM4HBNA8kyZJ3btbc91IyxIggpdVIIMcOiQtXmxr\nRS67LHQ1QGrKyZFuvNGuX3mFHYI4FcEKyCDz5km1tdLo0faOG0DLXHyx1KuXtGuXbQYBIghWQIbw\nXpozx66nTAlbC5DqsrKkz37Wrl95xd6wABLBCsgYa9fau+tu3ThwGYiH0lI7tWDv3uibFoBgBWSI\nyOGxl1/OonUgHpyTpk+361mz7PxNgJdXIAPs2yctW2aBit5VQPyMG2dnCB44IL37buhqkAwIVkAG\neO8967dTUiKdd17oaoD00XjU6rXX7GBzZDaCFZDmamqi6z+uuipsLUA6KiqShg6VDh+W3nwzdDUI\njWAFpLnycungQalfP3vxBxBfzkm33GLXb7xhU+/IXAQrII15L73zjl1fdZX9AgAQf0OHSuPH2wL2\nmTNDV4OQCFZAGtu0Sdq6VerUSZowIXQ1QHr73OesK/v8+VJlZehqEArBCkhjkdGqyZOl3NywtQDp\nrqAguo7x97+3EWNkHoIVkKb27LFzAbOy6LQOJMr119sI8YYN0tKloatBCAQrIE29/ba9Yy4rs27r\nANpehw7Ro25eeomjbjIRwQpIQ0eOSO+/b9fXXhu2FiDTXH651KePVF1N09BMRLAC0tDs2dLx49Ko\nUVJhYehqgMySlSXddptdv/KKdOhQ2HqQWDEFK+fcNOfcOudchXPuoTN8fqRzbr5z7rhz7jvxLxNA\nrGpro4vWGa0Cwhg92hqHHj0qvfxy6GqQSE0GK+dctqTHJF0vqUjSHc65otOe9omkb0l6JO4VAmiW\nDz6Q9u+X+va1F3YAieec9PnP2+jV7NnSli2hK0KixDJiVSapwnu/yXt/QtLzkqY3foL3fpf3fpGk\nmjaoEUCM6uvtvDJJmjaNhqBASH372qix99Kzz9q/T6S/WILVBZK2Nfq4quExAEnmww+lXbukHj1o\nCAokgxtvlLp3l7Zti07RI70ldPG6c+5u51y5c668uro6kd8aSHveS3/+s11Pm2ZTEADCysuT7rjD\nrl9+Wdq7N2w9aHuxvPRul9Sv0ceFDY81m/f+Se99qfe+tKCgoCVfAsBZrFghVVVJXbtKl1wSuhoA\nEWPH2jmCx49Lzz8fuhq0tViC1SJJw5xzg5xzeZJul8QeByCJeC/NmmXX115r55UBSB5f+ILUrp11\nY1+2LHQ1aEtNBivvfa2k+yS9LmmNpN9571c55+5xzt0jSc653s65KknflvSwc67KOdelLQsHELVq\nlbR5s9S5s50LCCC5dO0qTW/Y9vXb39roFdJTTO9rvfezJM067bHHG11/JJsiBJBg3kf75Eydau+K\nASSfK6+UFiyQtm61f7Of/3zoitAWWN4KpLgVK6TKSqlLFw5bBpJZVpb03/+7tUF5+207qBnph2AF\npLDGo1XTptkOJADJa8AA6frr7d/uM89Ix46FrgjxRrACUtiSJdYfp2tXO/gVQPK74QapXz9p927p\nxRdDV4N4I1gBKaquTpo5066vv17KzQ1bD4DY5ORIf/VXdj9njrRyZeiKEE8EKyBFzZ0rffyx1LMn\nOwGBVNO3r3TTTXb9q19Jhw+HrQfxQ7ACUtDx49J//Zdd33yzlJ0dth4AzXfttdKQIXZo+m9/G7oa\nxAvBCkhBb78tHTggDRxoHZ0BpJ6sLOkrX7EWKYsWSeXloStCPBCsgBRz4ID02mt2feuttnUbQGoq\nKJBuu82un33WFrQjtRGsgBTzxz/aVGBxsTR8eOhqALTW5MnSuHHS0aPSE09INTWhK0JrEKyAFFJZ\nKc2fb2uq6NoMpAfnpBkzpB49rCv7738fuiK0BsEKSBHeSy+8YPdXX227AQGkhw4dpLvvthYM770n\nLVwYuiK0FMEKSBGLFkkbN9pByzfcELoaAPE2YID0hS/Y9bPPSjt3hq0HLUOwAlLAkSPR6YFbbpHy\n88PWA6BtTJ4slZXZOsonnrB7pBaCFZAC/vQn2w04ZIh06aWhqwHQVpyzg5p797YRq+ees+l/pA6C\nFZDktmyxNRdZWdKdd9JeAUh37dpJ3/iGHar+wQfR9ipIDQQrIInV19taC++tS/MFF4SuCEAi9O0r\nffWr9kZq5kxp8eLQFSFWBCsgib3+urRtm9S9OwvWgUxTUmJNgCXp//0/afPmsPUgNgQrIEnt3Cm9\n8opdf+lLNj0AILNcc40taK+pkR57TNqzJ3RFaArBCkhC9fXSM89ItbXSZZdJo0aFrghACM5Jd9wh\njRwpHTwo/exn1qEdyYtgBSShN9+0RevdukXPEQOQmbKzbTF7797Sjh3Sk0/amy4kJ4IVkGQqK22x\nqmRTgO3bh60HQHgdOkj33y916iStXi099ZSNbCP5EKyAJHL8ePQF88orpdGjQ1cEIFn06CE98IC9\n2frwQ+mXv6THVTIiWAFJ5Pe/lz7+2LZaR3YDAUBE//42ctWunbRggfSb3xCukg3BCkgSixZJc+bY\nIaxf+5qUmxu6IgDJaMgQ6d577bVi9mzppZcIV8mEYAUkgZ07pV//2q4//3kagQI4t5EjpW9+0xa2\nv/lmtDULwiNYAYE1Pmy1rEyaMiV0RQBSwZgxNrrtnAUrRq6SA8EKCMh76Ve/shGrPn3s8FXOAgQQ\nq/HjLVxlZUlvvGEj3+wWDItgBQT06qtSebmUn299auiuDqC5Skul++6zdZlz51qfq5qa0FVlLoIV\nEMjixdJ//ZeNUH396zZiBQAtMXq09Dd/Y/2uliyxDu3HjoWuKjMRrIAANm+2Q1Ula6swZkzYegCk\nviFDpL/9W6lLF2ntWunf/13avz90VZmHYAUk2McfSz/9qQ3VT5pkh6wCQDwUFkp/93fWTLSyUvo/\n/8eOx0LiEKyABNq3T/rJT6TDh22U6s47WawOIL4KCqSHHpKGDrXXnEcekT74IHRVmYNgBSTIwYMW\nqvbskQYNku6+23rQAEC8de5sa64mT7bR8aeftnYM7BhsewQrIAEOH5Z+/GM7mb5PH9vBww5AAG0p\nJ8dGxb/4xWg7hsces9cjtB2CFdDGIqGqqkrq1Uv69rfthHoAaGvOWdPhv/kbe91ZuVL653+2xe1o\nGwQroA1F1jds3Sr17GmhqkuX0FUByDTDh0v/8A+2c3DfPnuz99JLUm1t6MrSD8EKaCO7d0s/+lF0\n+u9v/1bq2jV0VQAyVffu0ne+I910k41kvfGG9IMfSB99FLqy9EKwAtrAli3SD39o4WrAAHsxI1QB\nCC0rS7rhBunBB60lw7Zt0ve/bwc5s7A9PghWQJwtWWLTfwcO2An0rKkCkGwGD5b+5/+ULr7Ydg2+\n+KL0L/9izYvROs4HOgq7tLTUl5eXB/neQFvwXpo1y46p8d6af37xi7YzBwCS1cqV0m9+Y61gnJMu\nv1y6+WY7HgdRzrnF3vvSJp9HsAJa78gR6xOzYoW9MN18szR1Ks0/AaSGEyfsUPg33rApwS5d7HXs\nkkts+hAEKyBhNm6UnnrK3u116CB97Wt2ICoApJodO6Rnn7XXNck23txyi1RczBtFghXQxurr7R3e\nq6/a1N+AAdZNvUeP0JUBQMt5Ly1aJM2caW8YJWvTcOutdp+pCFZAG9q2TfrlL+3eOem662wLM+up\nAKSL2lpp9mx783jokD1WVGTLHEaMyLwRrFiDFb8GgGY4ftxeZCJbk7t3l+66y15kACCd5ORIV10l\nXXqprb166y1p9Wq7DRggTZsmlZSwBut0jFgBMfBeWrhQ+sMfrGuxc/aCM306Z/4ByAyHD0vvvSe9\n844dKi/ZiRJXXSVNnJj+uwiZCgTiwHtpzRpba1BZaY8NHCjdcYfdA0CmOXFCmjfPRu5377bHcnOl\niy6yVg2DB6fnNCHBCmgF7+2Q0lmzpPXr7bEuXWx3zCWXpOeLBgA0R329tHSprcNasyb6eJ8+Nn14\n0UW2XCJdEKyAFqittc7pb71lx9JINrw9bZp05ZVSXl7Q8gAgKVVXS++/L82dG50mlGz0qrTUQlaq\nH+tFsAKaobpamjPHXhQiu186d5auvlqaMiX91w4AQDzU1lqj5PJyadkyOy5HslH+wYOlMWOksWOl\nwsLUG/knWAFNOHpUWr5cmj//1GHswkJbJ3DJJYxQAUBLHT8eDVkrVljoiujSxUJWUZE0fLh03nnh\n6owVwQo4g4MH7V3UkiUWpurq7PHcXGnCBGnyZGnQoNR7JwUAyezYMVu3unKl3fbuPfXzPXtKw4ZZ\nyBo61NZmJdvrMH2sANkw9MaN0rp1dtu0yRamS/aPdvhwafz4zNgqDACh5Odbz6uSEnsN3rHDAta6\ndVJFhbRrl93mzrXnd+xovbIGDLAd2P37S926JV/YOhOCFdKG99ZjqrLSbhUVFqQaDz9nZ0sjR0oX\nXmj/wDt3DlcvAGQi56QLLrDb1Km2u3DrVmnDBtuFvXmzzS5EmpFG5Ofbn+nbN3rr1csWxSdT4CJY\nISUdOSJ99JHdPv7YjpbZuvXU3SiS/WPr1886o48YYUPN7duHqRkA8GlZWTYqNXCgdO219iZ5717b\nmR15o7x1qzUo3bgxekB0RF6ejWx95zsBij+DmIKVc26apJ9Iypb0C+/9D077vGv4/GckHZE0w3v/\nYZxrRYaoq7OdeQcP2j+uTz6J3vbsseHi0wNURIcO9g+sf39bKzV8uA0pAwBSg3PS+efbbfx4e8x7\ne93fsePUW+T3QWT3YTJoMlg557IlPSbpWklVkhY551723jcaoNP1koY13CZK+nnDPTKI9zakW1Nj\nt9pa69B77JjtwDt2LHqLfHzkiP2jaHw7cqTp75WXZ0PAvXpJvXvb8HD//sm54BEA0DrO2U7CLl1s\nOUdjR47E9nsjUWIZsSqTVOG93yRJzrnnJU2X1DhYTZf0K29bDBc457o65/p473fGveIYbdokrVp1\n6mONN0CeaTPk2TZItuTPtfZ7xVJLa+uLBKH6ehslilyf7bG6OvszdXV2i4SnSJCqqTl73c3hnNSp\nk926dYu+c4ncevZMvjl1AEAYHTok1+ajWILVBZK2Nfq4Sp8ejTrTcy6QFCxYbd4svfJKqO+eubKy\nrHVB41t+fvTWvv2p1+3b2wLyyK1TJ5u647R0AEAqSujidefc3ZLulqT+/fu36fcaNEi68cbG3/tM\n9cT++TM93tTnW/u92ro+5yzAZGfbfePb6Y9lZ9vzI49nZ0s5ORac8vKi1wQiAEAmiyVYbZfUr9HH\nhQ2PNfc58t4/KelJyRqENqvSZho82G4AAACJEsv4wiJJw5xzg5xzeZJul/Tyac95WdKXnblY0v6Q\n66sAAABCaHLEyntf65y7T9LrsnYLT3vvVznn7mn4/OOSZslaLVTI2i18pe1KBgAASE4xrbHy3s+S\nhafGjz3e6NpL+h/xLQ0AACC1sNQYAAAgTghWAAAAcUKwAgAAiBOCFQAAQJwQrAAAAOKEYAUAABAn\nBCsAAIA4IVgBAADECcEKAAAgTghWAAAAcUKwAgAAiBOCFQAAQJw4Oz85wDd2rlpSZQK+VQ9JuxPw\nfTIFP8/442caX/w844+faXzx84y/RPxMB3jvC5p6UrBglSjOuXLvfWnoOtIFP8/442caX/w844+f\naXzx84y/ZPqZMhUIAAAQJwQrAACAOMmEYPVk6ALSDD/P+ONnGl/8POOPn2l88fOMv6T5mab9GisA\nAIBEyYQRKwAAgIRI+2DlnPuec265c26pc+4N51zf0DWlOufcj5xzaxt+rn90znUNXVMqc8593jm3\nyjlX75xLil0tqco5N805t845V+Gceyh0PanOOfe0c26Xc25l6FrSgXOun3PuL8651Q3/5h8IXVMq\nc87lO+cWOueWNfw8/3fomqQMmAp0znXx3h9ouP6WpCLv/T2By0ppzrnrJL3jva91zv1Qkrz33w1c\nVspyzo2SVC/pCUnf8d6XBy4pJTnnsiWtl3StpCpJiyTd4b1fHbSwFOacu1zSIUm/8t6PCV1PqnPO\n9ZHUx3v/oXOus6TFkm7m/9GWcc45SR2994ecc7mS3pf0gPd+Qci60n7EKhKqGnSUlN5JMgG89294\n72sbPlwgqTBkPanOe7/Ge78udB1poExShfd+k/f+hKTnJU0PXFNK897PlvRJ6DrShfd+p/f+w4br\ng5LWSLogbFWpy5tDDR/mNtyC/45P+2AlSc65f3HObZN0p6R/Cl1PmvkrSX8OXQQg+wW1rdHHVeKX\nFpKUc26gpAslfRC2ktTmnMt2zi2VtEvSm9774D/PtAhWzrm3nHMrz3CbLkne+3/03veT9Jyk+8JW\nmxqa+pk2POcfJdXKfq44h1h+ngAyg3Ouk6SXJP31abMqaCbvfZ33vkQ2c1LmnAs+ZZ0TuoB48N5f\nE+NTn5M0S9L/asNy0kJTP1Pn3AxJN0q62qf7Qr04aMb/o2i57ZL6Nfq4sOExIGk0rAV6SdJz3vs/\nhK4nXXjv9znn/iJpmqSgmy3SYsTqXJxzwxp9OF3S2lC1pAvn3DRJfyfpJu/9kdD1AA0WSRrmnBvk\nnMuTdLuklwPXBJzUsNj6KUlrvPf/HrqeVOecK4jsSnfOtZdtXAn+Oz4TdgW+JGmEbNdVpaR7vPe8\ni20F51yFpHaS9jQ8tICdli3nnLtF0k8lFUjaJ2mp935q2KpSk3PuM5J+LClb0tPe+38JXFJKc879\nVtIVknpI+ljS//LePxW0qBTmnLtM0hxJK2S/kyTpH7z3s8JVlbqcc8WSfin7954l6Xfe+38OW1UG\nBCsAAIBESfupQAAAgEQhWAEAAMQJwQoAACBOCFYAAABxQrACAACIE4IVAABAnBCsAAAA4oRgBQAA\nECf/H1sD6/hHiLwdAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f042d9635f8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize=(10, 7))\n",
"ax.plot(ys, kernel, 'b-', lw=2, alpha=0.6, label='look ahead estimate')\n",
"ax.legend()\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment