Skip to content

Instantly share code, notes, and snippets.

@TomAugspurger
Created May 13, 2016 12:16
Show Gist options
  • Save TomAugspurger/7c87c26542db41f094673cb3b30174a1 to your computer and use it in GitHub Desktop.
Save TomAugspurger/7c87c26542db41f094673cb3b30174a1 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import requests\n",
"import statsmodels.api as sm\n",
"from io import BytesIO"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"wpi1 = requests.get('http://www.stata-press.com/data/r12/wpi1.dta').content\n",
"data = pd.read_stata(BytesIO(wpi1))\n",
"data.index = data.t"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Fit the model\n",
"mod = sm.tsa.statespace.SARIMAX(data['wpi'], trend='c', order=(1,1,1),\n",
" seasonal_order=(1, 1, 1, 12))\n",
"res = mod.fit()"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Fit the model\n",
"mod_s = sm.tsa.statespace.SARIMAX(data['wpi'], trend='c', order=(1,1,1),\n",
" seasonal_order=(1, 1, 1, 12),\n",
" simple_differencing=True)\n",
"res_s = mod_s.fit()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x116be5fd0>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAERCAYAAACEmDeEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8XNWd///XnaIy6pZGsmzZKi7HuGHkho1xwxTTCQRS\nSUISEsImbLJhk7BZknw3gWRJSCAJJDEt8Euy9Go6NriAsXGvx5Yt2bLVexlJ0+7vj5GMwEXyaKTR\nzHyej4cfD0+78/no2u85OnPvuYZpmgghhIgOlnAXIIQQInQk1IUQIopIqAshRBSRUBdCiCgioS6E\nEFFEQl0IIaKIrT9PUkrNBX6ttV6ilJoB3A94gS7gRq11rVLqm8DNgAf4ldZ65WAVLYQQ4uT6HKkr\npW4HVgDx3Xf9AbhVa70UeB74kVIqB/guMA+4BLhbKWUfnJKFEEKcSn+mX0qAa3rdvkFrvbP77zag\nE5gDrNNae7XWLcABYHpIKxVCCNGnPkNda/08gamWntvVAEqp+cCtwO+BVKC518vagLSQViqEEKJP\nQX1RqpS6AXgAuFRrXQ+0EAj2HilA08DLE0IIcSb69UVpb0qpLxH4QnSx1ronuDcCv1RKxQGJwCRg\nV1/b8np9ps1mPdMShBAi1hmneuCMQl0pZQHuAw4DzyulTOA9rfUvlFL3A+u63+wOrbW7r+01NrrO\n5O3DzulMoba2NdxlhFy09gXR25v0FXlC2ZvTmXLKx/oV6lrrw8D87puZp3jOw8DDZ1qcEEKI0JGT\nj4QQIopIqAshRBSRUBdCiCgioS6EEFFEQl0IIaKIhLoQQkQI0zRZ8fKe0z5HQr0fvvvdb9HS0tz3\nEwfRu+++wyOP/C2sNQghwutwdSsf7K467XMk1PvBME558taQGi51CCHCY7Ou7fM5Z7xMwFB6alUJ\nm/bVhHSbsydlc/3S8ad83Ov1cvfdv6Ci4hh+v8n1138B0zS5777f0tTUgNUax3/9188wTfjZz36C\naZq43W5++MOfMH78BJ599kneeusNDMNg2bKLuPbaG7jrrl/Q3NxES0sLY8fmM2NGMcuXX05DQz23\n3/7vPPzwE/z1r39mx45t+P0+brjhiyxefAE7d27n/vt/R0pKKna7nUmTJof0ZyGEiBymafKRriXO\nfvqx+LAO9XB48cXnSE8fwX//9//gcrm46aYvER8fx4033sSlly5jxYpHefzxR5k1azZpaen89Ke/\noLT0EJ2dHZSVlfLOO2/x4IMPY5om3//+rcyefS4AM2fO4frrP09ZWSm///3/snz55bz++qtcdtmV\nbNjwPhUVx/jzn1fgdrv51re+yqxZc/nd737DL3/5G/LyxvDXv/45zD8ZIUQ4VdS1U93gYqZynvZ5\nwzrUr186/rSj6sFw+HAps2fPBcDhcFBYWMimTR9yzjkzAZgyZRrvv7+e7373+5SXl/PjH/8Am83O\njTfexKFDB6mqquS2227BNE3a2lo5dqwcgLFj8wEoKCjE5/NRVVXFqlVvcd99D/Dii8+h9T6+971v\nY5omPp+PysoKGhrqycsbA8CMGcXs2dPnGmlCiCjVM/UybcKp130BmVM/QX5+Idu2bQXA5Wrn4MES\ncnNHsXt34Log27ZtZdy48WzZ8hGZmVnce++fuPHGm1ix4gHy8wsoKhrH/ff/hT/+8a9ccslljBs3\nAQCL5eMf9eWXX8WDD95HYWERSUnJjB1bwMyZs7j//r/whz88wNKlFzJ6dB5ZWU5KSw8BHH9/IURs\n2ry/Fpvdx8rGv5/2ecN6pB4OV131GX7zm1/yne98A7fbzU033cyrr77Mm2++xt//voKEhCR++tOf\n4/V6+dnP7uCFF57B7/fzta99k3HjxlNcPJtbbvk6brebKVOmkpV14q9KS5ZcwH33/Y7f/OZeABYs\nWMjWrZu59dZv0tHRwcKFi3E4HPzoRz/l7rv/Hw6Hg7S0NAoKiob6xyGEGAZqGl2U17SRP7mVGk/7\naZ9rmKY5RGWdqLa2NXxvHoRoXRY0WvuC6O1N+oo8A+nttQ8P8/Tqg4yZv516bxVP3vDAKQ+Fk+kX\nIYQYxrw+P2u2VWB1tFPnrWTSiAmnfb5MvwghxDC0q24vT+1/kYmW86hu9FA0q5FKYF7u7NO+Tkbq\nQggxzHj9Xp7a/yL1nQ180LaShMx6WuJKSbI7mO6cctrX9ivUlVJzlVKrP3XfvUqpm3vd/qZSapNS\n6n2l1GVBdSKEEDFidfk6/mfDb2nuOnGefV3Fh9R3NjDCyMMEjHGbaPe2M2dkMXbL6SdY+gx1pdTt\nwAogvvt2llLqVeCKXs/JAb4LzAMuAe5WStn7354QQsSOTm8Xr5a+RZWrhpWlbx6/v665gx//bS3P\n7Hkdi2mjdrvCfnQ2VsMKwPzcOX1uuz9z6iXANcAT3beTgZ8By3s9Zw6wTmvtBVqUUgeA6cDmfmxf\nCCFiygeVm3B5OzAweL9iI0vGLCDHkc0jK/fSkLAXu60L99HxeDvtfLF4Hrlji2nobGRU8sg+t91n\nqGutn1dK5fe6XQaUKaUu7fW0VKD3MoZtQFq/Oxxmqqoq+cpXPodSZ2GaJoZhUFw8i9tv/35Y6lmz\n5l2mTJlKZmZWWN5fCBE6Pr+PVeVrsVvs3DDxav6/fU/zQslKxrkvZH/jQRInHcZhT+LfLv0sXV0W\nxo1K7V7Mr7Bf2w/V0S8tBIK9RwrQ1NeLMjIc2GzWEJUQOm53CxMnTuRf//rHCY85nac/RXcwvPDC\nUxQXTx3U9w5HX0MlWnuTviKP05nCusObaOhs5KLxC7li+hK2NmxnV80+djQ1En9WNRgWvlr8WYoL\nxgb1HmcS6qdb93Uj8EulVByQCEwC+lyopLHRddrHnyt5ha01oT09/pzsaXxm/OWnfU5DQztut/eE\nEwUeeeQBNmzYiGEYXHjhxVx33ec+sQLjPffcxz/+8fcTVlvcvXsXf/zjvZimidPp5M47f8mePbt4\n9NEVmKZJR4eLn/3sV2Rn53DnnT+mvb2dzs5Obr75O3i9Hvbs2csPfvBDHnjgIWy20B+FKid8RB7p\nK/I4nSm8v6Wc+3c8CwkGyS0T2b6vmuSGacB+LOnVpNsyuXnGF8hPGnPan8PpPvjOJCE+ffbn8dta\n62ql1P3AOgLhf4fW2n0G2x52ysoOHV9gyzAMli+/nGPHjvG3vz2G1+vl1lu/SXHxLODjFRg3bHif\nysqKE1Zb/O1v7+IXv7ibsWPzWbnyJQ4fLqWs7BB33vk/ZGZm8cQTj7J69dssWLCI5uZmfve7P9LY\n2EB5+RHmzVvAxImK22+/Y1ACXQgxNHx+k0fWrIXcZnwNOfxj01HgKADxuVOYmJ/Evy24us+jW/rS\nr1drrQ8D8z913//71O2HgYcHVM2nfGb85X2OqgdLYWFgYa4e//znE8ycGVip0WazMXnyVEpLS4GP\nV2A8dKiEffv2nmS1xYbjz7nssisBqK6u4ve/vweHw0FtbQ3Tp8+gsLCIK6+8hp///A68Xh+f/ewN\nQGAd5XAu5yCEGLg3PzxMo2M3VuDrsy+nLNNKTVMH50zIonjiQhLiQjNok6HfKXw6RAsKCnnnnde4\n9NLP4PV62bVrO5deejkffvj+8RUYe1ZbvP32O/D5fDzxxKPHV1s8duwoo0fn8c9/PkFe3hjuuecu\nnnrqRRITE/nVr36OaZocOlSCy+Xif//3D9TX13HLLV9n3rwFWCwWCXUhIpir08Pj767HWtTA+NRx\nzM5XzM7v+3XBkFA/hU9fOm7+/AXs27eDb3/7JrxeL0uXXsiECeoTzznVaou33/4T7rrrF1gsFjIz\ns7j++s9z8cWXcuut3yArK5v8/ALq6moZMyafRx5ZwerVb2OaJt/4xi0ATJ06nV/+8k7uvffPpKRE\n75dIQkSrl9aX4c7YjxW4bNyyQX0vWaXxDETrlzjR2hdEb2/SV+Soa+7gjsffwj5lHYWp+fzHzO8M\n+HrDTmeKrNIohBDh8NqHRzBGHgTgkoKlg34BeQl1IYQYJDXNLaxvfBNbZiX5aaOZkjlp0N9T5tSF\nEGIQHG4p54+b/47V2UKaNYt/n/8NjK7BHaWDjNSFECLk6jsa+dO2h3CZLVjqxvHf825jdGrf67aE\ngoS6EEKEkMfv5aFdT+DyduApm8xl+ctJjIsfsveX6RchhAiRLo+PP7z/T474juKtG0WSaxyLZowa\n0hok1IUQIgRM0+Te11/naNIu6EhhcdYlXLK8kMT4oY1ZCXUhhAiBVVuOcMS6EYtp4T/nf538jKEd\nofeQOXUhhBig8po2nt65GktCB+fmzAlboIOEuhBCDEiXx8eDL23DMrIEu2Hn6okXhbUeCXUhhBiA\nZ949SF3cXow4NxfmLyIlLjms9cicuhBCBMFv+lmzfy/vVbyPfXQpSbYklo5dGO6yJNSFEOJM6YYS\n/rnvWeo667GPAQODaydcTqItIdylSagLIYavyvZq9tZr8lJGMTFjfLjLweVx8VzJSj6o3ASmgbd+\nFMU5U/jS/PNIsjvCXR7Qz1BXSs0Ffq21XqKUGgc8BviBXVrrW7uf803gZsAD/EprvXJwShZCRLuP\nqrbyatk7VLtqjt83P3cOn5lwGYm2xDPa1urydZS1HOErkz+HxQj+a8RqVy0PbHuEus566Eil8+AU\nCtPH8PUFxdisw+fryT5DXSl1O/BloK37rnsJXIN0rVLqQaXUVcAG4LtAMeAA1iml3tRaewapbiFE\nlNpTt5/H9vwfmBZG2ccxacQE9rVv5/3Kjexp0BSmjsVhd5ARn86M7KnkJuUcf63f9H8iuOs66nm+\nZCU+08eCUecyIaMoqJp2Ve/nod1P4KELT0Uh9tqzuHZeIctm5g2rQIf+jdRLgGuAJ7pvz9Rar+3+\n+2vARQRG7eu01l6gRSl1AJgObA5xvUKIKOU3TV7esJsnDj6OaQH3vlkcbMvgIGCzTSdrYiYt7GVr\n7c7jr3ml9A1yEkfiTBxBTWcNdR0NLBo9n2snXIFhGLx08HV8pg+ALTU7ggr1V3du4ZXqJwEwj05n\n2dhzufTafJIT7SHpO9T6DHWt9fNKqd5X0+u9dmQrkAqkAM297m8D0kJSoRAi6vn9Jn9buZPt5stY\nkt1M4DyuvnYJVfUuDle3oo80cXiPAcYosHowbB7sKW2QXkGVv5rqjirijHgctkRWH13H4aMeao85\naB29nbEpeTR0NrK1dgefnXjlGU3BfKSreOXIy5Bgcm7ClVz3ubk4EoZnmPcI5otSf6+/pwBNQAuB\ncP/0/aeVkeHAZrMGUUL4OJ3ReY3QaO0Lore3aOnL6/Nz5/89x37LBiyJ7czJncl/nP/FE64Q1NzW\nxc6DdZRWtFBa0UxDSydplukk+Ey27Kmg2WVjQlEibWmvc8j+If6sRCzAF6Z/ho+qtvL2wbXUmpVM\nze7fhSo27KpkxYZXsY1pZ5ZzDv+xdPmAex2KfRZMqG9RSi3UWq8BlgOrgE3Ar5RScUAiMAnY1deG\nGhtdQbx9+ETj9RMhevuC6O0tWvo61lzHHzY+istWi8VqsLjgPK7Kv5S6uraTPl+NSkWNSgXyPnF/\n/bkTeOJNzY6D9SSmzcGmPsAb34GvMZuVrzex5PyzeJu1vLN/AzmW0X3WdbiqlV//ax32qSUkWBK5\nftKlA/55h3Kfne7DIZhQ/yGwQillB/YCz2itTaXU/cA6AtMzd2it3cEUK4SIfg0tnbywaScfeV7G\niOskwZXHbQuuZ+aECUEFX2ZaArddN52DFS2MHOHgWMck3j6yhvqqcWw8UMOk/ImkxCWzrXYnN0y8\nGqvl1DMEHq+Ph1buwTJaY1i9XD3hymFzuGJ/GKZphu3Na2tbw/fmQYiW0dGnRWtfEL29RUpfbp+b\nvQ37URkTSLAFLhShjzRy/8r1mEUbMOxulP1cvjX3KuLjrCHvq665g58/sgmv38+8i+rYWLuRfzv7\nG5yVOfGUr3lqdQlv7t1C/KSPGJMymv+c9d0BHQrZI8Qj9VNeF09OPhJCDIr6jgb+tvNxjrZVkBGf\nzg3qalx1qTy65RUsE8owLH4+O+FqFo+ZP2g1ZKUl8vllE3h45V6qDqVDCmyq3nrKUN9f3sQbH5WQ\nePYuDMPCFyddF5JAH0oS6kKIoHX53BxuKceZmEl6fBqGYeD2udGNJTyx9ynaPS7GpxVxqLmMv+x4\nDNNnxTrSR5I1hS9MvoYZzqmDXuP8qSNZu6MSvbeRkfMz2FS9lQvGLmR0cu4nnlfd4OIvL+7CXrAX\n09bJ5YUXMyal7/n34UZCXQhxShurtrC9dheFafmMTy9kTPLo4/PRNa5a/rLj78fP+ow3EsC00kV7\n4MWmga1yOrs3jcZMyMWevwdrUisLnEu4dvIy4qxDc2igYRh8+aKJ/PzRTXSWTcKf9wFP73+R2875\n1vEjbKobXPzmn1tojTtMXGYlhaljuXDs4iGpL9Qk1IUQJ9Xc1cK/9HO4fW621QYOZkuJS2ZOTjEJ\n/gzeqngDt9lJhq+IxlYXHQnNGBYv/o4RmJ1JGE1jSDayGTk6noljxjK5YCGFo5JJsA/9cd6jnclc\nNHsMr31okj0yjwNNh9hSs51zss9GH25kxSt7aPbVkzx+DxaLnRsn33DaL1OHMwl1IcRJrSx9E7fP\nzVVFy8lISGd/40G21u7knfI1AJh+A0/ZVCrq8hiRGs+8KSOZVpRJenIcqUlxJMQNr3i58rxCtpXU\nUbW7iPhpFTy+8wX+VdJEQ5MfbG4yZ+7EZXr42qTPk+1whrvcoA2vn7oQYlioaKvi/YpNjIjLpKsi\nn46kBIpsuWzbk0OXUU76qCZmO2dSeF4BmWkJFOamYjFOeUDGsBAfZ+XnX5vDhj1VvHCghs6MfXgK\n36YQhSWliYrOFpYXLGPWyHPCXeqASKgLIU7w7P5XMDGp3JnP881lx+83gEvnzeWqBYXDbiGr/rDb\nLJw/fRRzJn+JJ/esZEfTNqq8O6ETznFO49LCZeEuccAk1IUQn7Du0E72Ne3H1zKCUfGFXHF1IW6P\nj1aXh/F5aYwfHfnLOsXb4rhx+jV4fJeztXYnVe01XFKwNOIOXzwZCXUhxHGmafLC3nchERZkLeZz\nV8yOyBF5f9mtduaMLA53GSEloS6EOG7LwUpc8ceI86XyxfPmnrColhj+ovcjWAhxRkzT5Okt6zAs\nfs7NLZZAj1AS6kIIADbrWpriDgFw4bhzw1yNCJaEuhACv9/k2Q92YU1tYGzSWDITR4S7JBEkCXUh\nBDsO1lNvCYzSF+TNDnM1YiAk1IUQrN52FGtmBVbDyjnZ08JdjhgACXUhYlxdUwd76jQWRxvTsibj\niKALQogTSagLEeNWbTuCLX8vBgbLCy4IdzligII6Tr37UnYPARMAN3Ab0AY8RuDC1Lu01reGqEYh\nxCDx+vysrVyHJcfF+aPOIy9lVLhLEgMU7Ej9m0CH1no+cDPwCHAvgWuTLgIsSqmrQlSjEGKQvLfn\nAD7nAeJMB1eOvzjc5YgQCDbUpwCvAWit9wOjgaVa67Xdj78GRP7KOEJEudfKX8ew+LksfzmJtoRw\nlyNCINhQ3wZcDqCUOhdwAr2/XWkFIn/VHyGi2Fqt6Ug4RoLHyQXj5oS7HBEiwa798ghwllJqDbAe\n2A9k9Xo8BWjqayMZGQ5stsi6uojTmRLuEgZFtPYF0dvbQPry+U1e0u9CMlw77SKys1NDV9gARev+\ngqHpLdhQnwOs0lr/QCk1E5gL7FdKLdJavwcsB1b1tZHGRleQbx8eTmcKtbWt4S4j5KK1L4je3gba\n1+rtZbQnlhHnT2LuyKnD5mcUrfsLQtvb6T4cgg11DTyplLoD6CDwxakFWNF9ZMxe4Jkgty2EGERu\nj48Xdq/BGOlj0eh5EXstTnFyQYW61roBuPAkDy0eUDVCiEH39uZy3OmHsGJlWeH8cJcjQkxOPhIi\nhvhNk7f3b8WS4KLYOZ2UuORwlyRCTEJdiBiyqeQwHc6tACwrOD/M1YjBIKEuRIxo87TzZNk/sMR3\nMj9rEWNT8sJdkhgEEupCxAC3z82ftjxMl7WZ+KbxfH7q8nCXJAaJhLoQMWBj1RbK24/ircvl4rxL\nsFjkv360kj0rRAzYXrsbALNyIudNyw1zNWIwSagLEeU6vZ3sayjB70qhuCCf5ER7uEsSg0hCXYgo\nt7N2H358+BuzuWj22HCXIwaZhLoQUe71fZsAKB45jaJRw2eNFzE4gl0mQAgRAcqqmqn0lGExE7hx\nsVxQOhbISF2IKGWaJg+tXodh8zB5xFkkxstceiyQUBciSh2pbqPWLANgSdHM8BYjhoyEuhBRatuB\nWizpNdiNOCZkFIW7HDFEJNSFiFJbjpRiSejgrBETsVnk67NYIaEuRBRqbndT4T4MwHTnWWGuRgwl\nCXUhotCOg3VY0+oAOCtzYpirEUMpqN/JlFIG8BCgAB+BKx/5gMcAP7BLa31riGoUQpyhbSXVWFIb\ncCZkkx4v14CPJcGO1C8CkrTWC4D/Ae4C7gXu0FovAixKqatCVKMQ4gx4vH721B/EsPiZnj0p3OWI\nIRZsqHcCad0j9jTAAxRrrdd2P/4asCwE9QkhztD+8ib8SdUATB6hwlyNGGrBfiW+DkgE9gGZwBVA\n78uotBIIeyHEENteUoclrQ6bYWdcWkG4yxFDLNiR+n8C67XWCjgbeByI6/V4CtA0wNqEEGfI5/ez\n6eBhLIntTMwowm6Vs0hjTbAj9WSgufvvTd3b2aqUWqS1fg9YDqzqayMZGQ5sNmuQJYSH05kS7hIG\nRbT2BdHb28n62ri7ijZ7BXHAnPzpEdl7JNbcX0PRW7Chfg/wqFJqbfc2fgxsBh5SStmBvcAzfW2k\nsdEV5NuHh9OZQm1ta7jLCLlo7Quit7dT9fXymoNY02sBGBtfEHG9R+v+gtD2droPh6BCXWvdBFxz\nkocWB7M9IcTANbZ2sbOqhLjJNYxNGU12Yla4SxJhICcfCREl1u2owJa/B4DrJlyFYRhhrkiEg4S6\nEFHAb5qsPvwBlqQWzsk6m3HpBeEuSYSJhLoQUWB7aSWdmXuwmDauU5eHuxwRRhLqQkSB5/a8g2F3\nM8+5QJYFiHES6kJEuLKqZuqsBzD8Nj4z+YJwlyPCTEJdiAj39KYPscR3olInk2CLD3c5Iswk1IWI\nYFUNLg527gLg0okLwlyNGA4k1IWIYC9t2I8lo4Y02wiK0vLDXY4YBiTUhYhQDS2dbK7ehmHxs3js\nuXJcugAk1IWIWO9uO4Yl6ygGBnNzZ4a7HDFMSKgLEYE8Xh/v7tuLJamFKSMmkRYfvYtgiTMjoS5E\nBFq3vYKulDIA5o+eE95ixLAioS5EBHpp3X6smRUk21OYmimXrBMfk1AXIsIcqmihtH0/hs3L/FGz\nsFoi65oEYnBJqAsRYVZtOYrVeRSAebmzw1yNGG4k1IWIIM3tbjYeOoQ1tYEJ6UVkO2TNdPFJEupC\nRJB3tx6DEYFR+vxR8gWpOFFQVz5SSn0F+CpgAokELj59PvAHwA/s0lrfGqIahRAEDmNctfUw9onH\ncNgTmeGcFu6SxDAU1Ehda/13rfUSrfVSAtcm/R5wJ3CH1noRYFFKXRXCOoWIeR/uqcGVfAjsXVw0\nfiFxVnu4SxLD0ICmX5RSs4DJWuuHgJla67XdD70GLBtocUKIANM0eeOjUuyjDhFviecKJf+9xMkN\ndE79J8DPT3J/KyAr9QsRIvuONFFt2Ythd7N07AJS4pPDXZIYpoIOdaVUGjBRa72m+y5/r4dTgKaB\nFCaECPD7TZ5dq7HlBkbpS8csDHdJYhgL6ovSbguBd3rd3qqUWtgd8suBVX1tICPDgc0WWSdOOJ3R\nucZGtPYFkd/bQ2+s52jCeqx2D1dNvoL8UdlA5Pd1KtHaFwxNbwMJdQUc6nX7h8AKpZQd2As809cG\nGhtdA3j7oed0plBb2xruMkIuWvuCyO6tsbOJB7b+nYqOY1hHQHZiNnNGzKa2tjWi+zqdaO0LQtvb\n6T4cgg51rfVvP3X7ALA42O0JIT5p5aG3qOg4hq8pi0snLOTyqXOwGHJqiTi9gYzUhRCDpN7VwobK\nLfi7HMxNvIwrp00Jd0kiQsjHvhDDTHunh9+++SKm4SPLM4kvXCCrMIr+k1AXYhhxe3z8+h+baU7Y\nj8W086NLriQ+LrIOJhDhJaEuxDDy1kflVPlLMOK6WDhmLsnxjnCXJCKMhLoQw0SLy82rG8qIyz2C\ngcGSMeeFuyQRgSTUhRgmXllfRldiJTiamJ41mazEzHCXJCKQhLoQw0B1o4vV24+QUKixGlauHHdJ\nuEsSEUpCXYgwM02Tp1aVYOSUYNpdLB1zPiOTcsJdlohQEupChNnaHZVsKz+MfVQZ6fFpXFJwQbhL\nEhFMQl2IIXa4pRy3zwPAsbp2/vn2XhKKdoPh57oJV5Jgiw9zhSKSyRmlQgyh9ys28Y99TzM6OZeb\nzvoyD754AAq2QHIDM5zTmOGcGu4SRYSTUBdiiLS523mhZCUWw8Kxtkp+teE+PBkpWNPrmJQxga9O\n/hyGYYS7TBHhZPpFiCHy4sFXafe6mJlyPsaxafgsbqzpdRSm5nPz9K9gl8vTiRCQkboQg6zL7eO3\nL6+iImMTflcyazbFY7OM4QI1EVtGDZcWLiPeGhfuMkWUkFAXYhCZpskjr+/iaPwHWIB873xGzxjD\n4nNGMyZbLkknQk9CXYhB9NamcrZ1vYMts435uXP54lK5YLQYXDKnLsQg2V/exLP6dWyZVeQn53O9\nuircJYkYEPRIXSn1Y+DK7m38CVgPPEbgAtS7tNa3hqJAISJRU1sXf1r1OrYxJaTa0rhlxlewW+QX\nYzH4ghqpK6UWAfO01vOBJcA44F7gDq31IsCilAxLRGzy+vzc88bLePO2YMPOd4u/TkqczJ+LoRHs\n9MvFwC6l1AvAS91/irXWa7sffw2QyUMRk/6w+gWaMjZiI45/n/ktRiWPDHdJIoYE+/tgFjAWuBwo\nIhDqvT8gWoG0gZUmRORw+9y8XbKJd498QLu1BsObwA9mf4uCtNHhLk3EmGBDvR7Yq7X2AvuVUp1A\nXq/HU4D8zO7/AAAU/ElEQVSmvjaSkeHAZousS3U5nSnhLmFQRGtfMPi9fVRawu82/BmfpROsYO/I\n4UdLv8b0sfmD+r7Rus+itS8Ymt6CDfV1wPeA3yulRgFJwDtKqUVa6/eA5cCqvjbS2OgK8u3Dw+lM\noba2NdxlhFy09gWD25vX5+f5DXtZ3fIkRnwnKW2TuGziQs5T47AYxqD+TKN1n0VrXxDa3k734RBU\nqGutVyqlzldKbQQM4BagDHhIKWUH9gLPBLNtIYYzt89DlasauhJ5+JX91GS+gyWpk5mpC7hp6ZXh\nLk+I4A9p1Fr/+CR3Lw6+FCGGt111e3ly/ws0dDYCYOYZWCwmc7JncuOUK8JcnRABcuCsEH3o9Hby\n+O6n2V6/E8M08NWPAquXjEwvU3KK+Ly6VlZXFMOGhLoQgMvTgcOeeML9ze2d3LNhBY1GOf62NNyl\nU0kxRnDzlVOYXDAiDJUKcXoS6iLmvXTwdd468i5fPut65owsBgILcb25qZyXDr4GOeVY251ckH4N\n067NojA3FZtVVtgQw5OEuohp+xoO8MbhwIFa/9z3DCOTshmbksdL68t4Ze964sYdJMmSzk8vvIXU\nBDkrVAx/MtwQMavN084Te5/CYli4vPBivH4ff9vxOM9s/IhXK14krmgHCdYEfjD7GxLoImLISF3E\nJNM0+de+52jqaubKoku4uGApJiYrS99kdddT2LIgJzGHL02+jpFJ2eEuV4h+k1AXMams5QjbandS\nmFrA/Jzz2FvWwKb3UvEmjsaW3MIV45dx0cS5WAz5ZVZEFgl1EZM+LN8NwL7N6dz29vrj98+dvIzP\nnFeEM/3EI2GEiAQS6iLm+E2TD4/sgQQ4K3MccVmJJMTZWDYrj8Lc1HCXJ8SASKiLmPPe9nK64upI\n8KXzH9fNDXc5QoSUTBiKmNLqcvPsps0YFj/FuZPCXY4QISehLmLKM+8exB1fC8C0nIlhrkaI0JNQ\nFzHD1enh/V1VJGY2YWAwPr0o3CUJEXIS6iJm7DzUgM/04k9sZHRyLkl2R7hLEiLkJNRFzNhWUocl\nuQk/PiZmjAt3OUIMCgl1ERO8Pj87DtaT5GwBkFAXUUtCXcSEA+VNdHR5cWQ2Y2AwLq0w3CUJMSiC\nPk5dKbUZaO6+WQrcBTwG+IFdWutbB1ydECGytaQOS0YVLUYVRWkFJ107XYhoENRIXSkVD6C1Xtr9\n5+vAvcAdWutFgEUpdVUI6xQiaKZpsqWsjPiiXcRZ7HxefSbcJQkxaIIdqZ8NJCml3gCswH8BxVrr\ntd2PvwZcCLw48BKFGJiymibacz7EYvXyOXUDo5JHhrskIQZNsHPqLuAerfXFwC3AP4DeF2lsBdIG\nWJsQIfHsvtexJLUyPnEac3NnhrscIQZVsCP1/UAJgNb6gFKqHiju9XgK0NTXRjIyHNhs1iBLCA+n\nMyXcJQyKaO2rsaWT0o59mFYbP1xyIyNSoudiF9G6z6K1Lxia3oIN9a8B04FblVKjgFTgTaXUIq31\ne8ByYFVfG2lsdAX59uHhdKZQW9sa7jJCLlr7Anjg9Q8hroPRtvH4Ok1qO6Ojz2jdZ9HaF4S2t9N9\nOAQb6g8Djyil1gAm8FWgHnhIKWUH9gLPBLltIUJi64FatlXtwj4Wlow7J9zlCDEkggp1rbUXuPEk\nDy0eUDVChIir08sTb2iseXUATMmSFRlFbJCTj0RUem/7MZpcLqypjYxJGU1afPTO0wrRm4S6iDqm\nafL+zirs6Q2Y+JmSKaN0ETsk1EXUOVLdxrG6drLGBL6UklAXsUQuZyeizvu7qgATj6OKZGsSBalj\nwl2SEEMmpkPd5XERb43HaomsY+XFyfn8Pp498ArrW0pxTDBo97Vy3uhZWAz5hVTEjpj9197maeen\n79/FG4f7PJxeRIhV5Wt579h6zPQKzIxjAMzNk0MZRWyJ2ZF6Q0cjXT43x9qqwl2KCIFaVz0rS9/C\n5k+gbddsbrl6KmOyk5mSV0hdXVu4yxNiyMRsqLd7A2ezujyRdVarOJHP7+fh7U/i8Xvwlk4hx+Fk\nZn4BhmFgGEbfGxAiisTs9EtPmPeEu4hMpmnyq1deoLyjDF+Tk2TPWD6/bIKEuYhZsTtS7wl1GalH\ntHW7jlEVtxWL38ZN0z/LrGvysUigixgWsyP1dk8HINMvkay908PTW9ZixHUxd+Qs5owvkEAXMS9m\nQ93VPe3i9nvw+DxhrkYE47k1h/BkHALgkqJFYa5GiOEhZkO997SLzKtHntLKFt7Tu7EkNzM18yyc\njsxwlyTEsCChDri6p2JEZDBNk3+9cwBrzmEAloxZEOaKhBg+YjbUXb1G5/JlaWTZeqCOkupqbJlV\njEzKQWWMD3dJQgwbMRvqMv0Smbw+P0+/exB7TjkYJovzzpPDF4XoRUIdOQImkqzdXkF1YysJuUdJ\nsjmYM7K47xcJEUMGdJy6Uiob+AhYBviAxwA/sEtrfeuAqxskpmni8n48jy7TL5Gho8vLi+tKiR9Z\nhdfo4oLRS4m3xoW7LCGGlaBH6kopG/AXoCcR7wXu0FovAixKqatCUN+g6PR14jf9JNkdgIR6pFj5\nwWFaXG5Sxh7FalhZmDcv3CUJMewMZPrlt8CDQAVgAMVa67Xdj71GYPQ+LPWceJSdmAV88ktTMTzV\nNLp4c9MR0nKbaTcbmZlzNunxaeEuS4hhJ6hQV0p9FajRWr9FINA/va1WYNj+j+uZQ8/qDvV2OaRx\n2HtyVQlen0nmuEoAlo45P8wVCTE8BTun/jXAr5S6EDgbeBxw9no8BWjqayMZGQ5stqG/QEWFzwQg\nPyuXTdXgMbpwOvt3YeL+Pi/SDOe+tu2vYduRw2RPq6LSfZgp2RMpLur/JeqGc28DIX1FnqHoLahQ\n7543B0AptQr4NnCPUmqh1noNsBzo8+oTjY3hmfaoqKsHwOqJI9GWQJOrldra1j5f53Sm9Ot5kWa4\n9lXT3MLLuz5kW90O4qdX02pARnw6l429uN/1DtfeBkr6ijyh7O10Hw6hXKXxh8AKpZQd2As8E8Jt\nh1TP9EuSPRGHzSFflA4jpmmy+sB23jy0nhZbOYbFDymQZmRz3eRlzHBOk8sPCnEaAw51rfXSXjcX\nD3R7Q6EnxB12B0n2RCrba8Jckejo8vLKzo2sq38Hr70V4sDmSWFi4mQumTSX8Zl54S5RiIgQk+up\n95xBmmRz4LA58Pg9uH0e4qz2MFcWexpbu3jmPc2W9jVYnEcwrQapXQVcXLSARROmYrHE7PlxQgQl\nJkO9ZwGvJLvj+LHqLq+LOOuwPWAnKrR52tnfeJDK9mqq2qspb6inprkN7J1YnF0kGxl8cfLnmD6q\nMNylChGxYjLU2z3tQM/0y8cnIMlxz6Hn9/tZV7qbNeUbqPQdBMP/iceNRAvx1njOyzufK4suwS6/\nLQkxIDEa6h1YDAsJ1ngcPSN1+bI0pOqaOnhn5wHWN76NL7kKAH9nEomufOI8GVg9KRRkZvPZxRNI\nT44Pc7VCRI+YDHWX10WSzYFhGB+P1L1yAlIwunxuypqP0OZpp93jornDxe6yekpr6rFmH8FI9pHk\nHcmC7EUsmTiNFIes1SLEYIrJUG/3uEiyJwGBL0sD97WHs6SI0uHuYl3ZdjbXbOdY1yH8hu+TT7CC\nLRfijASuGX8N5+fNluVxhRgiMRfqPSs0ZjsCSwQ47ImAXP2oL5X17by8eTe6dScdyaUYdjcA/i4H\nvsYczK4E8NpJtCcwb3Iu08c5KUwbc/znK4QYGjEX6j0rNDq6R+g9I3Y5AeljpmlS2V5NXXsTJZV1\nbDtWSq1ZisXRBhlg8ceR45vCtPSzmTqxgNysJGzdhx7abRYsFhmVCxEuMRfq7b0OZwRIsiV23y+h\n7vP72XHkKC+UvUSd/8jHDySCxbSQF1/EksLZFOdMl2P6hRimYjDUA3PnPaHu6HWceqxqbuvi+bUH\nWHVkPb5sjWH14WsZQbqZy6iMdFRuLguLppNoTwh3qUKIPsRcqPfMnfdMvzhieKRe1djGK5v2srlm\nO4azDCPXjc2MZ376RVw+7zySE+VIFSEiTcyF+vElArpH6FaLlURbQtSGumma1Hc2sq+uFF1bRkNn\nM22eVlrcLXQZbRg2E8sosBPPwrzFXFSwiOS4pHCXLYQIUuyFeq8VGns4bI5PXLN0uOj0drG/sYT8\n1DGkxaf2+3XNbV1sKS1nc90Wyr178VpPPFzT9Mdh92eQl5bNIjWV6WnTSbDJ9IoQkS7mQt3Va4XG\nHsNhpUa3z0NDZwPuLiu7DzZz2LOHEs9mOnwurIaV4uzpLB5zHvkpYz5xzPfRliqONddS195MZXMD\nB+qP0eKrx3C0YBhgYsXakssI60jykvPITckiJzmd7PRkxuYkYxhGVK9hLUSsiblQ771CY48ke1JY\nVmrs8vjQRxrZ13CADa1v0sUnR9Smz4qvbixmagObqreyqXorNncaia4CfF5wJZVCYssnN5oQOFIl\n3ZrDjBHnsLRoDpnJyUPWkxAivGIv1D2fnFOHj78sHYqVGk3TpLGtgze2HOD9/WW4k8ux5RzB9Bv4\nGnJJdlhJTTVwxuWS5ppEjdVHY3UnzRzDnVaKJ60Gb/r27o0ZJHaOJtXIJtmexIjEVGYVFnHWyNFy\nIQkhYlRQoa6UsgArAAX4CVzOrgt4rPv2Lq31rSGqMaROPv0yOCs1mqZJQ0sXm8r38V7tW7T468Aw\nP37CuMAOSLOO4MLcK5g6qwhn+unPwGxzt7O5ZjumaTIz52xS4mQULoT4WLAj9SsAU2u9QCm1CLgL\nMIA7tNZrlVIPKqWu0lq/2J+NuT0+rBYDt9lFq7uNeGsCyXbH8dGm1+fH7XPT6mmh1dNOjsPZZ5j5\nTT8uTwfxtnjslkCbpmnS5nEdX6Gxx+lWavSbfuo66qlx1TErdfIJj/v8PqpdtVS2Vwduu+NobDQ5\ncKyeg9X1uBIPY3MeC2yrPQ0LVuLtdpzJaUzIHkluipM5OcX9XnI2OS6JRXnz+/VcIUTsCfbC0y8q\npV7uvpkPNALLtNZru+97DbgQOG2o//ivH9Dc3oVZsBFLSiOG9VMLQ/mtgIlpmBi9R7hAhs1Jdtwo\nfF6DLo+PDm8nHb52unDht3bit3QdHxVb/HbwW/FbA/dZfAn86bmdpCXFkZoUR43NA8DjO15krHkO\nRkcGFd6DNFrL8MTVY1q8ABjbrYzwjsdJIfXeSpqtR3HbGz45+u4RD4wN/IBTjCyWOC9h1syJjEiN\nl8WthBCDJug5da21Xyn1KHAN8FkCId6jFehzHqPL48OZnogr0Y7Xn4LZmYC3yw4WL1jdYPVitVix\nWSx43Ra6XHGYXjuWpBYaUupp9NZ+vDFr4I/ps2J64jE96eCxg9WL3+7GsPowutLAG4+7PoetdXUf\nv9ZuIa4ok4a0Khp4LRDIPQP5zmToSMN0J+BPP0Z9vKYeDXYwTQOjIxU6UqEzBUeCjeQUkwSHH2dq\nMiNTU3E6sijOni5z3EKIITGgL0q11l9TSv0I2AT0ngxOAZr6ev3v/21B99/m9vlepmlS09SBPtJE\nl9uHz/TSTiOOBBspiXGkOxIZlTaCdEcSJtDm8tDW4cGRYCfFYcdm/fhal37TpL3DQ3O7m5buP66u\n6XRZGznQtZV2fwvTss7i3FEzyHJkftxUWjyvbF/HoebDTMgoZEbOJJLj5UQdIcTwYZjmSaYO+qCU\n+jKQp7W+WymVCmwDDgB3aa3fU0o9CKzSWj8d2nKFEEKcTrChnkjgSJeRBEb7dwP7gIcAO7AX+KbW\n+sw3LoQQImhBhboQQojhydL3U4QQQkQKCXUhhIgiEupCCBFFJNSFECKKxNyCXqeilJoL/FprvUQp\ndTbwIOAlcKjmt7XWHqXUcuBOwAQ2aa1v6/X6ScAGIFtr7R76Dk4u2L6UUunA40AG0E7gaKby8HRx\nor76AiYDfyDQkwGcC1yltX6z+/XDcn9B8L0BG4ngfdb9b/EW4CYCa0jdrbV+odfrI3afnaq3wfh/\nJiN1QCl1O4EFynrOI10BfF9rvRCoAL6jlEoG/he4TGs9HzimlMrqfn0K8Fugc8iLP40B9nUH8L7W\n+nzgHuCPQ97AKfSjr1u11tu11ku01kuBPwPP9Ar0Ybm/IOjenu7uLZL32XeUUknA7QQ+pC4m8MHV\n8/pI3men6y3k+0xCPaCEwHIHPcZorT/s/vt6YDEwH9gJ3KuUWgNUaq171hr4G/ATYLhdE28gfU0h\nsIZPz3MXDUnF/dNXXwt7HlBKOYBfALf1ev5w3V8wsN4ieZ8tIvCbh0ngjPRkoPdiUJG8z07XW8j3\nmYQ6oLV+nsCvSj0OKqXO7/77FYADyCIQgrcDy4HvK6XGK6V+Bryitd5J4FfhYWMgfQFbgSu7n3sV\nn1wGIqz60VfvtRu+DjyltW4AGM77C4LurbH7dkTvM621C/g/YA/wEXA/gFLq50T4PjtVbwzCPpNQ\nP7mbgDuUUm8B1UBd959NWutarXU7sAY4B/gi8HWl1GoCZ9i+Gaaa+6O/fc0Afg0UKqXeJbAS57CZ\nmz2J3n3VEOipxxcJnOnc40tEzv6CM+stUvdZNVCnlJoHzCNQ+1jgGqXUbCLr/xj0v7dZDMI+k1A/\nucuBL2itLyQwkn2TwCfqVKXUCKWUjcDc2G6t9USt9VKt9RKgik+uVjnc9LevPQR+Dfyb1noxgV8v\n14Wn5H7p3Vcm3f/pu9clitNaH+t5otZ6QgTtLziD3ojcfdbzbzEZcGmtPd1fhDYBaRG+z07XWzqD\nsM/k6JeT2w+8o5TqJLAC5eNaa1Mp9RMCO8gEntRa7/nU63qORhiu+t2XUqoL+LtSygAagK+Freq+\nndBX9/0TgbLTvG647y84s972EcH7rPvf4oVKqQ8JTGes01q//anXReQ+O1VvSqlxhHifydovQggR\nRWT6RQghooiEuhBCRBEJdSGEiCIS6kIIEUUk1IUQIopIqAshRBSRUBfiU5RSqUqp58NdhxDBkFAX\n4kQjgLPDXYQQwZBQF+JE9wGjlFLPhrsQIc6UhLoQJ/oeUKG1vjbchQhxpiTUhRAiikioCyFEFJFQ\nF+JEXmQFUxGhJNSFOFE1cEQp9U64CxHiTMnSu0IIEUVkpC6EEFFEQl0IIaKIhLoQQkQRCXUhhIgi\nEupCCBFFJNSFECKKSKgLIUQUkVAXQogo8v8DcePLEDkU/FgAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x116be5748>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pred = res.get_prediction()\n",
"ax = data.wpi.plot(label='observed')\n",
"pred.predicted_mean[1:].plot(ax=ax, label='Forecast')\n",
"plt.legend(loc=2)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x11752fc88>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAERCAYAAACEmDeEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8m+W5//GP9vK25QzbWY5zZyckYSRAEiCBhpUCBVo6\noYWW5nTQU85p6eR3SlvKaIGWDkYptJRCmGVvkrAhe90ZdoYdJ96WbW3p+f0hxU2IkziObVny9X69\neGFLjx5dVx75q1v3M2QyDAMhhBCZwZzqAoQQQvQeCXUhhMggEupCCJFBJNSFECKDSKgLIUQGkVAX\nQogMYu3ugkqpk4Ffa63PUEpNB+4EokAI+JLWul4pdTVwDRABbtJaP9cXRQshhOhat0bqSqnrgXsA\nR/Km3wFLtNZnAk8C/6uUGgJ8C5gNfAr4lVLK1vslCyGEOJzuTr9sAy464PfLtdbrkj9bgSBwErBC\nax3VWvuArcDUXqtUCCHEUXUr1LXWT5KYatn/+z4ApdQcYAnwWyAHaD3gYe1Abq9VKoQQ4qh6vKNU\nKXU5cDdwrta6EfCRCPb9soGW4ytPCCHEsej2jtIDKaW+QGKH6Hyt9f7g/gD4hVLKDriA8cD6I60n\nGo0ZVqulJyUIIcRgZjrcHccc6kopM3AHsBN4UillAG9prW9USt0JrEg+4Q1a6/CR1tXc7D/Wp08p\nrzeb+vq2VJfR6zK1L8jc3jK1L8jc3nqzL683+7D3dTvUtdY7gTnJXwsPs8x9wH3HUpwQQojeIycf\nCSFEBpFQF0KINBIKx454v4S6EEKkCcMw+O1ja464jIS6EEKkifVVTWzZfeQjxSXUhRAiDRiGwdMr\nqo66nIS6EEKkgXWVTVTu8TFznPeIy0mod9O3vvV1WlpSe4Lsm2++xv33/yWlNQgh+l9ilF4JwIWn\njT7ishLq3WQyHfYErn41UOoQQvSflVvqqaptY6byUlacdcRle3SZgP7y6Ovb+HBzXa+u88TxxVx2\n5tgjLhONRvnVr25kz54a4nGDyy67AsMwuOmmm6iu3oPT6eJHP/oZhgE/+9kPMQyDcDjM97//Q8aO\nreDxx//FK6+8hMlkYsGCs7nkksv55S9vpLW1BZ/Px4gRI5k+fQaLFp1PU1Mj11//Xe677yH+/Oc/\nsHbtauLxGJdf/nnmzz+LdevWcOedt5GdnYPNZmP8+Im9+u8hhBjYquvbuf/5TVgtJj59lFE6DPBQ\nT5Wnn36CvLwCfvKT/8Pv93PVVV/A4bBz0UX/RUXFFJ56aikPPvhXZs06kdzcPH784xupqqokGAyw\nY0cVr732Cn/8430YhsF11y3hxBNPAWDmzJO47LLPsWNHFb/97W9YtOh8Xnzxec4770Lee+8d9uyp\n4Q9/uIdwOMzXv/4VZs06mdtuu5lf/OJmSkvL+POf/5DifxkhRH9q7Qhzx2NrCYRiXHPhREq8Rx6l\nwwAP9cvOHHvUUXVf2LmzihNPPBkAt9vN6NGj+fDD9znxxBNpaQkyadIU3nnnbb71revYvXs3P/jB\n97BabXzpS1dRWbmdvXtr+c53rsUwDNrb26ip2Q3AiBEjARg1ajSxWIy9e/fy+uuvcMcdd/P000+g\n9Wa+/e1vYBgGsViM2to9NDU1UlpaBsD06TPYuPGI10gTQqSx9kCE+57dSE1DBwU5Tlo7wjT6gnz6\n9NGcMnFot9Yhc+pdGDlyNKtXrwLA7+9g+/ZtDBs2nDVrEgf9r169ivLysaxc+RGFhUXcfvvv+dKX\nruKee+5m5MhRjBlTzp13/om77vozn/rUeZSXVwBgNv/nn/v88xfzxz/ewejRY/B4shgxYhQzZ87i\nzjv/xO9+dzdnnrmQkpJSioq8VFUldpBs2LAOIURmamkPcfPDK1mzvZFgOMbW3S3sa/Jz+tRhXDBn\nVLfXM6BH6qmyePHF3HzzL/jmN79GOBzmqquu4fnn/80zzzzDli3byM7O4cc//jnRaJSf/ewGnnpq\nKfF4nCuvvJry8rHMmHEi1177VcLhMJMmTaao6NBDkM444yzuuOM2br75dgBOO20uq1Z9zJIlVxMI\nBJg7dz5ut5v//d8f86tf/T/cbje5ubmMGjWmv/85hBC9aPPOZsxmE6OHZWOzWohE41TuaeWvz2+m\nriXAWTNK+dzCCuJxA38wSo7HfkzrNxmG0UelH119fVvqnrwH5JKg6SdTe8vUviBze/N6s3n8Vc0D\nL2wGwGoxMaTAzb6mANFYHIDz54ziotNHH/UoN683u/eupy6EEOLYfbBhL397cTNZLhsnTxzCtppW\nauo7GF7kpqI0jyljCpla3uVVzY+JhLoQQvSxbTWt3PrIamxWM9+5dCrlw/vu65tlR6kQQvSh1vYQ\nv398LdFYnG9+enKfBjpIqAshRJ+Jxw3+8u+N+PwRrjx/ElPLi/r8OSXUhRCijzz37g427Wxm+tgi\nFs/tnyPXZE5dCCF6kWEYVNd3sGprPU+vqKIgx8FV503ot+s2Sah3Ye/eWr785c+i1AQMw8BkMjFj\nxiyuv/66lNSzbNmbTJo0mcLCvv/oJoTouW01rdzz7w3UtwQBsNvMfP3CSWS5bP1Wg4T6YYwenTgr\ndCB47LF/MmrUDRLqQgxgH+t6/vLvDcRiBidPHMLU8kKmjCns10CHAR7qT2x7llV1vXtq/AnFU7h4\n7PlHXa6rk7Juvvlm3nvvA0wmEwsXnsNnPvPZg66+eMstd/CPf/ztkCstbtiwnrvuuh3DMPB6vfz0\np79g48b1/PWv92AYBoGAn5/97CaKi4fw05/+gI6ODoLBINdc802i0Qhbt27hF7/4GXfffS9W64De\nZEIMSq9+tJt/vroVu83Cks9M6ZXjzXtKEuIwduyo7Ly4lslkYtGi86mpqeEvf3mAaDTKkiVXM2PG\nLOA/V1987713qK3dc8iVFm+99ZfceOOvGDFiJM899ww7d1axY0clP/3p/1FYWMRDD/2VN954ldNO\nm0drayu33XYXzc1N7N69i9mzT2PcOMX1198ggS7EALP/K+aeeXsHuR473710GiOHZqe0pm6nhFLq\nZODXWuszlFLlwANAHFivtV6SXOZq4BogAtyktX7ueIq7eOz53RpV94VPTr88/PBDzJw5EwCr1crE\niZOpqkp8X+D+qy9WVm5j8+ZNXVxpsalzmfPOuxCAffv28tvf3oLb7aa+vo6pU6czevQYLrzwIn7+\n8xuIRmNceunlQOKFk8rLOQghDhU3DB55dSuvflyNN8/J9z97At48V6rL6l6oK6WuB74ItCdvuh24\nQWu9XCn1R6XUYuA94FvADMANrFBKvay1jvRB3X3ukyE6atRoXnvtBc4992Ki0Sjr16/h3HPP5/33\n3+m8+uL+Ky1ef/0NxGIxHnror51XWqypqaakpJSHH36I0tIybrnllzz66NO4XC5uuunnGIZBZeU2\n/H4/v/nN72hsbODaa7/K7NmnYTabJdSFGCDaAxE+2lzH2+tr2V7jo6TIw39/djp5WY5UlwZ0f6S+\nDbgIeCj5+0yt9fLkzy8AZ5MYta/QWkcBn1JqKzAV+LgX6+03nzz8aM6c09i8eS3f+MZVRKNRzjxz\nIRUV6qBlDnelxeuv/yG//OWNmM1mCguLuOyyz3HOOeeyZMnXKCoqZuTIUTQ01FNWNpL777+HN954\nFcMw+NrXrgVg8uSp/OIXP+X22/9AdnZqP9oJ0dda2kK89nE1lXt85GbZKcxx4s1zUer1kJ/tSMlX\nOu5r9rNqSwNrtzewtbqVWDwxyJoyppCrL5jY7ztDj6TbV2lUSo0E/qm1nqOUqtFalyRvPwO4EngR\nmKK1/mHy9r8Bf9Nav364dcpVGgeGTO0LMre3TOyrur6dpW9uZ31VE/F419HgdlgpzndRmOMkN8tO\nJBonEIpis1qYM2UoE0bmYwJ217WzZnsjMyqKuvVtQYcTCsd4Ylklr360m/0VjRqazUkThnDShGIK\ncpzdXldvbrO+uEpj/ICfs4EWwAfkdHH7YeXnu7FaLT0sITW83swcKWdqX5C5vWVKX5FojEdf3crS\n17cQjRlUlOUxf0YpJ6hiOgIR6pr91NR3sLPWx45aHzUNHezYe2g4vrthL8OLPFgsJnbvS8wUv/j+\nLn745RM5QRUfc11rt9Vzx79WU9fkZ3iRh0vPqmDm+CHkH0OQf1J/bLOehvpKpdRcrfUyYBHwOvAh\ncJNSyg64gPHAEb97rbnZ38OnT41MHB1B5vYFmdtbuvcVicbZtLOJNdsbWb21gea2EAU5Dr54tmLB\n7NGdvTk9Ngo9uUwozYUThgOJ/V1t/ggt7SEcNgsuh5W6lgBvrqrhg02JL6qfpbyMGJLNM29XceO9\n7/GVReM5dcqwbtf33sa93PfsJgwDzj1lJBeeOgq7zUI0FKG+vme7CXt5pH7Y+3oa6t8H7lFK2YBN\nwFKttaGUuhNYAZhI7EgN93D9QogMZBgGq7Y28MhrW2loTZx16XJYWTCzlIvmjsHlOHokmUwmcjz2\ng74RKMdjZ2xJLl84exwATntiPePK8rjr8bXc99wmwpEYZ8woPer631hVw99f0jgdFr7zmWmMK8vr\nSaspI998dAzSfXR0OJnaF2Rub+nSVzQWxx+KsrfRT019Ox9vqWfjjmYsZhNzpw/npPHFlJfkYrX8\n59qCvd3bnoYOfvPwSnz+CF87fwJzJnc9Yg9FYjy1vJKXPthNttvG9y6b3qvHnA/0OXUhxCAXi8cJ\nhGK0+cNU7vGxtbqV6vp2AqEo/lCUQChKOBI/5HGTRhdwxYIKhhV6+qXO4UUe/vuzJ3DzP1Zy33Ob\nMJlMnDxhCGZzIhcNw2DTzmYefFFT1xKgOM/Fdy6d2m/19TYJdSFEt4UiMZav2cMrH+3uvGjVgSxm\nE26nFZfDSp7HgcuRmPMuzndR6s1ixJBsSr2efj8ssaw4i+sun8atj6zmnn9v5J+vbmXKmAIi0Thb\nq1tp7QhjMsGnThrB4tNH47Cl1wEcB5JQF0J0aW+TnzdX1bChqgmb1YzLYWV3XTvtgQh2q5nxI/Jw\nOay4nVZGDsmmojSP0mIPFvPA/JqG8uG5/OCKGbyxqpq12xt5d8M+AHKz7Jw0oZhzThrB6GE5R1nL\nwCehLoQ4SEcwwl+e2ci6ykYAHDYLhmEQjsbxOK1ceOoozpxZSo7bfpQ1DTwjh2bzlUWJS2rXNHRg\nt1nw5jpTckJTX5FQF0J0isXj/Omp9WzY0czY0lwWzCxlxjgvVouZaCyOycSAHYkfC5PJROlxnJQ0\nkEmoCyE6LX1zOxt2NDO1vJBvXzK1c2cicNARKmLgkq0khADgnfW1vPTBboYVurnmgkkHBbpIHxLq\nQggaWgM8+JLG5bDyrUum4nbKh/h0JaEuxCBnGAYPvbSFcCTOFxaOY2iBO9UlieMgoS7EIPfBpjrW\nVTYyaVQ+p0wakupyxHGSUBdiEGsPRPjnq1uwWc188RyVUYf2DVYS6kIMUnHD4MGXND5/hMWnjaY4\nX6ZdMoGEuhCD1FPLq/hocx1jS3M5+8SyVJcjeomEuhCD0Nvrann2nR1485x86+Ipcgx6BpEtKcQg\nU7nHxwMvbMbtsPLdS6eRnYan+4vDk1AXYhCJxeM8+OJmYnGDay+anLaXlxWHJ6EuxCDy+soadtW1\nc+qUoUwaVZDqckQfkFAXYpBoaQ/x5LJKPE4rl54xNtXliD4ioS7EIPHIa1sJhmNcMq88LS+bK7pH\nQl2IQeCp5ZV8sKmO0cNymDtteKrLEX1IQl2IDPfU8kqeeXsHRblOvvnpyXL1xQwnl2ITIoM9/95O\nnnk7cTz6/14xg4IcZ6pLEn1MRupCZChfR5inV1SRl2WXQB9EJNSFyFCvfVxNJBrnvNmjJNAHkR5N\nvyilTMC9gAJiwNXJ/z8AxIH1WuslvVSjEOIYBcNRXl9ZTZbLxmlThqW6HNGPejpSPxvwaK1PA/4P\n+CVwO3CD1noeYFZKLe6lGoUQx2j52lo6glHOmlmKw25JdTmiH/U01INAbnLEngtEgBla6+XJ+18A\nFvRCfUKIYxSNxXn5g13YrWbOnFGS6nJEP+vp0S8rABewGSgELgBOP+D+NhJhL4ToZx9s2kejL8RZ\nM0rlYl2DUE9D/X+At7XWP1JKlQBvAge+erKBlqOtJD/fjdWaXh8Nvd7sVJfQJzK1L8jc3rrqyx+M\n8MSyKmxWM59bNAFvmn7f6GDaZr2tp6GeBbQmf25JrmeVUmqe1votYBHw+tFW0tzs7+HTp4bXm019\nfVuqy+h1mdoXZG5vh+vrkde20uQLcuGpozDHYmnZ+2DbZj1d1+H0NNRvAf6qlFqeXMcPgI+Be5VS\nNmATsLSH6xZC9MDuunZe/aia4jwX580emepyRIr0KNS11i3ARV3cNf+4qhFC9EjcMHjoJU3cMPj8\n2eOwpdm0pug9cvKREBng1Q93s62mlZnKy5QxhakuR6SQhLoQaa5yj4/H3txOjsfOFxaOS3U5IsUk\n1IVIY/5ghD89vZ543ODqCyaSm+VIdUkixSTUhUhThmHwwIuahtYg580ZJV9PJwAJdSHS1oeb6/ho\ncx0VpbksPm1UqssRA4SEuhBpqLU9xN9f3oLdauar503AYpY/ZZEgrwQh0tCfn1xHeyDCxXPHUJyf\nnmeNir4hoS5EmvlY17N8dQ3lw3NYMKss1eWIAUZCXYg04vOHeeilzVgtZq48d4J836g4hIS6EGnC\nSJ416vNH+OKiCQwv8qS6JDEASagLkSbe37iPj3V94miXeeWpLkcMUBLqQqSBxtYgf395Cw6bJXm0\ni0y7iK5JqAsxwO1r9nPzwyvxh6JcfuZYOdpFHFFPL70rhOgH1XXt3Pav1bR2hLno9NHMmz481SWJ\nAU5CXYgBKBaP8+aqPTyxrJJAKMrnF47jrJmlqS5LpAEJdSEGmG3Vrfztpc3U1Hfgcli4+oKJzJ40\nNNVliTQhoS7EAFJV6+PWR1YRicY5feowLplXTo5HvjxadJ+EuhADRGNrkDuXriUSjbPk4inMGOdN\ndUkiDcnRL0IMAIFQlDuWrqG1I8zlZ1VIoIsek5G6ECnm84e5a+laqus7OGNGCQtnyQ5R0XMS6kKk\n0L5mP799dA11zQFmTxrCFQsqMJnkxCLRcxLqQqRAKBJjxdpanl5RRXsgwvlzRnLR6WMk0MVxk1AX\noh/sa/KzrrKRQCiKzx/h/Y37aA9EsFnNfPlTinnTS1JdosgQEupC9LGPNtdx73MbCUfinbe5HVbO\nnzOKBTNL5ZBF0at6HOpKqR8AFybX8XvgbeABIA6s11ov6Y0ChUhXccPg32/v4OkVVThsFj6/cBxD\nC9y4HFaGF7lx2mVMJXpfjw5pVErNA2ZrrecAZwDlwO3ADVrreYBZKbW498oUIr2EwjH++NR6nl5R\nRVGukxu+OJOzZpYyaXQBY4bnSKCLPtPTV9Y5wHql1FNANvA/wFVa6+XJ+18AFgJPH3+JQqSXhtYA\ndz2+jt117Ywry+ObF00mxy1TLKJ/9DTUi4ARwPnAGOAZDh71twG5x1eaEOkjEIqyrrKRtdsbWbW1\ngUAoyhknlPC5BRVYLXKOn+g/PQ31RmCT1joKbFFKBYEDz5jIBlqOtpL8fDdWq6WHJaSG15ud6hL6\nRKb2BX3b27bdLTz/ThVvraohHIkBkJ/t4MoLJrFo9qg+e16QbZaO+qOvnob6CuDbwG+VUsMBD/Ca\nUmqe1votYBHw+tFW0tzs7+HTp4bXm019fVuqy+h1mdoX9F1v22taefyt7WzelRi7FOU6OXXKCKaN\nLWTEkGzMJlOf/pvKNks/vdnXkd4cehTqWuvnlFKnK6U+AEzAtcAO4F6llA3YBCztybqFGIj8wQi7\n9rVTXd/Ohqom1mxvBGDy6AIWzCpj8pgCzHLikBgAerwLXmv9gy5unt/zUoQYWELhGCu31vP+xn1s\nqGoiFjc676sozeWSeeWMK8tLYYVCHEqOqxKiC5t3NnPvcxtp8oUAGDk0m4mj8in1ZlFWnEVJkUdO\n6RcDkoS6EAcIhWM883YVL76/C5PJxDknlTF32nCGFXpSXZoQ3SKhLgatxtYgNQ0dQOICW6u21rNq\nSwOhSIziPBdXXzCR8hI5MlekFwl1MSjpXc3c/ugaItH4Qbd785zMnlTGp04eIWd9irQkr1ox6FTV\n+rhj6VricYPz54zCabdgMsG4sjzGDMuRuXKR1iTUxaBSXd/O7f9aTSgS49rFk5k1vjjVJQnRq+T8\nZTFoGIbBX5/fREcwylXnTpBAFxlJQl0MGht3NlNV28ZM5eXUKcNSXY4QfUJCXQwaz72zA4DzZo9M\nbSFC9CEJdTEobK9pZfOuFiaNLmDU0JxUlyNEn5FQF4PCc+/uBOC8U2SULjKbhLrIeDv3trF6WwPl\nJTmoEXKtFpHZJNRFRtu118dvH10NwIWnjpZj0EXGk+PURcaqqW/ntn+tweeP8IWzxzFlTGGqSxKi\nz0moi4zk6whzyz9X4fNH+OLZ4zhjRunRHyREBpBQFxnpqeWV+PwRvnTuBOZPlWPSxeAhc+oi41TX\ntfPWmj0MK3Rz0fyxqS5HiH4loS4yimEYPPL6VgwDLj+zAqtFXuJicJFXvMgoq7c1sHFHM5PHFDC1\nXHaMisFH5tRF2qtt7OCF93extbqVfU1+zCYTl59ZkeqyhEgJCXWR1mobO7j5Hyvx+SO4HBYmjylg\n7tThlBTJ18+JwUlCXaSthpYAtz6yGp8/whULKjhzRilms5xcJAY3CXWRllrbQ9z6yGqa20JcekY5\nC2aVpbokIQYE2VEq0k4oEuPOx9dR1xLg/DmjWHSyXKRLiP2Oa6SulCoGPgIWADHgASAOrNdaLznu\n6oT4hLhhcO+zG6mq9XHq5KFcdProVJckxIDS45G6UsoK/AnwJ2+6HbhBaz0PMCulFvdCfUIc5PG3\ntvOxrkeV5fHlRePlAl1CfMLxTL/cCvwR2AOYgBla6+XJ+14gMXoXote89nE1L7y3iyH5LpZcPEVO\nLBKiCz36q1BKfQWo01q/QiLQP7muNiD3+EoT4j8+2lzHw69sIcdt47rLp5PlsqW6JCEGpJ7OqV8J\nxJVSC4FpwIOA94D7s4GWo60kP9+N1WrpYQmp4fVmp7qEPjFQ+2ryBVm5eR/3PLsRp8PCjdfMYWzZ\nsX3RxUDt7Xhlal+Qub31R189CvXkvDkASqnXgW8Atyil5mqtlwGLgNePtp7mZv/RFhlQvN5s6uvb\nUl1GrxtIfYUiMTbtbGbt9kY2VDVS3xIEwGI28c1LppLrtBxTrQOpt96UqX1B5vbWm30d6c2hN49T\n/z5wj1LKBmwClvbiukWGMgyDzbtaWF/VyNbqVnbU+ojGDABcDivTygupKMtjankhpd6sFFcrxMB3\n3KGutT7zgF/nH+/6RGYyDIOW9jD7mvzE4onQrqlv543Ve9jXlPjEZjLBiOJsJo7OZ1p5EeUlOVjM\nsjNUiGMhZ5SKPhEKx9ha08LW3a1sq2ll1742OoLRQ5azWszMnjSU2ZOHUD48F5dDXpJCHA/5CxI9\nZhgGHcEozW0h4snRd1NbkA821bFqaz3hSBxIHB5VnO9i/Ih8hhW5sSV3jmc5rcwaX0y2256qFoTI\nOBLqolua20Ks3d5AVW0bTb4gjb4gTb4QoUisy+WL813MVF5UWR7lJbl4nHIIohD9QUJdHKLNH6a6\nvoPq+nZq6tupqm1jd137Qct4nFaK810U5jjJz3FgS54IZLdZOKGiiFFDs+VsTyFSQEJd0Nga4K3V\nNazd3kjlHh+tHeGD7rdaTEwalc/U8iImjMynKM+J0y4vHSEGIvnLHESisTgdgQj+UBRfR5gNO5pY\nu62RXQeMwgtyHJ2HD5YWeygtymJooVtOyRciTUioZ5hoLM6OvW3sqPXR6AvS6At1zoH72sMYn1je\najExfZyXCWV5TB1byJB8d0rqFkL0Dgn1FGtuC+HrCBMIRWluD7GtupWt1S20ByLJ+Won2W4bbocV\np93S+c0+sZhBIBQlEI7hD0YIhGK0ByJU17cTicYPeg6L2URBjoNxZXlke+y4HRbcThsVJblMGJVP\nWUl+Rp7BJ8RgJKHeTwwjEcKNvhANrQG27G5h7fZGahsPvVSC3Wom221nx942tu/xdfs5zCYTJV4P\nFaW5lJfkdu7IzPHYMctOSyEGhZSHumEYRKJxDAPsNnOXR0wkRp5G5/HNA5VhGLT5IzT6guxpSBw9\nsqfBn5wGCRIKH3z4n91mZlp5Id48Fy6HlWy3jTHDcxkxJAurxUw8btDaEe6cBw+GYxhGYgLFbDbh\nclhxOay4HVZcDgsOm0WOOBFikEtpqF931wraA5HO08bNJhMuhwVLcqdcPG4QDMeIxuKYAG++i1Jv\nFjkeOyYSIRqMxAgEowRCUfyhGIFQFAMDlz0ReLF4HH8oRjAU7ZxPNpvoDMTO/+wWIrH4IeuKxP4z\nleG0W7BbLTgdFsKRxP2h5Ak2GAaBcOyQqQ8At8OKN9dFQY6DwlwnhTlOyoqzGD8i74hvVGazifxs\nB/nZjt745xZCDAIpDfVcj52iXGfi1HATBJNBuj/kTSaSo1gL0ZhBTX07K7fUd7kuE+BMjljNJhMt\n7SH2NHZgtZhx7Z+PTo5iY/E4bf4Idc2Bzuf6JJfD0hn4JsAAYnGDRl+AYCiG3WbpXGb/2Dg/JxHY\nBTkOhhW4KfFmUeL1yIk3Qoh+k9JQ//lVJx3T8oaRnI444BoiLnsiWB0HhPaByx9pOmL/1M/+HY62\nA98AzIc+bv+lM4+2XiGESJWUz6kfC5PJRF6Wg7ys7k1HHC14TSYTdpsFu81yTF/TJIEuhBio5IwS\nIYTIIBLqQgiRQSTUhRAig0ioCyFEBpFQF0KIDCKhLoQQGURCXQghMoiEuhBCZBAJdSGEyCAS6kII\nkUF6dJkApZQVuB8YBdiBm4CNwANAHFivtV7SOyUKIYTorp6O1L8ANGit5wKfAn4P3A7coLWeB5iV\nUot7qUYhhBDd1NNQfxT4SfJnCxAFZmitlydvewFYcJy1CSGEOEY9mn7RWvsBlFLZwGPAj4BbD1ik\nDY7pwocClFbuAAATx0lEQVRCCCF6QY8vvauUKgOeAH6vtX5EKfWbA+7OBlqOto78fDfWAf4VdZ/k\n9WanuoQ+kal9Qeb2lql9Qeb21h999XRH6RDgJWCJ1vqN5M2rlFJztdbLgEXA60dbT3PzoV+6PJDt\n/5KMTJOpfUHm9papfUHm9tabfR3pzaGnI/UfAnnAT5RSPyXxbW/fAe5SStmATcDSHq5bCCFED/V0\nTv27wHe7uGv+cVUjhBDiuMjJR0IIkUEk1IUQIoNIqAshRAaRUBdCiAwioS6EEBlEQl0IITKIhLoQ\nQmQQCXUhhMggEupCCJFBJNSFECKDSKgLIUQGkVAXQogMIqEuhBAZREJdCCEyiIS6EEJkEAl1IYTI\nIBLqQgiRQSTUhRAig0ioCyFEBpFQF0KIDCKhLoQQGURCXQiRsQzDoCXUSnXbHgzDSHU5/cKa6gKE\nEKK3dUT8/H3TY2xvqaIj6gdguGcoC0fOZ2bxNCxmS4or7DsyUhf9ZnfbHv6ln2RbS1WqS+l3jYEm\natprU11GRorFY7SF2zt/NwyDf2xeytqGDbisTqZ7J3OCdwp7/XX8beMj/PrDOwhGgymsuG/16khd\nKWUC7gamAUHga1rryt58DjEwBKJBmoLN+CMBwvEI5bkjcVqdnfe3hzuoDzSSY8/CbrHz4o7XeKv6\nHQwMltW8y3TvZBaXn0uxuyiFXRyfllArDYEm8h155DtzMQyDxmAzTcFmrGYr2TYP/miA13cvZ1Xd\nOgwMJhYoFpcvojR7eKrLH/Ai8SgWkxmzyYxhGITjEYLRIFk2T+dIe2/HPu5b/w/2+uu4pOIC5pXM\nYcWe91lTv56KvDF8+4RrMJsSY9fGQBNPbX+elXVreWLbc1wx/pJUttdnenv65dOAQ2s9Ryl1MnB7\n8jaRAToiflbWrWFN/QZ08zbiRrzzvjxHLp+puJCpRRN5o3oFz1e9QigWPujxxa4iFo6cz7u1H7K6\nfj3rGjZx9sgzOGfkGdgstv5uh3AsTNwwcFodR1gmwtt73ueN3cuJGwbDPEPIc+RQ5dtFbce+zuX2\nB8eB/yYHKssajsvqYmOTZlPTFmYNmc55o8/G6y48aLloPMoHe1diNVs5aeiMXuiy/8WNOG/uXgHA\nGWWnYzKZjunxjf5m7ln3T1bXrwPAYrJgYHT+22bZPJxQPJUiVwHPVb5MOB7BYbHz2Jan2dZcyfrG\nzXisbr488bOd2wWg0FXAlyZ+ln3+et7e8z5TiyYyuWhCL3U9cJh6c+eBUuo24H2t9aPJ36u11qWH\nW76+vu2gJ48bcToifiLxCIZhYAAOix2n1YkZE4FokI6on1g8BoDZZMJuseOyurCbbUTiUSLxCP5o\ngNaQD1+4DRPgsDiwW+ydzxGKhWgINNEQaMJqtjAyu5SynFLiRpzWkI+OiJ9su4dcew65jtzOP3qv\nN5u6Oh+BaAB/NEgoFiIaj+KyuvDY3Dgtjs4RRNyIE4wG8UeD2C02nBYnNrP1oBd4LB7rDAan1UHM\niLPTt5uq1p20htuwma1YzVbcVhc59mw8Ng8xI0YkFiaOgdvqwm1zk+/IZYi7GKfVQdyI0xBopD7Q\nRDQeIWbEMQwDc3LEU+DMY7hn6EFzil5vNjV7m6gPNBCORfDYXHhsHlxWJ2aTGX8kMdp8Y/dygrEQ\nACOySxiZMwKP1UUoFmZ5zbtEjRgem5uOiB+Pzc2sIScQiAZoC7dTkTeGM0fMxWa2YhhGcrT0LC2h\nVordRZxVNpcsexZOiwNfuI36QCPReJT5paeS68g54usubsTxRwJk2T2H3Of1ZlNf33bI7S2hVm75\n6Pe0hFopcOYz1FNMls2D3WLHZrYSi8eIxqNsaNxMa7gt8TqzOGkN+wCwmW1U5I+hxDOM5lALjYEm\nTCYTXlcRhc58okaM9nAHMSPGSUNnoPLHArC5aStPbX+e6vY9mE1m5gw7kXH55eQ6cmkINPJ81Ss0\nBpsBWDhiPovLF3UZiofrq7eEYmH2ddSx119HR8TPlKIJFLkSb0CGYVAXaMAwDAqc+dgPeEMORIP8\nbeM/WdewCYBPl5/LwpHzu/WcsXiMZTXv8mzVSwSjIYZ7huKxuYnGo4AJl82J3WxnW0sl7ZEOAFxW\nJ1eM/wyjc0bwl3UPsqutGoBrpnyJad7JXT5PTXstN394Jx6bmx+d/D2ybIe+bg5U1bqTf2xeyoIR\n8zhl2Kxu9dKV3txmXm/2Yd8pezvU7wGWaq1fSv6+Axijte5y+LI/1F/d9RbLq9+lJdRK1Ij1Wj29\nJcvmodBZgGGOU9fe0BlsXTGbzNjNNkKxMIm3pf+wW+wUu4rwuotoD7ezw7ebSDzSa3XmOXLpiHQQ\niUePuJzNbGW4Zxhmk5lIPEIgHqDR33zIciZMiTebeIxwPEKWzcNZZXOZOWQ6ha78g5bd56/nUf0U\nunkbp5WcwgVjzsFjcx+xjkA0yLOVL3VOyxyup69P+TIjcroeG1S37eHvmx6lpmMvX5/y5UNGXl39\nIUXjUe5Y9WcqW3cyIruUllArvnDXf2x2i535pacm33Q8+CMBmkMtFLu92Mw9+6AbN+KsrFvLs5Uv\nUR9oPOg+q8nCqSUns6lpC3X+BmYPO5HPqYsP2bHXV6EeN+K8Wf02z2x/8ZDX5vj8CoZ4itnQsImG\nYFPn7Tn2bLyuIoa4vVS27mCvv45x+WOp89fTEmrlyomfY9bQE474vJubtrJ06zPUduwjy+5h8ZhF\nnDJs1kEj7f1i8RhbWrazy1fNzCHTKXIVAIlPVc9XvUKW3cOCEfOO+Hwv73yDp7e/QLY9izG5oxiR\nXYrH5sJisuK2uSjPHUW2PYuP9q7ioc2PEYvHuHrKFw/7RtEd6RrqtwHvaq2XJn/fpbUecbjlo9GY\nYbVaeGDlo7yz+2MK3fkUuvJxWh2YTCYMDELRMMFokGg8hsfuJsvmxmqxggFxDILREP5IgFA0hN1i\nx2Gx47I5yXflkOfMxZQc4YeiYUwmMJssOCx2irMKGeIpIhANsr1pFztadmM328h35ZJl9+ALtdMc\naKHB30x9RyP1/iZsZivFnkKKPAVk2RMjWavZQkckQFu4g2AkSDgWIRQL47Q6yHZk4bY5CUcjBKIB\nWgI+atvrCMcimDBRljucisLR2MxWgtEQceKMzitjXNEYhmR5icajhGMR2kMdtAR9tIXasVms2C12\nzCYTHeHE89Z1NFDdWkttWx05jizKcoczLLsYh9WONRk8cSNONB6jtq2O7U072NW6B0h8EnLbXAzN\n9jIsewguq4P2sJ+2cAcdYT/+sJ+oEWPeqFP41Nh5OG3Ow21ODCOxPVxHWKYru1pqqGzeRXvYjz8S\nINeRzdBsL5VNu3hk3TPYLFa+OuOznDbyxM5pGl+wjRe3vcmTG18kZsSxmMzYLXZuWvA/lOYOo8a3\nl7+veYLpQydxTsXBf+D3fvxPXt62jDkjZvGdU67CZDIleo0ECEXDhGMRrGYLVouVfGfuMffTXdF4\njLV7N1LX0UhzoBWAheWnU+QpwBds45fLfk9l8y4mFyv+6+SvUODO65M69qtrb+DuDx5kY/1Wsh1Z\nnFo2i5KcoVjNFpbtfJ9N9duAxOh42tCJuG1O6v2N7G1voKGjqfON+bxxZ/GFaRdR49vLT16/lXAs\nwtdnfZ7TRp6ENfnm1OhvZkPdFna2VLOtaQeb6rdhwsSZY07lc1MuJMeZ3ae9xuNx7l/1Lz6oXk1L\n0NflMsOyiqltr8Nlc3Ld7K8xfdikPq3pGPVbqF8MnK+1vkopdQrwE631eYdb/pPTLwOZYRh4vdk0\nNLQffeEj2D/F47A4cNtcvVRdz+rYPwrq64/yx2Ndw0b+uuHhxBulxcnkovH4wu1sbd6OgUGeI5cr\nxn+GYDTA/RsepshVyLzSOQeNNM8oO42Lx55POBbm1V3LeGHHq5RkDeO/Zy7BkZyWG4iC0SAPJKcy\nPFY3l6lPU5ZdAoZBQYGHpmY/JkzkO3K73CcRi8d4bfcymoOtzC+dwxBPMbF4jA/2rWJN/TpG5Yzk\nxCEnYLfYeHnnGyyreZdoPMq0okl8bvwlZNuzDlrfvo46fOF2RueO6Bws7BeJRTo/dQzPGtp5u27a\nxt1r7iNqxMhz5DKjeCpVrTup8u066PEVeWO4pOICyrJL+vX12Hkce/seQtEQUSNGS8jH1ubtVLbu\nIN+Zx9VTvsQwz5Djfq50HanvP/plavKmK7XWWw63fDqFOgzs8DseA72vOn89K2reZ1X9OpqS882j\nckZwQvEUTh1+Ei5r4s3x39tf5MWdrwPgtrpYXL6I5bXvUu2rZUzuSGo76ghEA2TZPHx/5n8dspNy\nIDIMg+U17/HEtmcPO1Xntro4dfjJnF4yu3NarCHQyAMbHqHKtxNITKVNKhxPbcfezjn7/axmK9F4\nlAJnPheO+RSzhkw/5p2bR1Lnb+DN6hW8W/sR4VgYEyYq8suZUjSBsqwSSrKGHTTAGSivx1g8htlk\n7rV/i7QM9WMloT4wpEtfhmGwp2MvbquLfOehUxFxI84j+knaIx1cWnEh+c48XLlmfvXG3WxtqSTL\n5mF+6WnMLZ191Pn+gWZvxz6W1bxLJBbFZAKX004gGCYaj7GhcTPtkQ5MmPDY3HhsHlpCLYRiYWYW\nT2OqdxKv71rOzrbdWM1WTh1+EqeXzKaqdScf7F1Ja9jHvNJTOXX4yT3eT9Ad/oifytadjMwpO+RT\nwIHS5fV4rCTUByB5saUfrzeb2n3NbGupYkzuyM6joNLdgdssEo+yct8a3t/7MS2hVjoifiwmM58e\nex4nDjkhsX8q+YaYbc8ix96389XHK1Nfj/0V6nKZAJHxrGYr4wsqUl1Gn7GZrZw8bCYnD5t52GVM\nJhMlWcP6sSqRKnKZACGEyCAS6kIIkUEk1IUQIoNIqAshRAaRUBdCiAwioS6EEBlEQl0IITKIhLoQ\nQmQQCXUhhMggEupCCJFBJNSFECKDSKgLIUQGkVAXQogMIqEuhBAZREJdCCEyiIS6EEJkEAl1IYTI\nIBLqQgiRQSTUhRAig0ioCyFEBpFQF0KIDGLtyYOUUjnA34EcwAZ8T2v9vlLqFOB3QAR4RWv9/3qt\nUiGEEEfV05H694BXtdbzgSuBu5O3/xH4rNb6dOBkpdS04y9RCCFEd/VopA7cDoSSP9uAgFIqG7Br\nrXckb38JWACsOa4KhRBCdNtRQ10pdRVwHWAApuT/r9Raf6yUGgo8BHybxFSM74CHtgGje71iIYQQ\nh3XUUNda3w/c/8nblVJTgIeB/9Zar0iO1HMOWCQbaOmtQoUQQhydyTCMY36QUmoi8DhwmdZ63QG3\nrwQuAXYAzwI/11p/2DulCiGEOJqezqn/EnAAdyilTECL1voi4FoSo3cz8LIEuhBC9K8ejdSFEEIM\nTHLykRBCZBAJdSGEyCAS6kIIkUEk1IUQIoP09OiXjKOUOhn4tdb6jOTlDf4IRIGtwDe01hGl1CLg\npyROwPpQa/2dAx4/HngPKNZah/u/g671tC+lVB7wIJAPdABXa613p6aLrh2tN2AiiWsR7T9x7hRg\nsdb65eTj03KbcZi+gA9I822WfD1eC1wFxIFfaa2fOuDxabnNDtdXX/ydyUgdUEpdD9xD4jBNkj9f\np7WeC+wBvqmUygJ+A5yntZ4D1CilipKPzwZuBYL9XvwRHGdfNwDvJK/jcwtwV783cATd6G2J1nqN\n1voMrfWZwB+ApQcEerpus676eizZV7pvs28qpTzA9STeqM4h8ea1//Hpus2O1FevbzMJ9YRtwEUH\n/F6mtX4/+fPbwHxgDrAOuF0ptQyo1Vo3JJf5C/BDwN8/5Xbb8fQ1CXjhgGXn9UvF3Xe03ubuv0Mp\n5QZuBL5zwPLpus2O1Fe6b7N5JD59GCTOSM8CYgcsn67b7Eh99fo2k1AHtNZPkviotN92pdTpyZ8v\nANxAEYkQvB5YBFynlBqrlPoZ8GzyzFpT/1V9dMfTF7AKuDC57GLA1R81d1c3evMccN9XgUe11k0A\nab7NuuqrOfl72m8zrbUfeATYCHwE3AmglPo5abzNDtcXfbDNJNS7dhVwg1LqFWAf0JD870Otdb3W\nugNYBpwAfB74qlLqDWAo8HKKau6O7vY1Hfg1MFop9SYwEhhQc7NdOLC3OhJ97fd54N4Dfv8C6bnN\njtZXOm+zfUCDUmo2MJtE/SOAi5RSJ5Lmf2eH6WsWfbDNJNS7dj5whdZ6IYmR7Msk3lEnK6UKlFJW\nEnNjG7TW47TWZ2qtzwD2AgtTVvXRdbevjSQ+Bv4lec38bcCK1JTcbQf2Vkjyjz75hS52rXXN/gW1\n1hVpus2O2Bfpvc32vx6zAL/WOpLcEdoC5KbxNjtSX3n0wTaTo1+6tgV4TSkVBD4EHtRaG0qpH5LY\nQAbwL631xk88bv/RCANVt/tSSoWAvyWv7dNE4stQBrJDekvePo7EBeYOJ+22WfL2rvraTJpvs+Tr\ncaFS6n0SUxortNavfuJxabfNDteXUqqcXt5mcu0XIYTIIDL9IoQQGURCXQghMoiEuhBCZBAJdSGE\nyCAS6kIIkUEk1IUQIoNIqAvxCUqpHKXUk6muQ4iekFAX4lAFwLRUFyFET0ioC3GoO4DhSqnHU12I\nEMdKQl2IQ30b2KO1viTVhQhxrCTUhRAig0ioCyFEBpFQF+JQUeQKpiJNSagLcah9wC6l1GupLkSI\nYyWX3hVCiAwiI3UhhMggEupCCJFBJNSFECKDSKgLIUQGkVAXQogMIqEuhBAZREJdCCEyiIS6EEJk\nkP8PkViM2woVg0QAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x11752fcc0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pred_s = res_s.get_prediction()\n",
"\n",
"ax = data.wpi.plot(label='observed')\n",
"pred_s.predicted_mean[1:].plot(ax=ax, label='Forecast')\n",
"plt.legend(loc=2)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x117668978>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAERCAYAAAC0FCalAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnWeAJFd5rp8KnWemJ89smJ1N2l7tKiFpFUASkpAEEjJJ\nBsG1jUHXNgaML1yDTbDB2NiYi8EmC2Es22AbY5AVQIBAGeW4q029eXZ2dieHzl3x/jhVHXZ6cuzZ\nen7tTldXn+queuur93zfdyTbtvHw8PDwqC7kpR6Ah4eHh8fM8cTbw8PDowrxxNvDw8OjCvHE28PD\nw6MK8cTbw8PDowrxxNvDw8OjClHn8uZYLNYKPA9cF4/HD8zPkDw8PDw8pmLWkXcsFlOB24HM/A3H\nw8PDw2M6zMU2+XvgW8DJeRqLh4eHh8c0mZV4x2Kx9wD98Xj8l4A0ryPy8PDw8JgSaTbl8bFY7FHA\ncv57ARAH3hSPx/vncWweHh4eHhMwK/EuJRaLPQy8b6oJS8MwbVVV5vRZHh4eHmcgFd2NOWWbOExL\n/UdGqmdes6WlloGB5FIPY0FYqce2Uo8LVu6xecc1/f1VYs7iHY/Hr53rPjw8PDw8ZoZXpOPh4eFR\nhXji7eHhsSKwbZv/fvgQLx0cWOqhLAqeeHt4eKwIRlMaP3vmOF+/6xWe27/yE9888fbw8FgRGKbI\nXrZt+PY9e3ghvrIjcE+8PTw8VgSueG9YVYtPlbn9nt0c71t52Swunnh7eHisCAxTZC1vXBXlHdds\nwrRsjp5KLPGoFg5PvE/jQx96H6Ojo0s6hkceeZB//uc7lnQMHh7Vhht5q6pETdjv/G3lLrDuifdp\nSNLyaNWyXMbh4VEtFMRbkVEVcf3ohjXZW6qa+aiwnBd++NCheZ8h3rG1lXdcu3nC1w3D4POf/ywn\nT/ZgWTbveMf/wrZt/uZv/oYTJ04SDIb41Kc+g23DZz7zCWzbRtM0PvrRT7B581n8+Mf/xS9/+Qsk\nSeK6627glltu5W//9rOMjY2SSCRYt66TCy64kBtvvJnh4SE+9rEP893vfo9vf/sb7Nr1MpZlcuut\nv8XVV7+OV17ZyVe/+iVqa+vw+Xxs3bptXr8LD4+VjhtlK7KET5Gdv3nivSK55567qK9v5C/+4q/J\nZDLcdttvEwj4eetb/4izzjqXu+/+Ef/2b3dy8cU7iEbr+fM//yxHjx4hl8ty7NhRHnzwl3zrW9/F\ntm0+8pEPsmPHZQBcdNElvOMd7+LYsaP8wz/8P2688WZ+/vP7eeMb38TTTz/JyZM9fOMb30HTNN73\nvvdw8cWX8qUvfYHPfe4LrF3bwbe//Y0l/mY8PKoPV6h9qozqiffi8Y5rN08aJS8EXV1H2bHjUgDC\n4TAbNmzgueeeYceOHYyO5ti+/VyefPIJPvShj9Dd3c3HP/5/UVUf7373bRw5cpje3lP8n//zfmzb\nJpVK0tPTDcC6dZ0ArF+/AdM06e3t5aGHfslXvvJN7rnnLuLx/fzxH/8htm1jmianTp1keHiItWs7\nALjgggvZu3f3on4XHh7VjivUiiyjqkK89RUs3me0593ZuYGXX34JgEwmzeHDh1i1ajU7d+4E4OWX\nX2LTps28+OLzNDU18+Uvf513v/s2vvOdb9LZuZ6NGzfx1a/ezte+9m3e8IY3smnTWQDIcvFrvfnm\nN/Otb32FDRs2EonUsG7dei666GK++tXb+cd//CbXXns9a9aspbm5haNHjwCwZ88ri/xNeHhUP6Zj\nm/hUuWCbmCt4wnLZRN5LwZvf/Da+8IXP8YEP/B6apnHbbX/A/fffx7333suBA4eora3jz//8LzEM\ng8985pPcffePsCyL977399m0aTMXXriD97//f6NpGtu3n0Nzc8u4z7jmmtfxla98iS984csAXHHF\nVbz00gt88IO/Tzab5aqrriYcDvNnf/bnfP7zf0U4HCYajbJ+/cbF/jo8PKoaN8pWFKk4YbmCI+85\n9/OeLgMDyaq5Ba7UVpWwco9tpR4XrNxjm+/jenzXSe68fz+33XQ2Z3VE+cS3n+bK81bx3pvOnrfP\nmA4L0BK2YurZGW2beHh4rBwM0wbJQlE4I7JNPPH2WBB+9kwXz+7rW+pheJxBZIwMwQseYWf6iUK2\nib6CPW9PvD3mHdu2+fEjR7j/qa6lHorHGURP/jCST2NAP1nwvI0VXKTjibfHvKMbFpZtk84ZSz0U\njzOIHuMQACkjcUbkeXvi7THv5DQTgEzeE2+PxSGjZxiyTgCQNpNIsrBLPPH28JgBOV2IdzZvYFkr\n13P0WD68MrgPGwvbBhubhJZEkaUVnSo46zzvWCwmA98BYoAF/GE8Ht87XwNbLHp7T/G7v/tOYrGz\nsW0bSZK48srX8Pa3/86SjOexxx5h+/ZzaGpqXpLPnw9yJRF3Jm9QE/It4Wg8zgReGtgFgDXShtLY\nx3BuGFWRMYyVGzzMpUjnNwA7Ho9fEYvFXgv8LfCW+RnW4rJhg6iUdFnKvNr//u//ZP36T1a1eOed\nyBs88fZYeLJGjn1DBwjbjYyONTviPYqqSCvaNpm1eMfj8Xtisdh9zn/XAyNzGchdh37CS/3zWxb+\nqtZzedvmm6fcrlKh0te//o/s2vUykiRx/fWv5zd/851lHQO/+MWv8O///q/jugPu2bObr33ty9i2\nTUtLC5/+9OfYu3c3d975HWzbJpvN8JnP/A2trW18+tMfJ51Ok8vl+IM/+ACGoXPw4AE+97nP8M1v\n/hOqWp0FsK7nDZDJ6UBo6QbjseLZPbgPwzZpMzsZ0YQTPJwbQVWDnm0yEfF43IrFYncCbwV+c36G\ntPgcO3ak0ChKkiRuvfXt9Pae5I47/gXDMPjgB3+fCy+8GCh2DHz66Sc5derkuO6Af//3f8tnP/t5\n1q3r5Kc/vZeurqMcO3aET3/6r2lqauZ737uThx/+FVdc8VrGxsb40pe+xsjIMN3dx7n88ivYsiXG\nxz72yaoVboB8mXh7k5YeC8vLAyLoq9M7sfODgBBvn7LGi7wnIx6PvzcWi/0Z8GwsFjs7Ho9nZ7Of\nt22+eVpR8kJwum1y770/5LzzXgWAqqps23YOR48eBYodA48cOcT+/fsqdAccLmzzxje+CYC+vl7+\n4R++SDgcZmCgn/POu4ANGzbypje9lb/8y09iGCZvf/utgHgKWKyWBQtFVivxvD3x9lhgetP9RNQw\nqlGHraUBGMqNoCodzpPfymQuE5a/A6yNx+OfB3KAiZi4rEhDQxhVVWb7cQuGpiVQVZmWltrC3zZt\n2sRdd91FS0stuq6zf/9ufuu3bmXnzudoaIjQ0lLLueeezdBQH3/1V3+FaZrcfvvtnH/+Vtrb28hm\nR1i3bh3f/e536ezs5Itf/Ft+9atfEQ6H+fjHP0447Gdk5BSKYnHnnd9lYGCAd73rXbzpTTcSCPio\nrw+VjWe+Wch9A/j8RY9b9ikL/nkui/U5S8FKPbb5OC4Tg7A/iKIqYCnU+msY08cIBlSSGW1JvrvF\n+My5RN4/Av4lFos96uzn/8Tj8fxEG4+MZObwUQvH8HAaw7DKJiivvvpqHnnk19xyy9sxDINrr72e\n5ua15HI6Y2NZBgaSnHvuDh555Ne84x3vLHQHzGQsPvKRP+OjH/1TZFmmqamZm256G9dffyO33vpO\nmptb6excT1fXCWpqmnnssSe4996fYNs2733v+xgYSBKLbedP/uSjfPnL36C2dv5PgMWYjB0cThf+\n3TeYWpTJ35XavAlW7rHN13Fl9Tw1vgjpjAZAvb+eU+lTNGGhnXZtz4RHX+5hbUsNm9ZEZ/S+BWhM\nVfHvXlfBCqzUiwUW59h+/OhhfuqUxr/x8k5uee2mBf088H6zamS+juvDj3yKVZE2lMNXsvfYCJff\n1MPLg6+wqvdNHD2u8U9/ds2M14TNaQYf+PJjXLilhT9627kzeq/XVdCjasnlixOWXom8x0Ji2Ra6\npRNQ/IU1LBtD9QDY/iw2YM6iUCzrnMNuj5TliCfeHvNOTi+dsFy5E0YeS49uiXPNp/gwTAtVkWgK\nNYoXfcKqnc1qOjln0j3oX37zdC6eeHvMOzkvVdBjkdBM4XMHZD+GaaEoMk3BBgBMVYj3bHK93UKz\noH/5pux64u0x77h53pLk2SbVnva53HHF26/4MU0bVZZodMVbERPns8n1dq2/gM+LvD3OIHKaiSRB\nbdh/RncWfOC5bj7x7afRStoFeMwvmiVsOb/iRzctVFWmMSg8b8MV71n09HabqwUDnnh7nEHkNJOg\nXyUSVM9oz/vQiVH6R7MMJXJLPZQVSyHyln2YpoUqy4TUECE1hCYJ8Z6NbVLwvL3I2+NMIqcZBP0K\n4aBKJmecsdaB5kR8nu+/cORLbBPdtFFVIWmNwXo0KQXYhSyUGe1X8zxvjzOQvG4K8Q74MC0bTV+5\n/SUmw7VLznTffyEpet5O5O2k9jUGG7AkA1R9dp63I94BL9vE40wip5kEfAqRoIha0meodeJmLGTy\nZ+bxLwalnrdh2qiykDQ340TyZ9Fn4XkXI29PvD3OEEzLQjesgm0CZ65t4D5xnKnHvxicniqoqsXI\nG0AKZL3I28NjOpR6hQXxPkMzTgqRtyfeC4Yr3qosLDo38q711wAgzdY28fK8Pc40ciWPm+GA6C54\nptommifeC45rm/gkIbLuhKVf8YsNZBN9FkuhedkmHiueVFbnrscOF6LM0sfNM902yTu2yZl681oM\n3GwTGREoqLKwTfyy+L8km7OKvAtPkF6et8dK5dl9ffzkyS5ejA8ApWXFxQnLM1G8bdsuRt5nqG20\nGLi2ieJ0t1aV8ZH3nDxvL/L2WKlkHWEaS4uLyF05XlYNfH7xuHomRp6GaeE+rJ+JN6/FQjPFuSVz\num3iRN6KOcsiHRNVkQs3g+XI8nXjPaoCN0JJuOKtm4DNs+aP6RlcBWxYsZHnT586xhOv9AKij8vN\nr17P5dvbgaJlAmfmzWux0Aq2iSPeBdvEjbytWZXHu7UKy5nle1vxqApc8S5E3poJsknGTjCQ7wNW\nbuT58Es99I9kyeR0Tg1leG5ff+G10n4mK/X4lwOa5Yi3XTnyFrbJ7CYsPfH2WNG4EztJZwmqvGYi\n+cRqeEk9CdgrUrx0w2QkkWdLR5Qv/dFrgGKGAhS9f/DEeyFxI2/JFW8nVdCNvCV5drZJXjOXdY43\neOLtMUdcwUqURN6ueJu2iaTqVduc6m+/9wL/8rP9FV8bGM1hAy31IRRZxqfKZYJd2hIgmzewztD+\nLgvFU3t6+fBXHyeji3NNsoXQKk55fGnkbc5QvG3bdpqreeLtsYJxixnGMq54G+DTCq+HagzSVeh5\na7rJoZ4xXoj3V2ys1T+SBaC1IQSI7JrSRShKhdymOJHrMT+8cmSIREYnnc+hyiqWo88+xzZRZWc6\nbxaRt2HamJa9rHO8YZYTlrFYTAX+GVgP+IG/icfj983juDyqBFewkmkdy4lYJH++8HogrJNJVp9w\npbLiaSGdMxgay9FcHyp7vX9UiHdbQxgQKWU5zSRvajzc/Titxtll26dzBuGgbxFGfmbg3jw1S3dK\n48UNVnEmLGVJRpVULNnCmGGRTnEJtOWdzzHbyPu3gcF4PH4VcCPw9fkbkkc14Xreli287bxetE0A\nfEG9Kj3fZKZo9XT1jV8JvH9ELLHVUu9G3ip5zeT5vpe478gv2Du2GwDZWbW8Gr+D5Ywr3rql43M6\nCkIx8gbwyb5Z5Xnnq6CvCcxevH8I/EXJPqrT1PSYM6WTdGNprWzCEkAJ5snrsyuUWErcCViYSLwr\n2yaDmWEAEvoYANEaMXFWrb7/ciST0wtPRoazcrxrjSjyePGeqW2Sq4KOgjBL2yQej2cAYrFYLfDf\nwKfmc1Ae1UO+xOdNpLWyCUug8O9MzqAu4l/08c2Wssi7NzXu9f7RLHVhH6GAuISCfgXLthnKCvFO\nGUmgifqaACPJvNfTex4Yyo7w+ef+gRvaby78zbB1ZxUdYY24/bxBVFlKcg5jhv3k3Xmc5R55z9rU\nicViHcBdwNfj8fh/TbV9Q0MYVV3eX0YpLS21Sz2EBWM+jy1XcmFIioIFSD4NCQkbGzkgIthgJEBL\nS828fW4l5vO4bKWYs93dn6K5uQbJsUBM02JoLMdZHfWFz6yrDQAwaiQAyCFslfbmCEdPJZB96pzG\nt1LPx5kc19ETh8kaObrTPUAUsDHQiQRDhMIiMGioDxf2GfIHQRlBsZQZfU73sHiqamoIz/p7X4zf\na7YTlm3AL4APxuPxh6fznhHHI6wGWlpqGRgY/6i8EpjPY7Msu6wYpbt3jEQqj1STpzHYwJiWQLPE\nOoLdJ0fxs3DpcvP9m53qF/uqCfkYTeU5eHSIBkeg+0ezmJZNQ02g8JmSk5HSnxoEIKULEQ87GQv9\ng6lZj2+lno8zPa6Tg+K77U+MAFGQROAgWQrDo0JfshmtsE/FVkC2SGXyM/qcfmdbUzNn9b3P9+81\n0Y1gtp73J4B64C9isdjDsVjsoVgsFpjt4DyqE9cbDDvWQSKtkdMNJJ9GNFBL1F+H5kSg1TZh59om\n2zc0AnC8xPd2JyvbGooZKEGfCpLlFCZBzhY3rfpaERF6JfJzJ62L7z2lOYGgIs4/v+Iv2iZq0TYJ\nKD4kyUY3Z3burXTP+8PAh+d5LB5VhpvL3NIQoqs3KcTbzIJkU+evAySGcl1UY5WlO2G5fX0jz+zt\no6svyfmbm4HiZGVLiXgH/AqSP1v4v4kOskFjbRDwOgvOByld3BCzhvieJdkRb7k4YamWTFj6VXHj\n1M2Z3TirYRUd8Ip0POaAm2niRqDJjE7eFlFRnb+WaKAOsMGXr7psi2RWR5Jg23qxnFZXb2nkXZ5p\nAiJKkwLZsn1I/hz1jtVSbTev5Ygr3pqdE/ncrngrpROWRUkLOG1hdUtjJqz0PG+PM4Tn9/fzd99/\noSx1zsWNUOprAqiKxFhaK9gkdf5a6v11gMg4qbZsi1RGpybko6E2QF3Ef5ptUl6gA27knQNK1k/0\n56iv8WyT+SLtiLcla6xqChfEuzRVsNQ2cRdk0OyZffelPemXM554e0zKIy/3cODEGL94tnvca6Xe\nYG3Yz9BYFlSRGlgXqHEib5D8+bJy8WogmdGoDfuRJIn17bUMJfKFG1j/aJZQQC0sNgHlkffGaCcg\njrs27EeRJbJVdvNajqQczxtVZ21rpMQ2KRbplNomPifyNuwZet55T7w9qhzDtDjUI4pNHnzxRKEw\nwqV0seG6iJ9ERi/kdRdtE5B8uaoSb9OySOcMakMiclvXJmb7j/elsGybgdEsrQ2hQuogiO+gKN7r\nAXHcAZ9YUajanjyWI65tIskWq5uDZROWRmHCssTzdiJv05qh510led6eeHtMSFdvEk23iARF6fev\nni+PvoveoELUKcBx+5rU+UW2ifs3bYaFEktJKiuOqzYsLv5OR7yf3N3LaDKPblhlmSbgRN7+LCCx\nvq4DEMetKhKhoK/qPP/liGubANRFZVSfOKeEeDuRd1mRjvj99BlG3qVByXLGE2+PCTnQPQrAO67Z\nTE3Ixy+fP1E28ZbTS20Tp+lSSeRdHyh63ppRPZG3a4/UOIUf52xsZG1LDU/t6eU79+0FyicrQTSm\nkgJZAkRoCor0QiWYR5IkIkGVTN6o2J3QY3qYlknWyBX+H4lYqD7xffoVX1G8S7NNHNvEnGH3joId\nuMy7Cnri7TEhcUe8z93UxOsv6SCbN3jwxROF111vMOBXCqXvktMOttZf7nlXVeTt5Hi7tknAp/DR\nd13A6uZI4TtpOa3LoM8nJigDdg0RXxgsGdl5CgkHVQzTRpvFclwegrRRXuQXCJWIt1wSeVewTSxm\nGnkbSIDft7zlcXmPzmPJsCybgydGaWsIUV8T4NoL1xIKKDy+82Rhm+KsvEo07Ip3Hh9BVFklqAYJ\nKH4kX66sEnO5k3S8/cLTBFAX9vOxd15Ae6PIMFndFCl7j0YaSQLVjCBJEpIRBJ/wwN0iJi9dcPak\nNMcycR5edPKn2SaO5y0XbZNZT1g6q+iUzmksRzzx9qhId3+KbN5kS0c9AKGAyqqmCCPJfOHxv9Tz\nLkbeeYJyMYUu6q9zIu8qEm/HNqkNlzfSitYE+PhvX8gH3nIOG1fXlb2WsUQqoWyKY7e1ILaSx7RM\nIk4fb8/3nj2FNMG8eOLJ6BkUVYh3QJk88raZmWWVq4LFh8ETb48JcO2B2Lr6wt9qQz5Myybr2CX5\n0lTBiB8kE0k1CCvFBlTRQB2STyNnLI+oc1/XCDsPDU66jVsaXxp5u9SF/Vy8tXVcVJayRC8TSRPi\nbWkBkCChJQk7KYVexsnscdME7Zx44knrGWTFEWxJLUbeFYp0kC1MawbirZkElvlkJXji7TEB8eMj\nAIXI27IthqPPo3bsJ5ktWSke4QlHw/6C312jlos3QN4qZgosJd/7RZyv/GgXj77cM+E2ruddEyqK\nd6US6+OJE+wbPgDAmCZudnY+iGXZmHmny2B+rCDenm0ye9w0QTsnbo5pI4PkiLdkq4XIW5HHZ5vM\ndEGGfBWsXwmeeHtUwLJtDp4Yo6kuQHNUPKbec/hnDKpx1NbugrjlNBMpkOFHR/8bQ0kVMk1qfBXE\nm+XRVTLr9Bj5t5/H+fWuUxW3cW9Orm3Sm+7nTx77NI+eeLKwjWEZfGvXnXzj5e9yZOwYwzlxszNz\nIfK6ia2Jniaj+UTRNsl7tslsKdgmTuSd0TNITp43llIQ5/KVdNzI2yxE5lNh2TZ53Vz2mSbgibdH\nBU4OpklldbZ0iDLvJ3qe4VfHHwVAUkxGUkKIc5qB0tDHzqFXuO/4fYXsijp/sYWlWyKvS8tDvDXD\nJBrxEw6q3Hn/vooWyum2yf7hg5i2yQNdD2NaQjB2DuwmoSWxsfnXPT/gVLoPbNCzfjTDKhHvscKE\npWebzJ706baJkS1UWGLJ49awhGLkLSkm+jQzffJV0lEQPPH2qEDvkLhQOttrOTR6lB8c+B8ivjAt\n6hoABtLCIsjrJoqz2ML+kYMEVokinmigKN5u5G3I5U2blgpNt2ipD/FHbzsXG3h2X9+4bZIZjVBA\nLfinxxLiuEbzY7zUvwuAx3qeAuCi1vMZzA1zPHkCyQySz4vIDV3YJmP5hGebzIL+zAB7h+KF/xds\nk3wYCUlE4o5426aKaVqoilQ2F+Evi7ynJ97V0lEQPPH2qIBbBl8X9vF4z1NYtsVt23+L1SEh3sNZ\nZ7UYrSjeEhJ2jYhiG0LRwr6iAfFvS8kueZGKYYqJK79PLpS8VxLUZEYvm6w8nuzGJ/uQkHio+9f0\npE5xaPQoWxvO4t3bbmVtzWpApAnmdBOtzDYZK8k28cS7EpqpY9lFcT2Z6uXvn/8G39z5z8Ue3o54\nh5UIITVIRs9iO+JtmTK6aaEo5XI2G8+7WppSgSfeHhVw85xrQj5OpnrxK362NGyiPiSi6NG86HeS\n08yCVXJ959WF9zeFiml0tX7H//bll3wRYrdQyK8qIo+XogfuYts2qaxeKNDJ6Fn6MgNsjHZybvM2\nupLd/Mf+HwNw1dpXo8oq79n+LoJKgKDRhKZb5DQTW680Yel53qeT0tN88om/5vPP/iMPH9jJx//5\nQb7y4h2kjQw2trCjEHnetiXTWl9D2BcmrWewJQPbkjAMME27LMcbSmwT2Zq2bVIt7WDBE2+PCqQd\n8Q4GZXoz/ayOtCNLMk1hEUUnNLEgb04zQc0T8YW5cf111CpRsCVWRxsL+6pzxFvyaeSXuMrSLdH3\n+2RkSSIUUMctknDvoQdQtz1OJCwujeNJUVHaWdfBtR1XAHAscZyGQD3nNG0FYFWkjc+95lOsyl8M\nQDKtgS3jJ1Rum3gLMoyjPzNI1shxMt3Lj078O4m1D5IyUmyu3wBAryPeyXwaDB9t9WEiapiMIcQb\nSyGvixXiS3O84XTbZHpPffmSDKrljifeHuNwbZMso1i2xZqadgBaIiJtMGOksW2bvGZiqTnq/LX4\nFR9/suMPef/57yHiL1YfBpUgkq0gqUtfqON+vt+5MCuJ956hOHI4hVnTC0CX43evr+tgc/3GgkVy\nxZpLUeTiBR5SgwT9ItIbc4p8glKE0fyYaFoFjKbynBpK0z+SWXILabmQ0ERx0yVtF2Gno0iqQZt2\nPm/d/EaAYuStp7F1P60NISK+MLploNm5gnibplXW1wRAkRUk5BnZJll3wjKw/MV7+T8beCw6rniP\nmcLDXl2zCoDmiIi8s1Ya3bCwMLBlvdA9sCXcREu4qWxfkiThI4Tp05a8t4drmwRUcWGGgyqDY+UT\nqWOasIQS/mNAUbw76zqQJIlbzrqZXx1/jCvWXDZu/65PmkgL8Q4rNSSMQfJmjnBQ5eipJJ/6zjOA\naPb1hkvXzfMRVh9JR7zDeju5Pc1IgSz5mmbaw62ASNM0LAPd1rCNWlobQoxIIn01Z6exrTCabmGY\ndkXBVVAxZbOwWMNU5KukKRV44u1RgXRWR5ElBnL9AKyJCPGuc7JINDtDTjeLvbsDlVe3dvETIu8b\nJK8trW2QL7FNQPQcyeZNLMtGliVMyyRtCEtomG6yRpZjiW6i/jrqnYnXLQ2b2dKwueL+3UdtV7xr\n1FowRK73b98QI358hLxu8tSePnoGUgt6rNVCIi/Eu7vHACRqlajIdrJUGgL1nEr3FdMEDR+tDWFO\npIvtFzAV8pqIrEvbwbooqKCYGJ7nXU4sFrs0Fos9PF+D8VgeJLNiCbCetChicW2TsBoCW8KQc6JA\np0JedyWCUhhJtknmlzbXuzBhWWKbQPGCHdMS2NjYNliYPHbiKca0RKE/91ScHnnXqsUJ3ku3tfHu\nN2zl3a8XPvlYhWXlzkQSuriJHe7K0VQX4LLt7diI3jqrIm2MaQkGskNiY8NPa32IsFoUb9tS0Azh\naZ9umwAokoo0k2yTMyFVMBaLfQz4DhCYv+F4LAfSjnifTJ2iIVBP2CcuFlmSUawQtuosCeZE3tGp\nxNtpVDWaT0663UJT9LydyPu0/OuRnLBMrJE2AB7oEnHJummLt9ifK951fhGtu9k5IETB75ML21Q7\nJ1O9/NMXpKyxAAAgAElEQVQr3yvrtT0Tks45kc0oXLillc52cS519SZpjwjr5NDoEQAUO0Bt2Cda\n7ro4nrdhWmXrV7qokg9ka9q2Se4MSRU8BLx1vgbisTywLJtMziAYthjTkoWo28VPCMmXZ2AkW+hl\nUheoq7SrAmFFTGCO5RMLM+hp4nrufrU88nYnLUfcEvdEI23BVeRMcXOaaeQ95lRouotRjJSIN4jm\nVm4VZ7XzQv9OXhp4hcOjR2f1/oSWRLJlMHxcFGsprFrU1ZdkVUTcRA+OCPGu9Yt2u6eLd04zMa3K\nkbcqq2LC0pjeBPEZUaQTj8f/B2bY5dxj2ZPO6diAWiMeZ93JSpegHEaSLXrHEmXrVU5G2GlUldSW\n1ucdF3k74u3mersia2shzm86t/C+dbVrp7V/94JPOlF1Q1Bk54zmysU7GvGTSGscT5zg6y//U6EA\npRrJOH60e6ObKQkthW34idYE2Lw2SntjGL9Ppqs3VRDvI2NdANQHxXkUVosLYdiWUnhyOj1VEMAn\n+Zw87+lJVbUsgQaLOGHZ0BBGVZf/3cylpWVyQapmJju2XJ94jPXVpsGGras2lG0fDdYxkoeRfLIg\n3hva22mpm3ifLXUNkAVDyS/o9zrVvv1BkT3T0lhDS0strc3iiUAN+GhpqSV33JkY04LcuP0KHuh5\ngFW1rXSubp3W57c2i5uTG8lvWr0ajkLGTpWNrbkhzOGTCXaN7mHf8AEG7T42tJw/p2NbKoxD4gnC\nF5JmPEbbtkloSSwtxCXb2mlrFU8qm9bUEz8+wpbV6+EF0CxxM1zT2ERLSy1rpObiTkwFdy4yHPKN\nG0PIHwQdlOD0xmc7hT6r2+toaQhPsfXELMbvNR/iPa3lJkZGlkdjounQ0lLLwMDS+rMLxVTH1n1S\n9C3J2GKSqNaqL9s+KIkT+sTQYGHC0kjLDEziZ6uGmBYZSo8u2Pc6nd9saNiJErMaAwNJTCcS7+1P\nMjBQw8nRAQB8Zhi/EeTdZ99KNFA37TFr2XIf28hYhNQg/cnhsn0EnAixe1Bk8/QODTPgn/gzlvP5\nOJISVtjg6NiMxxipV9EtHVuvJ+JXCu9f1Rhi37FhDhwaoz4QLcwZRAMRBgaSaKVSYikMJ0S6p2VY\n48YgWSJgHE4kpzW+sYTw7jOpHAOzXHd1vn+viW4E8yHeXrXBCsLN8c7II6i2Qlu4pez1+kAt5GAk\nl0CK5lEklaASnHSfdQHxuJs1ltYe0CqkCkKJ550fBUshEhA3qEtXXTSj/Z/ewD/gU8rEx8Vddcid\nIJ3tZN9yID0H22QsJwTO1v1Ea4qrFhV8717he7vf3+p60eWy1POe0jZxSuTzxvTmGHJnSoVlPB7v\nisfjr56vwXgsPUK8bVLWEO2RtrIqQoBGp+lU2kgh+fJElJop1/pzs1GyS7wgw7hUQSfbJFvINhnF\n1oLUhWaXQHV6hoLfEe+MkSVvFqPyqCPeCV1ErRljeXRcnA2u5503Zi7e7lyArQeI1hS/80LGSV8x\n4wSgo1GId6nnLVkKaadnzOm9TQACzlJomjm97B5NN1FkqWxFnuXK8h+hx6KSzhpIwQwmJmtOm6wE\naHKqLG01Bz6NSMmqORNREwhjWzJ5e2lFqjBhqY6PvDVTJ6WnsfLBisufTYfx4i3TEBifLij2b5M2\nnTS5KhbvtDP2/Bwj7/qSyHt1cwRVkUTkHRaTlrap0FovRF2WZEKOgKuSOmnk7XeWQtMqrIRUCc2w\nCjf35Y4n3h5lJLMaUlBMvLmz/aW01YroRw6lkCR7ykwTEFkYtu4nby9xkY5R/khcKt6jeWcZM21+\nxFtVJBRZLlRmlmacRCN+UHUsxHiyenXaJpZtFW48s7FNRnNO6qgeIBopRt6qIrOmpYZjvUnu/B/R\nGEyxAsglT3gRR7wVyVfoGFgpVdDtLOhOek6FZliFm/typzpG6bFopLM6UkCISWOwYdzrrbUi/U0K\ni6hpqtJ4cGwKPYDO0vb0ntA2yRsF/1mIt7/yDqagNDfYzSWvrxB510X8SP6iYFdr5F1q98wm8i61\nTeoi5TfMmy7rJNZRz8ZG0UM+Gix/wnMLx3xy8X1KhfL4gCp+S92aXqqgpptlS6ktZ5Z/MqPHopLK\nGgVhaQhGx70eVoNgycgBd5vJC3RAZFfYuh9bEpFa2FeegpXK6gR8yoJdNE+efI6fHfsVjcY1QIUJ\ny5whJiuZW+StyDI+VUY3rIKQ1zvf4cgk4l2tnrfrdwPkZuV5iwAgotagnBY179jayo6twu++9/Do\nuInzSEG8izfaSuePu4K8Pk3PWzcswpHZ3bwXm+q4xXgsGqmsjuyKd6B+3OuSJCGbxeySxgoCfzp+\nn1JYnMBtAeqSyRl8/Pan+MGDB+cy7ElxFwgeVkUVoBsVq4qMX5XJ5A1GckXxLl01fqa4loz76O1G\n3mMl4h0OqCiBotjlqjTbJK3PNfIWtslU7RXetOkN4zJ/XPH2l0beFSYsg6o47wx7up539UTe1TFK\nj0UjldVRgnkkpAn9bNUuzvY3hKaOvH2qXCLe5VWWXb0JMnmDvV0jcxj15LiZBkl/F5JEWfe5UFAV\ntokbeedDs7ZNoOh7uyLuTliWRt6SJBGIFMUkU6XinTFKIu9ZiPdIZgzbkqmPRKbe+DTc5lSFpc6o\nHHkHfa5tMrV427aNpnue97IimzfYeWiwahrgHzwxyqmhpUmrS2V1JH+OaKBuXJqgS0Aqind0ir4m\nICJcySnUSZ4WeXf1CTHvG86MW5JsvnCjQt0/jD+SL0ttDAdEtkK55z37yNsV74Kvrobwy75xud6+\noBiTX/FXreedLrFNZpMqOJJNOJkmM0/NvKDlHLY2nEUNxWrL060XgKDiRt5Tn1uGaaGujTPQXB2N\nUs8I8f7JU8f4yo92sX8Bo7v54sndp/i777/I7ffsWfTPtm2bdFbDVrOFx/1KhOTi5FGdf2rxBlBt\nYbWMi7z7imLe3b8wvU9Kc6zVxlNlr4me3iLyVmw/WOocI2/ho7u+uiRJ1Aej4/qbuJ53W6iFrJEr\nW4C3WsiU2CYzjbxt2xY3cj1QliY4XWKNm/nQq36fsK9o4VXq5x3yiRvxdGwTzbBQokNovuEZj2cp\nOCPEe98xIdoL+Wg+Hzy7r4/v/nRfoZ+xW+24WOQ0E1PJgWQXHvcrEVGdx1xbdHqbDioiWj/d8+7q\nTVb893ySNzVkSQZbgvqTZa+FgiqmZTOSG0U1xaP4XCLvwGm2CUC9P0pST5VlPFhqFlv3E5Ij2Nhl\nN5hqIe3YJhISeVOb0Q0oa2QxbbFQc2ma4EwpzcmulOcd9Il9m9OIvDXdAlVDsaujy/WKF+9s3ihE\ndwe6R5d4NBOz6/Agd9y7l6Bf4aKYmFlf7PEKy0REUG5HvEq46YGyFRSiOA0Cjk9e2lkwmzfoG87Q\nUCsultIofD7Jm5qoyku2YgUT9Kb7C6+FAyoouogcjVBhYeLZcrptAsWME7clrm3b6HIGWwsiI6LO\narRO3GwTP+KmN90qRijexE8v0JkppTfJSnnebqqgeVoDVMuy+elTx+gr6bmkGSaSYqBWyRIFK168\nD/WM4VrdR08llnwR3FIMy+DF/l2k9Qy/eLYby7b5yDsu4HUXihakSyPeQkQms03c9EDVCk24zen4\nnYZWpZ53d38KG7g41krAryyYeGumRkDxYw2L3uQv9O8svBYOqAULw8oHqAn7yopBZoq79mHAV7y0\nTs/1zho5sf6nFkA2fYW/VRtutkkm6fQPmZV4l5fGz5SAv/g9V1qMwV1B/nTx3tc1wo8fPcIjL/UU\n/pbXDCTVwFcl4r3i87zjx0dBzdPUMcrQ0VaOnEywtXN88clS8HjP0/zo4L2E1BAa62mMbmLzmiia\nbqIqkhj7IuJOVsLkkXdTOAqj4LOnL94BNYBtyWW2iSvWne01dPTWcLhnjLxuzntToLyZJxqIkh9s\nIdyp8OuepzmveTsdtasJlYi3ngnQOAfLBEo875L2x6eXyI/mi5OjtiPepf5xteBmm9ia8J1zZp6p\nE0cF7tyHrQeon0Ne9VSRt5uNYlJuQbrnXjZfDOYSzjJ9frk6xHtFRN5PvHKKb9+7p+I6dQe6R/F3\n7ifT8hJyff+SWye6YfHNu3fz4oEB9g8fKPzdbNuHvvkR0noGv09hw6o6jvcnC30bFoNS8Z4s8l5V\n24xtyYSZ/k0woIoS+dIJy+O9SVB07h/9F+T2Q9g2nFiAhXk1UxP5wJZKY+Y8ElqSL73wdZ469TzD\n6iF8HXGxXcZP7RxyvKHoeZfZJqeJ90iJeFu62K5abRPJlsFJA51JxkmpbRKdg21SepOs1EzKjbyt\n0yJvd36l9Ek85ZybfmnyLpnLhaoXb9u2ufvxozyzt29cpJrXTY4ODKA09gKgRAeJL7F4HzgxyvP7\n+3nwxeMcHD1CW7iF93Z+AKN/LYaS4p7D9wMQW1ePbcOhnsUbb1nkPYl4b25rY93IG7m+43XT3rdb\nIp/UUoWUza6+JIG2U4xoI+gBsVDC8XmetDQsA8M2C5V4bfo5vO/c30WRVL6/74fs0h9CCqVZq27B\nHG2dU6YJlOZ5l9gmwfL+JqV9VHTNFe8qtE2MDLLtxzadhZxnkHHizn0EpTC+OSzSUtqSoKJt4kTe\ntlRul7qRd75EvJOauIFO1eJ4uVD14t3Vl2TIaaD+Qry/7LUjPWNIzcdBEmLhbxzmcM9YIUK3liDv\n273BdCVOkDc1zmrYxNCoid61jajSxBMnn+XQ6FG2dNSDqvHA8YcZdFfPXmBSGSHekxXogCiG+Pjb\nr+Q126e3PBgUqyxN2yRjZNENk5ODafxtovGQpIrH2vn2vd1JNFXyOeOQOa9lO3+240Nsb9rK2eEL\nye+8irW5K8H0zSnTBCaPvN2Ie7Qkp1zPiUtwqUvkNVPnqy/dwRM9z0z7PRk9i2z5wRHvmVRZJpzF\nO2qn0dhsMvwlN8kJ17CkPPLO5Az6R8T3XRp5p53l6ILy9O3ApaTqxfuF+EDh3y8eGMCyioK87/gQ\nams3PsnP1oazsHxpdCXJsVNJTvSn+Ng3n+Suxw4v6nhd20YLiBtNrGEzJwfTYMu8YfXNAPwgfhdm\nuJ/gOU9w1H6Or7/8T2V9JBaKVE6Id41aO2GBzmxx+5uAeGQ+MZCGmiEMn8jA0O08qiLWLpxP3Ek0\nBVe8xXG1hlv4wPm3cVXLddhaqHAxz6U0HiDkeN6lXmyNL4IiKQznRP5woWBHD5LLiktwqW2TXYN7\niI8c4rm+l6a1vWVbokjH8GGb4lhn0t9k1Mm8qQ/OTbwDU6QKypIMlgJyUaS7+4sBQl4vWq3uvENQ\n9SLvBce2bZ7f30+grYf2C/eS0JMcPFG0GXYO7EXy57mk7SIuaBULysrRQR7d2cP/+8GLJMJxXuxe\nuJ4ap6MbJkdOipNWrhMX8ln1GznpVFNe1LGF16y+lFPpPu7YcyeST8NKNDKQHeLOPf85aR5tzsgT\nHz7E2CTLkU1FKptH8uUn9btni9+nYGVFcc99h3/OsVNjKK3dgIiOMkaGjtYIJwZSFecuZosr3qor\n3qc9ortpgf1OythcbZNzNzZy1fmrOG9TU+FvsiSzIbqO48keDo8eK0TgYaWGTFY86i+1bfJM7wsA\n9GcGp7V9zshjY2MZYi4BZmabjOWS2KZC4yxK40spE+8KRToAkq2U2Sal9QSltol7A42os1+7cjFZ\nsmyTbN4gkzNoila+y+mGyVAiT3vjxF9kz0Ca/vQwoc17GZNMgtv7eDTeQGzd5eiGxYC6Dwm4tvM1\n+JzHJyU6yBOv9CLX9xHo3M+wdoyscQ0hdeYzzJmcQSav0xyd3mPWkZMJDNNiVUuQkdoRamii1l/D\nqcEM0Ro/kaCPt2y6kb1DcWRJZl3uSp7cnyP22oPsHY5z7+Gf85bNN5Xt87nu/Txx6mmOpg8USoDb\nI21c0Lydmze+vqwUvH8kQ31D+cWSyelk8yZN0SBj+RRSrU1TaOJMk9ni98mY/etYd3aOnYN76DHz\nKA19tARaqQuGOTLWxTltNRw9laRnIF1YTWWuuLaJIpVXPrq4nQUHx4R4ztU2idYEeM+NZ4/7+5s3\n3cSXXvgGPz50H3kjT1gNEQyHGU45ayYuYeSd0JLsHxZBzJiWIO+kVk6Gm2liamrB856RbaIl5pwm\nCFNnmwDItoopm9i2jSRJBWtOoly8M6b4DUKqZ5tMynfu28sn7niKV45U9nN//OgRPnnH0zz4wokJ\n9/F8vB/fuv3Yksm5TduQfHl2Svfyz7v/g089/gWk2mGi9mraI600hRppDTejREdANmjcKjrMSf4c\nPzn40JTjHcoOl0W+w4kcn/2XZ/n0d5+dVk+OvKmx97iwSs4/T0GSLdRsCznNYCiRY3WTENWwL8yn\nL/sof3n5n7JjXQyQkI6/ipZQM788/gjHEscL++wa6uPO+J0cTO1Fy/oxejupM1cznB3m510PlW17\noj/FJ+94hu//bB9D2RHuPnQ/GT3Dt+7ezSfueIrdR4dIOstyNS6EeKsK2DI3tr2V9kgbg8ohJNnm\n6nWvJuITVYbtLUIw5jPjxBUU2XbFuzzyDgedCS3HbZtr5D0RG6OdXNh6Hl2Jbnoz/dQHotSF/WSX\nQeT9fN/LWLZV6AMyUBJ9j+UTHBwZby26fU0MTYUZ2iaWbZExM3NOE4SpKywBZBQk2cS03InyFEG/\nQnN9sEy886b4DWr81RF5L4l4G6bF3q5hDNPm63e9wt5j5b0ELNvm6f09yPV9/OczT3L38y9XXMbo\n6eN7UBr76Kxdx/vO+102aNdg2zYv9L9Mykggp5v5ne1vK2x/duMWkA22XdlN2hpltRzD1vz8uvfX\nheq3Svyy6xE+/dTf8e1d/4plW4ym8nzxP19iYDRHTjPZeXjqR82vvnQHD+W+j1wzglInblipgTpO\nDYmLwBVvEM2KZEnmnI2NnLepif1H0/j7he3zXG/Rk/xZ/Bkk2aY1eyGvr303q/OX0PfCeXTkrwRg\nX0kq4jP7+rBsm0dePMFPjz7AL48/wnd2/Tt7j4nf4Ws/foXhrLCcJss0mS2FXh+2nw+c916RXmb6\nuGzVhYR9ItIJhMSFlJ6kLcDzfS/PyBpybRPJEe/AaRd4KFAu5nONvCfjzZtuQpWKfb7rIn4wnAUh\nllC8n+19EVmSuWqtWI62P1s8n+869BP+8aVv05XoLnuP6w/bhg/bmlnkndRS2NjY2nxE3qUTlpVt\nExkVZBPdsMhrJqeG0qxrqyXgU8smLHNO5L2ixTsWi0mxWOxbsVjsyVgs9lAsFts4k/d39SbRTI32\nzSNIq+J87fl/4549jxdeP9A9RL7jCQJbXiKw9Xl+mfgPPvX43zGQKUbpPYNJxqIvgQ23xt6MJElc\nt3kHuZ1Xk9t1BYH4TfzFFX/E2e3rCu85u3ELAEezcYJKgNe2vA695ywMW+cnR34xbpyHe8b4l53/\nw91O+t7uoX38OH4/X/zPl+gbyXLJ2aJZ/PPxXr6394f84thDFTsXJrUUxxLHseQ8ga3P8eLgS2BL\nJAdqC82yVjePP2EUWeaDbz2H7esbOLTfh2L5ebF/V+EJIJ7ci21LvPey63nLlRv56LsuYF1bDbt3\nSWBL7BsS4m3bNs87E7vDqRQv9O0S3/PYQZQ1h7hwSwu2bWMo4kZSP0mBzmxxIyRNN4n668ntvpxV\nQzcQVIPFBWVlIdpZrXIV7InkSe7c8x881P0YOc3g2X19U3aKLIi3VTnyDviUsorKueZ5T0ZzqJGr\nO64ARL+TurAfkPFJviWbsDyV7qM72cO2xi1sjHYC5ZH3MUe0Hzz+WNn73L4mtuErRN7TFe9CoDTL\nplSlTDVhCaXibdI9kMK2xQr1Ab9MXrMK55BmO5F3YG4+/GIx28j7LUDAWTn+E8CXZ/Lm+PERAluf\nY6zxGZRVh5EbT/JA333cf+RBbNvmPw/+CLlmjM7QZq5svRpG1pCxkvzji7czlB1mND/G7bv+FTmc\nYlPoHDrrOgDYvqGBGl+EOqWRP33XhbQ2lAviWfWbUJzI54bOa9jY1oI5sIagVc9Tp56nJ1XsOJdI\na/z949/nuaGniKqNfGLHh2kKNvHIycfotw/z+ks6eN+bttPWGGZP7ime7n2ee4/8nB8fvG+coBwa\nFRaNOdKCIimM5seok1vA9PHUnj4AVjVVPmF8qsIf3XIeWzsayQ+2ktCSHBo9wvGRPjTfMIFcK+ua\nxORYJOjjo+98FWsb6zFTUY4mjpPRs/QMpukbzhCN+FEaezFsnWs6rkA1I/jWHOZVF5l84C3nFhYI\nWIjIO1AQb4uxlIatB2mNiHaebmN9WxFCO5EN5WZppPUMj758ktvv2cP+KapQCyXbVmXPW5Kksug7\nsoDiDfCG9ddyaftFXL56B/W1Qrh8UoDsElVYPtv7IgCXtF9IS0j8Hm7knTVyhTTVlwZeYThXbOxW\nyH4yfMU872naJmOa0+NlHjzvMttkghXf/bIPSYL7nzlWmKxc11aDX1WwbLtgp2h2HtuG2pUceQNX\nAD8HiMfjzwAXz+TNL/btQa4ZY1PdJj50we+z3fgNrHyQnx77BX//wjcYlI5AuoEPXfQe3nnOTdyy\n4Rb07i2MamN8+cVv8ddPf4lhjiOnm7ntVUVbxKcqfPp3L+azt11SUQyDaoALWs6hPdzKNR1X0toQ\nQpEVwsPnYGPzQFexj+9dz+xCbu3CytQw8PwFnDguYx2+CNtUCG7azYZt4iRYvyWL3HaUOrWB9kgb\nD5/4NT88cHeZP35o9AgAxqkNvLHlf9FZ18EFDWJlENffXd088d0+4FP4vZu3YY+sBuD53pf5efxp\nALZGt5VtWxPycdsbz8Yaa8bG5sDIoUI65S2v3YS/rQdsuKzlcrL7zwdL5n+67qJ1lcF5W8VJO1lp\n/GxxG9znDZORpHOTcC5ct7G+KQuhzWmVxdst7Mib+UJu/1h6csEoNEuy3BVuxqdAhp21LCNBdUIB\nmC9Caoh3b7uVjdFONq0WN0nb9C2ZbbJ3KI5PVjm3eTvNoUYkpELk7QYzjcEGLNvikRNPFN6XLrFN\nZprn7eZ421pALMY8B0oj70or6QCsbhLf8wMvHOO+J0Qg1dleW3iv63vrdh5M37y3Z1goZnum1gGl\nDYqNWCw2rX2ZpsUp9WWw4Z1b38zWxrO49dIdGPsvQ9ZqhL2QD7GdGwj5xQ975XmriCS3Yp3cwmh+\nDN000Y5u5zfa30l9qHxh0ub6kPASJ+C2c36LT136f/ErPlRFprUhxHBPlFWRdl7s38VIbpREWuPZ\nQSGOb+i8DsUKcsd9eznZIxOzrkVRZP517w+4fdedHJQexbZk2hJX8OFXvY/2UBuP9TzFffuLj5kH\nR48g2QpWup5LN27mTy/+EDdsenXh9UhQndJrbawLcvn6bdhagOd6d7EvsQfbknhDbMe4bde317G+\nRjhZTx3fzQvxflRFoqPDhsgIZqKJR58bwUjXcUHgdWSNHN/Y+V0GtQFkSZ7WivAzxV8SeQ8nnSpO\np5ug63mbiIu/tN9EKUndFW+NVEZYLLkJtnVxBcV2xbvChRkOiO++ZoEmKydiw6o6VEVGz8tkjKVZ\nnDmhJakPRMX1IKs0BhsKkfeJpGife+P666jz1/JY99M8vU8kEBT6mpi+wo1xuqmC7hOUaocKFamz\nRZalwgo6Ey1fFg2KwKipGRIZHZ8qs6opXCioyjs2nUEe2/BVPEeWI7NNFUwApVe4HI/HJ03ObWgI\no6oK9+18EsJjNLOR8zecBUBLSy2vuyDGAy9KNG06yfCxVm545xZaWoof8fbXbeGf7tF4dUeMl17J\nUCvV8JvXb53zXXLDmihP7jrFdRuu5nu7f8AzQ8+i7VoPjSeoUaLcdvX1XLRmkL/7t+e48oI1fOCW\n8+nPXMm3n/s+u/v3AxAcOp+DJ21q3hwlF78Ye+39PHD0cV63+SpamlV6Uqcwk41sXtvIWRvEo2lz\ncw31tQFGk3k6V9XR2jr1ogbvfuM5PP3v7WhtXaDmCGZXcfG2DRW3/b0bruAvn3yIvcNx0gOt7NjW\nzv7sXgDMgbU8dEBchO+58gaeHYzwH7vuBqAp3EBb6/zbJi3DIlLz+VV0W0RI69fW09JSy2pT2D6+\nsBAvE8p+exfzhIiiLdkkZ4jTTVaVwraV3qOKzggE1CCQorWlZtx20doA9CVprAtW3MdCsnV9Awdy\nMkrQprbBT8hXOXV2IcYlCm3StNe2FPa/tr6Nnb37iNSrDBwV2VEXrT8bQ83xX7vv479feYTfuOr9\nmEecSWXDB0hgKpgY0xrnnl+LH2VdU8u0zvupCPpVdEOjva2uLDXW5VXJs3n61PPcdEMtD/1coaOt\nlva2KNFa8V1HasXvbkoaGDWsaq+bc9bRYpxHsxXvJ4CbgR/FYrHLgFemesPISAbbtrn3wM+xbbi0\n8QoGBopZA9e+ajW/evY4Q/s34PfJdDSFy16/aHMTPwz7ePIpHfDxG9d0kBide9Vho/vonl5Lnb+W\nXx56nFzvEaQWi+s7r2JoKM3axhBf+eMrUGSZoaEUCkHev/1/80zTCyS1FEOs5oGjJ/jjLz1MOmfS\n3LaWdLCbT975U7ZsDIEMQa2FP/yNbWXH1NFSw2gyT0s0WPb3iVCA7dHtxOkCYEvt1gnf11oTJGKs\nIhPsRgpkOLuzll8cuYuIL4yZase0YW1LDT5sXt10OSfW9PFYz1NEfXXTGstMyTr2xshYRjS9BxTb\nZmAgiVOVzEgqgSRFSCTzFcfQNyY811QuQ25M3AwGhtMMDCRpaamt+J6RhIjWM2mrMI7Tt3OzFII+\neUGOfTI2ttcSPyUuw+7egYqW1UTHNlfSegbTtghKocL+o4poNrav+xiHBrrwySpqLsTZofOxzZ+S\nrT1Ef3+C4ZRT4m84qZamSiqfmXKc9z5xlCN9/SgN8OFbLpuX4/KpEqoiMThYOcV0jU/Mie0fjPOX\n7wUPNkYAACAASURBVP1dZEliYCCJZYqIu7cviWTp2JKJbaokRjPkprDjJmO+f6+JbgSztU3+B8jH\nYrEngC8BH5nOm/YNH2DMGsAcbueSjZvKXmutD3HZ9jYAztvUPC6iDvgUXn+JyBypCfm4+lWrZzn0\nctwsj/6hPK9d+2rx6NdyBJUAV6y9pLDd6evjSZLEZasu5vrOq7k4Jsadzhlcd/Fa3r3jOgDytcfY\n1ScyPn7r8stprCuPqjrbheUz0WRlJd5+6Q6sXAjbkrkhdsmk217WcQ4AvjVHeDD1HyS0JFdvuJzz\nNogsmYudRR8kSeLtW97Mjeuv4/Wd1057LDOhdMKy4Hk7tok7YZnRs4T8KtkJPW+nmZCZJzlD28Q0\nxO9XaXFZd8JyIdMEJ2JLRz24bWEXOePE7aJX4yuef61h8WR4KtXHqXQfq2tWocgKhqZgjrQjBTLs\nHzpMWs8iIYGpIgG2qUw5YfnU7l7ufvwovqCGKql0NDdNuv10CfgUlEnmKhqC9bSFWzg4ehiwC9F5\nqeedKcmeqZbV42cVecfjcRt4/0zfF1JDkG6gZvScilWJb75iAwOjWW7Y0VHx/ddcuIZ9XSNcuq2t\n0Dd5rrj51SeHMty0dQf3HfolyBZXrrmU4DSrLjeuruOiLS20NYa55bUbsbGpD0RJt/eBFsKUZM5f\nvXnc+y7d1s6B46NceFZzhb1WZlVThFfX3MxIOs3GtslP/qs3X8BDAz9Dae5hJC9x/bqredd5b+ZF\nqY9EWuOK81YVtpUlmZs33jDtccwUN8tD08WEpSxJTqochVTBtJEhFFAmFOSC521oBfHOTFEg5Wab\nmLoj3pN43gtVoDMZm9ZEkZZoQYak88hT6y/OG7WExDn1yuBeTNtkbY0IksZSGubAGtTmk/z6xLNk\n9Iyz4oxEXY2fnKUW8qQnYo9TzxGuNQmotRUtjtlw6bY2RlOTLwQRaziLx3qe5GjiOJvrhdVYJt66\nOJ9kyz9v41poFrU83qc1kt1zKRc4EfbptNSH+MRvXzTh+4N+lf976wXzOqb2xjAScGowzbO7R9B7\nO4msOcV166+c9j5kSeKDbzu38H8JEZX//NiD4NPYWLcef4Vy4zXNET4+yfFOxO9cNX6SshJNoQbO\naTqbpJ7i1i1vobOuA7/iY9Oa6Kw+dy64WR55J/Kur/Uju3aFGkRCIqNnCAZURpOVI7hUSbaJ2/9k\nqurWaYm3k22ykDneExHwKTREakgAo9kUzH+iz4QUIm//+Mh795CYz+moFeI9ms5jJRux8iH2jOxB\nVVSxYDMia+ikqaBbOpZtTbg0nii+ssmYaVrClQO02fCm11Se9ylla+NmHut5kvjwwXHirelm4alH\nthb/Bj5bFvX54Ogpkd+5pWMRz9Ap8PsUWupDnBhI8Ytnj+MbOJtvvPHzc27OdPmqYvbkWfUzqmGa\nV95//nv504s/VMiFXyrcyDuvm4ym8gXLBETUH1ZDZAzHNsmb4zIvLNsqRIpCkMXrE6UVuripgoYu\nbhSVbJOCeM8xbW22rKoX51rXwOKuWu4+ydT6ipF3U1CkC+qWiERLI2+QMAfWoNs6WSOL5AhdQ21g\nWumC6byB7BMCHw0s7sTwloZNSEjsHyk2ois9J91y/2pZvxIWWbzP7mzg+os7uHRb5ch7qVjVFCad\nM0hmdF53UQfRyNxbQjaHmthSL3z9pRTv5YIb8Q6OZTEtm4ba8u847As5kbconNCM8uSlrJEr5M7b\n2CCJf0/HNvHJKrrhjmP8KX/ZtjZuvHQdF57VMqtjmyvrmsUk4YnhRV72ThM3w9LIW5EVQpKbASKx\npkZYa2OOLWEOFueaJLMo3sUS+Ynti3RWJ1Qjfog6/9yzTGZCSA2xvq6DY4nugj1VtE2sQuSt4kXe\nFWmOhnjXdWfNm189X7gFMgGfwvUXz1+E+ptb3sQb1r+OLQ2bpt54heNOAvU5KYMNp1XWhdUwaSNb\nyPvNnSbKpavOA6AIX/x02+TUUJpUSW+UvJknoAQK64KePvEMwut++zWby1ZlWUw2tAqfuXd0bIot\n55dKkTeAaoj/B61owe4bdbIvbC1Mo7xG/Fv3ociSmCuYRnOqdM4gEBG/VzSwuOINsLXxLCzbKjTa\nKoi3ZhYqRn1URy9vqPJ+3vNFR+v/b+/Mg2TJrvr85V5r791vm7dIb+blzGiZQRosBpBGQhZYLMYC\nB8FiWzEDOJCxARljW0R4AWOMbYzBNkEE2BNG7GAbImxsS4gQUggMQrKQhUZKjTSjWd7Sr7fX3bXn\n5j9uZlZVL6/7dVd1V1adL2JiXtfSfW9l1i9P/u6556iT9S2vuzDQRasLlXN8wyu/ZuCNDfKIrmnY\npp7tZuu1TUBF3kEU4CQP76xvslO8NV2JQO+GniCM+NFf/Bjv/d+fzR5rhx1sw6btR3vurhwFZooq\neNho1k60u9NekTdA1FI/G52udbjZsyA4H6jF99A3KTomBds4sCxsHMfUmz5OQf2e6SFsBDsId1bt\nK0kLttn2bs/b0vJjm4xWCHxKPPbgEpqm8bprh8/6EO4d2zIyO2Ruqv9LkqYLmvbeEXUaJWpoyjbp\nibxTf7zW9Gl3QrYa3ci7E3aYcqrUgnBPy2QUSOtHx7pPqx1kZWqHzfYeqYIArS0HyhDWugJ7p9ZG\n01TpXLt+H1//2q/hfb/nUyqY6o7lgPombV+VZDUK6thMnULk/YrpSxTNAh++/n+43VjlwdLrgFiN\nraMib0eXyDtXmIbOGx4+c6xGqMLB9Irnrsg7qW+iW2n+dr94p5kRc8kmFs0IKDoGYRTjJxeE1C7x\ne/zydqQi744fjuy252LSdkszAmqtg2vDD4qaX6doFrM+j6BEdvvGEsGty9Sun80ujJu1DmeSQm+1\nesjbX/FWWlsFFXlbxoFb5OtNNS/NUn7zzCmIt6mbvOu1T/HAzCv57Maz/M7138BYepG2H1LzU/HO\nRyMGEPEWTpBe22KX553UN9Et9SVv7WObzCd5yOghSzNKTNIovb5DvMMoVFaM4dAZYdskFW+M4K61\nzAfNdqdGdYdlsnKnCaGF/+JD+G2TrYavbIV2wNyUQ7lgslnvEIQRnSCi5KjI+yDbpN5S84pMJd7D\nqJ9zGK7OXOEHXvc9/P3H/g4AxuxtOn5EPYm889I5HkS8hROkN/Ke2RF5lxPrIO0iv3OX5VZimywU\n5tQDRsjSbDF5rRL6dONOmgOeZj44hk0nCPsK948Spm5iYKKZfiZywyaKI2p+ncqOxcq0EXPaD3Jt\ns8VmPfGpyw5TZZvthp9l+ZTSyPuABcv0ohTqLXRNz2yy0+Ly1EWWiovolTu0fJ9G0CCONBzz5HP9\nj8pons3CWJLaFlNle1fp1VJa01tPa3r3R96pbbJQVOJtGGFmvaSRd63VH3l3IvW7bN0iCOORtU0g\nWSgzgsxeGDYNv0lM3Le7ErrinZarXdtqZYuVMxWb6bJNrelnVR2LBRPbNg7splNP7KAODabs6r4b\neU6SK9XLaEbIdrSuOgMFFraZn2XA0/8EhYkhbUG20++G7oJltE9N7+1ODQ2NuYLKiS4U46x5cGOn\nbZJG3kkUaGoqg2ivDTqjgmMU0IyTi7xr/t6LlbfvKPF+6Ir6nFc3m9ypqc9xuuJk2VjLG8pm6Ebe\nd+8gr+YV044bp5ImuBf3z14BoKYt0wybxKE10ufITvIzUiH3pJHvTr8buhkXobZ35L3t1yhbpcwf\nth0oJuKdLm7uXLBsJ5G3qVl9f38UKRoFMIO+HPVhkq4h7I68lSg/dFmJd69tMlOxs1r5ab5+1/NO\ncqb3s01aAZg+ESHTJ7xBZz+uzlwBoGGu0ApbuarlDZIqKJwgmXhP7R95B2lDhl2Rd51ZZxoDJcSW\nHWfinUbe6a18Jt6BEh2DvVugjRJlq4TWjtlq1U/k76WlBnZF3htNpit2tnFtbbOVfc7T5a5431pX\nIl8smGrT3UELlk0fzVLPTZ3w1vj9WCotQGDRdpZV+mkgkbcg7En6xZjbwzZJUwXTJrC9qYJBFNAM\nmlTsSlba1bKjrJRrGqWnUWsQqqayqeedCv4oR1VVR0XAm+29a1IPmtoekXcQRqxttTgzU1R2iG3s\n8Ly7bcuW17u2iWN1Pe+72SapeM+MSOStazp6a47YTHaPhhbWCF/gd5KfkQq5J4u89xLvJFWwEyXi\n3ZMqWEtLl1pl/I7KgjDMcLdt0uMXB2Hc7Rwfq9c5I5oqCDCTiHe663HY7BV5r262iGNYnC2iaRrz\n0wVWN1vZ1vjpip3VPL/VI96moWHEB4h3M0CzRyvyBrDbPWWVA2tk00n3QsRbODHSSHmuujuX1tJN\nbN2iGSrx7t1h2evPdtrqlNV6xHu1vcpnVp6l1pOp4QdR5r/q8ejbJjNFJWh1/2TEe6/IO/W7l5LN\nOPNTBVqdkJurDRzboGCbmW2S+uClgommaaoGSnw3z7sbeY+K5w1QCLrFyOKc2SbieQsnxhOPXKDo\nmFy7tHdJ4JJVohk0sS29r7ZJr3gnjcfR9K54Pxv9IT/2oRWCxtuy9/hhlC1YEpnA6O6whK6INqPj\nt/Y7DN3Iuyvey0ma4Jkkf35hWl1k17Za2WPTO2r/pMegYJtJQ4b9FyytyuhF3uVono0YtfVfbBNB\n2Jv56QJvf8Nl9H06lZStEo2gQdE2+zzvTLytCq2kWUusB13bhDp+6Pd1ovGDkE6yYKmlneNHOKpK\n7Yt2dDKt0Lot0LqbZVYS8U43P81Pd++QppMMoZ01z0uZeBsQmftG3rVWB21qFduwOVsanZLQBatA\n3EjuBMQ2EYSjUTKLNIMWjmP0R95pTrJdodYKiEOdWAsoJTZMoCWibezwvJPIO07Fe4Qj77SyXydu\n7WpEMQy2/Tpls9RX8TLN8V6aScS7p+fqTEWJtmMZfaVz00YWtqXSBfeLvBvxHWK7zsNzLrYxOrsY\nHcsgqqk7wbzZJvkZqTD2pLssC4Vw78jbrqh0wMgkwMc0dAwdIi0RDHOH550KSTj6nncaecdme1cj\nimFQ69R2lYJd3mhSKVpZVcO+yLvcXWTutU7S2vwFyyAOjD1TBf0gIqzcBOCRxVcNbhIDwLEMgltX\nWPBdoq25kb7A7+RYZ7Pruu9wXfdXBjUYYbJJ65vYBVX0KK1Rkor3lF1hu+EThwZB7KNpGoVSBIkL\no/VE3mrBMom8Q3Waj3K2SSremukPvThVFEfU/Uaf3x1FMat3mpllArCwR+QNUC0rcS86RtaH1LEN\n4sDGjwJaOxopN1o+xuxtiDVeNf/gUOZ0VBzLIG6XmL7zeoiN3HSOh2OIt+u6Pw38c7KvjiAcj1JW\n07u/smBmm1gVthsdiAz8xBIplLr2imb62ZfPD8IszzsMR982cQwHLTbQzE5WB2RY1P1GUtekJ01w\nq0UYxX3iPVW2s89zuke8p5LIO11zAOV5x2313rXWRt/fu7m9gV7ZZCo+e+oFqXaSdtPZTi6Yk2Kb\n/CHwrkENRBBKSeRtZOKd5G93ali6hWPYbDd99NjKcritQo/QGX4WIfph1zaJss7xo/vF1DQNWyuA\n2Rl65J01YehJE7yxqrJPzs13BV3TNOaS6Hu6p6RBulGn1CPeKoJVx2+12d9I+c9XnwHgjHFwl/eT\nJj0n0s/cGuEL/E4OTBV0Xfcp4N2odt1a8v8nPc/7Ldd1nxjy+IQJIo28VU1vk1ayc3KrU6NqV9A0\nje2Gj6lZBLGq1W07XfHWzIDZQoGVOy2CIKYdJmVIgzTbZLS/mAW9SMvaGHpxqlrWu7Ir1DfXlHif\nn+/3wRemHJbXG8z0ZJmkxan6xNs2iNrq+K21+sX72W0PgEuF+wc1hYGRLr6mkbeTo8j7QPH2PO9p\n4Onj/qHZ2RLmiH95ellcHJ1c1EEzqnO71FkCD8xiGyjglGymZh0221s8tHg/c/MV6i2fOcMhAKoz\nNnapZ3HP8Dm7UOZzL92hWLaJ6gGWbmIknWLOLFVHdu4AFafCZrRKaOw+RoMc97NN9Zmdm1vIfu96\nsgX+1dcW+/7WE6+/SBDBww8sZbbThTPq+empQvbauZkScUdF3g1q2eMNv8mN1gtE9SleeeX8UOd1\nFOaThh7txKI7d3aacvH42TAnMa8T26SzsXEymw8GweJilZWV7dMexlAY5bmVQ1VDejNcBaa5ubxN\nI1JFgxadRb740jpxDFZSJfD67TVCrbs4ppk+xeQ2eG29Tq3dxNZttjbVaxq1FisjHFgVNCV+L6+s\nsLKylD0+6GN2fW1F/aNtZL/3uZc3MXQNM476/tZjDyzw2AMLbN7pfn/1JJXR1MheG/pBZpu8vHEr\ne/yTK58mIiK8s0joh32/exTOxc6Odntbmw0ateOdJIOe134XAtlhKYwMc4VZbMOmHq8DV2l1Aur1\nZQA+8JFN3rfyEUB1xgHVKUczu13NMYKs3KwfRHTCtH+lijRHecESkl2WDdhsDbc4VdY1Psk2ieOY\nm2t1zs6VMPSDhSvdIt+/YGlCYGFpNmvN7oLl9doNAKLaDJXi6MlNb866poGh5yf/4lifpud5HwI+\nNKCxCBOOrumcK5/hxa3roEU02wHLm+rL74QznLswjalrLCzMcHtDlR+NjbZahUFF3tM94t0O21Tt\nKu1A3RKPeibBdEH5zemC4rDYbG+pv5c0RdjYbtPqhJybP1wmyCvPT/HG157jK15zLntMLfxplPVp\nVlvrxHGMpmncSC6+UaNKuTA6m3NSelvj2ZaBts/u31Fk9C6FwkRzoXyWF7ZeQivUabZDnm9dB+Dx\nq/fz7V+lNnj8z+fXYEPV6w41Jd5xrGFYQRZJBWFMJ+zg6N3Ie9RzeGcLSkxrQy5Odae9CXQ7uN9c\nU5ZIWsP7IExD58mvfajvsXSzTlGrcidYoebXqdoVbtRuoccW+M5AvORB4/TcjY36xX0n+RqtMPac\nr6hoTi9u0+oE3G7eJmoXuHp2IXuNY6jouh228bWW6uLi22hmgJX0xuz4Pn4UqObDfoht6SMfVc2W\nlLfZDIa7PnSnvUnRLFBIuhKlaYKHFe+9SEXQidUcVpvr+KHPSnMVO5gBNMqF0YsV7RyL9+h9msJE\nc758FgCtVGO706Cl1YmbC1w+21206fW8O3GLOLAh0sHubtJpJXngtmHTCaKRTxMEmE5qereGXFlw\no73JjDOd/XxjbXeO971SSO547FjNYa21jqkbRHGE1prCMvWRXHPoi7xHcHx3I1+XGmHsOV9R4q0X\nt9noqKwIvTPVt/MvjbxbYZt23ATfIg4sYt3HNLTsOfVaFXk7I7xBJyXdNJN2ExoGraBNM2j2iffN\n1TqaBmfnind5591JRdAIkvZpzXVu1G4BEDUqIxl1Q//GrVG31XaSr9EKY0/VrlA2y2ilGuv+KgDz\n1kJfGdk08q51aoRxoCLv0AItJtZU6lda18QxnMQ2Gf2oKq1v4mvDE+/NxO+eTcQ7jmOur9ZZnCli\nHePuJI28db+7UedGXYl3Z7s0kouVsNPzHv1zpBcRb2HkOF85i+40WfFVpsnF6fN9z6eR92qyky8O\nbOJARXZhUmEwrWviGDbtnNgmJbMIMUR6OyvKdVxu1ZezBUpQlgmQRd7bDZ96K9i1s/JeyVLuko06\na82NLPJubpZGNvLuFW+JvAXhmNyXLFrGU7eIY3CXLvY975gq8l5LamjEfhJ5A36sRLvT63knC5aj\njqEbGLGNZnVo3GNxqiiOuF5TZVe/cGOTP/jEdaI44qc+/nP850//Wva6LNOkoMQ72xZ/jMVK6Iqg\n39Gp2hVWk8h7yp4iDq2RzDQBJdjpPZ2Tg7uzXkb/jBYmjjTjRDNC4naJq+dm+55PI+/1tHpdoDxv\nAJ/+yNvULOI4P4tRFkVVFvYe65v83+VP8uMf/bd8dv1ZfvePXuC97/N4fmWFetDgpe0bWYOHbpqg\nakDQLUh1vGp/qfC1/ZCFwhzrrQ3utDdZdNRO0VG1TTRNw07uGiTyFoRjcr7S0yarVd0lLKnnnVoA\nyvNWt+V+rMS7HauuMI6u0uHykgZm66qyYK1xb+K93FCLuy/XbtBItnx/6qWXAWiFLbY6arv2xg7P\n+8Y95njvh65r2KZOqxMyX5wjipXtM2upFM/yCO6uTEkvPHk5R1LyNVphIjhXPpvtmqxq87u2bKeR\ndyoQFoUs8k4zNdoklfNMFWHm5Za4aJTQNFhr3FttjHRjz3prIyuy9Pnby9nzz61d5yd//RO8kNQ1\nmXGmieKY52+q3ZbHjbxB+d5tP2S+MJc9VtXngdGNvKG7yzJP5WBBxFsYQRzDxghV2ty58u5mtWnk\nnf2sFYkTz7sVttAAX1PiXdZVfngePG+Asqki4PXG1j29r0+8fSXe1zdXsuf/5LnneOaLG7ywfhtT\nsygYDr/0Po/nbmzx4KWZbIfkcXAsI4m8uzaX0Va7OGd66oGPGhJ5C8IAKYRKAO6fv2/Xc6ZuYmjd\nKKlolCDJNmkGLSxTxzeUmBU0dRHIQ7YJdNMFN1v3GHknxabWml3xbkTd3/Hc+nU0DTSrjd+w+Q//\n7c/50J/d4NJShe/9ptcMZOwF26Dd6UbeGhq1TSXax7Vlhkm6HpKXC3xKvkYrTAzng9fRee7VvObC\n5T2f742+S0Ypi7wbQRPT0AkMVVsjD53je5lKdllud+6tvkk38r6TdSDSnG6++Fa4wbWLVTSrA36B\nP/v8KvctlvnBb310YJaGY6vIe66gLrxLpQVur6o1iEHYMsMijbyPk+d+GozuKoIw0Xzto6/ivi+c\n5dKZvWsZO4ZDI2iiazpf94YHWG4s8t/X/oiG38Q05/CNJnOFC3SSKDQvUdV04Xji3QpbdKI2juWA\n3UKPTRy9QL1Q58ELRV6swcMXzlMtnOcdb3xl1hVnECxOF/nC9S3ClsN9lfM8PO/y4Y83mK06feVj\nRw2xTQRhgFy7OMNfffPVvp2VvaSRd9ks8SUPLPLNb1JV7ppBE9PpgB4xV5jt1vLOSVQ1X0rKtDa3\n8F7c4LkbW4RRfNf3xHHcV4kwshq84lwVzWmCXyRuldGdFkvn1GdxeX6Jd/6lB7O63IMivdC+tNzg\nPX/hB/jq+97Gxnab8yMcdUP3wp6Xu7MUEW8hl6QZJ5WkA3rBctDQaARNjIKyC+YLs3SSWt55qG0C\nMF9WKXwrtU3+5a9+gh9778f44MdevOt7mkEry7wB0OwmjqPqm/tNm9qG+qyut74I0FfXZJCkxcNe\nXFaLxVnBqxH2uyG/kffo3ssIwl1wzES8kwU+XdMpmUWaQVNFnKjOPO1GapvkI6qaLSoBnD/XZGph\nhedur/PBl+7QLFziXPkMl6cu7npP2lC4YBRUto3TRLPVBSxuF4layor57Pqz6m8MS7zPqL/zwrJa\nKL25Opgc8mHT9bxFvAVh6KS2SaWnA3rRLNDwm2Cn4j3D7Zy0QEup2hVMzWAzvs2meRvrPDwbPM+z\nn/ljAH7k8X/IQnGu7z01X4nkpeoFPnfnC+hOk9hSn0HcKRA31Wd0M+lqM6zIu1SwWJwp8MKtbeI4\nziLv49ZNGTZpXZa8WGsp+brUCEJC5nnbXWEoWUUaQZPYVMI168xktkleboltw+JvP/pd/PWHvoV3\nvfZJAu9x5tbeyOuXHgFUoamd1JK2aZemVFqlZreILSXodlzh8uy5vtendU2GweUzVWpNn43t9kCa\nPJwEaQ76dGWwawDDRiJvIZeknne1L/Iu4kc+pqVu26fMaTr+bSA/kTfAA7NXeSD597yxydbNgEfe\ndJGP3/4kq0kxrl7SyPtsaQlTM4mcJr5ehxje+dZHcWev8qMff7/KgddNyubwFhAvn63yMW+FF25t\nc3OtTrVkURnRolQpTzx6ngfum+biUuW0h3JPHEm8XdedAn4ZmAIs4Ac9z/vjQQ5MEO5GFnn3iHfJ\nVOVIO+Yd4sDC1OzcpQruZH66wK31daaSbf6rrTU++pllojjmyx5WjStSz7til6mY02w4W0q8Q7gy\nd4apksOZ0hJf3HqRGWd6qO3gLicZJ89e32T1TotrF2eG9rcGhWno+6akjjJHPaP/LvABz/PeDDwJ\n/OzARiQIh2CvyLtkKfEO9TZxu4AfRFmqoJMzPzNlfkoV1sJX81xtrvPL7/8cv/L+z2WvSdMEK1aF\nsl5VVQlRFRdTi+RsSVX3G5bfnZKK4Ec/s0zM6FsmeeaotslPQVJ7U0XezcEMRxAOx4XKOXRN50K1\n26ihaHbbeMXtIn4Q0Q7yH3kDNGpqQfZ2fY1a8xIAzXZA0TGzrfEVq0xRU3nim+EKU3YVS1df8TPl\nRaBbCnZYTJVtZqsO61ujv7My7xwo3q7rPgW8G1XnTUv+/6TneR93Xfcs8EvA9w11lIKwg0cXX82/\nedM/wza6fmqpV7w7Rfww6rFN8hl5LyTivb7VZqEwx836CulXcW2rxX2LlW7kbZdxkgbAMXG2TR26\nBb5mh7hYmXL5TJWNbSXeEnkPjwPF2/O8p4Gndz7uuu5rgF9F+d0fOej3zM6WMHN067q4mD8P7LCM\n69yWZrtiFbWLVKpFtKSc7Pmz0yPbzeVuXL2kmko0/Ijzs2d4qXYDzA4EDgEai4tV2nELUze5eHaB\nqj2T3ROfnV7IjvWb5h9jPVrjLa94nPnScI//w6+c588+r/qPvvraEvPT+zc2Htdz8STmddQFy4eB\n3wS+xfO8Tx3mPRsbjaP8qVNhcbHKysq9VXXLC+M6t8XFKmGruxAXdwrcXt1mu66UbGuzQaOWP+vE\nTHZOvnRrizPzyhLRCw2imsNzL21wZbHMRmOLsllidbVG0Oimu5Wp9B3rJ5beSFSHlfpwj/9CVa1H\nFB2TsO2zsrJ3S7dxPhcHOa/9LgRHPZt/HHCAn3Fd94Ou6/72UQcmCIOitMPzDpIFS0PXMI38CTeo\nHGRDVxZJujkn3UG6tql2Udb9elYmwAi6NsVs4XQyPdJt8ufnS0PNbJl0jhR5e573VwY9EEE4Jd1j\n8gAACypJREFULmm2CSSed6A877xte+5F1zUWZoqsbrZYKKiMEbPUIlyD1c0WfhTQCtvZTtPIt4lN\nDU2PmTsl8Z6tOnzrWx/g4qL43cMkv2e1IOwgzTYxsCCwkmyTKLeLlSlLsyU2ax0cVEQ7Petn0Xjd\n72aaAHT8iLijPofTirwBvvpLL/LQlbmDXygcGRFvYWxIbRPV+kzLsk3ysjV+PxZn1byu3wiJY9AK\nTeamHNY2W1nd79Q2aXdC4lYJDa0v20QYP/J9VgtCDyWzyIwzzYKlanmktklemg/vx9KsypV+5vlN\n4k4BX68xP1Vgs97J2qWlkXfbD/Ffcnnnw9/WV7RLGD9EvIWxwdAN/tEb/h5fNv1WABV5B1FuN+ik\nLCWR958/t0bcLtGMasxOq+Wq5a07gNpdCSry1ttTfOnZR09nsMKJke+zWhB2UDAdCpZKl+v4EX4Q\n5a7U506W5lTkXW8FxC317/KUSr+7vb0J9NgmY3CnIRwOEW9h7EjTAhstH8jv7sqU1DYBcGK1aGkU\nVbrgenO3bZLWpxbGGxFvYexIUwPrTRWd5t02WZgpkmZLzzlqETK21aa3XZ53J8z9xUo4HPk+qwVh\nDzLxbieRd85tE8vUs0YBZ6uqwFRHV6LdW9cEoO1HFES8JwIRb2Hs2Bl556X58N1YSOqDXJlVG3Xq\n4SaaBq1QReBls0QUx0l2Tf7nKxyMHGVh7LASz7s+Jp43dEvDXlmcp2AUWG7cZrpi06FFySxi6Aa+\nHxEDtnjeE4G0QRPGjm7knYp3/mOUr3ztOaIo5uqFaR7avsYnbv8/qvd9npbepmKpglXtpPyt2CaT\ngYi3MHaYWeSdLFjm3PMGeNWVOV6VbDf/NvebeHHrZdZ4Bg1wdGWppOItqYKTQf5DEkHYQRZ5j5Ft\n0kvZKvE3X/M30JPYy0JZKu1O0nhCbJOJQMRbGDtS8Y5j9fM42CY7ua96nsfKaiepFSW7K8U2mSjG\n76wWJp6dJWDz2nz4IF6/+CitTz/OhfBLALFNJg0Rb2HsMHQdvacJwDhG3qAyUOL6NFtb6hYjs01E\nvCeC8TyrhYnHNHvFezzFbDZpN7ZRU63eMttEPO+JQMRbGEusnrZnea/nvR8F26TomFmn9pbYJhPF\neJ7VwsTT63uPa+QNMFd12NhS4t1JbBMpTDUZHLV7fAn4VWAWaAPv9Dzv5iAHJgjHYVLEe7bqcH21\nTrsTyoLlhHHUyPu7gY95nvcE8CvAPxjckATh+Fg9GSbOmNom0O97i20yWRy1e/zPuK6brghdAjYG\nNyRBOD6mMf4LltAj3lstOp0IENtkUjhQvF3XfQp4NxADWvL/Jz3P+7jruh8AXgO8baijFIR7pN82\nGf/Ie3273WObjO98hS4HirfneU8DT+/z3F90XdcFfhe4f8BjE4Qj059tMr6R6GxVbY3f2BbbZNI4\n6oLle4CXPc/7JaAOBAe9Z3a2hJmjL9HiYvW0hzA0xnVuvfMql1TzAsvUOXNm6rSGNDD2O2ZXA7VB\npxlEkGxMunB+hqKTj5pzk3AuDoujHuH/BPxiYqnowJMHvWFjo3HEP3XyLC5WWVnZPu1hDIVxndvO\necWREjXb1HM/37sds9hXcdPN2zWabfXvrc0GtZ4dpqPKpJyLg/h9e3HUBcvbwNuPMyBBGCap5z3O\ni5UA5YKJbepqo46mLlZ6DoRbOD6ysiGMJWm2ybjurkzRNI3ZqsPGdku1QJNMk4khH8aYINwjaZ73\nuEfeoDJOljeagCxWThLjHZYIE0uabTLOaYIpabrgVsMX8Z4gxv/MFiaSzPPOUYbTUUnTBUE26EwS\nIt7CWJKK9yREomnkDZMxX0Eh4i2MJdmC5QTYJnMi3hPJ+J/ZwkSSLVhOgm0y1SPeYptMDCLewljS\nzfMe/1N8ttIbeY//fAWFHGlhLOlmm4x/JFot2xh6ahON/3wFhYi3MJZ0s03G/xTXNY2ZJPqW/pWT\nw/if2cJEMknZJtD1vSdlvoKItzCmXLs4w1e8+iyvu7Z42kM5EVLfW2yTyUG2xwtjSaVo8Z1f//Bp\nD+PESHO9CyLeE4NE3oIwBizOFAEoF61THolwUkjkLQhjwFe+9hyGofHI/fOnPRThhBDxFoQxwLEM\n3vzohdMehnCCiG0iCIKQQ0S8BUEQcoiItyAIQg4R8RYEQcghx1qwdF33QeCPgSXP8zqDGZIgCIJw\nEEeOvF3XrQI/CbQGNxxBEAThMBzHNvl54D1AY0BjEQRBEA7JgbaJ67pPAe8G4p6HXwR+zfO8T7mu\nqw1rcIIgCMLeHCjenuc9DTzd+5jrup8DvtN13e8CzgLvB948jAEKgiAIu9HiOD74VXfBdd3ngWue\n5/mDGZIgCIJwEINIFYwBsU4EQRBOkGNH3oIgCMLJI5t0BEEQcoiItyAIQg4R8RYEQcghIt6CIAg5\nZOKaMbiu+wbgJzzPe4vruo8APwcEwLPA93ie57uu+3bgH6Myaf7U87zv73n/yNZzOercXNedAd4L\nzAJ14Ls9z3vpdGaxm4PmBTwM/DTdzKcvA77R87z3J+8fyWN21HkBH2WEjxcc+lx8F/AUEAH/wvO8\n3+l5fy6P2X7zGsZ3bKIib9d1fwj4BcBJHvoF4N2e570JuAH8Ldd1K8C/Ar7O87wvB667rruQvH9k\n67kcc24/DPyR53lvBP418O9PfAL7cIh5fa/neZ/0PO8tnud9FfCzwH/pEe6RPGZHnNdvJfMa2eMF\nhz4Xy8APoS5IX4O6SKXvz+sxu9u8Bn7MJkq8gc8D7+j5+aLneX+S/PsPUbtEvxz4FPBTrut+GLjp\ned5q8ppRrudynLm9CvhfPa994kRGfDgOmteb0idc1y0BPwJ8f8/rR/WYHWdeo3y84OC5PYG6m4iB\nKlABwp7X5/WY3W1eAz9mEyXenuf9NuoWJ+ULruu+Mfn3NwAlYAEldD8EvB14t+u697uu+0+A/+F5\n3qcYwU1Jx5kb8AngLyev/UageBJjPgyHmFe557nvBH7T87x1gFE+Zkec10by88geLzjc3DzPawC/\nDjwDfAz4dwCu6/5TcnzM9psXQzhmEyXee/AU8MOu6/4esAysJv/9qed5K57n1YEPA18CfAeqnssH\n6dZzGWUOO7dHgZ8AXuG67h8Al4GR8k930Duv26g5pXwH8B97fv5r5OeY3cu88nS8YI9z0XXdx4HH\nUeO/BLzDdd0vJV/fs8PO6zGGcMwmXby/Hvh2z/PehopK34+6Qr7add0513VNlHf1ac/zrnme91We\n570FuAW87dRGfTgOO7dnULdwP+953ptRt4YfOZ0hH4reec2TfLld150CbM/zrqcv9DzvgRwds0PP\ni3wdL9j7XKwADc/z/GRB8g4wneNjdrd5zTCEYzZx2SY7+Bzw+67rtoA/Bd7reV7suu57UAciBn7D\n87xndrwvD/VcDj0313XbwC8m5X3XgSdPbdQHs2teyePXgC/e5X2jfszuZV6fJT/HC/Y/F9/muu6f\noKyIj3ie94Ed78vdMdtvXq7rXmXAx0xqmwiCIOSQSbdNBEEQcomItyAIQg4R8RYEQcghIt6CIAg5\nRMRbEAQhh4h4C4Ig5BARb2GicV13ynXd3z7tcQjCvSLiLUw6c8Ajpz0IQbhXRLyFSedngPOu6/7X\n0x6IINwLIt7CpPN9wA3P8775tAciCPeCiLcgCEIOEfEWBEHIISLewqQTINU1hRwi4i1MOsvAi67r\n/v5pD0QQ7gUpCSsIgpBDJPIWBEHIISLegiAIOUTEWxAEIYeIeAuCIOQQEW9BEIQcIuItCIKQQ0S8\nBUEQcoiItyAIQg75/1RGwEF59mUiAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x116be50f0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pred_s = res_s.get_prediction()\n",
"\n",
"ax = data.diff(12).diff().shift(-12).wpi.plot(label='observed')\n",
"pred_s.predicted_mean[1:].plot(ax=ax, label='Forecast')\n",
"plt.legend(loc=2)"
]
},
{
"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.5.1"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment