{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# A Bayesian Test for Cointegration\n", "\n", "## R notebook implementation\n", "\n", "The code in this notebook is adapted from the original MATLAB implementation by Chris Bracegirdle for the paper [*Bayesian Conditional Cointegration*](http://icml.cc/2012/papers/570.pdf) presented at [*ICML 2012*](http://icml.cc/2012/)\n", "\n", "Contact me" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "The downloaded binary packages are in\n", "\t/var/folders/bc/dnbc_47x6j7036ft8dp3bckh0000gn/T//Rtmp9bibIM/downloaded_packages\n", "\n", "\n", "Bayesian Cointegration \n", "Implementation of a Bayesian test for cointegration\n", "\n", "Written by Chris Bracegirdle\n", "(c) Chris Bracegirdle 2015. All rights reserved." ] } ], "source": [ "install.packages(c(\"tseries\",\"pROC\"))\n", "options(repr.plot.width=4, repr.plot.height=4)\n", "\n", "cat(\"\n", "\n", "Bayesian Cointegration \n", "Implementation of a Bayesian test for cointegration\n", "\n", "Written by Chris Bracegirdle\n", "(c) Chris Bracegirdle 2015. All rights reserved.\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is a rather contrived function to randomly generate two time series, x and y, after making a random decision as to whether they are cointegrated, and generating according to the corresponding generating function" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "GenerateData <- function(T) {\n", " cointegrated <- runif(1) > 0.5\n", " phi <- if(cointegrated) runif(1) * 2 - 1 else 1\n", "\n", " std_eta <- exp(rnorm(1))\n", " std_x <- exp(rnorm(1))\n", " intercept <- exp(rnorm(1, sd=5))\n", " slope <- exp(rnorm(1, mean=1, sd=5))\n", "\n", " epsilon <- double(length=T)\n", " x <- double(length=T)\n", " y <- double(length=T)\n", " epsilon[1] <- rnorm(1, sd=std_eta)\n", " x[1] <- rnorm(1, sd=std_x)\n", " y[1] <- intercept + slope * x[1] + epsilon[1]\n", " for (t in 2:T){\n", " epsilon[t] <- phi * epsilon[t-1] + rnorm(1, sd=std_eta)\n", " x[t] <- x[t-1] + rnorm(1, sd=std_x)\n", " y[t] <- intercept + slope * x[t] + epsilon[t]\n", " }\n", " return(list(\"cointegrated\" = cointegrated, \"x\" = x, \"y\" = y))\n", "}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And here's what some randomly-generated data look like. Give it a whirl!" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "$cointegrated = FALSE" ], "text/latex": [ "\\textbf{\\$cointegrated} = FALSE" ], "text/markdown": [ "**$cointegrated** = FALSE" ], "text/plain": [ "$cointegrated\n", "[1] FALSE\n" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAYAAAB91L6VAAAEDWlDQ1BJQ0MgUHJvZmlsZQAA\nOI2NVV1oHFUUPrtzZyMkzlNsNIV0qD8NJQ2TVjShtLp/3d02bpZJNtoi6GT27s6Yyc44M7v9\noU9FUHwx6psUxL+3gCAo9Q/bPrQvlQol2tQgKD60+INQ6Ium65k7M5lpurHeZe58853vnnvu\nuWfvBei5qliWkRQBFpquLRcy4nOHj4g9K5CEh6AXBqFXUR0rXalMAjZPC3e1W99Dwntf2dXd\n/p+tt0YdFSBxH2Kz5qgLiI8B8KdVy3YBevqRHz/qWh72Yui3MUDEL3q44WPXw3M+fo1pZuQs\n4tOIBVVTaoiXEI/MxfhGDPsxsNZfoE1q66ro5aJim3XdoLFw72H+n23BaIXzbcOnz5mfPoTv\nYVz7KzUl5+FRxEuqkp9G/Ajia219thzg25abkRE/BpDc3pqvphHvRFys2weqvp+krbWKIX7n\nhDbzLOItiM8358pTwdirqpPFnMF2xLc1WvLyOwTAibpbmvHHcvttU57y5+XqNZrLe3lE/Pq8\neUj2fXKfOe3pfOjzhJYtB/yll5SDFcSDiH+hRkH25+L+sdxKEAMZahrlSX8ukqMOWy/jXW2m\n6M9LDBc31B9LFuv6gVKg/0Szi3KAr1kGq1GMjU/aLbnq6/lRxc4XfJ98hTargX++DbMJBSiY\nMIe9Ck1YAxFkKEAG3xbYaKmDDgYyFK0UGYpfoWYXG+fAPPI6tJnNwb7ClP7IyF+D+bjOtCpk\nhz6CFrIa/I6sFtNl8auFXGMTP34sNwI/JhkgEtmDz14ySfaRcTIBInmKPE32kxyyE2Tv+thK\nbEVePDfW/byMM1Kmm0XdObS7oGD/MypMXFPXrCwOtoYjyyn7BV29/MZfsVzpLDdRtuIZnbpX\nzvlf+ev8MvYr/Gqk4H/kV/G3csdazLuyTMPsbFhzd1UabQbjFvDRmcWJxR3zcfHkVw9GfpbJ\nmeev9F08WW8uDkaslwX6avlWGU6NRKz0g/SHtCy9J30o/ca9zX3Kfc19zn3BXQKRO8ud477h\nLnAfc1/G9mrzGlrfexZ5GLdn6ZZrrEohI2wVHhZywjbhUWEy8icMCGNCUdiBlq3r+xafL549\nHQ5jH+an+1y+LlYBifuxAvRN/lVVVOlwlCkdVm9NOL5BE4wkQ2SMlDZU97hX86EilU/lUmkQ\nUztTE6mx1EEPh7OmdqBtAvv8HdWpbrJS6tJj3n0CWdM6busNzRV3S9KTYhqvNiqWmuroiKgY\nhshMjmhTh9ptWhsF7970j/SbMrsPE1suR5z7DMC+P/Hs+y7ijrQAlhyAgccjbhjPygfeBTjz\nhNqy28EdkUh8C+DU9+z2v/oyeH791OncxHOs5y2AtTc7nb/f73TWPkD/qwBnjX8BoJ98VVBg\n/m8AAEAASURBVHgB7J0JgFxFnf/r9UxP7oMEQhIgmUy45BQ55UYQD0RXFAXXFcEDddG/9666\ngLieq7vioq4HIuKBB7KKgroilwqC3PeVZBKOBEIIuZPp6a7/51vvvZ6enp6ZnpmemZ7p3y/5\nTdWrV69e1fe9rm/96nrOmRgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKG\ngCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgY\nAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFg\nCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaA\nIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgC\nhoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAI\nGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAh\nYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKG\ngCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgY\nAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFg\nCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaA\nIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgC\nhoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAI\nGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAh\nYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKG\ngCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgY\nAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFg\nCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaA\nIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgC\nhoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAI\nGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAh\nYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKG\ngCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgChoAhYAgY\nAoaAIWAIGAKGgCFgCBgChoAhYAgYAoaAIWAIGAKGgCFgCBgCQ0AgGsK1dmltETiI5LK1TdJS\nMwQMAUNg3CHQQYnuGA+lMgKuj6co8v17fWTFcmEIGAKGQN0joDpzzJNwc93D3BgZTC3faRRX\nrTsTQ8AQMAQMgZ4ItBC0AZU75sUIuL4eocjXCLi+nonlxhAwBAyBYUEgMyypWqKGgCFgCBgC\nhoAh0CcCRsB9wmMnDQFDwBAwBAyB4UHACHh4cLVUDQFDwBAwBAyBPhEwAu4THjtpCBgChoAh\nYAgMDwJGwMODq6VqCBgChoAhYAj0iYARcJ/w2ElDwBAwBAwBQ2B4EDACHh5cGybVNt/84Tbv\ndmuYAltBDQFDwBCoEQJGwDUCsnGTyXw0ctkjGrf8VnJDwBAwBAaHgBHw4HCzqxIEIue1I824\n2JXGHqohYAgYAiOJgBHwSKI9Pu/VUjACHp9P1kplCBgCw4qAEfCwwtsQiU9IrOCGKKwV0hAw\nBAyBWiFgBFwrJBswncU+e6BzkXVBN+CztyIbAobA0BEwAh46hg2cQvT6pPA2BtzAb4EV3RAw\nBAaHgBHw4HCzq0DAOz9JQHiXST+naLgYAoaAIWAIVImAEXCVQDVatIW++bjFvmWffsodCDga\nwiSsVu8mcq/D+7mPnTYEDAFDYNwhYAQ87h7p0AvU5rMHNbnoaizcs/pJLRAwNvDgu6Cfm/8v\nzS66YfrNr5hFOuegi/u5p502BAwBQ2BcIGAEPC4eY20L4V20Z+QikevkvlLG8k0IeLDrgP2R\nW2468lwmcmVbHlx0Kve6KNG+bmvnDAFDwBAYFwgYAY+Lx1jbQrCsaJpSLCFYN8u76VjGd8/1\nboeuuwWSVrzBWsDHNM15tknpbb35UBGwZE7s2F9DwBAwBMY3AkbA4/v5Dqp0EPDU+MKYYOWf\n6dxcrOL9J7jsLowNP6Zx23QSFqcHS8DHNc9bFW5VWDPr+OBxbvvENccQMAQMgXGNgBHwuH68\ngyscs5oDAftuXdDZ6Uqt2fmZEG9bk/NzOQxd0J1rZzF+O1Dxe3LF8U1zn+nQlU3br2nH+Qxa\nYmFzZGIIGAKGwDhFwAh4nD7YoRUr7oJmQlQ6xktyfobShJy/hSWciVxTi7qoIWmfX7FgMF9D\n2lHJRVM2dird7S8++3ycS1GWNPk/oTvhNzEEDAFDYNwiYAQ8bh/t4AsGsaZjwCWTsGIC5lxC\ntn6CLzTt2Ll0ET3WTmQ6UFFX89qS8eMJjDovI+we9GWx+ikQ8dvwmxgChkBdIuDhEP8wemJd\nZq/OM2UEXOcPaDSyxyxoCNXnsG6LFnDkMqELOs3P1r8evjv267yOB/Zy0cRtgbDTc5Vdv4g0\nS7uq57hMYTWk25zEh4CD3Ju438S9Hv0B1+2fhJljCBgCdYGAfzm/yxeTldPQPdC3oyYDRMAI\neICANUZ0kaV7FHIsWsCFzpbZpWXf8tfDz46a8k0dD76okJm5lpnMPpk85bnGvxltR1vja8Kk\nrjvwXxEfh7+t0eTNK7uOfUrA7yPsE6jGoQ9Ozr87cc0xBAyB+kDgS2TjW+grUT6I5qyRDAgD\nFSPggSLWGPF3hnwfhEBnz/fxrOTcsoW7lxa9absXgjUbTdp2ZdN2a3Wqjfj/hHsb+lN0IaqW\nseQwdLvY9dcR7278H8rMWMc9OHJ+E05CwNE27v1FjtUIYCzYbUR3QU0MAUOgLhDwe5ONfdD9\nUP02r0H34pd8AK7JABBodAKeCFaadRsNALNGiDoNUnyKyVazsk8t+tAi33JqJvL7lha8aYfV\nAbMJL7nrmmhCmMj8E85fhurHmQpdVF7dyJ9En0LVpX0cqtZyduJRt1yJi0Tr6eJOLeA4yEXt\nhJ/AwRfQOUlgmePfQfofLwu0Q0PAEBg2BDTm625BmSwZfs8vxf0tqt/y5fwez8JF/BHo52O/\n/e0NgUYk4IMAQ12hMtu2oM8m7iO4X0OrGM8k1ihLy2G3ftJNW590+9YuMzt7NwnibVp/0TnH\nKNWtfzvk8IzzP8zuuuSlvtDVTmle1L7R55sKk468ia5qpLlzMX/zwR//0exmWbLHoiLdX6Dr\n0FQu3PHyUx9LDtbjlhFwGs2twbcHP+ZKJPzPnDu5GNM8hoAhMNwI0DsW6sjbcX+M6nd7M6rf\nN79T921+q/NwX40mZIzPpCICjUbAfL82dGuuxn0reii6NyqCOA9tQ+9DW9ChyiIS0AtZjS4Y\n6M3mX/26f9/uvM9+cKDX9Rcf1tTYq9t6x0uCxYt1C8lH+pG5rbcdIlJ1+fXTNk/Y/95no6bO\nNQUX5RQWteSewdGEDMlXUMaB3efQ36OSB1CFa7zoAq74EK56IOh+8M/zo522yGfftcCHZ6Dg\nVP6OR40iPbtEwoSut3OgPPIM01Z3ctocQ8AQGC4E0rFeEbAIlgZwpDpTdapEkypfj+6K7shv\nM/zG8ZsYAu7rYPCZfnD4A+fVehuK6OXzg9CqX9YFy3bbNueSd/5kKJmsdG2rd63sdOUnnXCt\nlzvvupevkiud9dlz/yC3rTDxWral3NzmWx5d6LMHKGzC236wO0XeE8UK9m/pSjuEQaLBjXDV\nAxGE6/cN6fns5bgPyA8JXxqfLf3r27nuzK4Q/28cQ+QejIPe2HWuUX1aJubVZW8yJAT8P4Pj\nG4aUxLi6WBvueKzaVPx3OaaO9CLaEvGMB4ff4nO4NLo9PYvhmF5Ff25JxKF6ZRypbn3pUBOq\nh+sbzQIWwW3tB/gVnFd36lDkcS5W9/C8KvWjyc2YTVydFJ6ftTmatnF+dbGrj9XkWqYqducT\nu6zzuazPTN5MKzaWaMK23wRfVNhMN7XGc+/jS0bBAt7+B2esoSX8MGGL0J+FePrjo2UQbXax\nn8APNuKHE6nlHPaWho2/LD9h7UQMmGdcxA++h6wkRFgi4Yd/Gh4uD13a8h9FOC1w/wq0z0bM\nDt5NheQP4ZrxJr+nQF24j7fSjUh5wsYzb+VWHxuR242Nm7ySbL6b3xVzZbx+c/+AXs7PL/SG\nlRRhOX7Ve+r1ejmquTWS96MfDj770wOBRiNgTRQQ2dFtErpKSgER8bwLVYX+p9ITg/RDSG5V\nlaox0AFJft2MtVE2N9SGQo975pN9oDuf3OkO7zM/bNnz0UCwPtfs1pz36f8tsPaPLmOVCwYs\n3J13HeGHmI0nZSh0xQI/YU9I9wuKs8g1HwNZ7593+ZDXxT77U/TAGS67iNOvYLIXv2rfznUT\nFJ/jspa1Qt3TKI2NUAHsi39vVGP3S9G7UFUMavCIhN5BvCPQT+JPxB/Lcbj/NNf8bl76P2LB\n/xc3U2t6jIp2K9NGJZKwzOswPMxMH8KnIUNajfonjFs+RenVOHsJOFZ6DxsInNBjdQkF/kxS\n6D1whY1+Z9clYSVOxPyOaDcCLkL/D30Q1e9Wv82Z4ImalCPQaAR8DQCIgC9Dt6Iat1yC0m0S\nrCm1fM9A9fLUtWSmr/97NHWTSKmmFQUTrqb4Qsb7TVPuiFq2bmSpENzKwO36Gd5tmbRyWZS7\nBzOWbigN5kbP0w/cIT9/Qjz52S/6KEj1wxBcExYteXQu45qCZY33VVz3Is5DIJLoGe8Kwj+V\nYjppAO4T6HGo4h2eHB+Nq2ep53cDmsqReP4VfV8aMPGYv3w9u+99H3ev+P2sbX877LPcczqN\ngg/N/OSX/tb80puPSeONMfeb5PdGnv92uLuU5H1Oid+81SPwYqIKS9WJegfPQkdJ1ND0e43S\nzdPbno3nTDR9t+h2DrOft/D7WZFG6unKMo5eRfg+qOpb3tEgSQ9WcmROQKCmlfcYwVStup+i\nsohaUVVYdF+6J9G7UTil/qVpxvq7M9M2vImctqGP1irHhY6WaVHHRCX3f9QCr0zTLWyYiiUc\nJdjkXtA8NUhsLb+2nGorXqQicWLFzmZZUcvOoTtas5eJWfzCkr4x7JlpnX7y0C+DkFmGFAtu\nMZ303rgi4LRCei3+R0gTlSvxp/DnNegnUOGh7q95hJNWlJv9xU8upkt9amf7Li9MPOy2SZwL\nMvvz5x6Q+cLH3/P8LcVKIj01FlxVcCINVXaL0M0o2IaJL/+Ay3BBtBzXpF8EwtKaU0uiMSkw\nzOa9HwxvLgkfAW/YA/0n3Ei9OHq2vM/Rt0fgxuW30LuVyk14pqP87t3cNLBvV8NN6lEMBsK2\n5LqHul/jZVX/CcWNNnU/1xhHjUjAerKqrO5LVMdjTpp3ffyRwroZeTIOIflv8gKrTEOW/GO7\nz8/MXiM+pPKO1Evg8s/P2uK3Tgx+HfPLCoQZufzaFghYYU0Qpz5RmHGZNsh1tsKyLkuXaCQy\npFUD4YYuX2096ScWXNO0jE4EC9htCF7+kHbxnWzzbsHSyK0g+Ofoe1AqI3c8ej5aIpGWlP2Q\nq5Xv/0anokqeytP/0WVf2tLcumzhluuOeXlh0+Q8W2dm2MVLcV3+qZ1PJs71eEk3bWDoTD1K\n2PrvDeSM7r4wu572j3sdOhO9GH0veiR6Ibof+m7UpH8E1Lg7syQaxOsOQj+E8g6NhIRGgHqL\n3o+qd0fyfygWp78cl995aJy+g/dUz3i4ZRE3WIXqd3IiDiQa3j8N/wxAZBFrlUNxTLj0Wrr6\n3U7oBah6sxpO4jqw4Yo99gsM+W1hghR8FSzDu3nJ31SLUnWumrOz3yxDyj1RcIUXsGaf8dta\nHsrMfGFJmv6y0BL2ubzLrOxgvpbCM64ly4SsUyDf/8cvdnuFQboi2wPj85mzYc9W+YkDMfvp\n8hNnHRZwkYAJChbwQu9ehHfZovChh+gJrtqNc7JSXob/c7iVRJWl3uln0PtQVaKfiLIdmczM\nda6wYdoBcPSGwpZJXb0cE7ZoHPVYdH+03kU9Ev+GfhlVOT+M6rlTQbrHUFV0x6ESKupK48Ee\nXMNM1mPiaOV/QyOmPHCcHWuinqeRUpywl1p1d1BQWWTC+XxU79xIiZ7hA+i7UVmMT6KfQCeh\nV6M/RL+CymocBil97uo5cjujZ6OHxeSrW0bUM9GP5RugPEf8UCd0XRcs/R8kx6d0hTeWTz9i\nkzGIQN5FW1w210TWX4y2ol9DE/H4u+9clZ7pzy1snDq3sHUipBptXeY6/3uby+3dPG/lY81z\nVxUJGAbdtt7lZrVHubupyYIFnHeF2VxzGp3NzJqOArnycjXTQkiJbVe6rc/S/YnzUUg4/UGu\n63Qd60vyFSzgZtfUSrxMxjWXtLijK7n6ppK4Zd7ocQJ0/lr0YPRnqD4W4dg6M9c0f+VDmamb\nnmfN8gsKl9AlfQMO5S2OdeGtW0mJQhlUBX1jktMVuOolUEVHhRmGVDK4O6Pl8hEC3on+gXck\n6n7Sq8K9pntYPR75V5P3XwwhZ3qn/h+qBp2Ed1YThiIabNp9jd9W6AGqiJ/iD4fIEtTvRvo2\nlEZm2JL1WfxHohIarhqnLn9u4dwQ/gTrm9+xPz1JhF4hdWo5GiRRexI2FEfvJfVDN1G9IKLn\nNx027Ol2slEO9CM1GZMIRFuoPaNo0pYDyL5eZCpnkZqWC4RurNcPqli5lh1crhnDFolcB/tR\nrsF3HURKhd0lbES5UUdbEwKGLHfjkErE70hctdq53GXRqYRdhiX9Av5/UThnCMsclfjXU/Px\nP534FqksWMYZdcch0a6xW/Xf9xLz41y3jaVPn93p1qPWRfSeN23/3NbZ//GvX+fctkzLtnVp\nalNPvUJWjxoX7yef3Nu/DL00PV8/bniue5bk5wnKeC/Hr0NPxK+KWs9KZPIrNI9+hrJANP4S\n9Ai0lbC3opIJKO9LN9mLo8XdQurz4FVki3JXsvCryrDITnJm7IQlbqsSf+o8jgfS8LukAbGr\new7Lt6rTZ0E+op+j30rue3vi3p+4IkaRdC3lJBJTuqpD5P4QfYY8PIVbC9F7mZYPr9dv/5so\n93Cf5z7fwW1IMQIeo4+d2cWblXUI2GX3fNjNvSr04uxLkFr3cF2xKxJv9eI7m2a7fLMsgKIs\niXLfWRrlvl8MKPHAnLSc1R+qBfuSSBXUobFf6yqjFsj3uxz/PQ6L/5LB/WJfYR25hSz8X3VM\neHPsRiIWJHoVY8EzYn81f0N3tX7YjCfnd594yG0zmuevZEeA8G3hibjbULDza2gorCD9N004\n4palRMfyCb0Jb8ZVhVRv8l9k6BVJpn6AqwoMia5CH4n9gYDlvQh9D/oP6BfRM9Hr0XNREXXA\nB3ceOAS8E/d4wmRl1buooaCGmhpNlMn/JzqQfKtngMevGfVhKORA/PeiJRLdxYGwOrokUF5h\n+SDXXVgWPtRDNZwlq2On+PfSxPcm3PWJv9yaTIIH43gwdFehwkO4qEE/GaUhWjPhd+ZauQVY\n+lfi/zW6EBX53oHbsGIEPEYffYfrCCQ54chbTp9+1g++Nfnkq/UDkuXz4qRIR/KyD/j5RlE0\ngz2eIaiqpVMx+YYwlXksWMOBMFmFNFMhdJdvwqEF77fFMRRfVkQgw8cUxvXqplItoIpVRJxa\nvv/A3K0TFVZRuLBiOIHcXxUJXxzujNSd7V1mJpEx2t0mLrsY93J0l/nXn7An7kr0IFQNAyrz\ngWPHdTWQMP5WKR3wCiIMz6NUv0iOSx1V3noeYB3K92H8E5MIwlWV+P+gdydhajTdTFmvwL0R\n3RvdjmMaUXUt6bv2JXL5GVTl1LOrVnYn4p/QDeg7UT3/+9ByeYKApCdGpwIux+KRBfo+jnnP\n/U/QkjicGZxsz2W3oGeVXS6yei/6MLojqmd8LNqP+M+Rr3/tO1KYn8G7FCb1PYArHL+P8k7o\na2g1k8dJaQ/0l+jvUN4xjbFH/43b0DLgCrqh0aqvwgcCnvvrU/LTP/aft0MskWvpKJDFM9AX\nUFk2+sEOSNjgcTobVokwqxMGc/mhd0JwJ+OqQiuRsEaVvq3cRkzcH5G5j0C8P1IE4os4v94e\n5VXBQNKFP3K9KqBAwLjTOO7QOW5RscXP1pWfbXPZ7ylOqdAHPkUbfpSMM4fTHM/Cs40bf4t8\nXFlwHReS+eWZbG5n15S/hnOnoqqM9btQ3BGWkD8w1N7Woev4EyUZSJ/la0BkRUl4qVeNmOc5\nT7GC3J64cnLoVHQZqopPjSyVVfoGFGswiMr+qcRfh07YylSNs6tQve+ppKScHvflqqFxLyoC\nvwBtQ59Cy2UVAaXkejrHpyWR9J62owo7Ce1H/Bd4pi/vGcnvQ/hhhO+A8v5Hd3aPE/EbUHe0\nnmlYlaDG027d41Q8OorQV1c80xVIuqFr+Mu4v02CyY+WsNVUHiK1PVB+0+EDOGtxP4s2vOjH\nZjIGEeDBqQKV+UfLP95lJpq+nh9wsICvT4pUWnkkQX05fm9mQO+MxQgJDEjUGNgFVUu3KHRL\nB0sY8t24PMrdySYe32DG82+pSVYqEjOgIYtYlked13PuQkhSDQfET8dK/YO6i9HZcVj3v4Sr\n0t25e6hYpvnFIPNmsNmr7Nxs0txGl/ql5OW29kjLLArfJm625YC7HiXucWhq/e1bdu1IHO7P\nTSag+6GvRM8FB6z4UEGrLN9A/4z2Jk9zQqSRyl143oieiH41CaT7OVJjQ0T8NVSVYiq34rkW\nVeVdr3IAGRP5/Tv6UvT96Bq0ync9NAqPJr7KeTWqRsl09Cm0XPSekm7Y4vTt+EWUT6IfR3mN\nwxK55bi8b6lofNN/MT7yB+Mnr+E9/Rhhf8B/FroHqjxIfohejMoCVjn6E1nAO/QXifMLUbBK\nf08Vr1BDZgX6dfT3qOoUGsJOxFxLuYfEhNcZqJ7TAt7B0sYTQY0pRsBj9Lnzq9kIAUFgngo7\nM1PFmPSSe76Dsxb9CboVTa0mvP1J6HI9JzNt45bmtqWP9Re7+/lIlQISPU6eij8syO5UhW4i\nr3Il7VHHz8jzdfIzjt2t0su4fI5faZafqpr7U7GKIY3oOqzlitYoBDuFZGYoLW1vyZeUtpOf\ny6nggiviKRVVcttKAzpcJwRMQq+9WhWr5A70r6jIeKTlJckN5+FSSYXJbLJiLkOfRS8Aj275\nJ6xUVJG/qisgWE2/5BpVqn9Hn0PvT87LXyp6Fl9DRc4VGzylkUfHH0jrLdybskS3ozyrSOSh\nckNwVW13eApx16PCROQpuQK9K/i6/1GDRu+ByOn76Nnoz7nnl3FXo5egv0Jb0VReh+cD5GUf\n3L+hJ6FfQZtQNZ70mxCBf484l+KqsbU7eiD6HNqf6D3Yoe9IYXnVTsSZipY0DkqvCqsQVG+8\nlvKsQ2+M/epaj14ojTl0f5ggeADp8ttX2lGxPhh62mM7BSPgsfr8IpeLaMWjk7Ea9WNzM/9w\ngiokyCq6gsO1ClJ4lSLSfk92r4ceYZtLXTsA0ZdPRHieSdNdZAtx7kfY3c9FXRttJImq8lL8\nlPRCMMxNV7ZrhiWncm2GdcVs+OF7LQdEO5myTo/Tcpc2ueazknTBgDOxFSBvEI7ZoSs0TNIg\n9dOr1e9mnPu5pThqwNyHrkTJxlBFlqt/BP1U95TC+ujyxoGiHIHS9ggzkdsUgMia2g39OyUK\nuCmwsoicI5FGBQnvxPySNFLsRUgihJdy7nLcNSiE4K/CrTf5ZzJ0LPpMWcY+z7EaX7uWhZcc\nqhclTAb8AoEigk5Uz/7N6Bn4aSf2EGEpshNWwp7Xp9jD8Db8P0KXoHrOpBFkL/5OQm9EVb8e\nj74JVd7VUADnUAblVdc8gGZRNZgfRPsT5aO/hvU+xBHh61nSSPC80/5Y/KXCMw7v2mNdgdGf\nKOLjXce19EX31DK18ZKWXhCTMYoANXUH+k7I6kwVYWIy4YYPIfxjZt7KDQRVScDq5gzjRX9u\n2ecBVTZU5AOR6FnINN/pct/lqm4/YEhSFV6ZxISdd50ivaIQN8dBlnIkpJpTGZSXCcVI3Tx+\nCjViiEvFMZtduObqNGEqTyXZmXwqzaK0kz5hmiWtSvO96PtRWSKz0SFIWDb07ySgiu4C6rqX\noBC835njE9E/4D8aNxH1ZLhj0FvRA9EWVFbpAvQ36FnoECX+clWSyG9x7wet/0W/jD6RhKvS\n5hHofaj1etPkDoN3DuFSkWaa1ySlaD0eEdgc8gxuYYvD5FzReTG+d6IiVDVqEglLfpRmJbmF\nwI3ojejJSYSVsRvx/CJhdRl6F/ppVLIHqiGZWWgn+ja0A/02+nNUVmkbiiUY5Gz+XoruT3p/\nCSF9/1HZd+k7Sjiv9JW3fdG90FPRUlH4Cu7ZW9lL45p/mBDIDFO6luzIIMCEojCZKdyNWVgT\n9H1eCOhHU9/wK5zwFZKTqJC+2rMyFeGE73q+nHgizWkTz7lYLfL9SFOVWdVCI0AV0LoVkXuQ\n8VXIw18EoT0fJ1BIKqzS5DL88N0Ty6MwDlk8wcuoLuhmatBAoNRazFZ2HeRHZFRJpnCvQMBY\nyrMZU1blimREphUk2oF47d1OEMDxNj4gAeloXCqMTYmAX0Xav+wWd2AHEFhY2rQB91kU6yJ8\nbevvuF9DJT/jHh+KvaEBJOL7QXKsPIiAs+gj5Etp1FCin5IYlX4PWZKENOHKqqwT8XuTkVb0\nEPSjaLmsJuCT6NfRSkSWvBthfeuT5RdXPg7rrE/i3OfB/1ZcPatru8dV923oQWjlWYp8F6If\nQ+9Glc9p6M+5Po9uw697K41rUEk74WeiKe4hsI8/+u3Q0Ox1trwu3R4VHp9FlQ/eM7cILZVj\nObihNMD8I4+AEfDIY17LO+oHXRRIbyJ8MkMBmakbRaIz0RPRD6IQgsbIwlaEap2rote5/0Pn\noffseNH7tEyHodRCt3Q516fwScJvw2IpkThI+AOkE4iX7SoTIu5KYqnr+Okm1yErr5uQBt3q\n7NLssrIS1D28CUd5adFxBZlCxaXv+56D28K1h833bhdwmJzG7WoIeFkjSKakyy0O4e8W7j2p\neBQ3KFRxngJe80vC+/BqDN23xhH8G3FFBpKfoKpw9SwkVJ7FrlL5X8Z1FDWQdTvu46hkFboc\nJWthDS9OraXSRJjoYe6ingzJ7Nipi78LyQXvdPQAWslqE15HoO9CISB9BEQSJrFNwJMS8G9C\ncNV/opu434Nx9OhC/LpPmYSZy3pW70PVcPkV8Q5Av4b/Rbhn4qZyPx41xihHkDWJW61zHxE3\noKVpll+rskLAEb895cPdg5Y/y0WEPYKajCICRsCjCH4Nbt2NKCHfCcw8DuTTtOOz+mEfjKYV\n/xfxX4beiaoSOhVNhYrFvZdrA/FFZeOkaaTeXCxZfaZQaRcF1uhEV2xxHcuLgamHjTdWRaGF\nnoYkbmx5NyWznttDd2OkbvaKBAzhhrKS31AW3D0mueyreakh03TNcTpuGpbufH6J67ii7KZ0\nWWvrwaZSAv4tcV5I4u1ZHr+X41cSTkUddhL6N/y7J/Eewr0p8WMFxRuXJMdyXoOKrN+JyuJ9\nDpWoAfMWlIozzFpW2EjJB5IbQWT9SfzO9BerBucXkYbw6U2WJSfoxQjLk45Pjr+N+1f0aPSP\nYCmCHg55lETfhm7iHiX5jMpJjkZdWP50Ay7XaJnRQESkGn7HSieR0O2udx7y9+cQ2IrKAk5F\ndUE5Ae9C2BNpBHNHBwHqKpOxi0D3Hy9ExU5PcfdtdtfH9ePaGxUB60e+D3oyOhE9HNU4lX6Y\nqrg+TEVwC2OoIhFMrtBVJu9QJEdK16ZbVlaTUHF3L5eRBdMJM4p8yXsYH+2RBNYtFrBE+0/H\nQthcrplEGZ5OQp6Vy5rfwyHfc0kzsYTjs8lfLXUqIeCA68s4J2z26haz20Go8OYlQbI0lIYq\ntp2SMFWWN6KXoe9G90XfjEruQ0UMEhGeGhmyjmStqEGhZ0Jew5eeOBxJCQ0h5T0tRy839+rC\n5h2Ke116iVSr4A+R0C19JPZYck6Np+XoecnxQtwD0TeiP0zChsPRu6Lf2jf6TlzvVpidfjPu\nHn3H7fXsvZzh96zlTeGznjq+BNXv+iJUjbkH0FT0LPlNpRJ+N+p9aU9DzB0dBIyARwf3Wt11\nW1dCHsLLzIFEAyllZqyXtTUNnY7elsSjCy90Pat1rB+/Kicsg8gvCl8dirtNowFawKTRQ0hD\nRCeSr1pyriN0LTKWq8oi+DOuQIVV2QImfCrESeWTjv3qVhEVS8T3hrXRgr+BgBUKpXZcDfkW\n5C8X0hAuJQSsGGErwrvx7KejXuR1hEOkweoVYUva0FnBF2a/hi/IqNL9LvoQ4degF6OqJEUq\nktmx4/6DOB69AsWaHlURbgv6yYEabGo4qMyIyMC/IfGfjj8b+4f61+9KCovRz/SR0n9ybndU\nRPtHdD4qmRM74X36feIfDkd5ux79yXAkXpbm/RzvhN6FfhXVb1nP6pOohNe/2OuiYzXqqAc0\nETDIi/mrOAo3GUUEMqN4b7v1EBFgrFbkpEHC/0Xv5hgejS1glhKJbGTt7oCq0pcl+FH08+jr\n+f1hIWvpRfQkx0ojISC/mi8g/VZhQxHSY0JVxbG6XpPlZdysk3QJ78C15C2IyqhKvrsQgYDJ\nxH4Ud7uuk+H7w/re4KOMRR9HOOVmR61eyDe5bgvd78KqXP5AwGvKA0uO5+EXeR6JHpWEq6dB\nvysaRKG7PwlOHfUuqBs0uo0QVaQiZYkqQz2jehHheijYfQJt6iVTByXhCxNX7xaNB21AEYjo\niCR8qM6LSEDvxuO9J6SegugxVLh/HxXhiGhaUclvOKeG5zCJup2jl6EjQWq3UojL0B3Rd6Bq\nbB+Ihh4s3PXon9BEwgQv8ldsTOr5PEJeN6QxzB0dBFRRmIxdBGS5QZ5h550ldOG24oeU+GlN\n2wAHBmnj7woUsoi+gVKxRn+NT3X9hc0m6IiLnn4y6nOsreuiPnyQKBVi4Zk+ovQ4xcu4RYHk\nZX+cUDnAs6pQexAwNQ/rnzUD3JNduUWZxlg451IyD+uIA07FGD09VFiZGT2DAymCW9ekrrI4\nqaX7VcKVz3ZUjRwJRNxf93H47N3pcXR1TVeaFJWcHXnn19zyTahI9Zhebr9vEv4hMBJJvj05\nlnUv+RThX0Z/hd6Azkd3R3s8zxA7/NHkqR7Lnw7h1CPgk77TXdEr+0Q2khtRNeS+g/4MHScS\neknOoDAiWomeUfj94t6F/jdYqQeqVNo5UF0g2RtVPJNRRoA6z2SsIgDhbqSy2pBzOVkGT1M7\nvRYykgVGR/SmlJQm6Vx/lbt3WcUT+d0kd6jCeOspy6JOdQtWLUsT0uSC10GityYXaqlVjwqb\nVsbU+Lx/pvQGXDcNHCDggiwm+pyjtRz3R8BYnums2dLU4u5rQlq7hXYdbJd4D8D9IkoRwhjw\n8yApC7IauY9IpxEfgqkrubkkNxW64UO3cCtxKGuY4PR+XDX+vo/KKpOcgJ6Nvg49BpXldQsq\nYq8g/hoCl6OKn0hYK/0JDh5KQ/p3w5pm5Ws6Sjki8qD1zuNO1Hulnq4/JyW7Dfc4ynpuclzq\nLONADaIX4c5Bu/1uSiOaf+QQMAIeOaxrfqclLndJweV2wwR8Hsv3aYgKwzDMxGRxDtzcJaGb\nueuwoi+0oOl+/peKZwcaGD4xOMCL2N2LCiK03BkH/q6uhjw7cHoQMJECAWMhd+tWJD4WsJsC\n8QYChohFiqp8ehWuWU06J7NLx6SySKs4XoJ+uCw8PRTxpnIFnseSA/mrFFm92qKv3iRqJ0ef\nRvU85qLlciYBemdkXUrejQprkUETmsq01IOLpRzGxyHEiuPDh3FeQwF7oKl8G8829LI0oEpX\nZH4v+oEq44/FaGqsvICqAS65g3dpXezt8XcFIarvD0F3QLv9bjg2GQUEjIBHAfSa3RLCWhbF\nLVksv2AhQMLhmTZN2XToxKNvkiVyOpr+QPu4dTwGClP3Zy32kUZNTgXipN7HKg1ScRlSc7JW\nGAOgrCKJZlMJzYCE1fXolkb5a/iWsSr2XgXsVhJ/LiyviqlEQpfwZwl4LYTxe7SEGMLa3aM5\ndz1KRRgmWJ2DfyF+CGY8SHQBpVDjYFH30niRot6rH6FYVaGrXqQry+ta9EJUvR/PoZJvor9A\nj0FpcIQxcyz/0pnnYab7dpzTNa/i3FRUx3uibwbTP+AORP6ZyCdz3bKBXDTG4v6d/K5FRa76\n3fYYWiIsFd7RIG/g7xw0/X3FofZ3VBAwAh4V2Gt/UxGN777F3t473fiylVRAP0U51bdEwZoJ\nk5X6jdt3SkM7iyUaiHNb0h1NN3sPC7jNZw+mofHvuhOuuhqLQjnmcbBzagEXT/Th4YMMVySn\nZX2Vy90EiJhfgV4NKXwdnYX/QFS/H1nHR6CIxt3CeuP4cHz8/QvFeCNl/kcUfILl+i7CRMpY\nmKGRospd7w2NlDC5T5bux9BWdAv6AHoPmkPViBFxyD0fTeXViedbuMeg70F1HxFLyYQijqoS\nbU857p5FecnVsHkv5czjHob+vDxCybEaQBej+6I7okbAgDDa0jzaGRjl+6vCnYaq1T2qxFML\nHOiGfhxCKiaFqZEtHvTrCVsxqrIbVSH3iQUcT8jKu/xWllaVE+ORxDuZ8hYil1EXHOK18cct\nlP8oHdEDnqajwz6F62RFIC0T4h7v+Cj5K/KgMg/drYtxZVnJyoOIg2B9dG8EJOHjxNG3aL3K\nLWv3bPQdaPqSrcIv+Ta6nODkWSgoNPo2ce3pHIjERRK/Jvxxwj6C/yb0BDSVT+NRY+vzqO53\nCvpS9KNcIxI36YFA9AxBf4yD+5t9rWfjf0Pcf0J5z4sflYgvt7+jgkAjWsAHgbQsHlW6+mGr\nJSj3EfRrqAh5TAqTsW5MMv5k4tKrWq00ieS2VRt7uOJBqilxBlcTzajt6Y4sleKSoS2droNK\nnqqFLme6mo/GXaljiDhNR4d9CmZr0vAoplsSP8zCFlmostNEniWorIn5KPcaz+RLCWNReSVq\n3NwcfPGfhIC1U1fUy1hrJNJdg0IA0f3JteoqvRzdhyfXlITtEocFsoXMi2uQb0/OmzN0BNRb\nJPKVpI2n+Mj+jgoCjUbA6jb8E7oafSt6KLo3ehx6HtqG3ocOgLiIXSfCZKzHIbDb0d+grOVo\nKlrAi33zx/lm7hm9Z7U+LGDyF6ydpxOXruQNVNzNrT5MzgnZ9y4zWR6IeYs08at7U2GBeLku\nELPC+hOYNSXgtHIqv+QWAiDeSFbZg2gzyq3C+4Mz7oWyB1nBX/027kb1GxokOYZu6w9xveRn\nkLB+g9uh31QAwuMvru1+IoTYn1ogQEMoiIyOlYnfnFFEQBVJI8mZFFZWrsi2XFTJUhnoM3Gh\na4xW/diTpS53CNTg+SThSTjFhgSkdTTHT1GiH1QqFRbjiyDtUW8Vi0DJxxaVQflklveGJup8\nGFI9E4Eo2TRjkviP8eLNtCADAeMPBExlTpzoifUuJ6KsTuiv5oOEeSKrF6CSfJnAbyQn1OW3\nB/ol3ScJG+eOLFj/Ogq5HF1AuenKHGqDTd2n/iTS+xX6BnQpeicqUSOYZ6xu/kjhJrVBQA0p\nke9F4Jr8XmqTsKUyOAQajYBVwSbWTq+AqZWvMaixKQlxQWQdhW5jwGHGqSq1HgJZf5HLsEj0\nXdjRFZEqFfOGNBe8oBvlpyUhAl4tP8ysjTYkW6id17c5fw95lwWlc1jF/tnno+ImBQquRngv\nKu85nVRW65JEZKX9mLBuk7+qucHYjhNdleQfrCUD+2RlfE35X30AXh+LD2O9NGyKm0fIwqbL\nujjOXn6hHQ8KAU1o9EdwqVm/g8Kv9hc1GgFfCYQ/QvWj/x3K8F9RNM54Onoa+tVi6Nj18Gk/\n18KejlOnuew78POhhvjrQWmR2nzz+fwgb4ewPkLl18za23r4YULArj3NIw9og/rRm1wWAi42\n2pOGRGwpM93nJqK8OrmGTwtWP/6b3gd8tkLevVnAaTTcMOO0wci3pPi198oq0zP/l66kw2Su\nWV3H5qsdApoEZ1IvCDQaAV8D8B9FL0NlUdG1FiysGbiyoFQZnIE+iI5poRtXM0phrew76V7+\nKuTSDtHuh7V71dKo47Vx4TJvgYDniXx1zKYVo94FTX42QIaPxPkLCxw3LHZeE7GOI+weTNxM\n10zneLyXsj5H7lN23kpcVegDlGgbNA8B5wd4nUUfIgI0AjWm3usGEkNM3i43BOoXgUYjYD2J\nS9CfoupmbkXnoOrafBJV1xc9t+NCsIA1BpxR+TRiyqxdfbTeacennSa47Gsgrp05tbfO49/K\nkp9Rb3hsLt+JS58P9KGrcqHySYZn4oQhAhoVyUQr/wDleljnKeNWLPlBELCGJnrrgo5Ttr/D\ngYCN8Q4Hqpbm2ECgEQlYT0YVtJZEqCu6XDTDVtYjvZ9jVyCkXDwG7FtYykNBuvZTzrrsxzn+\ngEKRIyEeyhpduDwayH678cW1/rsqisd5y9LV8wgTypqcU29FItpsIex29UssV1Tin6O86tkY\nkNAAYVMxd8CALrLIhoAhYAgMAYFGJOBzweuDKBNow84xH8MNFTmu5C7039Bf6GCQovGr/0CL\ny4D6SWe3fs4P5jQfMchAWj333KXgmohRKr+kW/oTpQH15Icc6R6O1y9mXHZ6V9586XMLwcwC\nfxcejOOBCd3at3KfQ+ii/zSXH8+a4qMGloLFNgQMAUNgYAg0GgHvAzwi3y+gsoLfh/4VPQat\n5cQaEYAsaJF8NVLzgUcyoElYWdwJ5RkgbOfE+g2nOK56zWx5WiN03EF+gwUMOULAce7J94Ye\n9x/MRyBIhG8pM84cTYGEDyf9BT3StQBDwBAwBGqMQKMR8Ong9z/oVxIcv497BXoNejxaKyJa\nS1rvQasVWW10BddS4j2URcIVUp1TGkacWpW7NNla+mUBT1WCjGX36IKuxY0g300Q+hTS0hyA\nhOxrkbKlYQgYAoZAZQSqtdAqXz32QtU1/FxJtrWJwxuTYxHxuGmQQKo5rDmR7wRcuKVLIJvY\nhEyCWLZT7wTMGHA6ht3E937DN4Cfgitf6CrV0HywLh9QjqYADHMAbDLW0NAcmavne7dLm3e7\nj8zd7C6GQO0RaDQCvhIIZZkeUgKlSPgkVN2OP0Z7dNkSNgYljI9iLYq4onZIi3Wu3T9enxaK\n5Ud1TcA0ECBgP2eRz17KDGc+NRhtgDBfnXedl6RlGKrLxC0sYN8GVpqIZRbwUAEdgesnuizr\n11suHIFb2S0MgWFBoNEI+HpQvA69AT0GTUWzZl+O7okuTAPHsgtprYSo5vlAJv7XbOl4AHoY\nRPZDlQuy4VQsY8ACVhf0IZDkGXRBzybvm5ZFHfeuiNKvGKUlGbxbcHlZwKFngD9GwIOHciSv\nnMjNaDR1CVbx5IU+a7PZuyAxXx0j0GgErIlRmngla1eznUvlaQ4ORN+BPlJ6Ymz6/UqIah5U\ny9rfaHN75B5G27Egk7J12xFnXT2XMeMKEHC6WUi0CwRZc4tdFnAXBlGWFkog464w89UhAjSU\nfLcG80TX8g8sVftdHebVsmQI9ECg0Qg4BUDjwKVLWA7iWOPDImh1a96LjmmBP9ZAvLMpRAsz\nfCGwWDpd5/ch5u9hSX6FyusOhTIF++7kdL06WgecyhmY7iVkmQYPzaU7e2lZCmYFlwFST4et\nvvmV5GcC7/hENpYp2eNcs+Sd1MQQqHsEGpWAyx/M1QQcVx44xo8hXVm/oTu1SGBPRO5p1ri+\nc0mU+w4krXJj6g3gy0GjAAr5LDYgVOGiNSdg9Q6A1+q0eJhVbW0++2OzhFNE6sfluRzc5DK/\n493eX7lipmGRcGlcTuX9mMRzwxA2MQTqGwEj4Pp+PoPOnXf5DixF6qbuBFyaIJYxcfwzkA/b\nMNa1FAlYuSTPNSfgON1oI6mHPaWbXfYQKvK3LHIt+9Y1Mg2YOd7rQ+Pn5ebL5YFNS2Fg85ng\npztLS8pMDIG6RsAIuK4fz+Azx4MVkWgjDrqgw5rgHolBMGsI1JacdS3qTk8zCPlq05I70uPa\nuoHYw3g49wmVO+mny9RqeytLbdAIMHyi4SJ6buJPULaE2dBpcj4QMOw7NQ2p5C722Z+2+abX\nVDpnYYbASCFgBBwj/S2cx0YK9JG4D6QVPkcYuuOSj9aX35eu6G+xdaPG0upcSr/SFD1Bvr86\nTBneDG5hdy0sqaN1Dyr7xcN0r2FLlvHRY8c3ucTEWwLg6amfbulAvHT99EnAPOeX8tvYK73O\nXENgNBAYNxtPDBG884d4fd1djuWLBRxhAfudqGi07WZl0deG6lzI4jNxT7oy2nP/51pln8pb\na4FZBw7tOn8iuJF0pKUu9Sk+DC/kQnZLcphxGX1ogzHQ/G9LgseNl7IFCzgtEMel3c37KDzq\n/v3oNGqJqx3VMqiJITB6CJgFPHrYD+udIeAw8QoSmcga194JeFhzUZvEmZ68jJSWKDUq22J3\ndG1S70oFqwicoqcIeRLc0t9GyQzbrrj14FvsWn622DV/qjwvYLQvDZXQFVt+bjwc85y2S8vB\ne85EQpao+eIGOgfF50q3LE1jJy4J4AOfMGO67KQdGgIjh0BayYzcHe1OI4JAU9cH6qlt+rCA\nRyQ3Q7wJ07S3uI79qDA1GSsQ8RBTrHg5di+bcXh1QRfHxTmuWwLGSt8BLtm+Z2H8HIhp3BIw\nz6RYNjB4TuXfWVZwIFbWcAfpnYAxn6clDazi7On4GvtrCIwsAtYFPbJ4j9jdZAGrmS+hlTWm\nLWCV4enIbV7gc8ev6LmBik7XRMBMy5Ay4FYkXcKK/prcpLaJqOu1exc5y29i8o2KJFXbW45e\nakycOjDvol3JwYI0FzyrtfIXIODWkvfcu8zMsMI9jVjiAljS9RxZF3QJLuYdeQSMgEce8xG5\nY97lcplkR8XOsW4BJ4ix9eRfhxM8JqT9P6Xf5lp+mjZeOKxjAtbHI7pb6Ds5NxPrTtnvcxKS\nyjnWhIbFOxnYfk/3fHsIONKbPoWG5nqdwyqGjx0EXFkUV2eIN+4aKZVLbKH1ioB1Qdfrkxli\nvmhZMQkrlozLjXkLOC3LsLr6lnD4nrBPNi5RlzebOtShMNP5GCiEGdrd88cMsmSC0ngkl67u\ndsgT4hWJRsvkZlyWDTjicWDcNeBSJGD2hy7rps+mvQalk7eUjIkhMKIIGAGPKNwjdzOshUAi\nVFR5WIQKyWQACCTYuReozOuSgBnj31PjmDznsvxlQ7cq4eOKXLTdJM+i5CML0TPQ713LXO7y\n+Ln6iXzWjPZHkCc5t4N8O3o3ha8mrWrzXRuq8JtIMRtXGMVFt79jCQHrgh5LT2sAeWVtUU6b\nGbPb1SmMn4aJKgO4vKGjUtEnFnCElRVX5PUGCB/VmM0HJPSEU2suZJHwKZCz/jXx8Jvx1P0y\ns2qwZWeyPYhXXJPNJzS1p/njKp9H2J0lizWREHC0grg7Kl26rGcJC34HoWGiMCTBLDICjvGw\nv6OEgFnAowT8cN+WCimQCJXP88N9r3GYfkLAfh2WZDeCq5ey8lzTbtVJ+jDBYt9ymvIG605O\n84gJWJd5T/M3QDcpS5gJr0lXtyyN8lcqDZohnWgLcx4SAvZMF/BzdW6CawlLlog/tc03H68w\nMAoWMHF2WOTHxudH+Rb2O8lraFSoDCbjAwEj4PHxHHuUAgbRkh0qKu1vbDIQBNLue6r2dVTs\nabfmQJIY9rjkS1+6EvnsicX7a7yfW+Bdm2YDpzdnPc44IuDY0mfMdwWW7z0513l7Wk49r4zz\nWQg1fMGK808RLxkLjy1fPt5wFLubXcu6LU1OS3CJ8GePTtOpZ5eK+qu0McZEXusZx3rLmxFw\nvT2RGuXnucht6HT+wHaXu6dGSTZQMsW9s5lVGzZ5oCez3iRiDasmi+mTk5GIpy3rWm5gOlLR\nAoacxxEBqyx+GyTLcEq0lJ1SGPItinosWnhIYYyY3oG1lD1piPgwGQtSPkmxYVzCm9IxYDVg\nijOhWeZ0M5bmu4qp1ouH1kQ8pl8cu66XnFk+hoiAEfAQAazny5dHuTupYfjtmgwQgWIXtK6j\n36/uiAxCgUjCWGdp0XbRmDXnwqz38UTAkOvBvMhb4SIIOC5fScFzBdeEBRxdoDB6fV4gTtIQ\nybQqjHP7y6WHYDL4TEIh8EDocHIspH8AmJ22k497F9Lw0XZ5/yaTf23NVmw4jHae7P61QcAI\nuDY4WirjCAF9ppHiMJO2EMYY6cqtn25odpgQ1BDFVEhkWTnsVNR7oWHcv2tMtDzWWDyOXk+p\nRZrXkfubS0tAefnwiLqf1U0tUs1vhEyDBcy5bt22bJ8tYuab0m6p0ikkny+kmSoC17emXzbB\nZT9Zmv5o+2OrnZLVYUNwtLEZ6/c3Ah7rT9DyX3MEvOv8Zd7lz2EnpceUOJV1XVjALKV5S5vL\nBvKhMta613blDyLGsCsKQ8E+EHC2bIZ0McYY9GD5iji38iWsC9FvlhYBLNj1LVpMmBpK67GA\n2VKUj2iExkrcBZ3GZ9No0okWcIwlHW2AuIMFTEBwFQ9c0wlu6WWj4mrpFd3iNxRctjXOgHVB\nj8qDGMabGgEPI7iW9NhEoD1yD7dHeU1sChPZ6oWAscznQw/zYlSjqbDufZDtn6DgC1KkIZ4T\n8K+Oj7svUUrjjEWXZxAIuJe869Ob54MN47nRepZnbVK8RS77VpziGK/CwHAupPthCP1vYKd9\nv8N51mOWxBv9LSpbvZuZdc0fojzHULaDlHfG960LOgZi3Pw1Ah43j9IKUmsEqPi2Kk3cuuiC\nZhYvHxFIxzbDGPAzWIMnQLrLy8q+CnLJj68uS1l/PjyPsrJymO5cpjN+XWey8xtE/AMC9gmh\nzj3KOXaHy7SCV1On6/gO4VjAmszGOibXUkLA3a1mnR9pybim4/ms5OeS++4nl+dpBDzSD2KY\n79foBKyuRe2YQx1rYgh0R4BdOJIKv7h1YfcII37kp1EJb8eaX9a5uplYc+uUBfbPTPLpw6Yb\nxFnD+RfogU2W4ox4Rmt+Q8rERKTKHxUhvLjtKjdejzVbXHoHwU6hMeLZuOOHpEF4xEQ1Ea5b\nhV/rvMMGHXRNvz/NNGGjjhuNrWKXOPl6m/LGUisj4PQhjRO3EQlY3TlXoNSvYSnDs4n7CO7X\n0JKWMEcmDYsAA6lpF3SdWMAa92WHK+cgkYiNJzqXxg8nn1qGzP5Va9KvgUSewt0pPj/G/9Jf\nTAmYuVz4r0oloawdXeHROrrm9ZsuChefucR1fgnsNoNejkj2AABAAElEQVTJ4RDyFn1di0l2\nmi0dCBjUXlG8wLldS/yj4qVMYRKZbk6+A/ESZgQ8Kk9j+G7aaAR8IFAyZhbGyDQ+dCi6N3oc\neh7aht6Hal2lSaMjELltVNaF0t2lRhMSKuBi41BWHey7LM5PxPIcTcSKgkUMqdyBHwKOxgUB\n7xxmLUcRM5af7AX/IgF7l/9Oe6TlSj40RhQ/79wjsJis5K1gcgRK20oivOKPNkDMRcLjBOu/\nR1fK8pNkJh3/H9282d1rh0CjEfCZQCcr973o1eht6IPoLejP0JNRWcInoCaGgEyvzdpfuT6g\n0M5NRdlM5kKX83rXyfsb/ROksyUQc5S/GvdZ8q3hlfEgwfKDOFmGVFGKBAzdPq4YxF0pFxy2\nshnNnYk/XE9DJkzSIpZIOljAhGmSVyqj3gBn+CB553zaGKDx4RekGTR3fCDQaASsMd+0u663\nJ7iCE1rSYGIIaOLLJr5xP5lZqXN38aP9XnSNAdIwKJLR85FbvzTq+AmPCwsv3QPcEZxuxzi2\nHyTrsIPln3cdxbHd0hLxjH4N0d6jMLoBAhkT9kQS517YOISBTRhSgMjCRiWQNeTmZ2i5D3ET\nAtZErdFfdoYFHBpblEM9chJ6Ndz84LM/4waBRiPgK3lyH0Vl6TaXPUW98NqG7jRU3dQmhgAI\nRJuZ/MKH77OfYKvHPy/07kWjAUubd7uRFzZp8gWUGc5RQiLdcqPGZUoyGt/crtvZMXvQEqxU\nupLpMu4py6LcRYR+R2cgrGTc3reD0zPg9J/pFZwLje+08cLkNcaEo0ktLvtVXII1du6vJf6E\n9JrRcskrFrBfzTP/RpyHwu24pVb6aGXN7ltDBBqNgK8BOxHwZah+jM+gS9DnUP24P4aegapb\n2sQQAAGPBewms6RlHjX0PHaX2md0YMmex/33JD+fpFJ+ljwULeCu/ARSDgQE8Wic83Asd64Z\n25JxhWABP9Xn2Kx/VKWEpJPyF36DFfmNZVHHz0tKn3ZBh8YLpJu4bt8kzqNc/1nCaXON+jyQ\n6eTp1wXX8Tvcq8jX7eRLm4vwGpiMFwTKrcDxUq6+ynEJJ3+Kqpu5FZ2DqrLSuJBav/RimRgC\nRQQ0c1bjcTtAfFoKkozNFc+PiIfKN1hl1L8baTUzgz/+LF/3m3uIOUos4AJLlZoyTa7lxfTK\nPtw93tg6wiqdSXm3AL+6hyvKVtd58yTXfC4DpWDDPpNR/rc40hJJsUmJNx1Tjg5XJBpa1AHx\nkq7tsYJplXeUXDzC3mg78vNke+Toyeh4HR+JOEQZaCVf7f0PozmTsYFAIxKwnoxavhpbScdX\nDsB/HVrLLrsdSY/hq6qECsakHhHA3MACDl8YUkNNrbNRIWBunb5L9Nyoa7J4rGwFIa/q0QkE\nvCzKX7XYZ/7GRKxxMJ+hCWuwOBkpKW13J15W1PnZ7qE9joIFTIMquDnXsZnvBSeR/C0Q/X+w\nDSn4tmh8SvNFNvRIYcQCPGuRfWhMxLeMGwa8fxqvVu+dyThAoNEI+HM8syMrPLephElvTM6d\nj3tD4h+MsysXPTaYC+2a+kJAXblU2PNwZ0FwGGHx5JiRziXWEB8L0P215EjL6Crt1uQh4NTK\nI4aL1tB1PuYbd2DOBiTpEquhIK+dtEgt2dADPNUQDwKx/WRZ1PlHfVNZAfQyiIBHRdp80ynk\ncyH5U89cEOaRhYZB1jUdS0f7/ybB5oxxBBqNgO/heX0AXYZeXvLs5uPfC/19Eraq5NxgvFoK\nsRBNrZb+0ngzEdQ4MKkzBNj84RYqwzOgsxlx5Z0uDxnpjOprPyKPPF/7yTwLEfd4twquk/HO\npkfSnIlgsOomp8dj1aUM+vCEhoiGJKQRLEfcjUoIN1jCSaLBDzmHOLDvqBDwDl6GQObH5G1i\n3hUYUoiFTG2R6cvz/Q8j4BiT8fC30QiYCsrdjv4IPRh9N/ocqi7ot6JfQGslKwaQkLoUTeoQ\nAchuJZXhbLKmb8hiDVfugm71rrWJmdJLotzZw1EM7hsIF5e6OPelzgr5WBa55VTOaCx0PzN+\n7Q5Z5JveSJf0FWn4WHMpA9tJdm0vOdj8M5HpYizb6/Mud73SwF8kYPUsKIx7JUTcAt+N/BDw\nVNfSRh4m8q5R5PzvlScJFnrIH94x36AKBbI/AQHewYaTpZT4KPTuRE9qOASswFUjgMWJteTD\nLFwuehqdUuliNs/fm+r73ZioTeE8faaLfMupfELwqkrxBxpGbZxYvBGf5HMrlkfuof7SoCLX\nMpuDaES8v7+49XyeMXhtwTlkC7g9yt2q9dJgt1LlbWfHLJ4tbRlJPhBvSsoQYGxwxidH7G86\n45tntyHOX3xrGg9Jw2B08jViADTYjRqRgPWIeZ/dZ9A3oF9F6dYxMQR6IsAYaunH3Z8iBl2E\nFSV0WbYmm0Yscs0n8OP6OZ3GL60Ye+CBdEGHF7c4btlfEiJgxcGt2Gjo7/r6Oe9rYgFXLk/h\nTUl4sDBFerH1OTpER8WUNPYiNgnpktWBgENjYVS6xrtyYr5aItCoBJxieCsedT8vQR9IA801\nBFIEMGT1ybrwO6Fipn6OKpIZhm+wmIjIjF2RXoaVLJoIVbMN9LGAfcdm13mv0q1G+DB9IGDK\nUDHP1aRRJ3H0RaMhW8CVylJw+QcVzqYcwcJM4mxludmodPXy3qS9Ld0ImBcqX3CF08mfEXDy\nkMaD0+gErGeoH/bF6Gt1YGIIlCJQcLmSpSjR8xBBL2TmQ8UYuZZAwJAlk7aC1KTCxIqFgKMv\nr46qHwvlE3zjxAKOajIJq/S5pn5wTYk3HWOF60LPwah0QWvGd5w3vy7NY+rSFf8kDcCI+QY1\neafSdM0dPQSMgGPsr8Y5bvQeg925XhFg54cwY1bWJzOiN1I599YFHSrszmTXJuKH5T9UmHxM\nKR2/HXwpsWJbmFS1fCApcO/EAu6t0TCQ1EYvbox5VNIQql1eWDQdLGuGGgJWccoRH7WIezRq\nd6eqUwoEzPN+vvwKcEgbC0bA5eCM0WMj4DH64CzbI4MAXcihIsRlT+iu8eDyu9NlmXZBJxZM\nprj+dlZtuqGz/Fh73QmqPD/JMfwyLsaA6U0o9LAIeynzgIL5ZiP7P+dfs8zlHkgv5Flvgez2\nSY9H1o1Yc+7/4l3HO3reNxes9BoOa/S8hYWMKAJGwCMKt91srCGgHZaoEDWbmM8Suk24Fbug\nqRQTqyTuisYCTrqiw/fuhtydCSFksYoGSsBJtyrbWPoeHx8ZQ48i0hrsYSFggcC2lVfT75zM\nhlaIdsry76lFz4VSG4jQBb2I+A+wpOyZ8uvMAi5HZOwfGwHHz/BbOLZz1dh/n4epBNFqCHYT\nFjDbUva2E1YmkKyIMsnEhDQzQ7FYtvdu2mKfvYj7TiadAS1MhbCL8Re7lgvS/IxBd1gJuBwP\nniFDDdGMNuday8+NwDGbcfkVle7DrKy1NAxykWveo9J5Cxt7CBgBx8/sfJyqZ5eOvcdsOR4a\nAtri0WMJ6/ux6TKR8hTTSVhNgYBLiBhDavAzoae65sOx/s6BELZjrm5xolD53Ssf50MXtM5R\nqX+4cpz6DyXvdOfnh80CLkeA5/VB7un5mIW2lB1RodE0m4bec5Vu+lzkNpC3O6m0D6h03sLG\nHgJGwGPvmVmORxyB6HoI9QUqRy1JmlypO5eKEYLUjkW+hwXc7LIvH2yWud9bdC2EsDnn8tcO\nMJ1kzNBfSzpj81N2dJ2Td7r9oxcGWPZBR18a5f7O836OiVjvHXQig77Q0+lRqEjASpJ80T2e\nKfauDPo2dmFdIGAEXBePwTJRzwhsdR2f6XA5iDDeHGFBBSsYkthBZUgtX8g6JeLfE/wVdsT6\nx8GUkfSS2bDROtagpLNgq0qKPKUW86O6oLUOPjRfVcZLIs1NGjadLkP368gJz+8vXUuCRui+\n8S5qM1m/vaa3O5KvDs4ZAfcG0BgLNwIeYw/MsjvyCGgi1hORe5o1weqC1l6TxQlWXbnxc+RP\nCRhvsnNV4UKs47tYk7tTV9zqfVyrXdskg1mGk3RBR6uVAD/2MVdxT3YtoeucWcEjSsD0ZdwG\nZCO63Gc+jQ0aTaxby/RKwOSJZ6oPc5iMBwSMgMfDU7QyjAgCTJMNBNy12UbpbaOw8xUzpVvo\nL8b61UYc/uftrvMGSJl1peGbwqUXVOXn2rC3NNZYsh65qstCJKylQMB0X4cuzbxr3rv6q+sj\nZjIrWCAkPQEjlS+tBR782P1gckkZeWf01YWOvrrb9UzrpyHFS9bq1bliMhgEjIAHg5pd05AI\nPFUk4EKwgFt99lC6lq8UGBBFWJ6Em21z2S9hyRxJ3fQQJ7bhMoFrcF+xgYCbNSGIWwzCAk4J\nOLaAm8NnFcfWowM7Zn/7Le3hwwkjl3dwp/uecfMRlBaXDWvHYdg+Jpz5DvJWNxZwq2s+mq+A\nPUqj07hkEO+KgTYI0OySBkUgcjlNhkKDpcKPZ28I9/gEjbSy1tjvCQpjMo3G6ySsHx703sKq\nbAsQ0SAIuCNYwHxhJ1jApPHKhb75uJCjMfPHT4ZwrhqF7Gr8PH2mPW6/yLs92nz2vkW++cQe\nJwcd4GfwbuU127mPJOrKAmabt1k0VLLMiwi/iT7ybacqIGAEXAEUCzIEekdAG0I0pZXNdEzT\nqUnc4jpgrN+wkT8/rrBxBuSrbQ4Htrm/Ns5gUg7pt0BAjAP7MI7be756nqEfcxUV+h8LrlMf\nG8EYdwuaXSY0DnrGrs8Q8jyJmeX3jnTumAilzTjCM610bz4/yXd7o31YMvTWSucHE8azkgUc\nhjn6uJ4tUevHAmZmQfgt8JyC9d5Hvu1UBQSaK4RZkCFgCPSOAAScfmjBT+frNRltlgG98Vvy\nWMiZLJVRsJywOIMFTIWprQ0HVEHRjf01sqAG8gSu76SyX9l7liqfeT5y6593uRN39G5K2kqg\nkh9YQ6By0iMZOomyD2j2d20yl+eZZXq1gBnrnxJbL1GFCXmDywHvznYQ+tq+ruZd2AYedTMG\nTH5m8G4j6j4P7c2+st/j3ELv5i137PrFLMUeJxsgwCzgBnjIVsRaIqCv1BQ/tBAqXwZ/wwQs\n7iLrpQWSSy2nYhe0xjIHmIuZVLStELss4E6swFUDvL4Y/Zkwnlk8TPNWDKhnj3Cj/CUfShip\n3IYlXH1gFU+q41nT+KqVRAt5t57oKzUIWu9U3YwBlzRGB9TAVBmZvDWR8eMHFrmmU/oq83g+\nZwQ8np+ulW0YEFAXdCa1eoJLJRIIGLJgo46wEUdiAcdbRxKuMeDUCK0yTxGVrN+RyFg7/oac\ny11Z5YU9o/EtWdIIex1DZgNtCPRMbwRDlF9MoxG3gLknX0liLbfvTnaLfNPJjP1eTr7S53ns\nIp89J4Wk1Te/gvPfS48H4vKOHMe7Esbt+7iO875uLGAwCl3QjJcMmICxml9MI3M70mjto7zj\n+pQR8Lh+vFa4WiNAxcvknNjCpfttL6VP71liAXtNnpF1klpO6RgwX9xx21NRv6vN9z1ZhQr8\nWCr9iPuQTgQBe7qgMw+wCQeTsAcvWGqvJC12d4rSvA0+sRG8knxDdPF3jUfwtnyZoXOF7teG\noVZ6X8Z+z+LZvB4sQ+MLAslwvGsaJ+MyxxJ2Ynpcrcu7cQhxX8rEvT/3dU0ysa+OCDidqBaP\nBfeV9/JzTA6crzDwmld+rlGOjYAb5UlbOWuFgCyUYOFS8R6qRCFINmwKsh6inUaFwimJXxm7\nERNb3e50H36n4JoPi8N6/qVLrpVNGK5vcy2ncxYC9ttznZY3pV3ZPS+qMmRp1PknqrpLid7N\nAtYYHPfds8pkhj0a1uM9oRHCnRb67AFgyeYU+duG/cZlN6C1w3+/jecFB5eKxv81BhvtnYby\njIqNGqzYw2g0zE7PVevSiNuFuOuWRZ2f7esa7qXZ2XVDwOSH9xSkBjjHYbF3bFyT+bGu5Rkv\nltuIYgTciE/dyjwEBHyoADWxiaqDiVeqfKJFcYLRWhFG7Pe5nOu8Q/5CWAYUT5yhgg4VVhyn\n+182/99XIVRmL6ISz5I+caPDSUGkP2QhXT6rWByfDuk1uZYPohcOOfEaJUCFvid6kJKjctoH\n5+mlkVtRo+SrTwageK5raRCVd62GxhcJvZEn9Tf0ZnANjZpkMp4s4ElqPCz2LT/Xy1HdTZv0\nrPttaJFcsQFYXbrDG4uuet7TIIdS1vB7qOaO9OzvCE4pluF5V3PdeItjBDzenqiVZ7gREAFP\nnOSaj4pvFDZGOICKmA0S/NMQ5yyFM2nqeLqNn4/jdD5C+D2oZkP3SsBQzgLFJx1N7EkrNrz+\nUYUPXQqsYY5aqShL0+ZetZxINPhctjIpR40OKuZgXeLuCmaPDz7FIV/JbljdGyykGEhD5AGW\nT0G+d5HnYAFjlhafbZNzXybuqa0uu381ueCZ69p+CZj3h/evSFzVJD2scdL3GTxOW+Syb09v\nxvrol9Ob8c30uNzNu0IYQwe/x/Kuwwi4HCA7NgQMgYoIqAI8nO7kLyVnH4PAtOvVUiqTjZyb\nq3C+7748vbo9cu18YefFhK5jDXGv3YeZhLy5juVNaWXuV3e6/HVpWkNxsVbuId1F9HUGolda\nHE9NLbihpF2La2kXJGPq7q3xhiF+V/K3shZpDyYN7q3JX8Xu5TiN0iVAXkuCinGwZsKzBc9V\nILub4tO1XMQ6vr7yX56NZrv3S8BcrQZgr+9Q5dSHM7T40RG9S8VuebDQ8Eyvs5tpoNCDJHz8\nxvZIeDWmNLoFrNasvmLDu2NiCFSDQESl6xZS8e6n2FS2kFqkbsplkDCznV2oeOknZBZtD+Ha\n3ruguT7pvo5KLeD7V0T6EHstpPNvSoXPIwbrQ37uiV+f+xt9YavMHyW5eAT/+eD6FvCtUdkH\nUz4PucbfeS65Wuuy/zc53opf3fqhCzolYI6Xct1OcZwm1THVCBawr2KoQV3QPfJUTfrDEof3\nvWj144/fK68GQnQ+N5zTvbelKwv0EKUWcKXfSVfEce5rRAJWd8cVqH7Yar0+m7iP4H4NVeVn\nYghURABLRRZIELon9cm6h5Ojp7F2iutVnwrWcBKxy1EFW6ywuoJjH5V56L6mIpuGP3QTQ5B6\nR2si7UmeqPxm0d27Z5LoFO4XCKQmNxlaIuG3R5lXUP6QJ0g46cYfWsKDuZo8iGDLLGA3kecc\nGjKkyfMs6KMNR7Mn+OngmFim0VPkGyNPUh1Zcq3ei34t4LzL14UFvMi3vLXNN/8reS5a4+AQ\nSBWTn2cXNYNBNN+5uUKhXNjkREMfWiVwW/m5RjpuNAI+kIfLbFC3GtUWcuomUbfJceh5aBt6\nH9prJck5k4ZGoKAKMJHCl9iycGV8ECrdQMAQc54+FZFtmfgOyG/73qwCurUTC1hLXDQBK0iR\n1MsSG/gh/eLkbSvW5SkZl71DCdCAoAs6/pDEwBOs9RUab9V6Zb+GlEP5mcA2agRMHrZCFFN2\n8m72Dt5NhWT3JW/0mEUvJCXnXYgY149Uj3IuG6xdlgo9kJzXRLKJXL9zety7q+Vm/RMw99P7\nV61V3fvthniG9+YMsHkX744+FpL0Uqg3JXy+K313Yefssems9tJb0uCYQkNrKUMzHykNbzR/\noxHwmTxgWbnvRa9G1fp6EL0F/Rl6MipL+ATUxBCohECRWGHZTVQiGxSJSuhZrOOULBkL7ilU\nOh2sE/1cm2v+aM+zCvHsLe3VJTeBuMECpqKrmQWsO5DuBvI8F9KYrFm7slrwz1zss1/RzF3F\nGREhE/FM8q67kZdJKDtoivhiAiZvo9gFHfaDPhsSWT3dZX+LtXsSCLbQ6EoIONrW6XJ/JM+P\n8pwg5tgaZA3xL7pKFe020bU8Md87LTPqS0Ra/VrA3EvvX9Hq7CvBWpxjmdoiGh7X8ILzSLqE\nA/VQqCHAuxo3QsEgEDDnivnj3CdpWP6q9HrWPF/KdfsTv+LvhHMNI41GwHphSiyYis95BaEN\nuy6tIiIWWIJAsEDCMRXLpiaXX6+DmCgKoWLG31vFklSw0a4lCZZ6tYb4eSpZkW+wIvCnpF4a\nb9B+6tGNCVmoBoU0ojlJYodjGd/JzNUfDzrxAVzIXtfvneqy/1d2ifZ9ZrOQ8AnHpAu+EPAt\nizdShzQEolkonBLthu4R3zh+zhDyVsbnacD733FuB3BNiWdlVwbDbmY80KwIulfBmqyKgElA\neeJWI9NL18Q6aAr/qh17DlMwCzx0r2ex+H+G/5s8t0DAvORFC5iMLia/M+Z2bdcqJn8b+gbK\n0tvvpFecxtuJRiPgK3mAH0Vl6ZavWdPL8y70NFTd1CaGQA8E2BLy91Q0v9GJvMtthNCCBcxm\nEWuXufwvOfdFVFZsDyE8JeCKXZKkhUXqsQAjuvUCCSuNmlrApHcnaR+lhKlcyYfHENXXhrQP\ncRANyQy7UAGzKVhxHNq1sgRJxEJe1nBzWXmBgMnXKBJwWPMN1/kXUPUahCECejp45l7PMmnM\nexpNYfONIwnXFlqKTzSOXHSM3P63ahSZ9bsNJSnFDcBZI9YNHS/D4nnJ4i2RsPSKZxbIVkR6\nP3nbuc03vTrjWtKGCMHxLOmJLhtPSmOqu54z4fTq99pQLbnP+PY2GgFfw+MUAV+G6sfzDLoE\nfQ5dh34MPQOlVWtiCPREgLW9jxdc/ns6wywbuqAzgSDyLrOWpj2Dv4WrsIx+3PPKEJJ0X3sM\ngoqi5Uca84R8fAuVuKyKSyvGHGQg+fsvKsDwu+deh8ekF/1BBBMnGc1p9dlDB5l8xcuYsHPq\nLO8Y1+4SyiVCm5V2Q8NWjP9SLTtP4ybMNA8EHJNd13Uj7LuF58D+3u6EBLODwelHy13nreR/\ns/Kp/MhPHKz3zIH4t6kQHIeuc9wFcZzMCeUYlJaFcuorWkkDrfRMd3/G5bhvmCk6pfuZ4Tmi\nPGG8mYfTjYAJD13QkOhcVJ9I1AoAnmnmXIqfEHBopISMNbnCPHmwpEutYxF3Q0ujEbAe9iWo\nWmMa73on+nn0LPRgdE9UVrKJIdArAuyO9KROMoWT8dSOQMCRy4UKd3nUeQvbPn6m8sXpTkfR\nTCw+vWtl4rGANevXy9Lajkr5h8ui3D1lkYZ42FkyphqdAqHcw/rjB0pIeV6Tiy4Y4k26XU6F\nfPlM13R0Gsh4KFtsRi/RMRV76JqlAg8E7F3nrwiGiMJHLSC3uIGTXjuS7lKX+y7P44glUe4O\n8qduZbIetUOwPJpgiSYWsMaKIwg4kFIgZSJcALZ3p/nl3CdmuJbXpsflLuchpi7CKj+fHlNh\nh+dHfvSOjIBkwnOhPN0ImLxquCBDvueRJ55XcQUAwygxyZLHJK+emeJNM5VZWlVF65jzFXuK\nRqBQdXMLsGtIUSvyPpTWtpM1cz36OMp7ZmII9I3AEpe7k1cl7HTV7txy7/InL2NHn76v0lmf\nkLVbwCzk+yHhUCmF62AaKjSsGv8Clfls/PTP5R7oP82BxZDVXnLFXvjvxyrmnl0CcZRVtl3n\nBuxje0LK0kTxioTBpKRzSOdFSotKOCHeuIKHvSCzwjbCwzKezqSBM+D71uICWlVLow7VE4j/\ncuwWyJ/EP8EMbZGyJCxXAjct6XpYATScLsK5Rf5UaIhMT/0VXFmaCaFXOJsE0V0XSC3rmot4\n9h67JmeCBcza8eI7wTyBD5JysSxUmlshkvBe0UBhKC+1gCN1xed5n58gbIZyw7hf0QLmGTe8\nBVw+DlqTJ1bHifD7cKeiJ6FnoAvR76PHoSJfWRtnoXejQ5GZXHwBGrrRqkgoVEZVxLMo9YAA\nNekS13ldyAr+pS7/22qyxZjxp/h04XoI6WyRUuSa9+NTvzfpWgZjU8tgnY6puDbzBYcN8tdS\nqOE3B8YjUfIwhZd+GRZv+X3SKEO+NV2OoVxU0IEwdvbqqvWf0t1j0dKdHLVyS7gnoWoca7MT\nLW/509PsBT3kTNQkAQ8Rg1iyLpvlM4ekyRKmPAcLmGVmP03DifM+yEqTpj6kMMqteSY9BEx2\npbz0Rrgf9zhZHkCjgF2qN2otd/mp4TjmuYXnwv2KBMx9vkCZAjHrnnQIPMeKgM2hxcS2ppRF\n1nwOV40FkTBuTMCEpe85S7RsDBh8G0qOo7TfRG9MSv1zXL1Ib0KPQf+C6lwYr8AdrKhhox/I\n7Cp1RMZzBlsYu642CCyP3DJSujhNjUpscernh5hUaOmkruFZfsN9Si1gbl94gWU13SwRKsma\nETA/hFAuKuFAwCzp2Ysqu6Thn+617N+bYEF3ZZRYwNFfqN2xoEZfGNtcleavPDfMCcAqVjnU\nDR2V4asejVgKLtODgLWhxQTXfDrXTqfc3Z5Del25S7znWc42IgRMuZLnl3RBx5Ooknc1zVnm\nWa0I0BHvDgQcr2kmnzSetGQrfENbRonOQ86x0INQhlV6pnHckh9CQxRaRPtf6CXoXPQQdG/0\nQVTyF5QF9e5E9AfoYOU5LvynAVz8LuIePID4FnXMIpCjGzqugyClYkMP6yeQHks6mPTDqGxx\nc4PaFvTJyG1p855JQpqJqgox2sxa1o3MXC3eiIqz1Nophg/GQwUTEibNQBh0d+/cFPcuh+Sw\nmlKyP1ABT9Kd2RbP/ObStLt3MHeu7TUF17FaGAFa0gVdmr4243Asy3Hqgt5ceiYmnzikkgVM\n2PvBZjc9DL6PW94T0T2p5IhrngsbulQ8W5tArRGf4rK7UaZTlbemxAKe7yp9zzr3bCdLkLLh\nUWvWcxPd09onWx8n0SQ0zWRPv5/cNUO6Z2OlNnkfS6k0mgWs8qYzDdU9oh9LeatzJWHFbhL8\nJoZAzRCgMitaRFgNagQGITxYFYSllbDez+GSIklQQW6C6cp/Aykp1uL+oVyULxAwH7QvswJT\nC9hp1vdWGI4sFVZxY2YFR3LrQpaxUoL8aYuzSuO0hEVsSOGn0ptQxFYZX+9yl1CgK+QHg2LZ\n2YzibfoWM2XU2u/QOwARlz8HXVZBomdppM2pcKJmQWw88o8MTfySvIV3lElUk1k3NmOSy369\n6yb6FKN/ZinfTqZiLZabsEO5bisW7h8p9TUiY64J7wH+Yt0KLsVrutJsLF+jEbC6nM9F34qq\na+urie6Jq27g96CvRq9FTQyBmiOwpISAqZwWdt0gHgOlEk4J+N6uc7X1UTlq3aref0ihsIkW\naVnF76drb+Na3BXSSSrcaHacXvdPH1KxB7InT4oXrEu6atlXO8p0ulHdhKN78eOu8L/QNwHf\nlEtMyhDKdhmX74ZlMo5/JXiv5fw0Xcla2ZPA5X/4DvPLOCxOZuJ8t2vL79J17FdDZDt0Hdfe\nx3uIsatNRGLCJP+TCi67iDfmzendci56d9751/MS5dudY1KaP4dzD/HstM/CumVR/iomsZ2P\nX42W8B50lnyMhHgN3wXdaAQsYv0A+nlUP/Y3oa9DH0LpJgnkfCZuhR8ZoSaGwFARiBxrJj0t\n/7ARPeOhsWSLm/bHO2thBf4uPVdrt8N1HEblh3UiAnab6fYNFT/5eiYOiyZjzX1ryPdlBnTG\nNW8fpxM2qpCX35dPGxnc37+n1btWwrWlYbCImKz2LHnxzC6HaOpHmFR1NPr38hyRX83yZdJR\nlGE8uEfdAQldzjWXcT7p2m/SWmz8/l8pZ7AwlSZhVREwRL2GuEmDpjw3tTmmoSHrXI2jMHaL\nyzahxd4KDZH4Fa7jIS27C3eEhJn5/Q3CwycvKdu6kpxgAcefceR9K451gFXDEzC9Tw0nP6DE\nP0EPQ3dG9QNQtyA/IvdnVN0lJobAsCFAxYY15FdCPppvEESVG+FUbJE2fmBK/vBZB09E7unF\nPp5hHe5Do6Dg3b9vc7n/ZpLUaXSDX0QlP+S6YbFrhmCij6iAJYSxO0fPEhSsQcJfxtjqERwX\nLWDlbyef2+GpyK3RtfUuyi/dsw/z9F60LN7Yp0eWKSebthRnEgerl7B9uiL6ixmL/33XcV++\niIlqaZduX/EGf47ntqPeQ55VWMkBIWMBd92Tcxt4qAR1F8pEe06iiVdFUZ0aLGCuKxIwca0L\nughRY3lyFFdkq5apuqHvQ+9EjXwBwWS4EfCPc4dbqaQmwLhNyd0YI/OdVErqmdESjWG1Dugu\nDJVfep9lUcd59Psyzpn7ke5PRTlkAobc2R85fCtZFlMY5yRplVMELFGjV+dkDWoMOJRdYWOF\nfJXXRCgTY9Z8cSoNKHUZD9VuWZRT0nM9MOd/s6L67z53QPZFIovTrPVfn1q+ej60zzKTmYgV\nhgs4zkPQxbkMpXemYZCQaqkFHHY2C/kl36HhpWtg72F9x0vzVa/+RuuC7u05XM2J43o7aeGG\nQC0RYCOP46h8LlaazEzS3ANVcqrctjAwGypwrKFhrZwgg1BRFpgBrfun0h6pN8jfWQsLmDRn\npOkGP40NCHkixPwX0r8C0vka92K2bPRGzkPMkRomY1VW8wyf6i3zlFF4H7TIN38Ef0LEiu1v\n0F8aQgOw9gtFi1LXDpMUCZj0Rbb6UtVkyqi9yq8l3xUJGFpNCDhs85tmjWkG6VhyWD71FGk9\nQlvl6TRCo7pDbuU2KnBWbkNg0AhgBmS8CLZF/bBTnmf+QcTsYColzUhW74zM4mEm4NjChu2L\n47FpeZhY82G227whPR6sS0VdtHbwNy2O18WTnL+CsdTblC6fQYR8/fF48xDytxU2FgWL8HEe\n69be8h43eKIWcEga+hoHj6bREPsJVtCxeddZNQHznmABp5Pbervj0MLJZ5GA8a/lGUHATfr4\ngnopVM6KBAwO2hNa1u2PunKgb2jHY8B0ZW9H/lfy/A/uOt+4PiPgxn32VvJRRACTYLMGxbB4\ngwVMBbcjldSGLa5jFV+OuQkTZ1gnIFEJqktUNakmH3YTyD9Y4ZjlfKGwcpdqtwvSA/op+dbx\nxyjHHPZR/hfKxFin7qJZu9EOBdcyV11uVOjFruY4HxG3dGxX2f9eyOmt6s1d5jrOC0XrJWMQ\n0ua47L41xiTS89WmFQ9BWr8DkCd6ubRCcNioZLi7oOm90LPj6cVr0rHa1UujZxd6LSoSMFeE\nhiMbdhXPpw0GfVuYV+RYrofQTYSAEXD8HnwL57HYa38NgeFHgB9eqKjYY3dqm88t4I76iky7\nxmExgo8Z/hwUNlPxF2CBkI/S+7HWtVOM2Er90J6ScWmECn4+ibPDZJddRgVNN6x/drHLatlN\nsOa5z11UzCeyycg8dbZiwMmCSqSLjAlQ1+rYFB4eGZdWFOYr0TUrVNOvI/nVNETagHfF0igs\nfax4XS+B4BRblL2cH1pwmJcQ0XhSgyiC6CM2cXOvJNETIGN6a7T3dTyHoPxGaeOKBkdJwy58\nLWlP3vXfEH9v9Gfl1zXqsX4NJs6dDwj3GhCGwEgh0J4QHxNbptAV/c9UXAwH9+wOHq78JF2F\nGzFyepAGlUJYI4ypExijmjywl/McyqCdoHYnPhpN5ThMvIJ8/qw02D5xe7mwR9EC5vbFBgDm\n9tglYBWsDwGXZGw07ZaPaGiF1k2PIYA+kgmneD4QY9ds4v7iD/T8Qtd8qK6BRJmkqk00Crdz\nuBNlWIi7ucPlPsHSq08qTrmwtCiUky6U4ixo8qtxfizoSOSrF65oHZdf32jHRsCN9sStvPWB\nAH10WBFYgp49gNP9doubcAx7HqlMIb7KH7un4g1d0LQMqu4hg2SLE66obOeUFODXW13npTqm\nnGpkSIoEzL2Cn3Os+82tjk+Pv79MPCtOdqOsbF3pw/GKQTW6QkMlLOsZDqRoCF0Vp5v7AhMG\n1U1+j47J903ouVp2hUm8stK9aWRpItoT2vK05HxZw6ryDOqS+A3jNQJumEdtBa0/BLQeuGk7\nyDCpTLs2qBjuvGLF/Jp9mf+p0n3YgjAQMLXmAAi4a0en0jRZ7vRLlhQ9SfWt7sxAzCwsLVq9\nVO4X6xyk/RBdsfDR+JRSAqaE61FtyLKJHoiA9cBKnefRdG3pOLBrq4odJs/xwJQ/GoqZ0DDi\nPf3xsqjzj32loG8nL4k6FpTGoeyJ9Z+Gdv/8ZRraiK4RcCM+dStzXSAA+WhLSFmFWgN7K9Yg\ny3JGRrBQnm+POm+odDe+wRu6oFkXVHUXNHsTM2YoCTt8xd7wt9itrG/Gvp2gx9ujMIs2nF0a\ndf4FYqY7trdlLSHamP/TVLLci8aGup01Dn7HYApGFy7cOExd0HztiOcRJnjhCVY6S9VCdzlE\nOqiu4/LtOUl/UOkMBqt6v8YIuN6fkOVv3CJARbpWyzIooCq8+9uj3K31UFhYN1hl1PJVW8AQ\naCBgyvS90jKwnWTofiRcXZJ8+tb9svR84uebseO7UqZFU+yCBqv1kNnXKftHKmBRRVDY2GJY\nuqBbi7PynWtP8kzGkyVSg7Ncy6x/Zn37B6soZENEGcAPrCHwsEIaAiOIgLqgHV3QkSzg4qSV\nEcxAxVtBFp2qGLCAqq4fkobEHZtd7rwpruU9lEcfotfmGsn4X/ErQmE8sfTGxNWylHG9NIU+\n943FQXLGfdl57N5SDAbox3rmm8oDXSZWxU2Ytj417vZg8lXSPa4PSszw7iu8Fj2eXRVJks2u\nuQ1LXMcBpJu8E9VcPb7jmAU8vp+vla6OEaALWmNj2gFLFrC6FetCqIBDFzRmcNUEzBznOVS0\n/7+9M4GWpKrPeNVbZmAGhkU2BxyGGWQVzBFiFI2GBBdiREVNRDkYYjBBc3I8RzCauBy3GKMR\nPRoVjGhM3BDBRIWooAYU3DjIwaBswwwCM+z79t6bV/m+6qr36vV0v+ru111d1fW7c75X+617\nf7en/vW/W92+JQhcra73ivC3zoxW0odtPPRIHvFc++98ZkPNrjTaBlh1uHP5zhqkeQadr6XX\nrw0mj+j8qs7O1A9Rfe/iclvwQqh23dNvChfMbtVZhDrL35tuxKlmBozvAm4Y4AU42IBAcQTU\n5mujtFxGSs+9sDQGWF5QXAUt6zu+bzRxVGdE4k/X3aEHrI23qlaj7ybXJQa4Mc+zvOI5Q5TG\nqzZnfbRh+kPp9kguGx8uSHoXL62znYYAuQ1ZqINT+81qTOPSkzjv6Ffc+h3FBlhl37LndL/u\nU8V4MMBVLDXSPBIE9J/PxknVtP4QQXk8YL0NxAZYfw6cCMZ+rC8nxb2Xt4Gu7rHzBjrcXSYh\nfmhrmsEzdO4NPl89quMXCxmL2ANWe+A2BlgfIbhmlHtAz3ObOsXrYhEb0Pn93a2lnaP0+3Ht\nSV+DpiGNe0Ar0r4ZYP2O5E1HM6r9OL2viR2ByDDAI1CIZKGyBGSA4xmNlo8Fs6XxgJWo2ADr\nM4FpdWSm+XKe9b7B5FNtoNdEwS56wNoYzI39lKFNvJ20DbhxTAY59obmY6nPmrg+qOp5fUko\n2yGr+/w7Hl+lqui+G2A1iyS92Xurbm6VG48JVg3HbhvCrd9udbzO+zDAdS598j5UAjJS9grl\nAcdtwHquliPIA47bgPWId9q0MVctuSCBmsXLQ6hsrV1tOampJjMvEdGmxslTcb5kLOJOVqpu\n3sYDbpw3+n8bk1OEFyunNy8ltw2DFqmMGuWzlLi2vXbcU1D6s5Rv2/ZY73tUw7GgTbn3mEbr\nyi46WYxWxskNBEpAQFP0BQfaCy5TFbQs5Iwtqt7ONX2gl3PVkk3IGt+MXR5/pi5YpvPjuZ8b\nJ03/RkNKfyWLrP5HDu7pHLonUrLd2Fu3vxuCqWOFQc3/Sw3hBjEdhAesKujwTnW4Sl6glppO\nrl+MAAZ4MTocg8AACajaWd7h2PrkFhnvcYA37SDq2xMDrFNjD1jVkmm74IKr1QS8s14gZE0m\n9bF2e8DzHck2xt8Vnj4svUDH1NM5eshDWtJ9tVz2xfjK9AazZ6lT05/1n6GnRs1+SKH/dyDG\neQJUQc+zYA0ChRJIekGn94w7KaUbQ17GVdBqr4w9rFl9q7hVemR8n+X9GpC6QsZANddZD3jh\nFbI7Gp6UtgsvPMZWLwRCzSfd505Y8SxYY89Vah7oJUVc0z0BPODumXEFBPpEYG5yCsW34LN8\nfYq/x2g0AYM+jaDa4vB1jqGdB6xDR/i4XiRW6E1+0Z7c08HM1yaD8et9PqEvBPTC1t9OWPsF\nk4fIqB+rF6Xz+5JCIsklgAHORcQJEBgYgbmOV/qM21wP4oHdrYuI9SC+Wl7tM3yJ1uPe0Nte\n3pgz2G2RqrbU6Jj2Q6n0BZ0b1Z1LIvSDgMrGv5e4iaAf8TXi8FhulXbmE5H9i5uYWhGoexW0\nf8C7S/7VESBQKIGxzBR9+gmWygDLC/p4CkNV0W0McOMrTmPBuDtrTcooZDphpVezHAQBjanV\ncK6oZdNA7/cL90quVe0HoQgCdTTARwrsuZKHRfih5wHnXl4rfUxq2eFE+wkQ6CuBmWBW/Z0a\nQW+ApTLAcmZ/nqZN33iNe0On2+lSaVYHaFdBRytVTe3RS6XpSJamcVSXetm5Ry9tq9Qbq4+1\nmJGdEQcMcIPDwP/WzQC7zepi6U7pROn3pEOlo6V3SuukqyV3KCFAYMAEZrakN9ADtembqemR\n4SzVrvtbeb5xm217D7jxTVoZYnti6gXdvhPWcHIxunedCWZkgINAXejjsdhLzel+0cSF+g0+\n2fGo9gMDvFSgHV7fx7enDu843NNO1u3t5drYNofLteOr0nekY6QLJAIEBkZAQ3du0+NO81iE\nEzMlq4LeqG/2ro6mj9o+mDxXBnbRKmgZaBlgf9FpfhjSwKARcUxAnlNsgNVd3QbYtXi9B40n\nUy3HC/UCdZ0jUY91DHDvNLu6sm4esNt884Z73Kxz9GJJgMBgCXh2IM2965c9hamSVUHr7SCM\nJ824Xg/mdgbYn1GcDYOx2APWw4Q24EZhDvzvxuTrUaqB1jSgSwuqd46bGNSMEHvTepGq7XSh\nSyPZ/dV1M8DnCdFp0oulZu/fDxFPlv4qydXUBAgUQGBmk42Y3OBSjr2U8X1YVZPbGmA9pe31\n6o/6UkQ76ZxxnUsbcAG/mPgWyWf9xH/JzWU7BpMHOU6V3y76Lf7i0WD6nKKyUff71M0Au1rZ\nBvgLkj1hd4Lx0Ii7JM9Verr0WukaiQCBgRNQVe9GOY67Jt7mwO/X/Q08hWRjysmma+MHv+zw\n3TICqeeEAW6CNNjNaEY96devj5bdsE/U26Qc66OJZ+nl6edOp1+itLhuSxj3kRls0ok9JtDs\nBdYBy9nK5FckVzOvlfaQ7pRukX4pqf8JAQLFESjzRPUysG5rPHJ9NPn6G8Pps1IqsrhxD2hV\nW96tc54jI+wHeOmq0dP0juJSHqur/I+S1uttaB8t405zXebV12dDXhNd9lzWl0igjgbYyB6R\n3NvZSoN/wH6O9Ct4OFOnfFsO8+hXQogHAr0SkFG1h7tSVZPPVRxzBlidKVIPeLOOP9PxTwdT\n/n9FKIiAuKvlIvgd3y4MJuRIzPRggEM3uWUDBjhLY8DrnRqIASej0OgP0d3+RlIHwuAfJbd/\nfE7aV/Lb/uulr0tLCfvr4uukfhr0paSHayHQE4HZYPbucX8XqWnaQ/2wZYP9A49uTX/m2ocB\n7olyzxdN6cqnNq6OZIB7CeFOTVfNzc7WtJ/NARComwFeJ4ZXSv8l7SV9WbLh/Yz0Leml0n9I\nV0gbpV7DDbrwUMnT83USjtdJ7+rkRM6BQJEEHghmLt85mPyRPOH4wwzpvVX/nGxHt6QGWA2I\nGOAUUAHLSC6vyiV+huvb0gfrluf3cNu4KSEdDqeyxAPuAWKvl9TNAJ8oUGdIb02Auc33Kun9\nybbX/0B6vjRX3ab1XsKvu7jIE4IQIFA6AvqG4AM7R9FFStjR2cTpwzn6ApLD2M3p/o0Y4BRF\nIUsZX7cBu/eymgHGntbLTdW0oJ7sqsfQTGwqzx01Bph2/F5A9niN65bqFFYrs7dlMrxF63qD\nXxDsvdImuwAJG3UmoAfzo3pAL/g/oe3YA46CqQu1rrHC0VY9x10lSiiIgLjHBljm8zcqIw+j\n7CUsV9kpquBumXKN/505s5dIuKY3AnXzgF3l7CrmB6R9pQOkw6XPS7+Q/Bb5cunfJAIEICAC\nqt58RA/4BVXQemj7C0iRe3CvjyKPInCTDqFYAqkHrFFs0T493nq5POBb1da/ScvNN4ZLnFWr\nx0TU9bK6GeBLVNAfkzwVpds6Xi0dLXkc3O3SntK/SD+VCBCAQIPAYzLAa/aNJo7eFM78wLvG\ng0gecWPYkQyx/i8xe1LRPxYZXRlglYzmMFAZ7N/j/f0d57fMBjM/vCkMNvcYB5f1SKBuBli/\ntdjA2sim4SdauVbyG7wNrztgESAAgYSAPCN3rlo1EYSnqrHwh24y1H+k7dV+lbQXhmo/jLJN\nO7ArgIDKQJ2wHML7VRkR90rv6raax1Jlq2KcvVnuL8a3K3j9OblubcDtqH1KB+wBY3zbEWJ/\nbQmEwda0Y84r92s024jF2PZ68Kf77QH/traAhpfxuApaBvRBJWGv/aLJU9KkrI/GX7Y2mnhB\nut1quWcymYomUmHoUStABezDABcAmVtAoMoE1AY893WcsWAy7uyjB4fbgONhR+o5+1O1CX+n\nynmsYtrl/SYGOLxfnuyeKpOPq4YidoqjYPx1mqbyxMXypTGSyxvHMcCLcRrksbpVQQ+SJXFD\nYCQJaG7Wh9I39dm47Tce+iIPOIo94BvDmQ+NZMZLnim9AE3PV0E7seFyDfPYR20Bqo2I9pRR\nXvT5rnHbsQFWT3Z6rw+prNP/V0O6fWlu+2ml5PrSpIaEQKBEBJI24DhFWv9wY+L/SAY499Oe\nJcrFSCZlxkOIZoOtGj7UCMuCyb3XRME6bWnGv0hfGmwfZLzjWbBkyKmCbo9poEcWfUMa6J3L\nFfm7ypUcUgOB8hAYC6ZVBb0sTpAM8NOXBeNHqw14O3W8UtsvYVgEVBaqgo7kvWZnr4q2mwjG\nD9GxFTKsuy2WtrFg4mU6Z+NGTeSx2HkcGxwBPODBsSVmCIwEAU2aPudhOUNqZjxYC/W6neuE\nNRL5rGAmPJ/BlB7ic230Gh62y2wwvrKRl2jV4nkK91dpXqZWYzzgxUEN7CgGeGBoiRgCo0FA\nXWznHvDOkaY9PEoP7p21igc81CKONPwomNILkXtBxyEKxr4+Fswe6A1VMSeGuHGsxd+18pRv\narGfXQURwAAXBJrbQKCqBPSJsIfV1qjJHhpVznqwH29pOx2GVNWsVT3dMsDh47PBtD3hOMig\nqmjCtfFGEE6quqLRdtDY0fx3lWbAuq95J9vFEcAAF8eaO0GgmgQ0EHhDMP0EPdw/IaP7vUYm\nQnfwwQMeaolG96kNd0oNwXMecFI2e6UvS2sW9YLDHVSmC5oXhpqdGt4cA1zDQifLEOiagIzw\njeHU6RqG9JHMtRjgDIyiVzU+297r4xp2tMAAq3PcXvKC73B6Zhb5SINc5R1UfY0BLrrgMvfD\nAGdgsAoBCCxOYGMw8z15V5cmZ1EFvTiuAR+d/ppmwXpjoxOV54VuBBlVGeAoNsCqf27bDqxy\n3DHKDGFKr2dZHAEMcHGsuRMEqk9AnrCqPTc4I6q+xAMeYonq4wm3bwhnLnYSZEzfoUUyHWi0\nh0rHn1oN0pnLvN4i7KDZsvCAW4ApahcGuCjS3AcCI0JAD41k2Aofby9LkcoQf1BmODbAejEa\nVzV0PLZXhtkTpmwT1kaBxnGH46rGxgBvQ6e4HRjg4lhzJwiMCoF46kJNUUkVdLlKdG48r6qh\n1Xk9rpZu+ZUk1VfHVdPjGOChliAGeKj4uTkEqkfAPW+danlQGOBSFV94mTzeuO1XbcPuoKUm\ngtYe8PLEAE8HU/EHNUqVjRolBgNco8ImqxDoD4G5r+dggPsDtC+xqJf629Wz+RuNyEKN23YN\nxXhLD1i1Fyt8ngwABrgv9HuLBAPcGzeugkCdCeABl7T0VfUszzd6aGsQXaU6isfG2njAUTAZ\nG2A1AC+Y5ayk2RrZZPExhpEtWjIGgcEQiILZKU1HGWwNtuIBDwbxEmKd/a6M8OZN4cxl66Jl\nLp+WHrDmjJYBDoO78ICXwHrpl+IBL50hMUCgbgRiD9geVt0yXvb8eljShnD6o410Ro+pGnr7\n/aKJ05JPFM4lX1XQ6oSlLylpWNncTlYKJ4ABLhw5N4RA1Qk02oDV3ogHXOKibJRP5OFGp00E\nky/MJtUesIwz7b9ZKENYxwAPATq3hECVCeihkXrA/hoPobwEVEMx5ipoG+G9FyZzbJX2MQZ4\nIZTCtzDAhSPnhhCoNoF0GFIUTN9e7ZyMdurVFvyIvOCV0naaw7vJAAfPUu6vGG0C5c8dBrj8\nZUQKIVAyAuHjGm86uzEI7ixZwkhOhoBmw3pI5aQPLgTLNOXkzplDmroy9oivz+5jvXgC9IIu\nnjl3hEDFCcxqsoex6+nAU+5ilOF9SN7vruosF2p9wZSUMs4rZZypgh5yEWKAh1wA3B4CVSOg\nnrYXKc0HVS3ddUuvjK87WR3XyPc2M2LpS0gBBnjIPwoM8JALgNtDAAIQGAQBVTPbA94tiXuB\nB6xjbhtmEo5BgO8iTtqAu4DFqRCAAASqQkATpsx5uDK2CwywttU2TC/oYZclHvCwS4D7QwAC\nEBgAAQ0zmvNwZWzjqSfT26j9V23A8wY63c+yWAJ4wMXy5m4QgAAECiEgLzc70cY2HnDWQBeS\nIG6yDQE84G2Q9GXHjorlbdJkh7Ed1uF5nAYBCECgIwL+XrM9LHm7nrFs3gBHgZ774YR6QjOT\nWUckB3cSBngwbP1jdy/RZR1G/8QOz+M0CEAAAh0RCINZecCadDIIzpEhfo0H/+pftPvcBxqY\ny7sjkAM8CQM8GLgaJxkc30XUp+jcs7o4n1MhAAEI5BGIPVzZ3A0atz2xl3pEb9HkKaqW8/SU\nwUwQ4gHnERzwcdqABwyY6CEAAQgMg8CspqJs3HfsBi+3D5bFNW1qG44NcBRM8TWrYRRM5p4Y\n4AwMViEAAQiMDoGGh6vezreoJVjTh86udt7mDXCAAR5yYVMFPeQC4PYQgAAEBkFADb5xFbPm\ngdZwpPBWGd59fJ9lwbK4Q5bahjHAgwDfRZx4wF3A4lQIQAACVSEQBlOxAd4aTD+sntC3yCC/\nd3UUPEnecFIFzfech12WGOBhlwD3hwAEIDAAAg8EwRYZ3ov0kN+s6G+TB7zX9sHEGu2LDfBt\neMADoN5dlFRBd8eLsyEAAQhUgsBdYfDgXcH085zYdVFwt5caG+wZsWSLo2n91SZhmAQwwMOk\nz70hAAEIFEBAQ4DvldX1pByxAdb23DSVBdyeW7QhgAFuA4bdEIAABEaFgMYC32PHdywYXyEj\nrFmwwgdHJW9VzgcGuMqlR9ohAAEIdEBA3q8McMMD1vpyrc59KamDyzllQATohDUgsEQLAQhA\noCwE7g2mz9ewo5uUnhXyhvUpwggPuASFgwEuQSGQBAhAAAKDJHBPGDwgw7tFWjkbjPljMXjA\ngwTeYdxUQXcIitMgAAEIVJlA0hFrN1VE0wZckoLEAy5JQZAMCEAAAoMkoLbfa9QR61At3RGL\nXtCDhN1h3BjgDkFxGgQgAIEqE5gNoo0yvPto8K8/kzpV5byMStoxwKNSkuQDAhCAwOIEHpcH\nvFwe8DIJA7w4q0KOYoALwcxNIAABCAydgI2uvN8QD3joRdFIAAa4JAVBMiAAAQgMkoC+ipQY\n4GCZhiThAQ8SdodxY4A7BMVpEIAABKpMIAq2TmkYkquf7QVjgEtQmBjgEhQCSYAABCBQAIHE\nA46ogi4Adie3wAB3QolzIAABCFScgMYBxwZY1c+qgp7FAy5BeWKAS1AIJAECEIDAoAm43dfG\nV/dZpgc/BnjQwDuIv+4G2B+m3l1SswgBAhCAwOgSsAHWV5BCaWXiDY9uZiuSszoa4CNVNudK\n90qPSncky2u1/JjkeVIJEIAABEaKwHjS8UqGeAdlDA+4BKVbt7mgjxDzi6UvSSdKd0qelHwn\naY3kfVdLB0j8QAWBEeDVegAAEHZJREFUAAEIjAYBe72Nqj5/DYnnWxlKtW4G+GRBt5f7zhbw\nL9e+r0rfkY6RLpAIEIAABEaCQBRMyalY5va2HTHA5SjSulVBu833sRz0N+v4+pxzOAwBCECg\nUgQ0B3RSqxdOhsHsI5VK/Igmtm4G+DyV42nSi6Vm79/tIqdIr5JcTU2AAAQgMDIE5Pm6z0sc\nVB39QLrOcngEmo3Q8FJSzJ1drWwD/AXJna3ultI24F20fqP0Wkmf7SJAAAIQGB0C40GQ8XrD\n+0cnZ9XNSd0MsEvqbOkrkquZ10p7SO6MdYv0S0k1NQQIQAACo0VggwzwfNvaNAa4BMVbRwNs\n7H4T/JXkHs/NYYV2uK1kpvkA2xCAAAQqSyAMpoMomta0B5PqhIUBLkFB1q0N2MjfIbnq+R7p\nTGmVlA1XauNl2R2sQwACEBgFAjK8aTX0faOQn6rnoW4e8FNUYG+SPiD5h/gG6cfScyUbZAIE\nIACBESYQ6rkXjW8Mc0eDjDCD8mStbh7wCUL/KenD0iel35Vuli6QVkoECEAAAiNLQD2hZYBD\nz/5HKAGBuhngXcX8rgx3d8t/RbJ9rpb9rBFwXPHE5x0s1UGRAAEIQGDQBCJ1OI0wwIPG3GH8\n/TQ4Hd5yqKd5HPDHpcuknyUpsRF+kXSJ9EVpubTU4M6Gnlu6W8NKD+ylkud6CECgLYHZIHpL\nGIw9qe0JHCiUQN0M8A9E9/vSD6Vjpf+VHNwp63nShdK+0lKDxxM/Q5rsMCK3TZ8lbe3wfE6D\nAAQg0DWBm8KZS7u+iAsg0GcCuym+5t7PvoVfSP5COtwbBYZn6l7qoBhXWRd4W24FAQhAoFIE\n3KznZ6WfmZUPdfOA0wLLtgOn+7z02N+zsztYhwAEIAABCAyCQN0M8EGC6E8P5oXrdQLDkvIo\ncRwCEIAABCDQIQF7t66+yNMrO4yvX6dRBd0vksQDAQiMMgGqoCtcuq9X2t3T+OmSjV67TxP6\nHAIEIAABCEBgYATqNg7YbbynShqPHrxZcq/jVrKHTIAABCAAAQgMjEDd2oANUpORBydIz/ZG\nyYKrVzoNHmNctxeoTtlwHgQgUB0Cdow6dXq6eUaWnkAdDbALxV9CssoS/FLg8GBjwV8IQAAC\nEFiEgL9YV/ngqlhCEBwpCBukYfZ8dho6nbjD554huTq9TuH9yuy3pctqlOkXKq9PlT5Yozx7\nXnb/vt8lba5Rvv3/2ZMCnVOjPHsSor+Sju4izza+V3RxPqeWnMDtSt/LS57GbPL+SBuutqlb\n8EvSyTXL9NuU3zq9cLh4d5NcJXmoN2oUzldeP1qj/Dqrfyw9XLM8z2WXNsQ5FKxAAAIQgAAE\niiOAAS6ONXeCAAQgAAEIzBHAADdQfFoLz35FgAAEIAABCBRCoK69oJvhurMHAQIQgAAEIFAY\nATzgwlBzIwhAAAIQgMA8AQzwPAvWIAABCEAAAoURwAAXhpobQQACEIAABOYJYIDnWbAGAQhA\nAAIQKIwABrgw1NwIAhCAAAQgME8AAzzPokprnjt6JOZC7RJ6HfNdxzynk/PX7Tdex7KuY567\nfOxxetkIeA7vdWVLVAHpWaN7dDpfdgHJKeQWK3SX1YXcqVw32b9cySkkNbvrLqsKuVN5bmIn\nsI7PsvKUACmBAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAA\nBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCA\nAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQ\ngAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCFSLwHi1kttRasdyzsrL\nc97xnOiHeniyzd3z8pR3vE20Q9+dl+6lHh96BlskoE6/7x1b5D/dtZSyzbs2vQdLCAyEwJsU\n66+lh6VvSrtKVQ8HKwNnSbdL10nvlnaT0nCUVr4nPSb9n/QcKRuqzuRYZSaS9shkapXWPyfd\nJZnLO6VsyDuePbdM6wcqMZdKLsurpDdK2ZBXlnnHs3GVZX1vJeQTksvxaunvpKwxzivLvN+/\noitVOEmp2dQiRXllt9jxPEYtbscuCPSXwDGKbov0Emkn6evSBVLVw8XKgI3NPtLh0s+l/5Yc\n/Cbt/8zvkZzn10l++UiNVdWZ+EVjs9RsgM/UPnM5QHqKdIt0spSGvOPpeWVaLldibpXOkPxA\n9YuHy/IwySGvLPOON2Ip399PKUn/I62WjpBsiI+T0rBYWeb9/tM4yrAMlYjXS365ajbAeWWX\nd3wxRmXIO2moAYEfKI/vzeRzvdZnpXWZfVVbddr9QPLDKQ0naMX/if0f+g3Sb6RsuEIbb0l2\nVJ3JN5SPs6WsAfbLhcs1NUxajb0mv5g45B1vnFW+v355+pU0nknaiVp/TrKdV5Z5xzPRlmr1\nXqXmFZkUfUbr5yTbeWWZ9/vPRDv01S8pBXYQnL9mA5xXdosdz2M09Iz3moBsNUivcXBdcQT2\n161+mrndjVq/Rzo0s69qqxuU4D2l2zIJf1qyHWnpl4xsnn3azyR7hQ5VZnKK0r+79EFnJBP8\nUuIXEFdXpiGb57zj6TVlW7oq1WW5RjpVeq3kF5BLJIe8ssw73oilfH8vU5JOkp4gHSC9QPqJ\n5JBXlnm//0Ys5fj7bSXDTQyuuWkOeWW32PE8Rs33qsz2WGVSSkKXCYG9RLcJZoMN8F7ZHRVf\n98uEDdP7k3zsp+XdyXq6SPNcZSZPVmbeK/nBvDXNWLJslWd7UdtJO0t5x5NoSrdwW6ibGezJ\nHyu5KvqXkqvh88oy77iiKG14lVJmw7RZulbyC8dHJIe8smx1PP39N2Ioz98vKin3S2FTkvLK\nLu94KwbZ/w9Nt6vOJga4OmU1oaS6vJof1rPaZ09xFILbx/z2/Hnps5KD/3O2y3NVmTjd/ym9\nQ7pRag7t8uzzXNZ5x5vjK8v2KiXkDyWX83GSjfFy6e+lvLLMO64oShmc7u9LLudjJFe5P116\nn+SQV5btjlfp/3xe2eUdb8fA/KrEweldEJxxQjUIPKJk+u2yudeztzdKVQ/PUwbOkz4kvSeT\nmdu03i7PVWXyl8rTQdI66QPSLpLDP0huG2yXZ7/1+zeQd1ynlDLYA7xY2pSkzuXnMrdBzivL\nvONJlKVbPFspeqr0JOn2JHV+aT5Lcg1IXlm2O75R11Yl5JVd3nF71K2eAen/h6pw2Cad9qgI\n1SFwnZLq/8xpcJWeq+82pDsqunyJ0u22wDdLWePr7FwvHe6VTDCDm5LtKjL5rdL+FckPFcvV\nyg42xNtLzvMOkqve0mAGaZ7zjqfXlG3pslrVlKjV2k4Ncl5Z5h1viroUm87fg9IdmdS4/F2+\nLu+8ssz7/WeiLfVqXtktdjyPUakzTuJGh8DJyoqN7RMlvxW6a/43pSoHt1/bq/uIZCOTlV8Q\n3UnJb8gvlhyeL90jmYHDKDDZX/lwVdoezlASLtDyc5IZeL/bDk+S0pB3PD2vTMs1Sszjkqth\nHX5fekg63hsKeWWZd7wRS7n+2gA/IJ0uuSz9u3VnJXeqS8NiZZn3+0/jKNPyBCVmU1OC8sou\n7/hijJpuxSYEBkNgUtGeKT0quTrr+5Lb0aoc3q7E2/i0kr0Eh1dK90lbpGul9IGt1WAUmLQy\nwOuUtyslv2w43x+TsiHvePbcMq3/iRLj/Dhf90pvldKQV5Z5x9N4yrY8RgmyQXJ+p6VLpP2k\nNOSV5WK//zSOMi1bGeC8sss7nseoTPknLSNOYKXy56rLOoUxZXbvRTI8qkxcQ+CHU7uQd7zd\ndcPe7xfH8TaJyCvLvONtoh36bv9+Xe3cLixWlnm//3Zxlm1/XtnlHV+MUdnySnogAAEIQAAC\nEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAA\nAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhA\nAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEI\nQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAAB\nCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAA\nAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg\nAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIDAHIHttRbObXW/4usJ\nEIAABCAAAQh0QWBHnRtJ7+7imuypT0+uPyK7k3UIQKA/BMb6Ew2xQAACEIAABCDQDYGJbk7m\nXAhAoNIE7MmulX4m/bU0KX1S2iilYU+t/K30oHR5ujOz3F3rvnZv6Tzpu5LDsdI66V+9kYTX\naXmbdGG6gyUEIAABCECgDgSaq6C/pEw/LD0u3StNSTPSoZKDjeis9Kj0mORzXYWdVkE/X+u+\nZlq6S/K535Qc/lTy9ke8ofA+ydsv9wYBAhCAAAQgUCcCrQywjeKfJxDszW6Vvphs2+u9TnKn\nrXHpGilrgG10NyfHtAhOlRyfDbfDtyTHd1yy/LKWBAhAoA0B2oDbgGE3BEaUgA3k55O83aml\nja4N8S7SSukcyUbX5/27lAYb812lW6QPSP8s7Sf53JMkh5dI90vfkGyoT5AIEIBAGwK0AbcB\nw24IjCgBVx9ngw2tX8QPl+z5XiKl4Xyt/FOycbCWPu523tXJPi+2SPck247L7csvkK5M9rGA\nAATaEMADbgOG3RCoGYEfKb/2Zl+dyfdrMutXJccv0tIdsFKdofVPSw4vldxOfKn0omRbCwIE\nIAABCECgXgRatQE/0oTgbm3bqDq4/ddV0kdLh0l3SNk24Ou17c5Z7nBlb/izktuAj5Rcfe24\nfyU5uP3Y295PgAAEIAABCNSKQLcGeCfRuVWyUbV+IWUN8Fptb5J8zNXN7iX9YcnBVc/Tkocx\nOdhD9rb3EyAAAQhAAAIQ6ICAjWdqSFud7g5bh7Q6wD4IQAACEIAABCAAAQhAAAIQgAAEIAAB\nCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAA\nAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg\nAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAE\nIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAA\nBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhCA\nAAQgAAEIQAACEIAABCAAAQiUh8D/Az5CKMlmo7+3AAAAAElFTkSuQmCC", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "gendata = GenerateData(1000)\n", "par(cex.axis=0.7, cex.lab=0.7, cex.main=0.7, cex.sub=0.7)\n", "plot(gendata$x, type='l', col='blue', ylim=range(gendata$x, gendata$y), ylab='')\n", "par(new=T); plot(gendata$y, type='l', col='green', axes=F, ylab=''); par(new=F)\n", "gendata[names(gendata)=='cointegrated']" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "LinearRegression <- function(x,y) {\n", " slope <- cor(x, y) * (sd(y) / sd(x))\n", " intercept <- mean(y) - (slope * mean(x))\n", " std_eta = sd( y - intercept - slope * x)\n", " return(list(\"slope\"=slope, \"intercept\"=intercept, \"std_eta\"=std_eta))\n", "}" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "