Skip to content

Instantly share code, notes, and snippets.

@ljmartin
Created April 27, 2020 05:12
Show Gist options
  • Select an option

  • Save ljmartin/c19617f4626fbfd4f64976aeedf04dca to your computer and use it in GitHub Desktop.

Select an option

Save ljmartin/c19617f4626fbfd4f64976aeedf04dca to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import pymc3 as pm\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from seaborn import kdeplot\n"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"-2.7112278332640463\n"
]
},
{
"data": {
"text/plain": [
"<matplotlib.lines.Line2D at 0x12b908150>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO2dd5zc1LXHf0czs9W79tped+O1jXHFjTVgGwjFNhgSQyC0kABphIRUEhLzCASSQEiBJISEUENIoYTy4GGqscEYDC649467171sm3LfH5qrudJIGs2MpJnV3u/ns5+d0Wikq9HVueeeewoxxiCRSCSSYKIUugESiUQi8Q4p5CUSiSTASCEvkUgkAUYKeYlEIgkwUshLJBJJgAkXugEiXbt2ZXV1dYVuhkQikbQpFi1atI8xVmv2WVEJ+bq6OixcuLDQzZBIJJI2BRFttfpMmmskEokkwEghL5FIJAFGCnmJRCIJMFLISyQSSYCRQl4ikUgCjBTyEolEEmCkkJdIJJIAI4W8RCKRuMiKHYexZNuhQjdDo6iCoSQSiaQtM3f9Pnzp8Y8BAFvuvajArVGRmrxEIpG4xE3/+aTQTUhDCnmJRCJxicNN0UI3IQ0p5CUSiSTASCEvkUgkAUYKeYlEIgkwUshLJBKJS4QU0l7vPNRUwJakkEJeIpFIAowU8hKJROIBr6/YXegmAJBCXiKRSFyDhNfF4k4phbxEIpG4RCSkoHNlCQDgpO4dCtwaFVeEPBE9QUR7iWiFsO1OItpBREuSfxe6cS6JJAh8sGEf4glW6GZIXCaeYJgwsAsAYP2eYwVujYpbmvyTAC4w2f4Hxtjo5N9rLp1LImnTvL++Adc89jEeendDoZsicZlYIoHKEjUl2J/eWV/g1qi4IuQZY3MAHHDjWBJJ0NlzpAUAsKnheIFbInGTRIIhwYCK0lChm6LDa5v8d4hoWdKcU+PxuSSSNgFfnJPGmmBxtCUGAKgoaT9C/iEAAwGMBrALwH1mOxHRDUS0kIgWNjQ0eNgciaQ4oKSUZ0yK+SDx85fVJcnVu44WuCV6PBPyjLE9jLE4YywB4FEAp1rs9whjrJ4xVl9bW+tVcySSooEo8z6Stgf3i2+NJQrcEj2eCXki6im8/TyAFVb7SiTtEanHB4uWpHDvUFpctZhcaQ0RPQ3gbABdiWg7gJ8DOJuIRkPty1sAfNONc0kkbZ1DjWqQjPSgDCYdygIo5BljV5tsftyNY0skQeOu/1sFQNrkg0r36lJMHtYdb6/aA8YYqMD2ORnxKpFIJC4wrGc1AOCC4T3x9qo9AICNDYUPiCqueYVE0o6QenywGNKjCkeaozi5T8dCN0WH1OQlkkIhpXygiCUYworeNFNoUw0ghbxEUjBmLN+FRVsPFroZEpeIJxjCIVWkXjqmNwAgFi/8SC6FvERSQC576MNCN0HiErFEQtPkp56sepBH44X3mZdCXiKRSFwgnmBa+b9ISP3fUgSBUVLISyQ+wwWAJFiINvmSsCpajzQXvnCIFPISic+EFCnkg4ioyfNc8r98dVUhmwRACnmJxHeKYTFO4j6xOENYUUXq6QPUwiGVJWEs234ID7+3ETsONRWkXdJPXiLxEcYYYjKfQSARNfluVaUAgOU7DmPagx8AAF78ZAfe/OFZvrdLavISiY+YCfjN+2TxkCAQSyQQTq63hEzWXQ40tvrdJABSyEskvmJmqnllyc4CtETiNjrvGiVdtBaqpq8U8hKJj7Sa+E0nZKKyQCB615gtrhdqvV0KeYnER2ImQl5mowwGsXhKkzemN1ApjJSXQl4i8REzm7xchw0GasSrKlIVqclLJO0THuZ+y/mDtW3SbT4YqLlrrG+mUqBkZVLISyQ+whdee3Ys07YVQ6ZCSf7EhIVXM6QmL5G0A7gmHwmlHr2ySKhQzZG4SNwk1bBIoQZzKeQlEh+JJjX5SIjwty+NBQB07VBSyCZJXGD+5gPYdbjZNqr1SFNh8ti4IuSJ6Aki2ktEK4RtnYnobSJan/xf48a5JJK2TCyhavJhRcHwXsVVQUiSO9c9MR8A8MGG/Zb7HG2J+dUcHW5p8k8CuMCwbTqAdxhjgwC8k3xfVDS1xjFvo/VNkUjchmvy4RBpHhjSg7Lt0xSNF7oJlrgi5BljcwAcMGy+GMA/kq//AeASN87lJre+uAxXP/oRth1oLHRTJO2EmGCT5+ZbGQwl8RIvbfLdGWO7ACD5v5uH58qJNbuPAgCOFWgaJWl/cD/5sEKaS50U8e2H5dsP+x78VvCFVyK6gYgWEtHChoYG3847f/MBTchLJH7B0xpEwooW/yg1+eBw/xWjbD//3INz0f/W13xqjYqXQn4PEfUEgOT/vWY7McYeYYzVM8bqa2trPWyOnisenqe9lm7KEr/gfvIRRdFc6mTEa3A4f3gP7fUL35qAi0f3KmBrVLwU8q8AuC75+joAL3t4rqxoLuJFEkmw4X2vLJKyycvcNcFBDIY6pV8NvnBKH9P9/LznbrlQPg1gHoDBRLSdiL4G4F4Ak4loPYDJyfdFgbG4LhUocZAkuGxsOIa66TMwe41+ApsS8iFNk3fzef/3x1vx5Aeb3TugxBGj+nYCUXpgmxj0JuJn2mFXKkMxxq62+Og8N47vNq1GIS9lvMRlFn96CADwf8t24pwhKZ+D5mTfKxU0eTfd7257SQ1VuX5if9eOKclMaUjBaf07p223KtoeSzCEfQp0LvjCayEw5vSWs2WJ22jukQaNrblV0OSTM8h7X1/ja9sk7tMST6DERGqHTYqHAP5q8u1TyBs0+QfeWY9zfv9uYRojCSQ8R42xEBQ315RHQiiNtMvHL5AcaYqixMQ0Y5WV0s86v+2ykLdRyM9YvqtALZEElZ++sBxA+gJbcyyOkEKIhBREQkBtVSkmDuxSiCZK8mTPkWZMe3Au/nrNKdi877hprd5isMm3S1XCKOQlEjd5ZM5G7bUxh3hzNIGycOqxqy4L+6rVSbLnww37cPC4WoT7uBA4+fryXdhzpAX3v73W8rtWSSl5DiM/aJ+afNxfF8qXFm/HD59dite/fyaWbT+EK8ed4Ov5Jf5yz2spG7vxIW+OxnUeGGFFKViBZ0lmWmJxfPGxjzGqT0cQEZZsO4Tnvjkep/bvrOUeysU7r81517Q1jC6UHMaYJzmfH5y1AQAw9U/vAwAuHdvHchonCRZGTb41ltDd+5BCnmjyhxuj6FgRcf247Q0+61+6/bC2bc66Bpzav7MmK3IRGTHjYo2HtEtJY2WuicYZ5m8+gD1Hml09n3GFvbFVBmO1F4xKQ5zpS8SFFPJEq7v+yfmuH7M9svtwuiw40qzmhed30U4xLC8x16OlTd5jrIT8pn3HcMXD8/DZP8919XxRg/1NRtwGgz1HmjMmtzOaa4zVg9zU5A81tmqvuZ++JD/eWrUnbRu/X+QgYrl3p3LT7UY3bi9pn0Le4gf+v6U7AQANR1tcPd+Og/pqMXKhLRicds87mPagXiEwamhGc42xDmhYIcRdWoRrjkqHArcZ2Se9sAuX6fxe8/e/uexk02OcOagragymsyl/mONeIzPQPoW8hSZ/vMUbDdtY99EYICNpu2xq0LvNGYX8swu36T+PM535zk1zDZNJi10naqoQqr/zHS+vTL5T35eEzcXpP792GhbfMcWT9jlBCnmB4x7llTc+etKbIrhkErRpmnzIPSEfjcl+5TbmssKwzpK8f8ZZW7HQPoW8hblmaM9qT85nzBcel3kUAkumWxtPJHQLrwq5Z5OfvdY0m7ckD8w88Ywzc17SUQr5IsJKk//Fq6s8Od/UET1175duk4tiQcUo5Mee0En33twm746Qf3D2Bu319RPqXDlme8dsfS6UJuRVeZJJyE8a2g03fmYghveqxri6GvcamYH2KeR9XNkG0jvFLc8v8/X8kvyZvWYvnl3wacb9jLO2zpWluvfxBEOIRO8axTWfaXGtZ/GnB1E3fQa2H5T1i/PhVzNWp20z3mOuNFpFt3Ieu24cpk8dgqqysK/pzdunkPc5rYFx8YZ3hqfnf4pVO4/42hZJbnzlyQX46QvLMxZ7MH5qVO681ORF4cODdxZtPejKsSUpjOY1vuDqNJAyrCi+pjWQQt6CQ42t2HesBWN/+TZW7jyccX87jJ2CRzze+uJyXPjA+3kdW+Iv05OJx6wwanlG7S6RMARDhci1B97MXGCcRUrypyXpqnrZWLXq094jqkmnJOzstw4plJad1EukkLfgo037MWddAw4cb8WjczblfC7GGGYs02e5/NzIwtd9lOTGswu3oW76DMsgKC7jv3x6P5RFFCzYotekVU1ecKEkcq3G6zfOGpC2rVgXA9syL3yyHUCqIMjuZIT8ih3OZuVuxkY4oX0KeQc2+dJwCOGkxu3U+6E1lkhzw/zk0/TpcucOJdJXvo2zzyJgjptz+netRHM0gQPHW3WfGyNew4p7mnxFiZr4TNTepZD3hl2Hm7DCMMOfNsqZ8hZSSOau8RonmnxpRMHPX1ZLqTm1mV75yDwM//mbum0Hj0fT9ovGEr4v/krcxUwwb9l3XHO5s7KSROMJnRAOKYS4Sw8876diyTkZk5Efg7p1MN0+/tez0jT3mooSR8d0MzbCCe1WyBt9XY2UhhUcbFQFtFNN3ixfyM3PLUk/fzwhUxu0cYwl+461xHD279/VPKcUi/4VSzBdBaFwyD0/eX4Y8fh3z/DGLbi90LvGPPeMGVZVoIyEfE4v7bmQJ6ItRLSciJYQ0UKvz+eElngCnSv1o24fw83ceSiVfS6fG3KkOWW+ef8n56BbVSlaYwnXtDeJu6zdfRQX/+WDjInHZq7WBx7xpHNz1jUAUGMiv35Gf1SW6Ot+RuMJnaadbVqDeIJZJrjjpiIxvH6nSRZFiXOyuTdOhXzYo/TSVvilyZ/DGBvNGKv36Xy2tMZUIT+ga6W2TSzkAADffXqx9to8f0X29O1cgYqSEBpb42iSmSiLkt++sQZLtx3CvI37s/qe0ca6+0gzFIXSopu37m/UDSChLCNeb/zXIgy5/Q0AwG/eWKMNKkBKIJnVGpVkz4cb9uH99fsc7x+xKNptxKv00la0y97QGkugNKzgT1eN0baVWiQXAqyLjFghPnhGKkvDONYSw/l/9C8LncQ5XBuLCQO7k0Vy4zrPI3M2QTF4zvBjirOAbKfubwupbx96dyOufSKVN54fRlRY+nZ2bm6Q6Mk2aNHKRGfEzcV2J/gh5BmAt4hoERHdYPyQiG4gooVEtLChwVo4uklrLIGSsILBPapQXaYm9Q/baD+NrdklLhMfPCORkIJZa/bicFP6gqyk8PB+EBUE75Mfbsn4PWNJyViCQSF9rnGusX/p9FT5x/ISBU3ReMYgKyOi186WZAFpPohUlKaEvDGlhtcs2nowMPUSjmf53DsliJr8RMbYWABTAdxERGeJHzLGHmGM1TPG6mtraz1vDGMM8zbtRzTOUBJW8NTXTuMfWH7Hzex+MjiluFm3+ygAvXB+OVlnwIr9x1rScrkzprovRuMMddNnYM66Bk3I9+ucMhN2Ki9BPMFwPMtqYc8JKYzP/v27eH7Rds2sWC5o8k0+ViHbdbgJlz30IYbc/ob06rEhcDZ5xtjO5P+9AF4CcKrX57SDuz0tSSYJ4yLX7jd3sqASc2i3lzK+uFm/91jatkwJ5e57e12aS+wPJg3STd//9dFWrN2t9j1xoC+NqI9gS5bar9G758f/XYq/vacG7YkLr68u25n1LCFXRGVo2fa2n4TPq59NcdFt1tH5vDw4EVUSURV/DWAKgBVenjMTxsT+PFjELg+4kxTEj76/WffemPb1/OHddeeTBIeIQmk2+YkndtUN6G+t2oPLHpoHQK808EVSN+Im+DHEhdeDjVG8tnx33sd2Qki4rvvfXoeNDekD5s5DTTjc2DZMlWKKil9cPBxXjeuLm84ZmPdxg6bJdwcwl4iWApgPYAZj7A2Pz2kLn0Z2LFfLcXGZy9dBzBZgnQRPvbdOL9Tvf2ud7v2wnmoZMWmuaRtEsvBQKQkraYvzCpEu26TI/M0H0s7jZtI8Y9u3Hjhusae7iAvU76/fh68+uSBtnwn3zsJ597/rS3vyRVzAvnZ8He69bKTOFJYrIUXxtaaEp0KeMbaJMTYq+TecMXa3l+dzAhfy910+CoAg5BnDml9egJe/MzHtO69ksMkC6Q/W8h36kGc+U5CafHHDb49CwNm/m4266TMyfoexdCGtkLW3hZj6gs8s3agrXFtVilP7d06brfpVMcoqBa+RfcdaTbcXG1eN65u2LeqCmSWskGPzrhu0OxdK7rrEp5ZiXueySAil4dxG6kyan50bnsxjUzx0ryoDoHrCbNmv5mK/sj79YRd5bO7mNIEWUigtzTDn5smDtdc88voLf5uXa5M1WqJxlIaVNCHvl7ueMYBslxCI9dzCbbj52fTo72KGK2RlkdTv+bUz++v2mTqiR/bHVVTXWr/WStqdkOeaPH+4+iT9iL8ysU63PVsymWH49GzuhvTgir8IFX0kxYHoHWIsxm1GS0y/cGpnrikXomCd+lY7oYkLeYPC4Yb26YQ/vL3e8rOfPL8MLy7e4Us73IIL4dk/PlvbVl0WwcybUw6CNZXO8tWI8H7hl27X7oQ87/BcKFeXRbDl3otw5bgTdNtFenfKLaBk+faUyaZHdZnlfnPW+xMfIHFOtlkC//SOXsCFFLI0zRmLhrgFdws2avL5prWNJ1haNk0zZq7ek3GftkSCqfeqZ0f9819eEtZe53L3+Bjsl5tpuxPyfOpqZV4xE/JO0hqYLdg2HEtNV685rZ/ld405xyWFJ5so5ynDumPrfn2ZPYXIUksXN2e7EN+1Q6nt56XhUNrgkq8nx2Pvb8LYX76NPUes8+AE0S8+wZipy7OYe4gIuHZ8P3z33BMdH5fXE/DrNwtn3iU4HDzeiv3JRR8rDSpXIX9C54q0beLD5ua0XOI9//OSfQUoEatn1eqWi/0i7DDfCae6PIx9x6wXac2UjXyFyavJojdz1+9Dz05lGHtCTVquJ7NztHUfgwQzL+knmsMIhF9cPCKr42qavE82+XYl5Mf88m3ttaUmb3JTnUzdjzSn+/5KT5r2gZk5hMhauIrPdraafKa9zeRGvjb5tXvUKOAf/XcpAODMQV3xTx4pnuTjzekJ3cpydGIoFuKJhIUmLwj5HB7x3YfVQfpocxQdSr0Xwe3OXMOxFPLCVOyNH5yJy0/pg2gGm2ZzNI5/ffSpbtvA2kop5NsguWhXVuYQq5QCoquh0/S0nEzFoo+3xtIWgfN11zN6Dr2/fl/atd32UnqMY1M0jqMmyk9bIRpnpnJCJ+RzOO4TH6iBk6/7FKTWboW81cMlavKDu1ehtqo0oyZvLPkHqMFWfmaak7hDLu6sZuY8AvDSEnNvEvEMbisCzdEE3lypXwD1wijw7X8vcrTf9Bedm72KjVgiYSHkRZt89vevawfVI6c6GZDpNe1GyG8w5CSxyv0sTp+JCOGQgliC2fq0mk3LYwmGW9twB2+v5PLgWSkBjS16bbcqOTUXNXm3zTUzV+/B2BNqsjpmLsxeq/cIs5J12ebkKSZicWa6dpeLYBf57RdGAgB6dbL2uHOTdiPkJ93/nu69lSZv9DHmnVSs8GTEzBMjGmdaMMi4Ou8fOok75DL7slqYH1BbqXt/2Sl9AKhuu5xsF0WdyBc/1/gXbT2IIbe/7sjFsq3xzIJt2JshEjkXed+xXNXk3UxlYUe7WngVsRLyRi8Ybj97Yu5m/HDySabfEZNLfe+8QVi67RB2HmrStt1wVv5JjST+YJUdcGBtJTY2mOeAWbr9sOn22y4aijteXonHrq1HazyBLpUl+NoZ/VFblXKDzNZeTjlYgcVBxW0efm8jmqOJtFTLnJmr9+LNlf7YngtBLveDK5LZxmLkSrvR5I04vTncZmpXOEQckW+efBKqyyO6xTjpPdl2sFpEPdwUw5Z7L8Ifrxzt6DgdyyMY3qsjXvjWBNRUlqB7dRnCIQV9Da622fqwZ9Ic59xyDvonZxB/v36co+94zY+TXjmcSJaLzYWid6dyXDq2t+lnP7toKIDcfltuKfYrE2W7FfJiPgo7bjpHDXII2fgzp2cgBDbvO66lKB7eq2OOrZT4jZX5hPumXzKmt86OPqRHlen+3WwinEXcdKFbf/dUnNClAj+aPBhPXF+Pc4Z0Q1VZuOCBSsblrLYSN9WSLBNqBr+mXIarsM/BUO1GyBtDvascTmE/O1Itn/a39zZa7mO0rb28RM1aeVL3DggphB4dUw+8k9z0ksLhRLsSH85PDzTa7JmZUX07YURvtU84SVhlJRieueF0zROkJKzg3CFq/YKQQmnZIf3GeP5CDzpOaInFse9YC/YeMbfJ86yyuWjyXEnwy/uu3Qh5q2RRmXASkWhlV315yc608z5/43jc+blhObVF4j3ZCqBGE1/4iSd2yeoYU4b1cHxuKz/+0weYn1MtJp6bUG2NJbBgywHTz/p3VU1Cb63KnK/G7PxLMlTbKjQfJBMJvrNmr+nng3uoA/PIPp2yPjb32PFr8G03C6/ZJIL601Wj0ZJcSKoqy/wTGbW/7557Iv48S80saRwjKkvD6O5wKi/xHycpLM4d0g2zLB7+FXedbznFt4I7AcQSDJmCRLMdhBQi5BoL9bV/LMD769OzpgLZeYaYLTAes/FWKwYyrdl95qRavHfL2ejXpdJ2PzM0TV4uvLqL6DVz/YQ6230vHt0bVyQLBtRUlqAkpOCikfqq9x9t2o9H56g1NY0P3kndU3Zas5nAiN7SRu82K3YcdqUQgxMh+tdrxppuv+/yUehQGs6qqhSQitnIZCr6dH+jLhGak8lpSMm9XoGVgC+LpFfCGpDU7FfcdX7a/n6WunMLJ/0gFwEPpAZ1aZN3GXGx7M5pw7P6bv+ulWkC5KpHPsLdr60GkN6JxVmD2QSib+cKbLn3Iq3uqyQ/1u85is/+eS5+9+bavI7DGHMkkMoiIbzynYl47pvjddu5H3y28Ic+mkE7Put3s3XvnbhGhohcT4SlEKHVmD9fISjkfCHZC4+f5xdtx0uLt7tyLL6GV9/P/RiXlE1eCnlXKctyCq37ronmwpm9Zi9mrdHbJcVZg11EozGTnyQ3DiYLQy/aml/KZv7Mmd2yE7t10L0f2acTTu3fWXtfWZL7veR1Q5scRodWlYYx5oROeOL6+oz7Kgq5XnmssTWeVng8kWBaX//hJPN4Eq/58X+X4ofPLs28owP4tdxy/uAMe2ZP4LxriOgCIlpLRBuIaLrX57NC9HDJltJwCM3CAyh6QXzlyQV4er6+clDYoZC/83OpGYVfpcCCCPe7dmJPt4N7O5gNvndZzP56JvvVcYtkZE7glaKcCvmjLTG89O2JOKVfZ8y8+Sy89r0zLffN1bsm08BgtMlHEwkt3P/7kwahS4aKScXuYcMHMaNXnhsESpMnohCAvwCYCmAYgKuJqCCuJfyH7VZlX3TBjNKIgo82HdBCtzOlbg3pzDXWQr6msgQ/SkbRFnunz8S0B+fijzPXFeTc3AbemudCFr8H5SZC3mox9fbP5t+dNU0+h4HixG5VGNbL2i1XIUIuP4tRUzeSYHqvsm0HmrJajPVCeLrJkSZ1dpjt+ooTuBKYb8Uup5CXGiQRjQdwJ2Ps/OT7WwGAMfZrs/079xvKJv/PE560ZeXOwzjWEseYvp2y7mArdx7BsZYYqkrDGNarGgnGTKs5jezTEeWREA43RbFmt5qDuyREGGOTMGrnoSZsO9iEcXU1bTo18cebVVe70wQThl8cb4lhxc4jeZ8/nmBYuPUgSkJKmpAb3qva1N584Hgr1u89hk4VEQzubh4YlQneX4b1rLb15uK/ce9O5ehT46wk5dLth1BREsYgg7kpE/y3sGNcXU3ac8B//zW7j+Jwk3Wa4cE9qtDJ5SyMbvZBfiz+TLsJlx81FREM6tYh74RnAPDcjRMWMcZM7XdeD6e9AYi2jO3JbRpEdAMRLSSihdGod7mnGVNDzXPRIPgtaE4uNlmNi9wnXnfPMtzANizXiwa31BR+nCyLNQHILfKRE8rSbzobxYyy3J/jpC12k89B3e0HlbVJJcgtvPI59+Lx5EL9YGPUMh+Sm3jtJ2/2G+nuBmPsEQCPAEB9fT171uCx4AaLth7EZQ99iMnDuuLRazMvVhm5+pGPMG/TftRUlODZb47HvmMtqP/VzLT9Hv7yKehWXYb5mw/giofnAVBttnbX9Pjczfjlq6vw2HXj0NGn/NJeUDd9BgDYXqtXLNhyAJf/bR4Uyu/8DUdbMO7umRjQtQNW7TqC4b2qkWDA6l1HcPfnT8bovumBLzOW7cJN//kE4wd2wV+vOSWn867ceRgXPTAXP5h0Ei4Y0cNyP/4bf25UL9x64VBHx66bPgNN0QT+/MUxWLD5YJorsBW7Djdh/K9npW0/uXdHnDe0G/44cz0eumYsTr3nHd3n4u/P22uFm33laHMUJ9/5lmvH5W3/65dOQe9OzmZNTmGMof+trwEADja2utLe5260/sxrTX47gL7C+z4Adnp8zjQue+hDAOlphJ1izFhpFcQQTh5fPE2mSNuUfa54bPIb9h41LYRSrPD7kW1udiP8HvD7fXLvjpqnjdV95BpkPlPuihJV12qKOvvNc+krVz38EW76zyemZSrNsLKvP3LtKeiVFHrZFDv3Gjefn3kbU6UMc42Ut0PsK3489l4L+QUABhFRfyIqAXAVgFc8PqcluWa/MwoPKy8Ovp94EzMV8PY7j4UTJt0/B9c+Mb/QzXAMf8DzfWD4PdAXjrH/Dhfy+aynpBZerfuA6N2Vi1fGpn2qWYA57GavrzBPDxxSSFuEzrQ46ydueqocakzlxg+COdVTIc8YiwH4DoA3AawG8BxjbKWX57Rjj0WyoUwYH2CrDsW1cnHvTLbCYtTkgex8zgvt/slr8OZrl+X3QKH0wZpZWP7HJ3PGfGViXc7n5S6UdumsxcpmefUVh0Lr3tfX6N6L3kV8RtxikUPeKVv2uWePFn+TYXe8gbrpM/D++gabb1gjPu5dO2TvjZcNfqQh99yPiTH2GpYekYsAACAASURBVGPsJMbYQMbY3V6fz455m9IryjvBeCOy0eQzPQirdqleIe+uza1Duk0uAvtYgU07vNBHvmMNH7z53QspqRwmVsfuVl2GLfdelFfJPa7JN9v4yYsDWD4RrLkOyL+6ZATCCqG6LOW8cOED79t+J9PAd/bv382pLWaIihdPGvflx+fnOCA6i3NxA2N9AS8obmfVosGZuYb71IrBUHYPLgDMT7pqvbPaPOGV3+TyTNz3VmH84zniA75DqMiVLVwg8DE6RISfXDAYXSpL0iJe3aQkrCCskGlGS44om62qVzkh11nAZWP7YMM9F6IsEnLsofbzzw33rfSl1W+Sb4Cc12zd34iVO80ri7mFFPIOMA7mVguvfD/RvJNpccoNH1k3ETXGfcdaMg5SAPDkh1s8bFFmxPWMiffOykpb3dhwDIeTaRH4feX3RFEIZw6qxaLbJ6PSxeIeZsQSzNZE9tzCbbp9cyXXr4prS9k4MHitCXOsZje5TFz8eCSvHd9Pe+11rVcp5B3Abzp/QKwWSblwEDv2cRs7KwBUJO2x2aRC9hJRyNf/aia+8dRCx9/164E2YtROX1vuvKboefe9h88++L7uOJq5xucBmAfgmPHvjz/VXueT2M4Nf/JIFrEmTuox5MqB46147P1NYIxZRo/mYtryY4npFxeP0F57EVUr0q6EfK6LKFwz5x0pU1oDMU9Opg5zW7JW5JQiyUhpbK9VulkzCiHiDx5v1ULQOYeaWi32NmfbAdXEwwdvzVxTJAOvkSnDrX3pM+HGAn82mryb4+Thxii+/8xi7DnSDAC45b9L8asZq7Fk2yHL2U0u1+t3JS2v+1nghbx4k62STGWCC3k+nc+U7L9jeQSzf3y2o2PzRE7Fkp8snw5eCMvTmF++jdtf1jtsNbbkliwspcmnL6AHBTcEmJW7qJmwclNgzly9By8v2Yk/vbMeADSf/2icWT6TuWTg7FShBiX+6SpnRdvzxetuFnghLy685Dor4jeB+wVzl70Km/SyTs0voSJzoczPO684hGKtIQnd4aYodh5qwh9nrkPd9Bnag//E3M26/bg2eMagrgDU6j9BI5dwjMnD9LNMKwuMWTpvJ/EDJ//8TVzwxzkZ9+MupnuTmnzK84lZenjlYq7hv1HPju5Gulrh9bMfyPJ/rbEE5m8+gDMGddVqNQLA3qP5+clzIcC1hrumDcfzi7ab2lIzBUFxuD0uWiTBUPloXsUSHGP0/pj6xznYebhZG3h3H2lGr07l+MWrq3T78YdtXF1nbLh7qhbB3JYZUFuJTUJ+lGzv79+vH4eJJ3bVbbMS3M9/a0LaNiezoaMtMS2hnx3/XaQWBJnJPdGShz7eGsNXnzRfO8pFk+e/kV+33+tHv+33YhN++8YafOnxj7H404NYtycVRNKlMjebPO+ncU3Iq3dlRO+OePab402z3jn18Cg2Td5pRCTH+BD5GRhl5ZVg3L7zsKr5ceE04d5ZaImlm3T4IB5SKBACHgD+96aJuvfTX1yWleA7tX/ntEGz1mJta2jP9JTHH+UYm2LGsu0pV8OjzVHN/fioTb3YXDyR3EhVkQ1eR7sHoycb2NCgCvaDja06e9eFJ+e+YCXCOw5Pk/DE9ePSaltyWde92n5g4fU9My3m+oVxeptpQmJ8iNbtOYZXlvqTnujvH2w23W4mwAG97dPMtZUvrBeLp5MbGEsEfrTpQMZYAnEQMHMdrakswW+/MNLR+Y0DrlvJvl5dtsvRfrkoT/wR8Cv1t9cLvYEU8vw3IxAG91BzfN81bXjOI3NNhbo4yrNEpnKcqD9fZWk4Ldc4X7y5sr4v7AglB4pcQ7DdxtjhMq38Gx+iLz3+Mb739GJLQesmvOwf571bzgZgHZsgPrRmMxa3Ep3lwsg+anF3My3b7dlRJrdYJ6bDTM4HVoQUwpdPT/mIbz/YaLO3NeJPYickc6qKpeUjyvqrOeG1lTOQQl6DUm59I3p3zPkwP71gCADVVgukOridxldVFsHqX1yAH2Sod8mP8e7aBry3rvCCPlshb5xqNiTXPf4j+HV7wbGWGJZtP6Tbxgdhq1QS4hhvtiBnzELpJ1OSi5tm7RJnS24Inkz2bycCXBzELz+lD7599kDT/TobygAayxGe8ZvZGc8FpA90YuCYXXtz0eSNOYy8RpprckC8rev2qB06nyl4eUkII/t01DpnzKEwKC8JZVyAFdvFvQYKiVHGZApmsXqIHn5vk1tNMuXGfy7Chxv19l6e6MvJAnDMZJ9oIvPg7RV8VmgmsESTxwfTz/W8LU6EfLMwkF44sid+klSEjBiVBKLcPLiM/eyFT7Zrr40mw1F9O+GqceoMOjdNXv3vm7lGavLZw0d9AnDPa2o2vXy1M4Uo5V2jCYP8fz7xIZi9tvD5azJp8it3Hsb6PSlN0GphqyzibdeauyE9SKskpIAIaLFIxSBeW9Sk3XGDGc5P+MBiptWJbsA9qnMvSO8Ubq755cXWcSVi0XG74KiRhhl0iAgvL9mRdZvsFlCNnykEnDlIdX/NxRTC5Ydf3SCfhHNOCKSQ5+jyuuc5KocV0uylXAt0Q+MT2/ja8t2YtWYPALWjzVy1J6MnxP1vq77fbnnnGA9jrNN50QNzMfkPc7R2WZ3XS++UT/eb23GJCIwBM5abL8qJ2qdZQisnZjiv4EqI2e8pavJ+eHxsTf6+TTZ5i04foJouzx/eHRMGdrHc755LT0aP6jI8+MUxAFSlwW5/K+w0cuOsjJByf8wt4lX9758mL4V8VjDGPFmtDimkaVRe2m73JnPev7R4B77+1EL86+Ottvv/7b2NANxLcuS0w/F6t/w3qe+nzzbopaDctO9Y2jZx4dtJ3UyzxcVC2uT572XmZcXNT049WvLl/rfXAgDeWrnHcp8JA7ti/d1T8fCX620Hnu7VZfjof87D1BE9cfWpJ+DBL47F8F7Zr4/ZafJGQa4QabOxtrDw6mbBEzMCJ+TH3f0OPtig2mrdvEclYUUTaFFN43P/5+MP9M6km9vuw/Z2ei4cnGSLdILTZ4IPKk3J9LjGtQfjgpubmAkVJ948Jwi5u800PNFP3m/4zMdOk8+1fOWphjiOMwd1tdhTf75MibOySawVUgi/vvTknFM2izMvY98yCsn6us6aJp+Pn7zXmjxPm+B1jEzghPy+Y6moVvEe5avcl4ZD2nRf86f2QOO7I5mHhc9AMwkcnoP83XXu2POdaj58MOLnn2+I+jUreu0WZr+IVZyBKPzFh8ksz3jcxbWWbLErA8mvzWkedyP//NqpmPvTc3D5KX0wpEeVdo/fX9+AGSb+5gu2qJ4r2WSb9BpRWB84rk9AJ85iH7u2Hj+ecpKQVLB4hTx375ZCPg/EXCqVpdZ5ZpxQGlE0gRH1wXYbd9jRuDAtCeV3fYCa5c9pkQX+YHFhbxyMvOy42RQZP3g8taYgXlvRafJcyNt41+SakrY0HEKfmgr87vJRqC6PaOf48uPzcdN/PrH8XlWZN1lPvjqxf9bfsetP4kA+pGcVwiFFu4c5mWuS3cRrIc/v+fq9mVM65EOghfz8zSkXu35dKvM6VmlY0YJs4gmGkEKeLoIlHAqcrh3UqWu+FXCi8QRG/eIt3PbSCst9eNY/QBU8R5uj+O5/FgPQ1wAFvBXyv5qxOm3btFG9TPfVeYEIbTRq/mJO8kIsvHKB8kAyw6JIa1y9hlw1eZEQUUbBd0V9HwDAj6cMzvt8ZnSsiOCik3tm9R2u9Izqk27PFwugVydjJUIuaPKeZ4dMnuCPM9e7Zm41PY9XByaiO4loBxEtSf5d6NW5rHhg1gbUdanAxaPNBUA2cHMNYwwPzt7gqhAz2khfX75LE06ZhDz/PN/kYFxbnL8lPdkaH3AOHEtNk1vjCTw1byt2J337jQLIiS00kWBoONqCg8db8dC7Gx1Hdl4yJv1+WmUEFafyI4QFP+P9S7DCavLPJ5Nvvbg43b2wNaZPo5EPIYUy9t2qsgg6lIbRv2t+ipEdThP4cbhNvjScfp9FDzCexoEf32jacYKW1sDjfiCaBfMtim6H15r8Hxhjo5N/r3l8rjSG9axGUzTuSuWV0rCClmjck5XwnxoCSb7170/weDINbiZtQhPyeXrX2F0XH0BEG21rLKGbPVx+Sh/dd5wMgg/MWo9xd8/EV/+xAL95Yw0+3nwAG/ame84YsUqQJTJv43589ckFOg2JCWFyd7y8QreQmWBMEySF0OQnnGjtVsh/f+NsKRdCCiFTrFMsnvDcw8jsJ7Yb5PlahZkSIgZGcQ4lU158+9/W5igrtCyUHqvy4tKPlxlcA22uOdwUxZ4jLZZ+1dlQElYQTSQ8cc+0O2YmuyB3FctbyNt0Mi4oE4aFS1GQjxd8n7t2KHUU4PFm0kVv8z7V5fGhdzdi0v3vpaUrSGurgwHkmsc+wqw1ezFHSBUhtnfN7qO6ByueYAXV5E+tS89kyonmaZMXaYrGsXTbIWw7YP1MRBPM84HOrF/bKQbZzpzzmWnHfXKhFB+RP89KN9O5hddC/jtEtIyIniAi07LtRHQDES0kooUNDe7mbuHZ9txY2FCIkEh4U8GpX2fraXEmbSLlX+2dJv/MArWItPjgxBN6L5deQnbBsEKmwUZGVu86AiAlvLhw50LfCicPMN/lvrfXadvsMn0y5s9aixWdKqxdTvlg5IZNnntBnfnb2Zb7xOPMcw8jMyFv6wvv4METKznlo4zxZnjdD8Q2PjXPPh4mH/K6k0Q0k4hWmPxdDOAhAAMBjAawC8B9ZsdgjD3CGKtnjNXX1uZXicdqumfMVpgLIUXtaF5o8h0rIthy70Wmn2XqZ/xhyVeTtxskeKbAmE7IMzz1UapjhhUFM753Bj669bykScD578TNJtxXPNO1GIXBpKHd8N3zBgGw9wG3GxzW7D6CjzfvL1hef+4/blaNit+bXP3ksyWa8N5c88mnB9O22f323CPox1NSCf/OG9JNt484UNo9p4s/PYhDjda2euaTJu9XJbW8fKQYY5Oc7EdEjwJ4NZ9zOYEXZPaCEKkLVlc+/JFn5zA9b4aeZixNmCt2Wi4vgyY+OIwxze4JqIuCJ3ZTFzatFvfiCYZ4gllqpHxWcsgwKB843orG1hj61KjBTE9+uAUAsOzOKYgoipaYDADOOLEr3l+/Dxee3AOvLd+tO45dtr9vPLVIF2NRKN5b14AjzVEcbY6hc0UJyktCmleXG+YaJ8TizPNzmc3W/vruBtxyvnmiM96fqoT8+OUlIYSVVE6piYLJ0Kp0H2MMn//rhxjWsxp//8o4dKksSUvBwc2SXrtQ9u2sb2PD0Za00pVu4KV3jegj9XkA1r55LtAcjeOs31lPQfOF276X7zicYU93ydTReIfkQj6RYLjvrbXYe1QfKTt3/T5bDdnOJv+7N9cm9xE0eYOmJD4oVkJ+yO2v46SfvY61hlS3fPDgA83dr+ldJMf/+h1dSlqeUriqNKwT8PzcQMojRcROUywGAc95dekuTLx3Fr7y5HwAgibvU3BSLJEoyOLzX2ZvtPyMC3Ix8V1ZJKQtRvfsWKbrg6f274zhvarTomN5H1i16whOu+edtBKQgJC7xuPfwGgOuvFfizw5j5e95rdEtJyIlgE4B8APPTwX3lntbQbHQlWDy9TRuLDlAnzxtoP486wN+MEzS7R9Pvn0IL70+Mf4/VtrLY/jZCYgavJGgSkKBSshz4X4Wyv1GjY/rlVqAmMRkNF9O6FXxzJTmykX8sZBDlAFhdfZMd1gaE81EvKjTar9PBWJ64/gjcZZ0ZU/jGtCPjWoR0KKNvCZKUPj6jqnKS/GpGsvmbis+p27htOQYw3qTHh2JxljX2aMncwYG8kYm8YYc1avK0fsIvcmDe1u+ZlTvB7VAWCiiRtdptPyzs+F/GUPzQMAzX8dSHWeTTaJu5zkEBdt4Uabp07Ik70vtqKQbv1kTzIpmyjM7UoIJhJMq6hlhLdDrAfKiSdYmp/11BH6kpA/ucCbAKBsMK45uFnE4jeXnZz5/PGEKz75bsJ/g3JByCuU8ps3WydWkllJRYz9ysw/3a+0BgBw8+TUGoNXg3hxDdce0acm/7qSXvvMAsADV41J25ZpHTBl6tB3VrG4sZOFJCfVaYzeNSJGc00mN0ezNLaiOel7Ty+2TF8QSzDL+2GXC161Neu/d/Zg/UJnpAB5a4yIv8Pv3lyj5QdyQ9FwYvON+eBCmS2akBfMc0QpE5bZepxCemVk24FGLNyiX/A1m8H6mWpYfKa8UiS9SU5RZIjZB3PFD99pM7e1TFGgXAM32tsTOq1b/W93DWY2bCN2Cb5E4amaa6wHDYUIzybdMu3Ye7QF/U0KSccZs7wWOytDLJFI+57xNy8GDVYM4PrL7I2oSv4GbnRBJwuq0XiiIEna7OBKiKjJE8g2QExRSKck2bmNiviV1kA8F+CdIllcdzIPeA4XM66fUJf38Y3C4Rc2VXNyPoeJgMnk0sc7iVEj2S+EczuZfmaryRuFvCgUjFGVt764HHXTZ+jaY1eQgrN1v7l5KR63E/LWXXrdnmOaaUhrt+E3L4bMi6t2HtG9P5qc0bihaDgR8rE4K0hOfcBaqeF9T2z/Z06qtV2MJnLmX29cC2IOlCK3EIW8V5p84Xu0S0wb1dt0e+fKEld+POMNdyPE3IjZFNlKyH/hoQ9x3n3vCjZ5686cCu6wPrdRaD96bT2++ZkBlm0xzhzsNPmn5+uLepdFQo6Kk/z8lZVp56+bPgNvrNxtKcyzNTMY76Nfbop2iAFcIm6YD5x46EQT/i68isnKjIMwh/c98TmcNKy71reNayuAqhk7yYf08SZ9qgS/XCgBvTnWqzGl8D3aJRgy38x8MN5wcZXfq3MA1hG2C7cexMaG45rG3BpPWHZoJ51W9JMPKYTJw7rj1qlD0btTOT6XzPAoakXGmYPo6ZIpCVYiwdDdQa3SrYZ0FOLAYiWDsh3QxxoqWvkVcJQLbggds0Rua3brZw6xeAIRj7XYGd87A2GF8OH0c/GXa8ZqVa+OtZgHLlpV7eJR02ZZHBUiROMsY0GZa5+Yrz+Xj941YjqLMSd4U4OheHt0llhphrlkoTPDqCGaZcNz+xxA5ukmv+5oLGEpWJ24hOl84IXjlJeEhHquKSEbtfG5D2cQ8i2xOLrYmNesEEv2ZaPJ2wnujuUR/HBSysOhGDR5K9wQOuUmysnVj3ykRYDOWrMHK3ce8dxcM7xXR2y450ItHUbnZLTq8ZaUQJ6/+QBW7lS9pPhCflghXDq2N/5w5Sjd8WavTU+Jwvv9tD9/kFXb/EprAADnDOmGpXdMARFQY5PaIh+Kt0dniSgMzxzUFTedM9DV4xs1RGMQjhfnADLn4NDMNfGEpUeLk+AOq7QGISLNXi/uYudXzzX5yfe/h0v+kv6AtcQSeNWkIpEZPP8QoNfWrGSQmR2V5xg38sntk1EaDukEWiGSk3FEdzoz3Gib2TEONka1HP1ffXIhAG8LsZvB3Wd/88YabdsVD8/DRQ/MBaB3I73/itH4/Jg+6Qcx8Nd31eCqtXsy564SlRLGmK8+8h0rIqbunm4RGO8aUeac0q9Gi4p0C+PKt5lG5AWZbNdiMJTVgOBk4dVKyCsKab+tfuHVul1cyK+3SBv81so9jh48QK1WxRF9mq28P8w0+Y7lYdOIVi7cFSoOIT9+YBfgbevP3dAsrfpA2hqLz78Dzwn/4cb9pp+nNHlvBp9oPIGQoj7TCRvvLa8wunu6emxPjloARHt0PMHyTthlxHjTvYqcvP+KUboF0kzrkzyroDH1b6+OKZu3Mz958xOFhUVUsRPa/b6q9m/dcKcCHtAv6IpJrayedbPZilWGR+4Tr4/Wddw017ETK24JHavjGGW/35o8N99ZXaZWKc0jM5Ko5CSYP6Ya/fmZZylTAiPk9YE6TJv+XTY287TOCUbh0bcmf997My4d2wd1QqlCp1kRWw02ef4yFk/gpy8sB5CHJs/S2/LacmtzS2s8gZUGN8BcESMUvy+kashGk6+2qFWqafLCd/zwqLDCrmCKW4ql1eUt3HJQd36/g6F49s3vn2dustJy/VtcgF32USeIM9OEz+Yazvvr93ly3OAIeUHL7FQRwdwN6g/mVpFcY6evqfRmkQRIz/bohJZYQhflyoW2uCBlp51Y5c0IK6RpUWK7jPlkRNzsrH+etcF0u1VNTDNN1WqRnN9TY96dQmG3ZuLW4GMlJHccatKZtPz2k+f3wOon4M+31f35+eeGpW17+hunOz6/mOMmkWAFHezdJjBCXrRdf2Vif3RLhm9/dmR2BYOt8POmi3LdaV72aDyBqX96X3vPF0bFWXeJzYMrml90Pu/Cwqso5Os8rP/phIVb0/ORA+YaviiwqoQIWj7oicLVj/QVVtid22tzDaDv435HvPJ2RRMMj8zZqFuLAYB43LzI+oa7p+K1752JE7tVpR3ztP7W1bYANdEdR3Qk+HDjfi2VhJ8M7VntyXGDs/CalD8fTj8XkZCidZpuVZn9sZ2wwscUw6LPv5m1Zs+R9AyLrfEEjrWka/IlISGhk80DLtrQjdGrZt41blWXt/KpnzKsO95atSeH45mfg2PmFVQsmrzdgO4kQtgJdrM58bfz+3cgIigEfLBhHxZtPZiWYI73T2MfDocUDOtlLhwzxUxMHtYdS7ap1chEc41bpsZsmHPLOeicg1uxEwKjyR9viWFA10rN79ZJlGc2HLNIluUF4rNu5l3z2zfSUwYb9+OdVtTK7XzARUErar6lEQXN0VSueo6dkDfLpmmFlTkq19zpov/8NaedAECvoZqtPWxqSNmiCynkpyWDzsxwy/HC7vrsUkL4QVhR0JTUoLm3DYfPInNdK/it4JrJEQMa8y2fmS8ndKlAB5M8TW4QCCGfSDDMWrMXm4RqM6kkQ+48tH4VbAD0D/THmw9g24FGrBc8Usw8e4xjAa/CJF6/nVupKOTF4KGKkpCmRYqa5rGk/f+68f3wzA1626fdYHLt+H6691ayK9eHWfyeloZW+A3unJaec2ipSVriQuBFFLURcVAd0kNv4vDKhc8pimLdhnyLrHOfeRHxOXLbG6+YCIaQN+sYWrpQd87hZ7i7+CAu2XYIk//wHib/YY62zWzcEn8DHrq+81CTLvGYnaeOKMDPG5qqnVkeCWvalfh9njTr5smDcfoAveb+E4sSbgBQXaYfaKzkSq5apSgE+EMs3rprx9elfUfMX1NYMQe8+YOzPD0+V1amjuiB31+ujxp9/P3Nnp47E2FF0fqYsV/wwvDZCvlJQ7tZflYWLh5N3ksCIeTNbJmaJu9Ssdywzuzh7ZTeKIubDYUNjJ93ry7VCWC+aHTPa6s1jVv9nrkI+9qTC/CckPpXrLNZURLCjkNNmPbgXPzsf1MVHHmud7PIXysbKWA96Bq9I3LV5EUhwGcUmRbNRe+bQvtUDO6RvoDoJqXhEJbdOQUPfnFs2oxwhuAWWwil/lhLTAugMwpdzbsmy5m53aAgzpwy1T9oywRCyJt1yJF91JVzNwqGAMCV4/pqr1f94gJXjmlFpmRrxrwxCWb+G7y+YjdueX6Z9r4lltC0cpF31uzVdXLRNMWFuHEhrLE1jrBCWZuxrMxnxtnAZwzFPBwfX3jNB4pMAsuLjKLFTHVZBCGFitpN8OPN+syQqukx+wR0dl5C4iAnPlMDayt1mTHbOoHo3WZmiG+eNQBv/OBMjOrrTmY3sfq710msRmdoc7NByIeIdLOZr5/RX3stLmA99O5GDL3jjYznN6YNtiIXG7LxO1zOhBXS2fZP6FyBIT2qUFORe3oK7st/rNV+0dzP9ZZiwu+oznzItVqV3aAgpiZpNUS8FnIB3m3y6t1EdDkRrSSiBBHVGz67lYg2ENFaIjo/v2baY2aGUBTCkB7e+J16zR+uHI2Xb5po+fm63foAr0vG9Nb9Bl06ZC7xZofOhdJGEORitrJ6dkIK4fQBXYQqSIQ1u4/iYKN56lmr+AexJzw1bwsAYEYyGZqVMB8/0Lk3UJBwGmhXCEQngYajLWpd3xwEr10XLdV516R+C7MKYm2ZfFWYFQAuBTBH3EhEwwBcBWA4gAsA/JWIPHMd4Iq8Vy5IflNREracgRxqbNXlfikJKSgNKzqTRDbmBzMXTVF422lC2cxouDujEX50bWBJbsj0kE0YmF0Y+ye3T8bCn00y/ezyU9xJfeEl3RzUZs0WOxHvdX2GTJw/vLv2+oqH59nW9bXDrj69zlwjaPKtsURR1xXIlryuhDG2mjGW7rQNXAzgGcZYC2NsM4ANAE7N51x2cEH1oyn2qVrbOokEw4NCmP/Xz+iPF789IU0gZpOr3bhoff7w7rppvJ2szUbImxWrAFImA2PiqUyHdvK8i79L58qSNM8eYxuKlT415fjLNWNdP6546wcYIpgLreSL59+87zjiOWrydllcSy28a5pa456kEi8UXg1XvQGIlZq3J7elQUQ3ENFCIlrY0JCe+N8J3FQRpCkWAPyvwWQTZwyrhSo+t144FCN6d9QEMU/SdN7Q7nCKWCxkSI8qPPxlndXNNXPNFfV9oRAwdYS5mYXbW/kRjYuComYHAP0t0iqc0LkCJ3SuwBPX1xf1wmI2zP3puRhXZx+inwvFbK4xyuYjzdGMGVnNj2Mn5M395JtjCZR6lGW2EGS0bxDRTADpBRSB2xhjL1t9zWSb6a/NGHsEwCMAUF9fn1Ov49posWtk2WKcMsYTTLeND2rcpBJWCCd175Ax133d9BlYedf5qCwN6/zozQZJM3PNOYNrMXttQ1bpaAd1r8KmX19k+Tk/t6bZG8PXhXWCx6+rT/PG4ZRFQpjzk3Ms254J7pVVDCy/c4qnfVpMwVzoQCgjxpJ9L36yI6fjWF3XpWN7o1Iw78YSTEvX3RpL+FYvwg8yCnnGmLkh057tAPoK7/sA2GmxWPyRrQAAEMBJREFUb95o1dUDJuSN3l+xBDMtdMw11lgye56TGc2nBxoxtGc1PtiQKtJg5r1gPFa/LhWamcYuNa4dlSUhHDe4ckYMF2vUwvt2TqV2tvPD1x/DfPusH33GUssvJk+bKgvzklvUVpXiw+nn4t7X12DxNvOEb4XCLQ82q7Hr/itG60qDRuMJDLrtde29H9HHfuFVj34FwFVEVEpE/QEMAjA/w3dyJlUazKszFAbjoBWPM6zalZ48ie+3dX8jjjSle6OYCW+u4cxZnzKRmWm+xjbc8dlhOSUOE+kuFDThR+c2eX464+By4cmpyaRTs5yVIB9Q26HgWTSLhV6dyhEJKUgYAj4Lbcpx6/x2Jh7xHMa0BmVFNNjnS74ulJ8nou0AxgOYQURvAgBjbCWA5wCsAvAGgJsYY57l7tTK2wVMyhun6l9/aoHpfvy6Pz3QiJ2H0zNUmuUG5/07LtjkzWZC+dS27dvZPBDt0WtTdn/RTx5ICWajIBcFdmnIWRu4ScfrCGW3ubK+L07pV+Pb+Yyl5y46uSe+e94g385vhl0N4WywGyysXCgBb2o4F4q8fA4ZYy8BeMnis7sB3J3P8Z3CtZCgLLRxjPbEBVvMp9SZxraIouAb5w7QFeDgx45msMkb5aPTPOMb77kQBGD/8da0EPWBtR201x3LS7DvWIsm7PlimDHMXLy3VRaVnox8ZWIdXlq8w1GqgMqSkG6GUUh+84WRvp4vpJCur3nhyZMtYj3fblWl2GtR1CYTdvOBDqVhzLz5M5h0/3tpfVSaa4qMlHdNgRviMk4z45kJ5xe/PUF7HQ4RfjRlsO7z11fsBqD3rjEV8oZtjUL0qJ2wDSkERSHUVpVq6Z/NeP7G8bhr2nDNna0+6UVi1L7FscWLGdvyO8/HzB9+xvXjtgWIKCfPFa/o2bEMzcLC66DuHWz2tifTgvLA2koQqak9RKSQLzK4d43Xmnx9vxqM6O1fFK3dlFWsXWvmgTFGCKYy84J5KJl6tbNQxtCJd83W/Y3aa6uyetlQ17US102o097/7gsj8b83TUS3qjL8VdAo4wmGx6+rx08vsM5waSSqVRPK3M2V5KDUHlFI9ScnAr537okFa0fvpDKgEOk0ead1js3IlF2SiBAJKdhgKC4fJCEfiBBR5pOQf/5bEzLv5CKDu1ubGX520VDttZktnYhApNreIzbCa2jP1DnMarMaj91JyCWTT14ZK8oiIS13jzgAtcYSOG9o96xiAHj2TrP8+5IUIUXNfcRYYd2QX/3uGdh7tAXfeGqhroawcVE4G0TvMStKQkqay2a7cqFsC8QDapOvLA2jpiJimr9FvFZRhl8vaMVKMnGZMZpUpNUu7hvptvFJgpC9TRhosuXD6edmTL8gmmyMC2NOGNRNneZfZ5JDXpJCoVQJxkI+QjWVJaipLEFIIV3lMad1ju2ws+tHQoRjLfpzBEkxCMSVBNUmDwAvfnsizjZJu0sWNmoxfQDfbPRBF4llmM6uNrhsip46+Uxpe3Uqz5hITSwckktRh27VZdhy70WYGqC0sV5AJLohF15RUgg6TT4fcw3nfy4cir9fP870MzOf/CCZawIhFlNaSOE7qNv071qJa07rl7ZdrGojTi3FDst/DzMXSk6mqvRvrtytex9RFJzaX10ctcsL4gaif399nX8uhe2NEFGqUHYRPEIhhTTzCVF+0bg8sdvpA7rgnCHmVaLMhHyQzDWBEPJBjXjlmAUziZGZYoZCXQZJzQfd+jbvP5aK+psyLN3e3a1K71aoKKSdw2uPDD441VREXFnklZijKKR5chWDoqQQaea5kpCi8wDLFq6QVJRm13+ClLsmEFeiedcE4mrSMXq9fNGQtldcoBS1Ej715kL5S6enp/sVzSClJtrLFfWqF88vLxmB//vOGQCAwd1VD6NOHiy8ivDBLWiJ54oNcUZWBDIeIYW0Sk0lYSUvTf73l4/CWz88yzIDKQDsONSUti1ImnwgFl4H1lbiP18/DUN6ts0iIZkwavLG96L9MGwi5Pm2oSa/z+y1Kf9gM/PLdRPqcM6QbujXJZUGYPrUIZg0rBtG9O6YzWVkDbfJF4OdOMjsF3K4uFUTOR9CCmnuwyWhVHHv7+cQhVsWCeEkGy81u+8FhUAI+aqyCCacmF0RibaE0c/dqNmKIdglpuYa9YW4gMVr324/mNJizBa4iEgn4AFVu8q2aEcu8HbnWtRb4oxjLakAt2JIAaEQaQuvkZAirBf41zavS3z6SXCuJMDY5XEBbDR5hZtr1G2iG6Io3DnFVrGez9Lba5CSX4iR1cWQhVPs7yVhRWtfgOSurwRCkw86Rk124VZ9DhsxY56ogRxK+tfzBcz6DEmvii2neFxzjZVC3kvEdZliKHsnOlBEQimfeTcWhf989RgMqLXPQDppqLkXTlul8HdUkpE0IWcQxqL2bjbd5t41o/p2wvq7p2KgRScvNiHPyxh++fR0F1KJe4iafDGYKUQHCtFc48Zg/7lRvTC8l/Va0rRRvfDYdeb+9G0Vqcm3AYx+7nZTajNXw5DhoTm1fxccboql7edG0ImbVJdFsOVe62pSEncQcyQVo7mGe9r44SL9wNVjPD+H3xT+jkoyku5dY33bzPJgf7TpgO79ocZW7DvWoltwA4pPk5f4Q52wsF4UmrzOXKNo6bClk1VuFP6OShxgnXbXiJl/72FDtSieZvjFT7brtucTdCJpu3xXyDxZDN41oiYfCZGw8Fr4trVFpJBvE1gX0DBSaRLZd4JQH1XkjpdX6t5LTb59MkjwIy8GTV40y4QVRYuslkI+Nwp/RyUZyUb2VpakL7NkmuaedZKaAK3ITPKSAlAMQl50mRVfOy2iky1WjghBofB3VJIRo+y10+TFHB2ZUvlqx+euitLo2e4pBmVZbINoPXpmwTZPzlcMA5uX5FvI+3IiWklECSKqF7bXEVETES1J/v0t/6a2X/p31WsadrJY1OSfv1EtcpLpueVVeIKa+0eSGZ57vxgC4kSzjKjQNGXImCoxJ9/HegWASwHMMflsI2NsdPLvxjzP064xahp2mry48Gplw7zlfH29V17AmpttJO2PHsk+YFdy0i/E/v3hxlRlJ+bRmlExZN70kryEPGNsNWNsrVuNkTjDbkqtGDwTzJg6oofufd+acsy79Vx86zMDXWmfpO3RvVoV8sUQ8SrWERYznXo1yZgwsIs3By4SvLyj/YloMRG9R0RnWu1ERDcQ0UIiWtjQ0OBhc9o2j15bL+SNd6Z5cA3FqKkYZwLxBEPPjuWB12gk1tw5bTjumja8KAQed/mdNqoXugrVw7zy/rr61PQU3EEio5AnoplEtMLk72Kbr+0CcAJjbAyAmwH8h4hM8wAzxh5hjNUzxupra6W5wIrJw7rjwS+OBQD062LuEmmEy2yj6DaacYrBDispLB1Kw7huQl1RDPSXjO4FAOhYHtHNWsW6Ce4S7P6fMa0BY2xStgdljLUAaEm+XkREGwGcBGBh1i2UaJzavzMevbYeZ53kLM2vleJjfI6LLZ2BpH0j1vYVB52BycVht2nxyDWzWPAkdw0R1QI4wBiLE9EAAIMAbPLiXO2NySYl+gBg2Z1TTIw4SeFt+CBdkw92J5e0LZqSWScPNrbqNPnbLhzqyfm88r8vFvJ1ofw8EW0HMB7ADCJ6M/nRWQCWEdFSAM8DuJExdsDqOJL8qS6LoMpQ4qymQp3eThmmX2g1+sNLTV5STDyz4FMAwKvLdukUkl6dyj05X21yravSJO9TEMhLk2eMvQTgJZPtLwB4IZ9jS/KnS4dSLLhtUpot02h3jcqcNZIiQlQ6/ChH2KemAvP/5zzdIm+QKLy/lMRTaqtK08wzxvdSk5cUEzyOo6rMv0zo3arLAluBTAr5dog010iKmQuScRyl4ZDmNnn7Z4cVskltGink2yFkuOveuaZJJNlTGlJt4/FEQnNuHNXHupqTxB5ZGaodImryv/vCSHx2ZK8CtkYi0VNdHsbXz+iPS8b0xoOzNgDQF6uXZIcU8u0Q0SZ/eX3fArZEIkmHiPCzpHnmN18YiXOHdMOI3lKTzxVprmmHFEFQo0TiiI7lEVwxTioi+SCFfDtE5o2XSNoPUsi3Q+xSFUskkmAhhXw7hPsDhwPqFyyRSFLIhdd2ys8uGoozB8msnxJJ0JFCvp3y9TMHFLoJEonEB6S5RiKRSAKMFPISiUQSYKSQl0gkkgAjhbxEIpEEGCnkJRKJJMBIIS+RSCQBRgp5iUQiCTBSyEskEkmAIcaKpyoQETUA2JrHIboC2OdScwpJUK4DkNdSrATlWoJyHUB+19KPMWYawl5UQj5fiGghY6y+0O3Il6BcByCvpVgJyrUE5ToA765FmmskEokkwEghL5FIJAEmaEL+kUI3wCWCch2AvJZiJSjXEpTrADy6lkDZ5CUSiUSiJ2iavEQikUgEpJCXSCSSABMIIU9EFxDRWiLaQETTC90eJxDRFiJaTkRLiGhhcltnInqbiNYn/9cI+9+avL61RHR+Adv9BBHtJaIVwras201EpySvfwMRPUDkf+FZi2u5k4h2JO/LEiK6sI1cS18imk1Eq4loJRF9P7m9Td0bm+toc/eFiMqIaD4RLU1ey13J7f7eE8ZYm/4DEAKwEcAAACUAlgIYVuh2OWj3FgBdDdt+C2B68vV0AL9Jvh6WvK5SAP2T1xsqULvPAjAWwIp82g1gPoDxAAjA6wCmFsm13Angxyb7Fvu19AQwNvm6CsC6ZJvb1L2xuY42d1+S5+2QfB0B8DGA0/2+J0HQ5E8FsIExtokx1grgGQAXF7hNuXIxgH8kX/8DwCXC9mcYYy2Msc0ANkC9bt9hjM0BcMCwOat2E1FPANWMsXlM7cFPCd/xDYtrsaLYr2UXY+yT5OujAFYD6I02dm9srsOKorwOAGAqx5JvI8k/Bp/vSRCEfG8A24T322HfKYoFBuAtIlpERDckt3VnjO0C1M4OoFtye7FfY7bt7p18bdxeLHyHiJYlzTl8Kt1mroWI6gCMgao5ttl7Y7gOoA3eFyIKEdESAHsBvM0Y8/2eBEHIm9mm2oJf6ETG2FgAUwHcRERn2ezbVq/Rqt3FfD0PARgIYDSAXQDuS25vE9dCRB0AvADgB4yxI3a7mmwrmusxuY42eV8YY3HG2GgAfaBq5SNsdvfkWoIg5LcD6Cu87wNgZ4Ha4hjG2M7k/70AXoJqftmTnJoh+X9vcvdiv8Zs2709+dq4veAwxvYkH8wEgEeRMosV/bUQUQSqYPw3Y+zF5OY2d2/MrqMt3xcAYIwdAvAugAvg8z0JgpBfAGAQEfUnohIAVwF4pcBtsoWIKomoir8GMAXACqjtvi6523UAXk6+fgXAVURUSkT9AQyCuhBTLGTV7uQU9SgRnZ70ErhW+E5B4Q9fks9DvS9AkV9L8tyPA1jNGLtf+KhN3Rur62iL94WIaomoU/J1OYBJANbA73vi52qzV38ALoS6Cr8RwG2Fbo+D9g6Auoq+FMBK3mYAXQC8A2B98n9n4Tu3Ja9vLQrgvSG042mo0+UoVA3ja7m0G0A91Ad1I4AHkYy+LoJr+SeA5QCWJR+6nm3kWs6AOoVfBmBJ8u/CtnZvbK6jzd0XACMBLE62eQWAO5Lbfb0nMq2BRCKRBJggmGskEolEYoEU8hKJRBJgpJCXSCSSACOFvEQikQQYKeQlEokkwEghL5FIJAFGCnmJRCIJMP8PhhGibfhgcI8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"####Generating correlated curves:\n",
"np.random.seed(2) \n",
"def gen_correlated_curve(rho, num=100):\n",
" num_with_runup = num+5000\n",
" y = np.zeros((num_with_runup,))\n",
" for i in range(1,num_with_runup):\n",
" y[i] = rho * y[i-1] + np.random.normal()\n",
" return y[-num:]\n",
"\n",
"timeseries = gen_correlated_curve(0.99, 3000)\n",
"print(timeseries.mean())\n",
"plt.plot(timeseries)\n",
"plt.axhline(0)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x12b761e50>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD4CAYAAADiry33AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXxV1b338c8vJxMhCSETIRMJIcwzIQJSUKFlqEqrrYqdJ2odqn06XGvvffq09/b23j5Prdpap2qrrYpaJ2pxQAQcGEyYA0kgBELmhMzzcM56/sjBG0MwJ+Qk++Sc3/v1ysucvdc5+5ctfNlZe621xRiDUkop7+VndQFKKaWGlwa9Ukp5OQ16pZTychr0Sinl5TTolVLKy/lbXUB/oqOjTUpKitVlKKXUqLF///5zxpiY/vZ5ZNCnpKSQnZ1tdRlKKTVqiEjRxfZp141SSnk5DXqllPJyLgW9iKwVkXwRKRCRu/vZLyLygHP/ERFZ2GtfhIj8XUTyRCRXRJa68wdQSin1yQYMehGxAQ8C64CZwEYRmdmn2Tog3fm1CXio1777gTeMMdOBeUCuG+pWSinlIleu6DOBAmNMoTGmE9gMbOjTZgPwlOmxF4gQkYkiEg6sAB4HMMZ0GmPq3Vi/UkqpAbgS9AlAca/XJc5trrSZDFQDfxaRgyLyJxEZO4R6lVJKDZIrQS/9bOu75OXF2vgDC4GHjDELgBbggj5+ABHZJCLZIpJdXV3tQllKKaVc4UrQlwBJvV4nAmUutikBSowx+5zb/05P8F/AGPOoMSbDGJMRE9PvmH81Cuky2EpZz5UJU1lAuoikAqXATcDNfdpsAW4Xkc3AZUCDMaYcQESKRWSaMSYfWAUcd1v1yiPtOVXDC/uLOXS2nuK6Vq6YFsuNGUlcMS0Gf5uO6FVqpA0Y9MaYbhG5HXgTsAFPGGOOicgtzv0PA1uB9UAB0Ap8o9dH3AE8LSKBQGGffcrLPJd1lnteziE82J9FkyJZNiWKN3Iq2Xa8kstSI3nsaxmEBwdYXaZSPkU88VfrjIwMo0sgjD4PbD/JvdtOsGJqDH/80kJCg3quI7rsDl46UMK/vpLDlNgwnvzmYmLDgi2uVinvIiL7jTEZ/e3T36OVW7x1rIJ7t53guoUJPP61jI9CHiDA5seNi5N5/GuLKapp4cZH9tLQ1mVhtUr5Fg16NWTNHd38fMsxpseF8d/XzyXgIv3wK6bG8OQ3MymubeWnLx3RG7VKjRANejVkv9t2gorGdn71+TkXDfnzFqdE8uM109h6tIKn950doQqV8m0a9GpIckob+PMHp7k5M5lFk8a79J7vfGoyK6bG8MvXjpNf0TTMFSqlNOjVkPz+nZOEjwngJ2unu/wePz/h3hvmERrkz7+9kqNdOEoNMw16dclK6lrZdrySjZnJjBszuCGT0aFB/HjNND48U8uWw33n3yml3EmDXl2yv+4tQkT48pJJl/T+GzKSmJ0Qzq+35tHS0e3m6pRS52nQq0vS1mnnuaxiPjNzAgkRYy7pM2x+wi+unU1FYzsP7ihwc4VKqfM06NUlefVQKfWtXXx9WcqQPmfRpPF8fkECj79/mqrGdvcUp5T6GA16dUme2lPE9LgwMlMjh/xZP1g9FbvD6FW9UsNEg14N2qnqZo6XN3JDRhIi/a1QPTjJUSHcsDiJZz48S0ldqxsqVEr1pkGvBu2NnAoA1s2Jc9tn3nHVFESE32/Xq3ql3E2DXg3a1qPlLEyOYOK4S7sJ25+J48bwpcuS+fuBEs6ca3Hb5yqlNOjVIBXVtHCsrJH1cya6/bO/d0Ua/n7CQztPuf2zlfJlGvRqUF53dtusne2+bpvzYsOC2ZiZzIsHSiitb3P75yvlqzTo1aC8nlPB3MRxJI4PGZbP37RiMiLwyC69qlfKXTTolctK69s4XFzPutnu77Y5Lz5iDF9YlMjmrGIdV6+Um2jQK5e9k1cFwJpZE4b1ON9bOQW7w/Dou4XDehylfIUGvXLZ+yerSYgYQ2r02GE9TnJUCBvmxfP0vrPUNHcM67GU8gUa9Moldodh96kalk+JdsskqYHceuUU2rvtPPHB6WE/llLeToNeueRIST1N7d0sT48ekeNNiQ1l/ZyJPLm7iIZWfb6sUkOhQa9c8v7Jc4jA5VNGJugBbr9yCs0d3fxl95kRO6ZS3kiDXrnkvYJzzIoPJ3Js4Igdc8bEcFbPmMATH5ymWderV+qSadCrAbV0dHPwbN2IXs2fd8dVU2ho6+Jve4tG/NhKeQsNejWgD0/X0mU3fGpKzIgfe15SBCumxvCn9wpp67SP+PGV8gYa9GpA7508R5C/Hxkp4y05/h1XTeFccyebs85acnylRjuXgl5E1opIvogUiMjd/ewXEXnAuf+IiCzste+MiBwVkUMiku3O4tXI2He6hoXJ4wkOsFly/MUpkVyWGskjuwrp6NareqUGa8CgFxEb8CCwDpgJbBSRmX2arQPSnV+bgIf67L/SGDPfGJMx9JLVSGru6Ca3vJHFbniS1FDccVU6FY3tPJ9VbGkdSo1GrlzRZwIFxphCY0wnsBnY0KfNBuAp02MvECEiw7cgihoxB8/W4TCw2KJum/MunxJFZkokD7xTQGunjsBRajBcCfoEoPdlVIlzm6ttDPCWiOwXkU2XWqiyRtaZOvwEFiRbG/Qiwo/XTqO6qUPH1Ss1SK4EfX/z3c0g2lxujFlIT/fObSKyot+DiGwSkWwRya6urnahLDUSss/UMmNiOKFB/laXwuKUSK6cFsPDO0/R0KazZZVylStBXwIk9XqdCJS52sYYc/6/VcDL9HQFXcAY86gxJsMYkxETM/LD+NSFuuwODhXXszjF2v753n60ZhqN7d26Xr1Sg+BK0GcB6SKSKiKBwE3Alj5ttgBfdY6+WQI0GGPKRWSsiIQBiMhY4DNAjhvrV8Mot7yR1k47iyZZ223T26z4cWyYH8/j75+mpK7V6nKUGhUGDHpjTDdwO/AmkAs8b4w5JiK3iMgtzmZbgUKgAHgMuNW5fQLwvogcBj4E/mmMecPNP4MaJlln6gAsGz9/MT9ZOx0R+PXreVaXotSo4FLHqzFmKz1h3nvbw72+N8Bt/byvEJg3xBqVRfYX1ZI4fgwTx42xupSPSYgYw3dXpHH/9pN8bWktmRYP/VTK0+nMWNUvYwxZZ+rI8KBum95uWZnGxHHB/OIfx7A7+o4NUEr1pkGv+lVS10Z1UweLPOhGbG9jAm3cs34Gx8oadcEzpQagQa/6dai4HoAFSREWV3JxV8+dyIqpMfzmjTzK6tusLkcpj6VBr/p1uLieQH8/psWFWV3KRYkIv/rcbBwG/u2VHHpuFSml+tKgV/06XFLP7PhwAmye/UckKTKEH35mKtvzqnjtSLnV5SjlkTz7b7GyRLfdwdHSBuYneeaN2L6+viyFeUkR/NurOVQ2tltdjlIeR4NeXSC/son2LgfzksZZXYpL/G1+/O6GeXR0OfjRC4dx6CgcpT5Gg15d4HBxAwDzPfhGbF+TY0L516tn8N7Jczy554zV5SjlUTTo1QUOF9cTERJAcmSI1aUMys2ZyayaHsuvX88jp7TB6nJcUtnYzt/3l3D3i0d4PruYbrvD6pKUF9KgVxc4XFLPvMQIRPpblNRziQi/+cJcIkMCue2ZAx69wqUxhnvfyuey/9zOj144zKuHyvjJ34/wmfve5Z28SqvLU15Gg159TEtHNycqm0ZVt01vUaFBPPilBZTWtfHjFw575JBLh8Pwi38c54F3CrhuYQJbv/8pjv1iDY98ZRF+Imx6aj8Hz9ZZXabyIhr06mNyShtwmNHVP9/XokmR3L1uOm8dr+Tx909bXc4Ffr7lGH/ZfYZvL0/lt1+cx8z4cPz8hDWz4njxe8uYEB7MHc8epLHdc38jUaOLBr36mMMlPTNi5yaOjhE3F/Ot5amsmTWB/3o9j/1FtVaX85E3csr5694ivr08lZ99dsYF3WPjxgTwwMYFlDe0c89LRz3yNxI1+mjQq485WtpIQsQYokKDrC5lSESE//vFeSSMH8NtTx+kprnD6pKoamrnpy8dZXZCOP+ybvpF74EsmjSe//Xpqbx2pJy3c6tGuErljTTo1cfklDYwOyHc6jLcIjw4gD9+aSG1rZ3c9dwhS1e5NMbw0xeP0tJp53c3zB9wxvF3V0wmcfwYHtYnaSk30KBXH2lq7+L0uRbmJIzubpveZsWP45fXzuK9k+e47+0TltWx9WgF2/Oq+Je100mfMPD6Qf42P769PJX9RXUe1fWkRicNevWRY2WNAMzyoqAHuCkzmRsyEvn9OwVszx35oYvtXXb+641cpseF8fVlKS6/74bFSUSEBPDIrsLhK075BA169ZHzk4xmx3tX0AP8csNsZieEc9dzhyiqaRnRY/9l9xmKa9v418/OxObn+tyEkEB/vrJkEttyKymsbh7GCpW306BXH8kpbSAuPJiYsNF9I7Y/wQE2HvpSzzj1W/52gLZO+4gc91xzBw++U8Cq6bEsT48e9Pu/ujSFAJsfT3zgecNE1eihQa8+crS0gdle1m3TW1JkCPfdNJ+8ikZ+9srIDF287+0TtHXZ+en6GZf0/piwINbPjuO1I+V0duvyCOrSaNArAJo7uik81+I1I24u5sppsdy5Kp2XDpTy9L6zw3qsE5VNPLPvLF9eMokpsaGX/DnXzIunvrWLDwrOubE65Us06BUAueWNGINXjbi5mO9flc4V02L4xT+ODetSA7/6Zy6hQf7cuSp9SJ/zqfQYxo0JYMvhMjdVpnyNBr0C4GhJz41YXwh6Pz/hvhvnMyE8mFufPjAsk6l25lex60Q131+VzvixgUP6rEB/P9bNjuOtYxW0d43MvQXlXTToFQA5ZQ3EhAURGx5sdSkjIiIkkIe/vIialk6+v/mgWydTddsd/OfWXCZFhfCVpZPc8pnXzIunpdPOjjydKasGT4NeAT0jbnzhar632Qnj+I8Ns/mgoIbfvpXvts99ck8RJyqbuWf9DIL8bW75zCWTo4gODeIfR7T7Rg2eBr2itbObgqpmZsd7943Y/tywOImNmUn8cecp3jpWMeTPq2ps53fbTrByagyfmTnBDRX2sPkJn50Tx/bcKlo6ut32uco3uBT0IrJWRPJFpEBE7u5nv4jIA879R0RkYZ/9NhE5KCKvuatw5T655U04DF49tPKT/PyaWcxJGMcPnz/M6XNDm0z169fz6Ox28H+uneX2B7esmR1HR7eD3adq3Pq5yvsNGPQiYgMeBNYBM4GNIjKzT7N1QLrzaxPwUJ/9dwK5Q65WDYvzM2LnjPKliS9VcICNh768EJtN+N7f9tPaeWlXzLsLzvHywVI2rZhMavRYN1cJGZMiCQ3yZ0e+9tOrwXHlij4TKDDGFBpjOoHNwIY+bTYAT5kee4EIEZkIICKJwGeBP7mxbuVGOaUNRI0NJM5HbsT2J3F8CPfftID8yiZ+9nLOoCdTnWvu4K7nDpEaPZbbrpwyLDUG+vtx+ZQoduZV6Tr1alBcCfoEoLjX6xLnNlfb3Af8BNBpfR7q/IzY0faMWHdbOTWGH6yeyssHS/nzB2dcfp/DYfjBc4eob+viwZsXMibQPTdg+3PltFjKGto5Ualr3yjXuRL0/f3t73s50W8bEbkaqDLG7B/wICKbRCRbRLKrq6tdKEu5Q3uXnZNVzV4/I9ZVt185hU/PnMAvXzvOn95zbdXIB945yXsnz/Hza2Yyc5hvaF8xLRZAu2/UoLgS9CVAUq/XiUDfMV4Xa3M5cK2InKGny+cqEflbfwcxxjxqjMkwxmTExMS4WL4aqryKJuwO43NDKy/Gz0948OaFrJ8Tx3/8M5d7t53AcZEx9g6H4ddbc7nv7ZN8fkECN2cmD3t9ceOCmTExXMfTq0FxJeizgHQRSRWRQOAmYEufNluArzpH3ywBGowx5caYnxpjEo0xKc73vWOM+bI7fwA1NEedN2JneeHSxJcq0N+PB25awBcXJfLA9pNc99BujjifpXteSV0rdz53iEfeLeTLS5L5f1+cN2JdX1dOiyG7qE4fHq5c5j9QA2NMt4jcDrwJ2IAnjDHHROQW5/6Hga3AeqAAaAW+MXwlK3fKKWkgIiSAxPFjrC7Fo/jb/PjNF+ayZHIUv349jw0PfsD0uHBSokKob+1iT2ENIvDjNdO49Yq0Eb2/ceX0WP648xTvnTjHZ+dOHLHjqtFrwKAHMMZspSfMe297uNf3BrhtgM/YCewcdIVqWOWU9cyI9fUbsf0REa5flMinZ03gifdPc6SkgfyKJkTgB6unct3CBJIiQ0a8rgVJEYQF+/N+QbUGvXKJS0GvvFNHt50TlU18a/lkq0vxaOHBAdy1eqrVZXzE3+bHZamROnFKuUyXQPBhJyqa6bLrjdjRaGlaNEU1rZTWt1ldihoFNOh92JHSnhuMOrRy9FmWFgXAHr2qVy7QoPdhOaUNjBsTQLIF/cxqaKZNCCNybCC7T+lTp9TANOh92JESvRE7Wvn5CUsnR7HnVI0uh6AGpEHvo9q7em7E+uqKld5gaVoU5Q3tnKlptboU5eE06H1UfkUTXXbDXB9dsdIbnO+n1+4bNRANeh91fkasjrgZvVKjxxIXHqzDLNWANOh91FGdETvqiQhL06LYq/30agAa9D7qaKneiPUGS9OiqGnp1GWL1SfSoPdB52/EarfN6Kf99MoVGvQ+KK+iiW6H3oj1BonjQ0iODNF+evWJNOh90NGS8zNiNei9wbK0KPYW1mC/yLr5SmnQ+6CjpQ1Ejg0kIUJvxHqDpWlRNLV3c6yswepSlIfSoPdBR0r0GbHeZOlH/fTafaP6p0HvY84/I3audtt4jdiwYNJjQ3WBM3VRGvQ+5nh5I3aH0f55L7MsLYqsM7V0djusLkV5IA16H5PjnBGrI268y9K0KFo77Rzu82xbpUCD3uccKWkgamwgE8cFW12KcqMlk6MQgd0F2n2jLqRB72NyShuYk6g3Yr1NREggs+LDdeKU6pcGvQ9p69QZsd5sWVo0B8/W09Zpt7oU5WE06H3I8fJGHEZXrPRWy9Ki6LQ7yC6qtboU5WE06H3I+RmxcxMjLK5EDYfFKZH4+wkfaD+96kOD3occLW0kOjSICeFBVpeihsHYIH8WJEewR/vpVR8a9D7kaGk9c/VGrFdbmhbN0dIGGtq6rC5FeRANeh/R2tlNQVWzTpTycpenReEwsK9Qu2/U/9Cg9xFHShpwGFiQpP3z3mx+cgTBAX667o36GJeCXkTWiki+iBSIyN397BcRecC5/4iILHRuDxaRD0XksIgcE5FfuPsHUK45eLbnRux8DXqvFuRvY3FKpI6nVx8zYNCLiA14EFgHzAQ2isjMPs3WAenOr03AQ87tHcBVxph5wHxgrYgscVPtahAOnq0jNXos48cGWl2KGmbL0qI5UdlMdVOH1aUoD+HKFX0mUGCMKTTGdAKbgQ192mwAnjI99gIRIjLR+fr8wywDnF/6dIQRZozhYHG9dtv4CH28oOrLlaBPAIp7vS5xbnOpjYjYROQQUAVsM8bs6+8gIrJJRLJFJLu6utrV+pULSuraqG7qYEGyBr0vmJ0wjrBgf122WH3ElaDvbyxe36vyi7YxxtiNMfOBRCBTRGb3dxBjzKPGmAxjTEZMTIwLZSlXHSzu6Z9fkDze4krUSLD5CUsmR/GBXtErJ1eCvgRI6vU6ESgbbBtjTD2wE1g76CrVkBw8W0dwgB/T48KsLkWNkMvToiiubaO4ttXqUpQHcCXos4B0EUkVkUDgJmBLnzZbgK86R98sARqMMeUiEiMiEQAiMgZYDeS5sX7lgoNn65mbGIG/TUfT+oplU6IB7adXPQb8m2+M6QZuB94EcoHnjTHHROQWEbnF2WwrUAgUAI8Btzq3TwR2iMgRev7B2GaMec3NP4P6BB3ddo6XNWr/vI9Jjw0lOjRIx9MrAPxdaWSM2UpPmPfe9nCv7w1wWz/vOwIsGGKNagiOlTXSaXewIEn7532JiLAsLYrdp2owxuiyFz5Of5f3cgeK6gD0it4HLZ8STXVTB/mVTVaXoiymQe/lss7UkhQ5hgnh+uhAX7Nias/otV35OlzZ12nQezFjDNln6shMibK6FGWBuHHBTI8LY6cGvc/ToPdip6pbqGnpJDNV++d91cqpMWQX1dLc0W11KcpCGvRe7MPTPY+UW5wSaXElyiorp8XQZTc6S9bHadB7sawztUSHBpIaPdbqUpRFMiZFEhJoY2d+ldWlKAtp0HuxD0/XsjglUofW+bBAfz+WpUWz60Q1PaOglS/SoPdSZfVtlNa3abeNYuW0GErq2ig812J1KcoiGvReKutMT/98ZqoGva+7wjnMckeedt/4Kg16L7XvdC2hQf7MmBhudSnKYkmRIaTHhvKOBr3P0qD3Ulmna1k0aTw2P+2fV7B65gQ+PF1LQ1uX1aUoC2jQe6GqxnZOVjWzNE0nSqkeq2fE0u0w7Dqhk6d8kQa9Fzr/wInlzqVqlZqfNJ6osYG8fbzS6lKUBTTovdD7J2uICAlgpvbPKyebn3Dl9Fh25lfRZXdYXY4aYRr0XsYYwwcF57g8LRo/7Z9XvayeMYHG9u6PRmQp36FB72VOVbdQ0djO5dpto/r4VHo0gf5+bM/V0Te+RoPey3xQoP3zqn9jg/xZlhbFtuOVOkvWx2jQe5kPCs6RFDmG5KgQq0tRHmjNrDjO1raSW64PI/ElGvRepNvuYE9hjV7Nq4v69MwJ+Am8kVNudSlqBGnQe5EjpQ00tXezLE2DXvUvOjSIxSmRvJ5TYXUpagRp0HuRHXlV+EnPTTelLmbd7DhOVjVTUNVsdSlqhGjQe5HtuVVkTIokIiTQ6lKUB1s7eyIAbx7Tq3pfoUHvJcob2jhe3shVM2KtLkV5uLhxwSxIjuB17af3GRr0XuL82OjVGvTKBWtnxZFT2khxbavVpagRoEHvJd7JqyI5MoS0mFCrS1GjwPo5Pd03rx3Rq3pfoEHvBdo67XxQcI6rpsfqYwOVS5IiQ1iYHMGrh0qtLkWNAJeCXkTWiki+iBSIyN397BcRecC5/4iILHRuTxKRHSKSKyLHROROd/8ACnafOkdHt4NV2m2jBmHD/ATyKprIr9DJU95uwKAXERvwILAOmAlsFJGZfZqtA9KdX5uAh5zbu4EfGmNmAEuA2/p5rxqit3OrGBto08cGqkFZP2ciNj9hy2G9qvd2rlzRZwIFxphCY0wnsBnY0KfNBuAp02MvECEiE40x5caYAwDGmCYgF0hwY/0+r9vu4K1jFVw5PZYgf5vV5ahRJCYsiMunRPPqoTJd+8bLuRL0CUBxr9clXBjWA7YRkRRgAbCvv4OIyCYRyRaR7OpqfQqOq/adrqWmpZOr5060uhQ1Cm2YF09JXRsHztZbXYoaRq4EfX939/r+8/+JbUQkFHgRuMsY09jfQYwxjxpjMowxGTExMS6UpaBn1ERIoI0rpmn/vBq8z8yaQJC/H1v0pqxXcyXoS4CkXq8TgTJX24hIAD0h/7Qx5qVLL1X11W138EZOOatnTCA4QLtt1OCFBQewesYEXjtSTrc+ecpruRL0WUC6iKSKSCBwE7ClT5stwFedo2+WAA3GmHLpGev3OJBrjLnXrZUr9hTWUNfaxWe120YNwbXz46lp6eSDUzVWl6KGyYBBb4zpBm4H3qTnZurzxphjInKLiNzibLYVKAQKgMeAW53bLwe+AlwlIoecX+vd/UP4qn8eKWdsoI2VU7WrS126K6bFEBbsr2PqvZi/K42MMVvpCfPe2x7u9b0Bbuvnfe/Tf/+9GqIuu4M3jlXw6ZnabaOGJsjfxrrZcfzzSDntn7frnycvpDNjR6l38qqob+3i2vnxVpeivMCG+Qm0dNp5J0+fJ+uNNOhHqb/vLyEmLIgV6dpto4ZuyeQoYsKCtPvGS2nQj0LnmjvYkVfFdQsS8Lfp/0I1dDY/4Zq58ezIq6a+tdPqcpSbaUqMQq8eKqPbYbh+UaLVpSgvct3CBDrtDv5xuO/oaTXaadCPMsYYXsguZl7iOKZOCLO6HOVFZieMY8bEcF7YX2J1KcrNNOhHmWNljeRVNPGFjKSBGys1SF9YlMiRkgZd0dLLaNCPMpuzzhLo78c1OklKDYPPzY/H3094Ibt44MZq1NCgH0Wa2rt4+UAp18yN1weAq2ERFRrEqhmxvHKolC5dEsFraNCPIq8cLKWl085Xlk6yuhTlxb64KIlzzZ3s0DH1XkODfpQwxvDUniLmJo5jflKE1eUoL7ZyWgyxYUFsztLuG2+hQT9K7Dtdy8mqZr68RK/m1fAKsPlxQ0YSO/OrKK1vs7oc5QYa9KPEX/cUMW5MANfM1SUP1PC7cXESBnhOr+q9ggb9KFBc28obxyq4cXESYwJ1wSk1/JIiQ1iRHsNzWWd1nXovoEE/Cjz+/mkE+PqyFKtLUT7k5suSqWzs0IXOvIAGvYerb+3kuaxirp0XT3zEGKvLUT5k1fRYYsOCeHrfWatLUUOkQe/h/ra3iLYuO99ZMdnqUpSP8bf5sTEzmV0nqjl9rsXqctQQaNB7sPYuO3/ZXcSKqTHMmBhudTnKB31pSTIBNuHJ3WesLkUNgQa9B3v5YCnnmjv4rl7NK4vEhgVz9dx4XsguprG9y+py1CXSoPdQDofhsfcKmRUfzrK0KKvLUT7sm5en0tJp54VsXdVytNKg91Bv51ZSWN3Cd1emIaKP3VXWmZM4joxJ43ly9xnsDmN1OeoSaNB7qEffLSRx/BjWz46zuhSl+MblqZytbWXb8QqrS1GXQIPeA+0vqiW7qI5vLU/VRwUqj7B2dhyp0WP5w44CjNGr+tFGU8QDPbyrkHFjArhBHy6iPITNT/jeyjRyShvZeaLa6nLUIGnQe5jc8ka2Ha/kG5enMDbI3+pylPrI5xYkkBAxhj+8o1f1o40GvYf5w44CQoP8dbkD5XEC/f347srJ7C+qY29hrdXlqEHQoPcgBVXNbD1azleWTtInSCmPdENGEjFhQdy7LV+v6kcRl4JeRNaKSL6IFIjI3f3sFxF5wLn/iIgs7LXvCRGpEpEcdxbujf64o4AgfwXHb5UAAAzzSURBVD++vTzV6lKU6ldwgI07V6WTdaaOt45XWl2OctGAQS8iNuBBYB0wE9goIjP7NFsHpDu/NgEP9dr3F2CtO4r1ZkU1Lbx6uIwvXTaJqNAgq8tR6qJuWpxEWsxY/vv1PH2u7CjhyhV9JlBgjCk0xnQCm4ENfdpsAJ4yPfYCESIyEcAY8y6gHXoDeGjnKWx+wiZd7kB5OH+bH3evm0HhuRY2f6grW44GrgR9AtD7MTMlzm2DbfOJRGSTiGSLSHZ1tW8N3yqtb+PFAyXcmJHEhPBgq8tRakCrZ8SSmRrJfW+fpL610+py1ABcCfr+5t/3vQvjSptPZIx51BiTYYzJiImJGcxbR71Hd53CGPjuSr2aV6ODiPDza2ZS39bFv7+Wa3U5agCuBH0J0HvmTiJQdgltVD+qGtt5NquY6xcmkjg+xOpylHLZrPhx3LJyMi8eKGGXTqLyaK4EfRaQLiKpIhII3ARs6dNmC/BV5+ibJUCDMabczbV6pYd3FdJtd/C9K9KsLkWpQbvjqnQmx4zlnpeO0tzRbXU56iIGDHpjTDdwO/AmkAs8b4w5JiK3iMgtzmZbgUKgAHgMuPX8+0XkWWAPME1ESkTkW27+GUatioZ2/raviOsXJpISPdbqcpQatOAAG7+5fi5lDW382ys5OrbeQ7k0x94Ys5WeMO+97eFe3xvgtou8d+NQCvRmf9hxEofD8P1V6VaXotQly0iJ5Aerp3LvthPMT4rgazqr2+PozFiLlNS18lxWMTcsTiIpUvvm1eh2+5VTWD0jln9/7TjZZ3Q0tafRoLfI77cXICLccdUUq0tRasj8/IR7b5xPUmQIm/66nxOVTVaXpHrRoLfAmXMt/P1ACTdnJjNx3Biry1HKLcKDA/jz1xfj7yfc/Ng+CqqarS5JOWnQW+D+7ScJsAm3XqkjbZR3SYkeyzPfWQLAzY/t5XhZo8UVKdCgH3EFVU28cqiUry5NITZMZ8Eq7zMlNpRnvnMZfiJc/9BuXj+qI62tpkE/wn739klCAmx8V9e0UV5s6oQwttx+OdMnhvG9pw/w6625tHfZrS7LZ2nQj6BjZQ3880g537g8VVeoVF4vNjyYZ7+zhJsvS+aRdwtZf/97fHhaR+RYQYN+hBhj+NU/cxkfEsB39Gpe+YjgABv/+fk5/O1bl9Fpd3DDI3u49en9FFbrjdqRpEE/QrbnVrH7VA13rZ7KuDEBVpej1Ihanh7Nm3et4M5V6ezMr+bTv3uXuzYfJKe0werSfIJ44pTljIwMk52dbXUZbtNld7DmvncBePOuFQTY9N9X5buqmzp4aOcpnss6S0unnctSI/nW8lRWzZiAza+/hXCVK0RkvzEmo799mjgj4Om9RRRWt3DPuhka8srnxYQF8b+vmcmee1bxs/UzKKlrY9Nf97Pqtzt57N1Calt0fXt30yv6YVbV2M6q3+5iXlIEf/1WJiJ6xaJUb912B28cq+DPH5xhf1EdgTY/1s2J4+bMZDJTI/XvjIs+6YrepUXN1KX75WvH6bA7+PfPzdY/sEr1w9/mx9Vz47l6bjz5FU08s6+Ilw6W8uqhMqbEhnJzZjLXL0xkXIje27pUekU/jHbmV/H1P2fxg9VTuXO1rlCplKtaO7t57Ug5T+87y+HieoL8/fjs3Il86bJkFiaP14umfnzSFb0G/TBp7uhm3f3vEmDz4/U7P0WQv83qkpQalY6VNfDMvrO8crCUlk476bGhXLcwkc8tiNe1onrRoLfAj144zIsHSnhu01IyUyOtLkepUa+lo5sth8v4+/4S9hfVIQJLJ0fx+QUJrJkdR3iwb3ftaNCPsH8cLuOOZw9y+5VT+NGaaVaXo5TXOXOuhZcPlvLKoVKKalrx9xOWTI5i9YxYVs2Y4JPPeNCgH0Elda2su/890mJCeeGWpTqcUqlhZIzhwNl63jpewdvHKzlV3QLA9LgwrpgWy8qpMSyaNJ5Af+//e6hBP0Ka2rv44sN7KK1r47XvL2dSlD4HVqmRVFjdzPbcKt7OrWR/UR3dDsPYQBvLpkSzcmoMK6fGeO3Vvg6vHAHddge3P3OQk1XN/OUbizXklbLA5JhQJseE8p0Vk2lq72L3qRp2nahmV341245X9rSJHsunZ07gmnnxzIoP94kRPHpF7wYOh+Gel4+yOauYX183h42ZyVaXpJTqxRhD4bkWduVXsyO/ij2nauh2GCbHjOWaufFcOz+etJhQq8scEu26GUZddgc/euEwrx4q05uvSo0StS2dvJFTwZbDpew7XYsxsCA5go2ZyVw9dyIhgaOvs0ODfpg0d3Rz57MH2Z5XxY/XTOPWK9J84tdApbxJZWM7Ww6VsTnrLKeqWwgL8ufa+fFszExmdsI4q8tzmQb9MDhUXM+dmw9SXNvKLzfM5stLJlldklJqCIwxZBfV8ey+s/zzaDkd3Q7mJIzjpswkrp0XT5iHj9PXoHejhrYuHtl1ikffLWRCeDC/u3G+TohSyss0tHbx8sESNmcVk1fRREigjWvmxnNTZhLzkyI88jd3DXo3qGpq54XsEh59t5CGti6uW5DAz6+ZpQstKeXFjDEcKq7n2Q/P8o/D5bR12ZkeF8ZNi5NYN2ciE8KDrS7xI0MOehFZC9wP2IA/GWP+q89+ce5fD7QCXzfGHHDlvf3xhKA3xlBU08ruUzVsz61k54lq7A7DFdNi+PGaacyKHz19d0qpoWtq72LL4TKe/fAsOaWNAMyKD2d5ejQZkyJZkBxBtIXPgh5S0IuIDTgBfBooAbKAjcaY473arAfuoCfoLwPuN8Zc5sp7++OuoDfGYHcY7MZgDB99b7cbmju6aensprm9m+aObupbu6hobKesvo2Tlc2cqGyixvkAhInjgvncggSuX5jIlNjRPQRLKTV0eRWNvJNXxY68Kg4V19Nl78nR6NBApk4II3H8GGLDgokNDyI2LIiYsCDGBvkzJsDGmAAbwYE2gvz9sIlg8xO3dAUNdcJUJlBgjCl0fthmYAPQO6w3AE+Znn819opIhIhMBFJceK/bzP/lW7R0dONwhvqlCAvyJy02lFUzYpmTGMGytCgmR4/1yD45pZQ1pseFMz0unFuvmEJ7l52jpQ0cLq4nv6KJE1XN7Myv5lxzB67GkAjYRIgNC2L3T1e5vV5Xgj4BKO71uoSeq/aB2iS4+F4ARGQTsMn5sllE8l2ozV2igXPnX+SM4IE93MfOiwL0nFyMnpcLDfqcnALknks+3kWH/rkS9P1dyvb9d+pibVx5b89GYx4FHnWhHrcTkeyL/crjy/S8XEjPSf/0vFzIk86JK0FfAiT1ep0IlLnYJtCF9yqllBpGrqzdmQWki0iqiAQCNwFb+rTZAnxVeiwBGowx5S6+Vyml1DAa8IreGNMtIrcDb9IzRPIJY8wxEbnFuf9hYCs9I24K6Ble+Y1Peu+w/CRDY0mX0Sig5+VCek76p+flQh5zTjxywpRSSin38f7HriillI/ToFdKKS/n00EvIv9XRPJE5IiIvCwiEb32/VRECkQkX0TWWFnnSBKRL4rIMRFxiEhGn30+eU7OE5G1zp+9QETutroeq4jIEyJSJSI5vbZFisg2ETnp/O94K2scaSKSJCI7RCTX+ffnTud2jzgvPh30wDZgtjFmLj1LNfwUQERm0jNCaBawFvijczkHX5ADXAe823ujj5+T80uBPAisA2YCG53nxBf9hZ4/A73dDWw3xqQD252vfUk38ENjzAxgCXCb88+HR5wXnw56Y8xbxphu58u99Izzh55lGjYbYzqMMafpGU2UaUWNI80Yk2uM6W9Wss+eE6ePlgIxxnQC55fz8DnGmHeB2j6bNwBPOr9/EvjciBZlMWNM+fmFHI0xTUAuPSsDeMR58emg7+ObwOvO7y+2pIMv8/Vz4us//0AmOOfO4PxvrMX1WEZEUoAFwD485LyMvgcjDpKIvA3E9bPrZ8aYV51tfkbPr15Pn39bP+29ZhyqK+ekv7f1s81rzokLfP3nVy4QkVDgReAuY0yjpyyG6PVBb4xZ/Un7ReRrwNXAKvM/kwpcWfZh1BronFyEV58TF/j6zz+QShGZaIwpd65cW2V1QSNNRALoCfmnjTEvOTd7xHnx6a4b50NR/gW41hjT2mvXFuAmEQkSkVQgHfjQiho9iK+fE13O45NtAb7m/P5rwMV+M/RKzocvPQ7kGmPu7bXLI86LT8+MFZECIAiocW7aa4y5xbnvZ/T023fT82vY6/1/incRkc8DvwdigHrgkDFmjXOfT56T85wP2LmP/1nO41cWl2QJEXkWuIKeZXgrgZ8DrwDPA8nAWeCLxpi+N2y9logsB94DjgIO5+Z76Omnt/y8+HTQK6WUL/DprhullPIFGvRKKeXlNOiVUsrLadArpZSX06BXSikvp0GvlFJeToNeKaW83P8Hc97cRi8quEsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"kdeplot(timeseries)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Auto-assigning NUTS sampler...\n",
"Initializing NUTS using jitter+adapt_diag...\n",
"Multiprocess sampling (2 chains in 2 jobs)\n",
"NUTS: [center, tau, rho]\n",
"Sampling 2 chains, 0 divergences: 100%|██████████| 2000/2000 [00:02<00:00, 987.20draws/s] \n"
]
}
],
"source": [
"with pm.Model() as ar1:\n",
" rho_ = pm.Uniform('rho',-1,1) #we assume process is stationary, so -1<k_<1\n",
" tau_ = pm.Gamma('tau',mu=1,sd=1)\n",
" center = pm.Normal('center', mu=timeseries.mean(), sigma=5) #set the prior for the true mean initially centred on the population mean\n",
" likelihood = pm.AR1('likelihood', k=rho_, tau_e=tau_, observed=timeseries-center)\n",
" trace = pm.sample(target_accept=0.9)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Analysis:\n"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([-13.37825194, -9.30304861])"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pm.stats.hpd(trace['center'])"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x140f55990>]"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO19d7wdRdn/9zm3pId0CKGEEhJ6CyAgSA9Fpb8i+Co21Pf1FQsqiP7EjhXFhoAir6+CKKJo6NJ7AiQkIQQIBBJSCOnJTe695575/bFnzpmdnZmd2XLOuZv5fj7J3bM7OzM7O/PsM08lxhg8PDw8PIqJUrM74OHh4eGRHzyR9/Dw8CgwPJH38PDwKDA8kffw8PAoMDyR9/Dw8Cgw2pvdARFjxoxhEydObHY3PDw8PPoVnnnmmbcZY2NV11qKyE+cOBEzZ85sdjc8PDw8+hWI6HXdNS+u8fDw8CgwPJH38PDwKDA8kffw8PAoMDyR9/Dw8CgwPJH38PDwKDA8kffw8PAoMDyR9/Dw8CgwPJH38ADAGMMtMxeju9zX7K5khnVdvXh84dvN7oZHk+GJvIcHgAcWvIUv/fV5/PCuBc3uSmZ4/3VP4vzrnkJvX6XZXfFoIjyR9/AA0FMOkue8vrorszrXdfVi/ZbezOpzxQvL1gMAKltBYqDn3liDNZt6mt2NloQn8v0U85auw63PLGl2NwqDoQOCCB+busuZ1bn/N+/Bflfck1l9SbEV0Hic+avHcd61Tza7Gy2Jlopd42GP065+FABw9sE7NLknxcDAjoDf2dRTHJk8x9ZA5AFgwYoNze5CS8Jz8h4eACpVQtiVISffKtgaxDUeengi7+EBoK9K5Tf3FpCTb3YHcgbzHzEjPJH38ECxud0iPxtQ34V5qOGJvIcH6pw8UZM7khEqAuUrOI0v/EcsLXIj8kS0PxE9QURziOifRDQ8r7Y8+idaaZvd10J9yQJbRKeuYj1aBAV7dZkjT07+egCXMsb2BXAbgC/m2JZHP8SuX7kD5137RLO7AaDO+RKKwcpvFqyEis7pFv350iJPIj8ZwMPV43sBnJ1jWx79EIwBT766utndAFAX1xQF4s6kWE8WhafxZuRJ5OcCeG/1+FwAO6oKEdFFRDSTiGauXLkyx+54tBIqLUZUOTdYFJm8SNmLzumywn/G0iEVkSei+4horuLf6QA+AuC/iegZAMMAKH2OGWPXMsamMsamjh2rTDbuUUCs3dw8d38VeHiXAtL4wnO6LcYvtBxSebwyxk6IKXISABDRHgBOS9OWR7GwoYkxXVQomuJVfJxWUnDngax2Kr984BXsv8MIvHPSmEzqaxXkaV0zrvq3BOCrAK7Jqy2P/odW475qiteCyGtEEQY/+sQfZuI7019oTodyBMsoyOYP716AD/z2Kevyl/1tDj5/yyzr8v+evwJzlqxL0rVUyFMm/34iegnAiwCWArghx7ZaGo8vfLvlZNDNRqtxl0VTvIY5+eDv3fNW4LpHXmtOh3JEnjL5F5aux9w31YT5pqffwN+efdO6ro/eOBPv+cWjWXXNGrkRecbYzxhje1T/XcpabVU3CHfNXY7zr3sK//fU6w1p760NW3DPvOUNaatI4OKaVufj731hBf49f0VsOXGxFV3xmuf3+dSrH8G7f54tYV7b1YN/zLL/OKSF93jNGa++vREA8ObazQ1p7/zrnsJFf3gGPeXWThTRamSnttNqcSr/8f+diY/eODO2HNuKTCgb8RG77bnswnp/5uZZuPjmWVj09qbM6jTBE/mcsaU3ILaDOtoa0t5r1YnT6mZlrcZcFlnxWnRRYV5E/u2N3bXj22ctzazeZVWGr7tBjJgn8inBGDPKl7dUoxoObBCR5xO+9WlWa3Ww7vHauujqsQ+DLL//cpFTAOY0lX50dz0VZJbrt9G6fU/kU+KIK+/H1G/fp71eI/LtjRlqvriLLofNGo1QvM5ftj4Vse3ts++juJOrMIZv/qt4VjUceb06kbAPyGH9Nmq37Yl8SixbtwWrDLklG83Jc7T6Dr3VvkGcfuZlQvnqyo045WeP4Pt3vRg6v2pjt7Vs1sV2Qbauuf/Ft6zv7W/Ii6EZPaSzdjygPUNOvrpfbNQa8EQ+Z9Rk8p2NJvL6GVSpMLywdH0DexNFi9F49FWC95TXTnrlhkC+O3tx2Bzv8CvvxzE/etCqDheiwDTHRUReRH7b4QNrxzw9pA1ufHwRnl+yVnvdi2sKBi5HbS81dqhNyrZvT5+PU69+pGHafRVajpPnYQ1yWoC1x5Xqd7GCchkykeuvMNZy450l8no28eNhuxNfsX4Lvn77PLz3F4/FlvWcfEHARTkq+dv6Lb248IansXzdlszbNYlruM1+udLayri75i7DSw1KzlwLUJYTL88XdJraXThWOXZNkd1U8no00eKqrWT35k67Ot6m/sXljU047ol8jmCM4bk31laPo9f/8dybeHDBSvz8/pdTtaGCiSBw7rGZ61786D3+yttYJZircXzy/57FSVc9HDmfB/JWvPLnTbNTcBLXSLFrTLd2l/vw4vLmiu/SIC9xjTglbN/b24p5rINXvBYAS9bUHaBUE5FPolKKlc8YMPfNdbju4VeluuMnUDN5O7F751//FC643j5mSB7I/YOXQf1u3HjYGcp06+W3zcXJP30Eb23IfkfZCOT16sTx3m6bQTm1kj88kc8RcYSWTyLLnaC6DgDv/vmj+M4d88NtW0hi4vr3xMJVeGNVV/LOGSA3/fJbG3Npxxacq8qLK+S1phEHucnk68cVxoxc49OvBYlbxGxSKnzsxhk46gf3O/SiMcjrnYm7u7YclDVeJl8AyAtNRs2TPhUn7y6uqd9rvv7+657E0T98IEm3nNEqTki5Efnau05eRzqZvL4sJ2ZxO8r75r+FxauThed47JW3sT6n8NJ56Rta3QzZFp7I5wh5ocnIIhuRbh5mQeTzhMxZNjvCLx+LvMakuTJ58y6AE/n2NrvOLVzptutas6kHF1z/FP77j8863WeLLIix+KF4fdUm/OahhVgnJLZp9TAhJqRKGuJhhkhoTRMxrUze5XyoTBMnrty/ZifQ5t3JnZNvlLhG8ng1PVa5OjltRRIfuP4pPHHZ8dZ94TFa8rKUyuKVievzXT98EAAwcfTgTNtoFjwnnyPisvNwgpJOJq+efTbWIi01cVtEXpPXFr0mk08jrnHoXFQno7+XO4LZ1s4du1zRCHv2pFCtz8GddR64lZaKKzyRzxWChYNSXBP8TSeTV5/vb7Frmk7jq+PlYkq5uacPDyywCxfQaDt1F06eP7NtF8uOX8K8RXGZEPmM6p287TD7Nr3itXXw0Esr8aDlYhYhrgUVx52FTN6mbR3ymGTrNvfi3Gsex+LVZquciLim6VQ+gAsx/urf5+LDN8zAK2/FiyHqnHzyB3VSvBpk8owxLF7dhd6qm2+NyOfEr+ZNzLIR10QrCX3MLBuZOnEkAGDCiHiTS3G8f/XgK/jbs9nFrBexVRL5e19YgYmXTq/FXo/Dh373NC68YYZzO2HrGv31fGTy8ZMyD27/zjnLMGPRGnzk9+bxiiheW0Ymb38PV0Bu2GIRArgmk0+OpK+rIoXDfnPtZhz1gwfwo3uCULrcszMvURWvP8vqL7j+SfzpqTeCejNRvEbPiRFDbZtIaub615lL8MCClQ5322OrJPL/nB0kADAFEcoCIhGVJ9GLy9fX4tqkWvg6mXyTnaHi7N5bjZPn/UmSPKRRyb+TEhAm3fvyiuDdzF4czH/OyT/3xhon5scWeSQteeyVVfjKbXOC+nPaKriEduZI8yHOaxZtldY1XNFpmhxrNvWgrY0wfGBH4nZ0dvKbe/pw8k8fEfqTg0zewhkqDzlx0hpbRFrjFs7Xpd4MTCjd7OTDDIZ464r1gWfruGFBlEUuluBJqR9+aSV2GTMkeUcluMr8XZEFkVeLawRO3rEJm3nEpOO8eIWtkpPnRNWUv+HAb92LA75xT6p2QgtNOC9HHkzr8apCy4c1kH5nzQ3fPntprF4g3B93xSuHTc/rJpTJkSp2jXCCmzQOroa/rvsIsNDftLj1mSVY19Wbe2rFbOzko+fEuSCOyRW3z8OfZ7yhq0n43639NMyeCVsnka9S1ThCmHby6Ewov/GveeFyqdpQ3523M9Qbq7pwxPf+7URIw23LMvls8ZmbnsO7fx4fEbDen+Cv7p3PWrw2+qwOsua6x2v23s3KstKx+Jt3QRZHZEmKX16xAV/4y2x87pZZDcgxmw8nrxPX/P7xRfjyrXPUPWlBcU0qIk9E5xLRPCKqENFU6dplRPQKES0gomnpupktOOect1mbbOHAwbfFHGkiIMpWExx2VSZv9/bZb2Lpui34vydfD9doWWWkmDTD497Nxu4yfvnAK7XMWyqIHou20LV7xi8fw1E/eCBctvrXZdeUTv/iUFacCxUWupkTXV2o6SxWBd8tLF+3ReDk81lvmXDyinOJFK8s/NdcNixSy0u3k5aTnwvgLACheLBEtBeA8wDsDeBkAL8iosamRjKAb4vEyfHi8vW44vZ5iQn/F26Zjd2/ckfonGyrrEMqIq/5kOTNyU+symxlL8akZnjy9I4bkjueX4Yf3r0AP7n3pci1JO+wRoYSdL/cx/DsG2tw55xl+vozMJdNHLtG+s0Z1LKGU82CaPLnZMg/jLO4U5BzM9wyYzGeeX11bB2qoe2tMON1ZT01cY3bMzPGWlMmzxibzxhboLh0OoCbGWPdjLHXALwC4NA0bWUJoqi45gPXP43fP74osTffrc8uiTiJ2JrZujqXhMCUh5YmlIZqY+7nOS8fW7jKuo7Fq7sw8dLpeHDBW5HxKEmKiTiCNrKaf1PlKp/m4+XGLQd/y5UKzvrV4/iUITZLvd7slexxZWUTSn7cKymlSLqeBvU8pqxmBJDXxlms9ulFYYL+pVufx9m/fiK+DpXi1ZGTL/dVjHq+SJvScUuKawyYAGCx8HtJ9VwERHQREc0kopkrV+ZjJyqjbl1TP1fntLIb6tDCMpRLJ65R32sz2UyLOa5L/F5ZiSxWKT/XvKVBftM/PvUG5BGRRz2OINTel+JaEmuLukw+em932RyCV8cRq+pPM702dVvY49dbDB1WJKIPRJmLLGlwjZNnycxSXSC+s7czDLnguiwP/Oa9uLXq0LRhSxn7XnE3Hnvlbat7K4w1T/FKRPcR0VzFv9NNtynOKYeMMXYtY2wqY2zq2LFjbfsdwfJ1W7DRchHwwZTzYAbXstteqj4iKqRJw6dT7qa1rslCKS0TkW0GBdz3jEWr8eSr5i10bCz+lH2L1lfdZivu7eqOEvmr7n0Jc94MPlp280X/UbLFOdfEc6QcYaKuviZz8roh7+2rKDN3mVASdstZi2vkteQqplRBd5+L/m6DQH+6evqwYUsZP70vKk7kkEWteYlrYu3kGWMnJKh3CYAdhd87AFiaoB5rvON7/8auY4fg/i8cE1u2xslXRKLIr5HxxbjBTqaXmeJVOE4rk7dNeGKCXAc32Vvb1Ysf3h2W8rlyMWn6boLq3k09UebhZ/+up2yUiaUKcZx8IJMNLq7c0I3HF76N0w9Qbn6NmHT5HdhvhxH48slT6nUjnDSEP6Ou3/IQfO7Ps/Cv5/X6BhVKgkyet5f0rby+ahNmLFqDcw7eQdm/vGLXAEB7qYQei/er2+25GCK0quJVh9sBnEdEA4hoFwCTADydU1s1vLrSzlOPD2Zf6EvKOXnCg4J78QtLk+e+tOUwGAtMzlzjdAf3qj8kVtp9w7KLu1/3XRJP8/y2qmsy5Pkdv3CT991UnerWrpiMSVYRP6t/deEbxD5/5PczcPHNs7CmmgTetg0gMPt75vU1kXkR4uy5dY2lCaUrgQfq7zMLTv49P38Ul/xldu23XJtunbmYburmTMmSQlqFtoi2KrTfoopXIjqTiJYAOBzAdCK6GwAYY/MA3ALgBQB3Afhvxph5pTQQKnENq6/CkHPSqVc/Ahf8x2/qW+pwgDI9GANOvOphHP/jh5zaMsFqYVlyw6rFoiXCwnk5b6sLxxWvEzDdm0AmX6s3em/cWPY6EBMtJy8cv7k2yL4kPofNbkFXX4Wx0ImauEbT70wDlbH0YQ3WVwkofzfy+9WFD7HhwOv3qfvYVqMV5vt1z2gi3DJTlsYp0oS01jW3McZ2YIwNYIxtyxibJlz7DmNsN8bYZMbYnem7mh1UYQ1qxyxq6eECni8TsJeRy4vq+kdetY4fEhbX2LWnuleGOGdVirM0yk0V5K1qnDio9k1WrKI023fVrXHVlTMQ16jmovhszkReYjBU9cv9rnu8OjVlbJ9BCFCWsmL+sTUSeeG8E5HXnOe0IO7Dp93ZWj5y4AzVv8Q1LY3aixNewKbqlpyBZZa0N0SADS9bvNbVU8a3p88P7QiMbWhENDaTi5vWxXHqfRWG2YvXhpJ6axl5Y1/1V13t5M2cvPledX1RArJmUw+29PbFLnAbE9ha7BoLcQ2HWNLGgkfVXlA3C32oeXd1O5Qs+HheNctQ8crHWR4r0cJLfH+y5ZcJuvnUpqAVKsRZEKnmPpOOW1Jc0+pgjKG3r4Jl68LJh+vywuDv2q6e0PU0nLyI8DbSxMkLx9UfGy1lfDoCZMPN9pQrOOeaJ7Cr5MQFAExYHxXGcPovHwsl9U4jErEqG2vdw4mm/t4ki0Zs9cBv3Yvzr3vSgpO3IPK1rYf6uiliKQD0ulpgSR/8MANg5tgz4eTBP5rpFa8cFU2/u8vqQGIuux/dfLZl+OJEUnFjGohr8qHyhY5Cec1Dr+L7d70IAPjxufvj7Kp2viQ5Q4mml4zZv9hYWHLWKsWRtVxUU8yGe/r0n57DZk1YgLj8tEkIhGkhyBM8dlEYrtUybpmrULYXEMS6pcuzb8SHo+6zIMAxNB6rN/Vg+2qiCRVRdOfk64iaUKrFHvV701P5mgMUmJODkApEwXvRiWt6QkQ+uPYf1zyB5xavsW4jXlxjRhKTX53COGsUmpMXI8XNEDzhZNtXUfbJUN+ipYXJVllESJbu6B2oEwnZTBodgZfvV30wtATC0K5pDFyta8zt1LnI6x5+1VhPrb5Q3Va3KO/VllHMNREX3iAYn1UrFMfAmcgrOHcOTnT179CpKSV43es3l5Xz5665y7B+i11sIT5i9Q9HGKLs/cXlgQf004tW1wKM2fBsOgbEVvEax1TF6iNy5OQLTeTFKHLiGMuxa+ShzUpcEw41bKTyNfDF0V2uYI+v3om5VYcb7a3iYra0y1fXE75BZXKna9e6DcMYuMrka/cpXpVIvL5zx3y7igTITdtstWXo7KZ1M+sVIckKry4kX3a2rhHnnnRN+AjmBd71dZt7a85+nBC+unIjPvl/z+JLf3neqq66yXOUk69UWIiT/9fzy/Bd6Z132NpBKmDL8Ol3ReG/oWuSMtzL5BMgFPRfjB5SHc2f3PsSKpXw4DLG0JaVtMaak69DVOD0lCv4wxOvR28QsFTSN6jqtEEk7KzIySsmcNykVl4zcvKSuCbmCWzFX7Yw7YLkvkQ9LqMNHvvDBxP3SbT04nD1ijZx8qqPiHxPWohjxok7l53z6KDi3O3qKeP1VWqLMk5na8lHhKHoYyyiYL1W2r3Z0Hjds/N7k1rXxNVfu47+F7umJRAO+l8/L36c31y7OWJulhUq4ZWmLRcJC+uAs371uLIJ1wUrc55xnHwSLjBpImrldeFNBe8wzBW5wrQLivutGoulUjTEuMxQIZFhTVxTv56lTF4n2663n34VqPrOiTG3kuloq5OfC2+YgXdJH0aOGievSDjeV2HO5qXq/uYsrlFQFvEMY9lJEGQUmsiHxDXCeVGxOqC9lIgoLFnTheN//KCxjGmh6SBzzS5KsHB7bs+0pTe8UEIyeamu/7jmiVp+TReYOfnw73iZfPB33tL1OPLK+0Nx7dOKIeLj9qg5YxP4LUvXqnde4uPXLVPqNbvbyet7VdNZaKrkt3aX+zA9gber2AYQJYD8WTqELbPoXyKjxskrxEwVBScvw8b+fJFuFyER3l8/uFBZLl4mb26/ZZOGtDpCoUJFTl54cRUWnuy2tPGPT72BhXFhFEI7BAMnLxzLk8WFVttGvVQhysnrF6kczvU701/ALTMXV/ugb8NEPJes2YxFggNY/KII/i6rcsxiCAV559FXYdjvirtrfQQCG/hdL5uOJ6qhko3iH+l35ENs8ZJ4kRmL1sSWr9mYC+dck0qbdnV1SyI1Q8HP/viel/Dff9KHT1ZhwfINkZ2V7EfAn0Xk5DlUu0ZOpPm1mwWDitWbevD0otXobE9Hyi68YYbyvGxpx631ZCTR28gShP4Wu6YlUA6Ja7hcsC8UgKzCWCLOzAZxZoj1vgn3pN95AgD++OTreHG5fdwdmZN36dN1j7yGL/01UKKJj3nYLqO0darw4d/XF5qrdc2QAe3CtXDZjd1lrN9Sxjdur6ddnLV4LSoMuOahKGcWlVUn//CqoCLYSgWyMGlsvGpFmJT+/Pn6mNoRLm7XYcK0nz6MI6+8PzRGfZWwiWNvlfPuVBB5lf5Hlsn/4K56cLuL/veZ2Iim8th+8S+z8V9/fMZ4D0fdGUr/0m+fvVTrDFXzSVCKa8L0ySteE0DliXjj44siYhzxBWUZt8NWRm4Ss7j0Riz71GurcfJP7ePuyGn0dJYdcbHVRew5frjUP/PTiERFNV7zl9U/WvJlkcjLY3jetU8CCO/gTHl+I5yvdN1VYclY+Mn5GGqJaI3Trp9yTSwj3rtsbVg/IJqYqpys6vqD5FQnFBZD+LGltyKIaxREvubVyvD4wrfBWD3OuoqQvvxWYDLp4t36l2eW4I45y63K2ni8fuam52LFacr7JfrgxTUJEFK8Vv9u7pFkzxWJm7FcSzYvRKzK6PFqkH+7iWvsy8qQRTLiT3EBHfbdf8f0Ibp7Avg4m/ug81wEAm7plJ89Uk+xJ10fOqCeXVImwvzjIJrDyfGLXBS3Mr21sQQS6+fPKYqnRLkxb/+eF5Zj4qXT8fqqTc46FrG4bEYqhhxQKXQZA/b5+t345+zk0cFDNv4V8dn7akxDu8KMjc/Dm55ejPOvewp3zFle43BVcu9xwwYm7qMNbE0o12vyCfMPpertXfznWZh46fTadW8nnxI6Lz/Gwov2pJ8+HAo1rIM5uly0LdMSDXHykWw9LlQ+vsgzr6u9AE1Bn0TufW2Xe3Lsu+Yux65fucM66JqqPy9VnVxertqTy+Mi61lUEOWrbZLFhgj5jExfTXqTV1dujF6X6uREXh8FMsCfZwQ6hPnL1mfqoMXnWIXJIs16GdsEPNo2NDqd3r76h8Ukrpm9ONCxrN/SG7GuETGoMz51dFdPXy1kuCnxuwolA5EWoUsazxRMBIeYatTbyWcAk7JJnDwmImZr4VDnlIRzlobyeaZKe+b1NTj7148rr0W4U1Fc47AVDtVR/XvT04GiTHT4sb1Xe10qUAkRK/Xd4oegludX4UXJIo8rE3W9SO24Hz+Eq6Tk4nKo3+4qoQnJ2SVfDaCe7m9wZ3sCTl5fXmR4VMYJmcxAjaipt88srrml+mFbtSkggKOGdEZk8iJsrQ4fXxik4Xv0Zbt0fBy2AcrWxDA//HZVfxljgbjGc/J6yBN67pvratugWhnpL4cYQCkOJ//0Yaty3HHF1tpF5TjC8bdn37RqU65HBTkQW+jeCCdfP+5OSuSrdfAFO2JwR+w93JRO9074MpCvioTEipOXZPKhD3KMyEwXC4bjFw+8ErlffDfH/fghvL5qU0g3JC5vXj/npocMaHM2CzVy8tWLfRWmlvVnwGjoZPIikf/zzMWRHKjfnh6IllZVE6YMG9Be/yArlbJ2hJHfytu2pae2cazWGdaW2H67wjOLi6+8TN4A+d3f+8IKRRm1uEZlXaNDrMlkFXWnjXA7OmRlXRP3GKJyUoaJcCUN2coJ2+qNPdXf8eAhlnXmfap2gDAh0Tq2CGwUb0e1c4rrp6uZKwOLlHlwwcqQF6uKlnAi315y9+UwFRc/bKrdaRacvE4mHxD5+m85sQwHn3MMUeuaNHCtwdbjdX1M1Nha5jkFxd1cDXPuxTUGxJm8AQInHxFLZB/Do6/CsGpjdyiGhtG6xsBFuiDuzoEdevll5OMnrP20nDyXPbuY5KneExBOKyfChsirFljdJl1/v1ybq0eoJK2p1RHnxcr7VmHRj4RFq7UjmXgwod5yyNIsuqtJirBMvj5/evvc7NdEMYbqvdqKOLp6+tDVU7ZeX/zDYiuuWb3JjpNX7Qx4ekmveDVAHP8/PPk6rr7/FW0hlezZhUPg20uTF11fheHrt8/DkjWiSaC+jbWbhVyeOcrkbeS0qt86pZKyjZB9dgA+Ui65QnWvRGetoFMghupUfKvmL12Pz9z0HLp7w7JpXe5cVd/iiH5fJUqkGdx0PEk+LBwy8RDDGqjENVmYEYs1yJy8i/MYA6sR3B/d/RLmLAkH7LMli1fd9xL2veIe6w/YqCEDANgT3ts1lkgbu8v4/l0v1kQyqtAFXdVE8XmJawoRT158cV/7+1xlGV2cdtcFdMH1T2HRlacZy5QrUS7N1ILozJEmH6atJ6XymkRvxKrkpCv2/Qn+JmFQTNz0U6+uwrf+9ULoui47kAjVx7ynr4LbZy/FYMFKgzEWGqtIqN4YmX2kXcaU4iaR+JmYBrk/NhCLlwgQbUrEsAaqwGdZ8Bkhs+C+MJF3q6c+Nk+8ugrv+cWjifvUV2HaHb2MbQa14+2N3dbWNTq8tGIjXlqxEcMGBqS2XUnkq5y8j12jh80WrG5dI51HdunJav2psIh9re02MU1fthjEKm+u3YzvTH9Be5337+UVG7CuqzfU35+rdkYOSJK70jRc77v2yYieIBRxVHOvaZck+1QYdz0KE0kTKhpO3taL1cU4oFa/UFwef36tTxbXaHa7SSB+O8SPWU85qgMzmTUymC1o3KNz2j0cV5DW1nHKLx834VTZ3W/qdjPrdEUhOHkbMDD8+sGFuO6RcBjSSsVdRNJTrkQsKER0lytaOWgc0ohrzrs2mheWZzn67M3PGbMc8YV34lUPY7exQ/Cjc/dP1Idw92633cAAACAASURBVIMfrpz86k09zkTtrfXd2NLbh4Edbdp7Tbsk8ZZ/Pb8MY4Z21q8Zyqp+R9rVyORF2lbTNSgjftobB9TqN8jkRRNKFWdtk+kqDjo7+XKlEhmvd3wv6mBX57iZUe7uGtPHFpyrLlE1M1XK+rioUiX+6XG0+HFFIYi8zfxnTB1cKIlS68EFbxmvb+wuR77Ytk2kWV+LV0fFKj+6ZwG+OG0KemIWgzgGC1duSszNibfVxDWOdcxeshZjqjJRuS4d7nlhBT7y+xn408ffoe27iVCKH1dZFCTeNvHS6Ti3mkaS46r7wnbxkboDmWCkzrJgkN/V04dpVz2MDYpsSUnmqDgGMmERrWtUMvksCKdJJi/D5JsS1xNXHw7bceTLt0RJ9qF6qDh5/gwtqXglonOJaB4RVYhoqnB+NBE9QEQbiegX6btphpW4xnDNVURiMkUEAicWWYtuu038wG+jJmVp4nvfMnOJVTkbCyVX1GXybpP35RUbEin/Hq9FlNRw8iadhOF55b785Rm7MRXrjnDyYBFiumDFhkgceoBb57iKa+rl5TSP4jiIRJef5orANAjJ5AXOpads9yR8xvColjpwLvjuzx5t1S91rKLoufYaJ29nXWMLFZHnHuWtaic/F8BZAGQvoS0AvgbgkpT1W8Fm/LUWFwm2wiZTRADY0F2OKFFuenox3ljV5dQORxYy0rgJFLWTT9aOOJRJrTRWrO/WxoexytepadYk283RqEltXcPsE4FUmD7uz5OvrsKsxfHJxkP1SWEG6n0Kjm0Zg6VrN2vndNgZqn4cWNfY9/XKO9WhfTl6yhXsOX44Jm83zKo+VdtKT1pO5EsBk5JV4MJ+x8kzxuYzxhYozm9ijD2KgNjnDjvFq567cyXyceU3dZeVWnRZH2ALeRK+vGKDcx2xykHpmbizUBokta7p7asYPF7jK9PdO2W74firIxcOAOdfp3bYsUWFsVAETaCqeLWUzclzVBzP8659Emf88rHIPaYpGopLJHD5Lh/279/1Io648n4c/cMHsLG7HBGb6Ozky5VKiGDK4ag5bLvS21eBIjqCFqp6VXowvhPPXFyjWAzcD6WwzlBEdBERzSSimStXxgcGU8FKJq85b+KSdIjjwDZuiXLyQGDDnwQy0TrnmvQEONJGRhZGpgTStugpK7g9h8p0RP6JV1fhkr/MzuxZrftTAW6esTh0jrGAi+5sK2G74eFIigfsOCJ8vySTt4uAqn9G8fG/KegfXJgdMUPSPl+/G+dKTIFYVSjkRCV8TRWJ0iVRSU+5gjaHRN0qZk+1/jnHTUQoVxh++cBCpb7EFSq6kDQ2lHWbcQWI6D4imqv4d3oWHWCMXcsYm8oYmzp27NiElVi1oznvbtESJ8Pv6atYx7xI0l6SSeEqrkmKkLgmISffY3CYsakr7nX2ONpqp4VqfpX7Kij3VdDeRhExkrwLZJJIUaXjMCmLZYh1iQ57aebAbElkpLOuYQgzFCoC65JysFxhUHwntFA9ompHVSPywjk5QXgSqOhC3uKaWOsaxtgJubScIdIoXhlTZ8cxIW6brbKTTwOZSLS3EWDJVNj2IuAW01H67nIffnh3RHrnPHl7ypUIwZmpCZGswuLVZt1Hd2+DiXyFYdSQzpDre29fBeUKQ3uJIkRb5vYqlfD8VY3mbx99DV+cNrn22/QqtfqpDHc4Wk6ehR27nnh1leZ++744rTVFtSZOPnE7Gqg4+ZritajimixgMx24V5mMQN7p1l4cJ19h2X6V5cXnUje/M+6WZPFRwjhXEiPVlKWO9ahc33kcfJu6Lv2bOcm4S3arLFCuVLBGilLY0xfYqHe0lSLvRub2ZOMA3bv82I0za8emVykzRftXxUNZOgVqI6sm0IHFwYX4qsRYSpm8os6hMVZ1NlDpD1pa8UpEZxLREgCHA5hORHcL1xYB+AmAC4loCRHtlaqnBthMGl02eMaYs7gmzo64wpiTMigO0aiH7osk7haWweJ7XoorklQor+Lks0SjxTWbusuR8e/tq6Dcx9DeRlgn2YlHvaXDjIhOFfioELbXNEe4qSnHsCrxypL4itzx/S/W/Uqy2DHKcCHyqnml+rjxD61Idwd3ZkDkm6B4TdVrxthtAG7TXJuYpm63fiS/N0nwJxsikWUcijQpAVdu6I5Ydqggb6OzQH0X4Siu6asYFYdxOGinEUbv3rwVXTJUbus95Qp6KxW0l0rY0N0duibPnRXrt4QJgI1ewqF/vL0saa/ug8GgJrS7jxuKru6y0k8gDi4c8K0K6yre17HDBtSyNamiTw4ZEJ+FKg5qO3kfTz4WaQiCaxRKIJ5IVBjLVPEaCR7meP93ps+3ENekG0cV+MfTZSSIgN4y0yptbYZ1bUzUzLxc4XX44O+eDv0eNrC9xsl3KLSG8qmv3z4PK9bXiZ/VeCoe8fJT91QW5e1lGQFVV5Vux7jjyEE4d+qOidpSmSvroNLt8PUv1pOlTk2E0k6+FtagBcU1rYI0c1PHWZgQT+SznSRRTt6tw0TAm2vMkSSzkMnLqFXnMBRtRHh60Wr86ak3XG+tYWNMAgfXSIhZY/jADvT0VVCuVNCukOup5o6otLWyMFJQ+YljhijLchPETBWvGoZBt2MkosQy6bRrjRN5sf1SzYSyXi6L9aF6Rm4I4BWvBqQT17hb1/TEKO4qjCnzVyaF3L8ka/GtDd3G65WKu+dvHGrcuMM93BJj+pzAjE7mbmyePa5Ms4l8Z3sJvX1BWAMVF6omdiQcxY+o2rNT/dy8D2l1FSLzoXsHDGompUT2+Vqj96ajjnzeix+LdoW4JgvFtImTb0nFa6sgDXFK4vEatxgYy1a+1lcJREp8B7H7uKFO99vIoJNYGdnCZRsqF5UXls1Ci3ufjRbXyOhoI/SWAzt5FTOgIgTiuNgMp+pDpntu3t6SmN1eHMRhv2NO2Nb9/YfuVCukfj+UmGtWOVS5gE8pcdxV4tYsxFlKmXxva8euaQmkGXqddY1Jznf3vGgOWREBV5yiUxJWbuzG52+ZhT2+eifWbOrBnDfXxd8kwCZ9X5IYPnGoKV4d7lERvefeqMtR7T5Y5uf4+yz75Oh5oKOtVLeTVxAoldLeUe8aCUoG6HcwWRkJiOMuc6X8teqYCaLkMumkHPAfnwo80Gu6I6EalbgmC3GWUlxT9LAGWSCNSRZj6pe30+jB2nueiXHMSbI7MOG8a5/EP2YF6cVej3H0UcGGyDPGwDKWYvD34jISnQoiP3dp3TroZ/9+ObaOuMXI5f3NQkdbCT19lcBOXuGSb8oDCtgRw809Kk7eLK5JC3HUH3opHKKkHs1RzUzwuO1JkFQmf9PTwTyocfJCB5Timpg1bdN/1VjXnaG8uEaL1DJ5xf2u3MGe44eH6tSHUXDvrC5J9cjBHVb3x+kQgnqzd1LhtbkQERVna5tBCQiyDGWtQM4anZyT71Nz8iqiddtz9d2HzWhuUbzzOHFNWvy/f6hTbwICkYd6vRLCilcXS6+k/ef94GsyTvEaxzzYzHNVBNseb0IZj1TWNUwt53UdcPH9mhys0u74xI/EoZoIfjIWrtwUWyYPcQ1qi8j+FqUsNGbQdh1btxqZ8rW7sKE7fTz0PFFTvFaYtXVNCBaTc7PCw1v3sczK3Pempxdrr3ECqmMmKIXiVdX/JJEOxGp4neK7iJuHNh+bAe3R9+0VrxZIY9+tI26uAx764huUmK45KWWI6zRb55XsFa/T5yxz9ihWLZTpc8wBq9pLhCmW8cRbAR1tVOXkK+hQxUjJYLEHqRDDy/vAnUYqyzoEcUyMeqpU9XojSpYLGFDPmXaHh6ooOHl+KL6LOH29zXvrVBB5b0JpgTTEycRZuECcICauOK0ZlovowgWMscydoYDA+cRlh6BasM8ZvFeBgDjkJc/MAx1tJfSU9eKaOPpkq3iV3fD33n64sqxt8pI04O/1pRUboJrCRMmta1RzxobGz6vqejjfJd6jMquME9ek5eQ9kTcgjeJ1bVcPvntHNPuMK9EQS4spPb9z5j6hcqqcmjJk2d4ZB2xfO+7Nyc6xr8LwiT88k3m9GxVxW0xIIl9Ns9VvBjraA5l8r4MzlAg7xWsfBknyX919jfAb4DLuW2YuwW3PRUMLEJKLK5QRIy3r2tLbV2NuxPb7aoQ/zLyZoHqXMlSc/JaqJZQX1xiQhu7prDVciYa4gHj44iGdbbjgsJ1D5R556W351ghk7u7vVcsaQO/QkhYVFs8xJ6rX0ckqjsANHxgNt1SSuMCkGGGpyE6LQPHKsKm7jMEKRVwcgbKZm719lQhB0d2nIjxZQ2x7iyLUs8zJu9Vtp7xWQUxHKK9hQBbXpOfkVWO9Yn3gqDhu2MDItSxQCCKfxlJeZ3GQTibPtOGGbbLeqMzqOPLaWqvsqrMAd+SyRdxCUVpClfShpOOw06i6qezIwZ2J6tBBpxjvaCN0l/uwamMPxg4bELkeN1o2nHwfi46l7r69xg/H4M70wbdMiFtPgcdrhpy8JZHvLlcEmXz9fBJxjY11TWebfpx30YSdSItCEPk0CsisOGNZXFNhLDFnYvLgyytM7tUW9ucyVM8nB8FyjYkTx8WqFNdptrk3fuTQmvlr1ptllc0/EMjkV27oRrnCMGZolMiPHhI9J8Kmn32VirXpaqlEOOugCVZlkyLuHRGkODEOdauec0C73UerW0g1GdarRc/FLb2knDzH8EHpQxmrUAgin0Zco+MyXTlbCk2QwE4+qSehSbbX6KxGJqiebvyI8Jazr+Jmfx+3UJTmroZky6I+QwUClJEgs4Bs3cLR0VaqzdkdR0Wd7sYOG4DLTpmirZco4CpnLV6r7XufxptWhRJRprGWVLDRM7jqwXiVqnVmGxb4yCvvx+8eey3oo7SGgXCSD5W45qPv3AUH7zyyWjYdkc8LhSDyaaxCdEReTNor2mHrELaTN2eHUmnYRZi4WRvv1bzwiXftGvotK/aAaN/7HO3vVeILESrFdYn0nJ+KUw7fS3VOMGNar+MmxYV+3JRxij7Fz5HfPvoazvjlY1pxY1+FWZsRloi0uw4TLjxionXZWCKPqK9JHPjzqdaLy0eLe5OrHJ/EJOEqcc2gjrba/ElL5JOakMahEEQ+jcRFZ+0ihqsdF0N4gCgnX2FMr+iKmYCmuZJV6jpbpkncCsvb4h+cs7+i3nCZSoXFbnNFfPP0fYzXVWu/RHoqH2fxQGRnFZEEA7ScfDBG2wzq0AYj6zSKGwgvLt9gbHttV6+9uIbciGIS6OYbX1vBbqxeyMaBj9Nf9Ri6E0yVuCaOkwfq4lUbi54BhnH2JpQGpHWGUtfphrBMPlC86iZanBmlSczjwsl/9TR1kggAGBvD4XKIHKXMGW47fEDEIkVecH0Vt3RvKuuZOBD075ETOt0CIqoT3azXmI6T5wR1mOZZicjM8VmYjL781kZr5WOJKJFo0UUXoivbJrwfXRfEkCHDFHlWeR13XnwU/vMdO0euW/dRZScvEn6lqNDN8cqLaxIineJVF2OmfqzaRh26yyjc8OFDsFtVlCN7vDIDJ99X0V8DzNs+USYf99jv3k8vj7ZNSjxAEMnInKFKFi533VVck0SJWiLSjgUfy30nbKO8LnKQi1Z1ZZrsRSeT57FKdKF94/QEtsG8rGXyJcrdz0DH5daIPPTv/v+9u54e+pEvH1s75u+N17Hn+OHYo+r5nMR3JhQ7hyteY+zkxXdlY0UW9/HOA1s9kU+qtJ04ejCOnTyuRlzCE4SLa3ScfMVITMwyeXtxjYrBGDWkE/d9/l3WdYQ4eWmrqXoE+ZkrFWblAFa7PwG16eotaz8kdU5RXW+JUBN99FUYZlx+Ak7bb7xzH1TQiUDi0hMSkVEmLwfz0sFeJp8slML6LebnEKGrXuTkdVsp8YOnYrjEtaSbPnuNV3v7hvuoULzGWdcQ1dq3Cd9hEtV6mbwBebjji5NSNUHroVOjZSoVs+I1Lj2gab3FEQgRqoW7zaAOp6QjnSFxTbg+1fPJpxat6sKGmHR84Tqti9Yw98312g+9yCmqQFKyilFDOrHDiEHunVDWrQbX98je0Bw2MnIrTj5ncY2L6FArrqmeL5H+wxViLhRFxHmuI5S/+c+DY/sobnyOnRwoxI/aY2ztnI6T5x9TG07etLvynLwBeYSWjauTakQ+6hI9d+m6WDt5E+dk+gC4xELPwk06zMnHE3m579c8tDD0++yDdsC3z9ArV5P2OY7I6wheSSELfvYNc74AW+h2D5uqUTJ1epFA8Rq3NC04eUtxDVmKf2S4iOFsOHldF8R2lAyX8AJ17Wy3zcDYHaw496ZOHIlFV56GA3YcUTunI+LtDuIa0/zOS2KWisgT0blENI+IKkQ0VTh/IhE9Q0Rzqn+PS99VPbbbJh93YBP4vFK913lL12NzT5/xhW4yeGjaErq4NZZFxh9ReRgR1yhmT6zTSwxBSU7kzYpX7YeTotyfTfYpFd4tiXl0j7JXNVAYT6wtf/RsTBrFx9G1YyuuaSslS6LtkilJV3v9vZA2qFhXt5AwRXFd/IDrd2zxZqniR1n1gd6i8J0hqj+DDZHPi1s3IS0nPxfAWQAels6/DeA9jLF9AXwIwB9StmPEtsMbR+RDMkTULSRkTnFTT1+miYkvOWkP53qUpmXVv7bLU3TKkXcfNuIaVfsm2WMSMcR3ztxHq1vhfTRxtfIwJVW+2prtXXz8JNx58VHYY9tASXjavlEdQIdEkESrI4awvke3K3RxhnJ95B1GDnLk5HU7KS6u0X/gD99ttLEekZmRxajifXG7o9CHU3FdFTqDQLUwJDbjYZ77LSiTZ4zNZ4wtUJx/jjHGo2rNAzCQiOxs9loQ4rvjAaX4ZLrug1Px1dP2xISRYTnuxi29mcbi2Hm03iFr6s7qOOF88W+X4iN42K6jax8y+XFKCi/FtJy8jVNONPAWafUyNsRLfoakdvNyU2u7enDOwTtEyrW3lUJmgW0SMSYFJz9Oeodil3Vjbjv7ksSN+eE5+zsZLWh3G9x81STa1HDq/J2HPnIUviaejtNzqOLJi9ioSEZDBHz6uN1xyMSROFXxsVaV116LvTsZGiGTPxvAc4yxbtVFIrqIiGYS0cyVK1eqijQd4oQZ1Bkm8tuPGISPHbVr5J4NW8qJs8irxCxJnFVKJWDOFSfhwS8ek6gfQEAAOKc5RDK7tLGukUEGXubpy4+3EjHJ2+4S6a2k+GmdZygQFTsldWKTH31wZzu+fLI+PAFHVKEdfUa5S+I4H7jTCKhwx5zlsW0DwUclSUA+J3FNDCdvazEkFuHRLEcPrQeW04prKF7PIc4D1Szt6lEQeQShKf7yySOwzaD4KKZGIp8TlY+lHER0HxHNVfw73eLevQF8H8AndGUYY9cyxqYyxqaOHTtWV6zhEAm7yMnziRL3Ql5+a2NyTl5xm8luWqvUIsKwgR3hvJKOXSIifOuMffD05cdHHJVUCzduy2ri5G3HS3YyIiKtfoKf1yZbYdEFrZOtDoqJ1CjWctkpU/A/x+1utZOQn5tgjiUjP6suFIRtMLu2BOKaEpGbuEbXtqh4TUjkjtx9TL0d3a7GSs9h5uQ3dZvNl2263wxxTaxHDGPshCQVE9EOAG4D8EHG2MK48q0MPpen7jwSKzcGGxIbghTneq6DMp1ZAvvarJIQdLSVMG7YQLzSvlGqP1o2TvmUxXZVDhdA0Cte+dm9t98GO44aXItTEqTgi4aYBfRc/+DONqyzNGE9dd/xGDG4E2s29cSWVXHycVynOI62fdKhVNIr6ccM7cTbG6PPQIbdkwpif9tLVPOd4HO9RGSnuFTMEjGGkmkOxQWiiyOySk6edD90bcQWyRy5iGuIaASA6QAuY4w9lkcbeUOcTJxjuWTaZGF7GUaWZpzqnJXhcxMsbLnTWtf87b+OCP2WOaES1UeJp5aLzzJEsUrSr5waFXHsv0PdY1XVD93wc+Lf3kb42XkHCr0QuTY9J/+pY3arHR+rCCYmQqkUtFjV0bjvZvEck+oVzfySwBTF85EvhQ3juI9FICJz4eTrLYjiQ1Gevt7iY6UazpAzVHXYGIsqtOOIeJzFkioybcgix1i7fZmskdaE8kwiWgLgcADTieju6qVPA9gdwNeIaFb1n3mFtBhEcY2YQIC/0yzME3U4YMeROGWf7ULy3LCiDank7LYYPSScREMmPCKh4R8hG05eV4bXxx1RdJA5+VIpnuDIxHbvCcFHqbO9FFl44ofq4uMn1Y4vOmpX/Pmid2jbUDJ1FtMkSnyiliBRpXf9WOxjEpQMshJRr0QkEmVHcY1Q/Q4j6+GVRXFN2h1J0Kt6Q1e97wCne0NOVYrxiPediW+jGbmIU0WpZ4zdhkAkI5//NoBvp6m7lVAPXSAoinJ8V20l4NcfCDz0vn9XkH9WnIAdbaUqweXutvZ1p+l2hMgLP7k4yaTg5FDZGwOodU61kxFrlWXyJSKtEpBphuhXFxyE11d1YcTgzshF/hE6ca9tQ/oMImAXU9jpECfIxRD64tpqyGzTzaTgd2mjaJoc80yKTJfor7pxEMU1a7uScfKq6wzuwcDiuHLVFAsb9mxF4pqiQcwBWUtUkOPb2ndCdPttDIOQoA2bAE5ykc52WXZc/23rEELQh5Hlj7jrWHPYBVm2SiZxTfWK/LqGD+zAO3YdXeuTCC4vvlRK3EGgsBJbQkgEVP2bdJ7EhjVIVKsa4tSSLXVkzlNMem0Tq6VWj6bHYtiJd02ON7zIK74LEC+uiVszdkp2x05lgMIQ+R+cs1/ie7975r647oNTI+cPrxIB/nIJYeeNPDBhxCCcvM92kfN5iod0kKe07EEpJtA+Yc9AxBIXF4dI7fzD69P2RehMRIat6qx0n8meX26XpxiUc+0SAQMNcd5JQSSS0HjGmOIZox+2rECGsZAdhMSdkYrofeFEjdOeprvizvioSWPx8/cfqC5Y66vxcqpxCVnXKDqsmmJh3Y6+7p+ddwB+cf6BaIZUvjBEXjSjcsXQge3K7fGNHzkUz19xUijfY30bLnM42WDHUWqFqm2UwI8cuYvxussikBdxNHYNah+kCw7bGXOuOAmTq6Fete2D8M5JYzBycNSm2NQ1UUcSJUR6+TA/G5F6G7i23ao7CZXJpG2qQL74E4ULiJ1MdvmDvyaE6DVh7LC67kWeZ1FOPkAw5tG69tlBE9JZ03Y93n/w1yaDlAn89kShhkU7eUVD6tST9eMV65WuQACA0w+YgHfvt31TxDX5ZI5tAtIwugPbS8rB72wvobO9JGRzr4trGv2ubKNWHrrLSDzz+mrMXrIu8z7IIgQiwhXv2RufPWGPiKOUDjXFdUJLFCA6FqZ3X9EI5cWmJm07DC+tqJuHXv3+A/H84nVK+3PTR1J1JcmiFj9YormhCJv5fvyUcfjWv16ILbfdNnXGwjb/RWBCGe2XLhCcbtzkMCEu6/j5K04Ck/QCacQ54dg10esqvY9o6bV8/Zb4NpJ1LRUKw8mnkZEHFhbxogLRCEE3aT97QjpLBx3M0evq1xgD/vyJwzHj8kTuDSHIU1pFXNvbSrF5VEXwGlTjZ+Tkhc6okpfowj7UabxeXPP9s8OivuEDO/DOSdGdoa3STzxOw8k/+7UT8YePHhapmxnCWOv6o8K51ZALYgKZOE66JrqkeKIX6oumvlJNJh/m6HUQrw8f2IFtpB2hy3DL1mNxIYtVOghxt/eDs+NFxs2wrikMkU8zdnGecLUtP8Vvw10IngpxCqpQvxSx7BmAgR1tWi9Il2GKmFAqZPJx4Mm/D5k4stpXPn7RsnH1cVPBqHwd+O5Z+xrvjZog1k8MifFkrdcRJ0oQiUT4rws4IR01pFObXcqm3rjx/P7Z++Hl75zidE9d/KUW1+gsWnQ7BDk9Y5xYMu65+XUbYY2s5xJ/quanauciOmKNHNKJX11wUC1bnLLNJrDyhSHyaTj5DoW4RnyfdcVrvZDpZT1+aXxkZU70bGH6Dol9d7FdljFYIHZPf+X4wLxQwFAprIGNMviyU/bEoitPwyETR4XOqyKHmmpjrC4PZwx45EvH1u8jwgE7jsDzV5wUuY8TkYioSTy2nDtxpd7aIGzXU3HyYvx0xcc9uOBcr4xSKRo+Iba/XD+l8U0QzVvFRB1xXtk1q7UYimSteLVYBjJDFhdqWGUyKn+ET913PPY3OKflaR2kgyfyCAhAHIEBwhPQROAGW3CGFx/vFjo4qxjzJuwqcCCDNQmTjxHM3Fy4Etna5Lcfmor9JCVd3DO21V4Aw46jBtc4cN6PoZ31Pv/vRw7FNR84CB86YiI+9s5d8Ml3hYPIpaGR791/eyX3/8CCeoC9uvjBvX6RmOjG2MpcLwHbaHsLQR2GQCR6YtJt3Tjwb4ytojpenGO8XMO9nzsaO0qRY8PWNVGoxDVKk1rDGowq+Q1+FxnBK14RbzEhKl51CMmMLZxTXBe/2U5ekMm7VRtCWXBk0inQxHXt8mGtEb3q33HDB+Kcg3fA84KC2GxdU+8TH2vR6gkIE7WjhbRtX1VYmSSRjfJbrn7/gfjto68ZlZpxuhsTQpy8hv2wcrxxbtlCJl/9W6IgprwcnykkrjFUddHRu2LYgHYsXRfsfuJ0GD8774Ca1ZMJts88adthRgW+qh+qj9ogg9+EVT/GDdP6jWSFwnDyusU0akgn/vcjhxrv7WgrGWdHpSauiZ+MRPG5Na94z17a5vRx0Q0fmFDETLsAYWKpcVX5vehurs1PK0x0JyKvKCovGhNBZIzVTDj5XTpHp7ygkrnryyaH+ApVz8YYqxGkWz91RLRAFUl2t7EyeUHx+uP/iIYNEMU1pvhAXzl1T/zP8ZMETj6A7iNz+gETsM8EtXmmCJePasQUNy5riAIm5zgV5O41Yu4Whsjr6Oq79xsf4upUiPMurCte7WTycdzQhUfuop1EunCmpjrTCHqBXgAAGtdJREFUBkebvN0wXPfBqfjJ+/aPbU8kzC4TlBcV7zF5x/7z0++MnKtz8qz6Nzifp/exCLGZuJ1jGiuKuJymQH1O8sBtH33nLoo+uLdty8kTkTJ+usjJ24wXV7RS7a99X1VwuV0WZ7m8Xw4VkTctx6jjXYaRDXVt5t5Cg5AmO05HmymVhdqKxURY4jj5oF/qMis3qB0q4uK5cCRRvA7qaMOJe22L4QPri1b3CGL9TsSVL2Lh1BkHTsA+E4Yri++rcKrhMnn5CZtglaYk4roMRq4Q173KY5ch+oH74OE7K/ro3rat+EFXtd6EUrfzDRP3tB/sWhTK6u/vnrlvzRtbhiylDceTt+uHe3yc8G9tnoMMUUgiv8uYujLD5mV1tqmdoTiYQiZvKi+3ec/njlaUUd8bstAQILZ96C5hSxXZhtoVKkWxbtzCRN6+DRUnP2boAPzrf46yroOLazgRFL0vGwGxFVWTHW128yMONh/qCguLqpT5dhN8agZ1tuGprxyP2V8/CQ8pIp3G7Z5CnLzYFx0nL9nJJ82vW28zvNs7/7CdcP2HDlGWlZ8hbCdvBxuGTtU/Ds/JO0B8Xw9ccozTvXHiGv4ewpNWfrn6l7W9Iva7bmro3rk4+a/5wMGha0zB+emw1/go5yybSpogilhcFmRaOsxYVFwDxQ4rV8SIbEU/AlcCK3LAcXoVxkTZeNCOypJGPHV1TEwYEdsOH4htBnUoHczi9CCh3YwF0eTn64H/rLtprtCmqPQQtjt1EapxN70/udr4/AvpUUgir/ptQocipjgLHdc5eXlSRvpRLbH7uKH4+nv20paVbc51GFV1SBK5DFMIgbjYMVdKXnmn7Tsel0ybbNUXIPwRcpE7Z2EfLH9UbKyeskToGVR21CGNqVvdd322vqNRhrSVfssJLFRzTHw/J+65rVuHYuAaVlf3imoerxmJa1zulnm7uLAGKrhz8mGULUJzp0VhiHw0T2b1rzSq75u6Y+Te9hIZCRa3W3b50t/3+Xfhw9VgYeKCOKgaynXv7beJcOQqPPKlYzH7/51kdBLhtOWcg3eItUCQFUVfnDY55NYeB07IbvnE4db3ABlw8mC1vnPup9HiGhGqFvsMClOTgwwQDq+sEtfI9V33yGuh3+pYQMJxRivdVhxYInmM1O+oHoWSQr+TwoXxkJmGtgRE3lZXpqvXi2scECHyNUVf+Pz3FSGJY61rBPknr3fVRrWCVGUCKXbtb/91ZO1YFVJYxpAB7dhmcIfZuqbapk22+LTgBMjVPjgLMlwn8rLpZQaVWyCOMw1ZHgnnH7/0ONz08cOs2zn9gAnG67YRFsMe2ukHSVT4mqq74cOH4IFLjrGyVpHPp1a8Vm+3GSFVeIzaNcsZ68zJR6xrvLjGGvJY6zh5FcS0fiocNSkwwRw2oAOfODrwnIzKK5NNzhGKkLsqZM2tJjW7VO1qbFCzz3do96aPvwPnH7ZT7b6BVaVeT7nKydfENW59SQqxGVtnGSDQyQzutN8tjRpiryMxIoGM2QTbd3fs5HHV9aGWz4vg/eL+F+lNKB04eYOdfBpO3gQvrkkBnRLF9ApO2288zjpoQmy5b52xDx7+4rHYZnAHTtl3POZccVLtPru+6a/d9/l3WdVhM5kaQets5OB7KpS7SezGD99tNL522l4YOqAdX5w2uRbxTxbXNCqyn0t6uKz6pCKstmQhLuBWGtg8XmjnoynD+yV/H1VzyKVfNh8kVQL1er/sBswUOFAFeV54xasjtt9mIL51xj4A6l9007v65fkH4ScKrz0Zne0l7DS67g06bGCHk3WNacKMGToAJ1goxVTR+ZIkRkgLnvlpmEFxfOfF9maRcRjU2Ya535iGk/beriau6eFEvsHOUCLimsy6R0naiwu4lQQuc05sMa75rLyX426/4cJD8JVTpyjbEkUv4qW7PnuU0g8ByICTb4BMvjCxawDg8cuOr/+QRnPKdsMicTZEZEUnVNvFuKovPWUy7pu/wljGFGzKRk6qg+s9PzhnP5x/6E7YcdTg+MIp2lGBZ++SuR+XdfbrCw7C4wtXJWo/RLRi3moe350nLjsOX751Dp57fY1V+TzEWC67J9la5YQ9x2E3KT1kb5XIyXqxxAxMTSavvv/YKeNw7JTAOUpmnMQ+iJembDccO2nmu5wSM2jb0D1Z8doAcU2hiLyIukw+OJr+maNinEzy4wbjF0T9+o/P3d9QLm07UbiupcGd7TgiQarFLEwo+SLsLYc77cLJn7LveJyiyTEbBwpTeXPZjOaTSKzGbzMIu48dihmvrY6UU8nxXfpwxgHb48EX38JnjjcnvdEk2lJC/iiqnJJ6q/oV7ieQdtxc5sKFR07E9Y/WrZTaQ0ReFuXodAr6+t+5+xg8+srbxv61vLiGiM4lonlEVCGiqcL5Q4loVvXfbCI6M31XXftW/Vv93aaInd2wvliW23XsEJxdzdbjgqnV2PSH7zZaW+avnzwcv/+w2vOvEciCs+XelKYt7p8+fpgy7k0WcApQlrUMXGhRZZHR3lbC504Ih6926cOwgR347YWHKOP8K/vjKpPXlOeiN9vcuXHgRNSGeRED8gX3mupVnzcxVudOjV/L3eX8iXxaTn4ugLMA/EZxfipjrExE4wHMJqJ/MsbKKduzRkmm8jHIU6ybV918Hh+880i88M1pRguOqVLSjlFDOvHG6q5Icu6/fvJwvLl2c9ZdzYSvHT2kEx975y4466Dw4hG5oyN2S57Q3QV5KHun7b0t3rGr/kPNIZuQcqRxCLSFTgwSB914cU62o/oBT2KFJYLLyE3B73Qw7SJcdgguHq/rNvda15sUqYg8Y2w+EH2BjLEu4edApAtznghdPYFHoO7FnXVg2DomUsqxxzYa9ViLgcRiSHIy0QOAa//zYNw7f0WEm5k6cRSmau7JAmkIDxEpY8Nn5egT3wHlobpoguf8zX9GR16eV6Z6bWzOzztkR9w8YzFujAm/rYMuZ64KNjufnnJYJi93+eCd3TKocbGP7kOYFJnp7KSKdho1GG+s7tKUzga5yeSJ6DAAvwOwM4D/1HHxRHQRgIsAYKeddsqhH9Fzi648LfN2TO0BwK2fOhy7jolPehAHF+9UE8YNH4gLDlNbDOSBPHdKTQk1LH1YLjhsJ7x3/+3xvmufDMpmJJPnIiobf4qIHFlR5sqz94uEtnBBXfEaX9ZFXDNAiubIwPDEZcdhxCA3n4GO9qChJE5GpmdyCuFhKCteufdzR2PkkE5M/fZ91nUnQSzFIKL7AKhcMy9njP1Ddx9j7CkAexPRngBuJKI7GWOREIuMsWsBXAsAU6dOzezzO6C9hO5yxXqppd1+v2f/7XHzjMWRXKYcB++sPh/uhPnyLZ84HDuOigY7a0UcMnEknntjbUPaaoYzlEzEv3DS5JDyM6vvzpTthuGbp++Nd++3faQPkf5ZcPJZwbVm3UePK147FIrX8du4z3VeT9ZWKy5zzFZcM2lbc5yprBBL5BljJ6RpgDE2n4g2AdgHwMw0dbngzAMn4OYZixvVHI7cfUyuOwQgGmK4lfGXT6ozFmVp2t9eoqoStkFUXkBE/i1fz6wdwgcPn2hXVo4Wk4dM3uH9idy5jkhOGT8Md81bjgmKSK1JUBfXJODkDdeSfDBVjGMzfDpyEdcQ0S4AFlcVrzsDmAxgUR5t6SBHt4tD44degcb7NjUMeWSpH9TZhg1byo3j5E3b8IjSM59OmUUKjeiDvQ2lGGJblQQbAP7nuEk4dvK4WgC3sdVUlKcmNHOtcfIZOxm5jGRzDLX1SGtCeSYRLQFwOIDpRHR39dI7EVjUzAJwG4D/Yoy9rasnD9TSilkOaxM+sC3Rdn8GT3bSnKQhsvw7Xh6eNxrxsXv/oYHebNiAeB2BGBJ7S6+as24rUShC56ghnZhzxUn4zHFme30duClmIusao1I7ASdvfTJfpLWuuQ0BEZfP/wHAH9LUnRZtjpy8jCN2jzdl82guAoui7oa4hgOSIjFyUV822z6YlHr5U5DPn7gHPnP8JGufk9FDOrFqUw+29KpzF6swbGDyaKodjun4bJGnR3zeKFTsGhG1ONWW5cXBv/Ejh+LXF8THevdoLni4YxcCkhXibNIbETTtyydPMfYhDxCFnQr/9PHDcMYB22vLv2tyEMG1Ue9Il2M2LZw4eYHn+N5Z++LWT9VzLzRj117YsAa1d205qmKx3cYOqUU8bCQKLJLPBcdOGYsXlq1vSBx9IMwI6JLUNBK7jQ2Hu25UNE4RR+w2BkfsNgZ/n7VUeZ0HldvSAM9OID43hAlGZ6gE1RLVxVu1eoqieG0F1BSvTe6HDZL2sQlBKFsKXzhxMs49eEfnYGlJYRLXNCzcsXAseyu34lznH+BGRUxNmwhcB9UH4B//faSiZFw9jUdhibwqNK8tmsERpUI/625WKJUIE8dEk003AnEmlI1Am8ReNsrKyAX/c9zuYAx43yHRtJstBwfLpTMPnKBN6WgK/eDFNRminjvSrrxRqdai2MoZ+YahREFSi/Bcyt8mXQmhHTn1XCsyJ4M723HpKVPiC2aMrD94yaxrovd4xWuGMMVfj0Oz1krSLW3rLe1iwSbVH1+8nz5294b0CYiKJlRT/r37b4/LT90zk/ZGZ5WWMGdc98Gp+PcXjrEq+6sLDqoduzhDmdZqkyKaa1FgTj74a2tdFw6m5MlmXuiPuw/VR/+VtzYqy1wybTIumTY5v74IczOSRFrR0avff2Bmbdumqmw2TtwrPtMah22awSQ8o2reNEOkVlxOvhZX2jKzPamPG4lW3G5nhTQhYJsNbiEiTqVT990OwzIKFpcUkRylObc3sp9w8i4wObiFymVlJ++ta7JDnZNPEle6OWhGztZGod3giXj1+w/Euq6eRnfJGn/71BG4e97yGrEHguBSc74xDVO+die29FaaYukkp55rhnlef4ftmEXCqRvKGsOOW7WWLQpL5PlLsY1T5JLaLWsUmYPn4ItJ5Z363v31zjStgEnbDtNGDDxmj3G4a97yhsW0F6dKhJMv/jTKHLYGF0ms9VR3eOuaDMEXQH/gjkcPDbbB5x2afTz9VgGXH1daWFyTxKnqp+cdgOXrtmBAe+Od5/qDnXyrw5bo7rm9neweAAZ0BF98lfFHM/R9hSXynLbbkpRmKl6HD+zAa987taFtNhr8o9uoODOueOGb0xKJOwZ2tDXUVl/sYdS6xpN5V4hjZhq+CSMG4Ufn7o9lazfjx/e+ZBTJXPGevTF+m0E4fsq4yDXPyWeIHikLfByarXgtusimLpNvjHu7K1zTJ7YCotY1zekHANx58VFY25V/vtKs4bLszjl4B/xj1pux5UYO6dT6BnginyF6+oKASJ0JotIVm9w2B9w7s1U5+f4Ck0y+mZy8rSliq0F0mmzE8OUVQM2E4hL5sjp3pA62plQeydDej00oWxWydY2fte7gYxbZFWlwyj7j8eShq/H5E/dI1p43ocwONXGNLZHvj2EN+oFSmYNznb0Z597cmuGta9KDE90SkZUurrO9hO+dtW/e3coUhXWG6q5x8u5WD/1lsXz3zH1x5oETcMRuY5rdlVh0tLhMvr/A5PEqims62rKdxP0lpIEr+BAWWWntOfka+t9L3nHUYFz1vgOa3Q0rmOzkPZKhTTahrP58z/7b4+cZhjNIanmUJ4YPzIZ0cU6+vUQNY+4mjRuKTd3lxjSGAhP54VWb55GD7Wyfw+Ka1prQRQCXH3uZfDqI87RDlsk7hvKwRatZHj152fGZJfXhwynatOdN7O/53NENlc231tvLEF8+eQombTsUxylsVWPhaXzm6M+xa1oVutg1/UhVkwjbbTMws7p42BNxLPNe/o1WvhaWyA/qbMMFh+1sXT5sXZN9f7Z2mGLXeNhDnJrRePLBX1PSCo8w+HwMO0UViwAUVvHqigkjB9WOi/WKWwOt7vHa30AUdZuvR15tRo/6J4ZUI4meffCE2rmiMSKpiDwRnUtE84ioQkRTFdd3IqKNRHRJmnYaAdEKp2hf8lYAV5RNGDEopmSxMWbogEzqUdl18zNJIq9urRgyoB3zvjENX57W+OxVjUJacc1cAGcB+I3m+lUA7kzZRsPhSXz22H3cMFzzgYNx5O6jm92VpmHeN6alTzRdZUBUpsHkOflE4Nx8uaA+HKmIPGNsPqDmfInoDACvAtiUpg2P4uDkfbZrdheaiiEZJhkZOSRqNcaXYcGkDQ1D0cQ0HLnI5IloCIAvA/iGRdmLiGgmEc1cuXJlHt1xhpfWeLQq+NQcOTjqnFSftsUkVnmjqPqiWCJPRPcR0VzFv9MNt30DwFWMsY2GMgAAxti1jLGpjLGpY8eOdel7bvB28h6tjhEKIu8Vr+lQVG/s2P0jY+yEBPUeBuAcIvoBgBEAKkS0hTH2iwR1NRyek/doVfC5qcovWxfXeCqfBJyT/9Dh9qbX/QG52Mkzxo7ix0R0BYCN/YXAe3j0B6hi09Q4+UZ3piDgMvm2RuVybBDSmlCeSURLABwOYDoR3Z1Nt5oLz8l7tDraVXHJveI1FTiRl9Mq9nekta65DcBtMWWuSNNGM+Bl8h6tCj43TXby/SkEdSuhXImGOCgCirUv8fDYSqDiNlstUmR/Q01cU7Bx9ERegYK9Y48Cgc9NOSsUAOw0ajAA4KhJrZ9foBVRVE6+sAHK0qBYr9ijiFApXieOGYKnLz8eYzMKnbC1gZtQ2qYC7C/wRF4BH7vGo1XBZ6bOAmTcsOzC8G5t6KuaycvJWPo7vLhGgWK9Yo8ioa+qVO0sGCFqBZT7isnJeyKvgGfkPVoVdTM/v3SzRtnbyW898OIaj1ZFb18xlYOtAO4p7Dl5Dw+PpoGLFFSKV490KKp1jSfyHh79CJwQqUwoPdKhr89z8h4eHk1GueI5+bwwbZ9tAQAH7TyyyT3JFt6E0sOjH6HcV0zlYCvguCnbYtGVpzW83QcvOSZXb2VP5D08+hG44rVoQbS2ZkwcMyTX+j074OHRj+DFNR6u8ETew6MfodznFa8ebvAzxcOjH8Fz8h6u8ETew6MfwStePVzhZ4qHRz9Cb0GzF3nkB0/kBRy9x9hmd8HDw4iax6vn5D0s4U0oBVz3wYOxcUu52d3w8NCi7Dl5D0d4Ii9gQHsbBgxta3Y3PDy0KGo4XI/84Pd8Hh79CNwZqsOHGvawhJ8pHh79CF5c4+GKVESeiM4lonlEVCGiqcL5iUS0mYhmVf9dk76rHh4edXGN58887JBWJj8XwFkAfqO4tpAxdkDK+j08PAR4Tt7DFamIPGNsPuAzKXl4NAq93oTSwxF5zpRdiOg5InqIiI7SFSKii4hoJhHNXLlyZY7d8fDo/+hsD5ZsR7tnrDzsEMvJE9F9ALZTXLqcMfYPzW3LAOzEGFtFRAcD+DsR7c0YWy8XZIxdC+BaAJg6dSqz77qHx9aHGy48BP+cvRTbDR/Y7K549BPEEnnG2AmulTLGugF0V4+fIaKFAPYAMNO5hx4eHjXsPHoIPn3cpGZ3w6MfIRdxDRGNJaK26vGuACYBeDWPtjw8PDw89EhrQnkmES0BcDiA6UR0d/XS0QCeJ6LZAP4K4JOMsdXpuurh4eHh4Yq01jW3AbhNcf5WALemqdvDw8PDIz28HZaHh4dHgeGJvIeHh0eB4Ym8h4eHR4HhibyHh4dHgeGJvIeHh0eBQYy1jpMpEa0E8HqKKsYAeDuj7vQHbG3PC/hn3lrgn9kNOzPGlPlLW4rIpwURzWSMTY0vWQxsbc8L+GfeWuCfOTt4cY2Hh4dHgeGJvIeHh0eBUTQif22zO9BgbG3PC/hn3lrgnzkjFEom7+Hh4eERRtE4eQ8PDw8PAZ7Ie3h4eBQYhSDyRHQyES0goleI6NJm9ycrENGORPQAEc0nonlEdHH1/CgiupeIXq7+HSncc1l1HBYQ0bTm9T45iKitmjryX9XfhX5eACCiEUT0VyJ6sfq+Dy/ycxPR56pzei4R3UREA4v4vET0OyJ6i4jmCuecn5OIDiaiOdVrV5NLYm3GWL/+B6ANwEIAuwLoBDAbwF7N7ldGzzYewEHV42EAXgKwF4AfALi0ev5SAN+vHu9Vff4BAHapjktbs58jwXN/HsCfAPyr+rvQz1t9lhsBfKx63AlgRFGfG8AEAK8BGFT9fQuAC4v4vAhyaxwEYK5wzvk5ATyNIG8HAbgTwCm2fSgCJ38ogFcYY68yxnoA3Azg9Cb3KRMwxpYxxp6tHm8AMB/BAjkdAVFA9e8Z1ePTAdzMGOtmjL0G4BUE49NvQEQ7ADgNwPXC6cI+LwAQ0XAExOC3AMAY62GMrUWxn7sdwCAiagcwGMBSFPB5GWMPA5ATJjk9JxGNBzCcMfYECyj+/wr3xKIIRH4CgMXC7yXVc4UCEU0EcCCApwBsyxhbBgQfAgDjqsWKMBY/BfAlABXhXJGfFwh2oSsB3FAVU11PRENQ0OdmjL0J4EcA3gCwDMA6xtg9KOjzKuD6nBOqx/J5KxSByKtkU4WyCyWioQgybX2WMbbeVFRxrt+MBRG9G8BbjLFnbG9RnOs3zyugHcGW/teMsQMBbEKwjdehXz93VQZ9OgKRxPYAhhDRB0y3KM71m+d1gO45Uz1/EYj8EgA7Cr93QLD1KwSIqAMBgf8jY+xv1dMrqls4VP++VT3f38fiSADvJaJFCMRuxxHR/6G4z8uxBMASxthT1d9/RUD0i/rcJwB4jTG2kjHWC+BvAI5AcZ9XhutzLqkey+etUAQiPwPAJCLahYg6AZwH4PYm9ykTVDXovwUwnzH2E+HS7QA+VD3+EIB/COfPI6IBRLQLgEkIFDb9AoyxyxhjOzDGJiJ4j/czxj6Agj4vB2NsOYDFRDS5eup4AC+guM/9BoB3ENHg6hw/HoG+qajPK8PpOasinQ1E9I7qeH1QuCcezdY+Z6TBPhWB5clCAJc3uz8ZPtc7EWzLngcwq/rvVACjAfwbwMvVv6OEey6vjsMCOGjgW+0fgGNQt67ZGp73AAAzq+/67wBGFvm5AXwDwIsA5gL4AwKLksI9L4CbEOgdehFw5B9N8pwAplbHaiGAX6AarcDmnw9r4OHh4VFgFEFc4+Hh4eGhgSfyHh4eHgWGJ/IeHh4eBYYn8h4eHh4FhifyHh4eHgWGJ/IeHh4eBYYn8h4eHh4Fxv8HDzXiGo3i9QsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(trace['center'])"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x1417881d0>"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD4CAYAAADiry33AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXzU1b3/8dcnk42E7AuE7EBYgqwJIAi4VURai1oXbNWqtZRetbW3/lq97e1tb9vb/XbTFqlLq94WpRZFxYJVERBUEgg7hBCWhASykYQkZD+/P2Zox5hlkkzyneXzfDzyYGa+52TeicnHb873fM8RYwxKKaV8V4DVAZRSSg0tLfRKKeXjtNArpZSP00KvlFI+Tgu9Ukr5uECrA3QnPj7eZGRkWB1DKaW8Rn5+fpUxJqG7Yx5Z6DMyMsjLy7M6hlJKeQ0ROdnTMR26UUopH6eFXimlfJxLhV5ElojIEREpEpFHemk3W0Q6ROTm/vZVSik1NPos9CJiAx4HrgOygdtFJLuHdj8BNva3r1JKqaHjyhn9HKDIGFNsjGkF1gDLumn3IPASUDGAvkoppYaIK4U+GShxel7qeO2fRCQZuBFY1d++SimlhpYrhV66ea3rkpe/Ar5pjOkYQF97Q5EVIpInInmVlZUuxFJKKeUKV+bRlwKpTs9TgLIubXKBNSICEA8sFZF2F/sCYIxZDawGyM3N1bWTlVfr6DQECDh+J5SylCuFfieQJSKZwGlgOfBZ5wbGmMyLj0Xkj8BrxpiXRSSwr75K+Yrqhhaee/8kO45Vs7uklsjQIKanRHHZ+Hhun5PGiGCb1RGVn+qz0Btj2kXkAeyzaWzA08aYAyKy0nG867h8n33dE10pz1FQUsuXn8/nbH0zU8ZE8bm5adRdaGNPSS1vHa7giS3H+OrVE7htdiq2AD3LV8NLPHGHqdzcXKNLIChv8VJ+KY/+bR+JkSGsuiOHS5KjPnL8g+JqfrrxCPknzzFvbBy/Wj6DUZGhFqVVvkpE8o0xud0d0ztjlRqE/JM1fOOlveRmxPDagws+VuQB5o6N468r5/HTm6dRUFLL0l9v5Z0jFd18NqWGhhZ6pQboXGMrD/55N8nRI1h1Zw7RYcE9thURbs1N5dUHLyMhIoR7ntnJ/2w4RGt75zAmVv5KC71SA2CM4eG1e6hqaOXxz84iMjTIpX7jEyN4+f7LuOPSNFZvKeaWJ3ZQUtM0xGmVv9NCr9QAvLq3nLcOV/Do0klMTfn4cE1vQoNs/OCGqfzuc7Mormzgk7/ZysYDZ4YoqVJa6JXqt9b2Tn6+8QiTkyL5/LyMAX+epVOTeP3BhWTEh/Ol5/L571cP6lCOGhJa6JXqp798eIpTNU18c8lEAgY5VTItLoy1K+dx9/wMnn7vOLes2q5DOcrttNAr1Q8NLe385q2jzBsbx+UTut21rd9CAm1899NTWHXHLIqrGrn+sW18eLzGLZ9bKdBCr1S/PL3tONWNrTxy3SS3L2+w5JIkXntwAbHhwdzx5Ae8UnDarZ9f+S8t9Eq5qLmtg2d3nODqSYlMT40ekvdIjwvnb1+ez4zUaL66poC1eSV9d1KqD1rolXLRq3vKqGpo5d4FmX03HoTosGCeu28OC8bH8+jf9rH1qK7mqgZHC71SLjDG8Mx7J5g4KoL54+KG/P1CAm387o5ZjE8cyZef38Wh8vohf0/lu7TQK+WCD47XcLC8nnsuyxi2pYcjQ4N45p7ZhIfYeODPu2hp77rdg1Ku0UKvlAueee84MWFB3DBzeDdIS4oawU8+M41jlY2s2lw8rO+tfIcWeqX6cKaumTcPnmX5nDRCg4Z/TfkrJiZy/fQxPP5OEcWVDcP+/sr7aaFXqg/rdp+m08Btual9Nx4i//mpyYQEBfCtdfvxxKXFlWfTQq9UL4wxrM0vYXZGDBnx4ZblSIwI5RvXTmRHcTXbiqosy6G8kxZ6pXqxu6SW4spGbs5JsToKt85OJSEihNVbdKxe9Y8WeqV68df8UkKDAlg6NcnqKIQE2rh7fgZbj1ZxsEynWyrXuVToRWSJiBwRkSIReaSb48tEZK+IFIhInogscDp2QkT2XTzmzvBKDaXmtg5e3VPG0kuSiHBxvfmhdsfcdMKCbfxhq57VK9f1WehFxAY8DlwHZAO3i0h2l2ZvAdONMTOAe4Enuxy/0hgzo6f9DJXyRG8ePMv55nY+4wHDNhdFhQWxfHYar+4po6z2gtVxlJdw5Yx+DlBkjCk2xrQCa4Blzg2MMQ3mX1MBwgGdFqC83qt7yhgVGcK8sUN/J2x/3LsgAwP83wcnrY6ivIQrhT4ZcF5ZqdTx2keIyI0ichh4HftZ/UUG2CQi+SKyoqc3EZEVjmGfvMpKXdtDWet8cxubCytZOjVp0GvOu1tKTBiXjY/nlYIynWqpXOJKoe/up/xjP13GmHXGmEnADcD3nQ5dZoyZhX3o534RWdTdmxhjVhtjco0xuQkJ7lnnW6mBevPgWVrbO/nUNOsvwnZn2fQxlJ67wK5T56yOoryAK4W+FHC+UyQFKOupsTFmCzBOROIdz8sc/1YA67APBSnl0V7fW86YqFBmpsZYHaVbi6eMIiQwgFcKevxVVOqfXCn0O4EsEckUkWBgObDeuYGIjBfHSk8iMgsIBqpFJFxEIhyvhwOLgf3u/AKUcre6pja2HK3kk9M8b9jmoojQID4xeRSv7y2nrUP3mVW967PQG2PagQeAjcAh4EVjzAERWSkiKx3NPgPsF5EC7DN0bnNcnB0FbBORPcCHwOvGmL8PxReilLtsPHiGtg7Dp6aNsTpKrz49YwzVja28p3fKqj4EutLIGLMB2NDltVVOj38C/KSbfsXA9EFmVGpYvb63nNTYEUxLibI6Sq+umJhAZGggrxSUccXERKvjKA+md8Yq5aS+uY3tx6q47pKkYVt3fqBCAm0suWQ0bx48q8M3qlda6JVysvlIJW0dhmunjLI6ikuunjyKhpZ28k7o7BvVMy30SjnZdOAM8SNDmOGhs226umx8PEE2YfORCqujKA+mhV4ph5b2DjYfqeSa7ERsHjrbpquRIYHMyYzlHS30qhda6JVy2H6smoaWdhZnj7Y6Sr9cOTGRwrMNlJ5rsjqK8lBa6JVy2HTgLOHBNuaP96y1bfpyccbN5iO6dIjqnhZ6pYDOTsObB89yxaREQgKHf1/YwRiXEE5q7Agdp1c90kKvFFBQWktVQwuLs71jto0zEeHKiYm8V1RNc1uH1XGUB9JCrxTw9qEKbAHCFRO888ajKycmcqGtg50naqyOojyQFnqlgLcOV5CTHkNUmGfsJNVfczJjCQwQdhyrtjqK8kBa6JXfK6u9wKHyeq6e5J1n8wDhIYFMS4liR7EWevVxWuiV37s4B/3qyd5b6AHmjYtjb2kdDS3tVkdRHkYLvfJ7bx+qIC02jHEJI62OMijzxsbT0WnI03F61YUWeuXXLrR2sK2oiqsmJXr8ImZ9yUmPIcgmOnyjPkYLvfJrO4qraGnv5CovHp+/aESwjRmp0byvF2RVF1rolV97+3AFYcE25o6NtTqKW8wbG8e+03Wcb26zOoryIFrolV/bUljF/HFxXnc3bE8uHRdHp0Hn06uPcKnQi8gSETkiIkUi8kg3x5eJyF4RKRCRPBFZ4GpfpaxyoqqRUzVNLJqQYHUUt5mVFkOwLUDn06uP6LPQi4gN+z6w1wHZwO0ikt2l2VvAdGPMDOBe4Ml+9FXKEluO2hcBW5TlO4U+NMg+Tr9TNyJRTlw5o58DFBljio0xrcAaYJlzA2NMg2MzcIBwwLjaVymrbCmsJC02jIz4cKujuFVORgwHyup03Rv1T64U+mSgxOl5qeO1jxCRG0XkMPA69rN6l/s6+q9wDPvkVVbqcqtqaLW2d7LjWDWLJsRbHcXtctJiaOsw7C2tszqK8hCuFPruJhebj71gzDpjzCTgBuD7/enr6L/aGJNrjMlNSPCdP6WVZ8o/eY7G1g6fGra5aFa6fRvE/JM6fKPsXCn0pUCq0/MUoKynxsaYLcA4EYnvb1+lhsuWo5UEBgjzxnnXJiOuiA0PZmxCOPkndeaNsnOl0O8EskQkU0SCgeXAeucGIjJeHLcVisgsIBiodqWvUlbYUljJrPQYIkK9c7XKvuSkxZB/8hz/unSm/Fmfhd4Y0w48AGwEDgEvGmMOiMhKEVnpaPYZYL+IFGCfZXObseu271B8IUq5qvJ8CwfK6rnch6ZVdpWbEcO5pjaOVzVaHUV5gEBXGhljNgAbury2yunxT4CfuNpXKSttK/K9aZVd5TjG6fNOnmOsly/WpgZP74xVfmdLYRVx4cFMGRNpdZQhMzZ+JFEjgtilF2QVWuiVn+nsNGw9WsmCrHgCArx7tcreBAQIOekxOvNGAVrolZ85WF5PVUOrTw/bXJSTHsPRigZqm1qtjqIspoVe+ZWLyx4s9MEbpbqalWYfp999qtbiJMpqWuiVX9lSWMnkpEgSI0KtjjLkZqRGYwsQHb5RWuiV/2hsaSf/5DmfXPagOyOCbUwZE0me3jjl97TQK7+x41g1bR2Gy/1gfP6iWWkx7Cmpo62j0+ooykJa6JXf2HK0khFBNnIyYqyOMmxyM2K40NbB4fLzVkdRFtJCr/zGlsJK5vnQblKu+NeNUzp848+00Cu/cKq6iRPVTSzK8o/x+YuSokYwJipUL8j6OS30yi+8e3E3KR9e36Yns9Jj9A5ZP6eFXvmFLYWVpMSMINPHdpNyRW56DGV1zZTVXrA6irKIFnrl81rbO9leVMWiCQk4VtP2KznpsYB9gTPln7TQK5+365Tv7iblislJEYwIsunwjR/TQq983pbCSmwBwvzxvreblCsCbQHMSI1m1ykt9P5KC73yeVuOVjIrLZpIH91NyhU56TEcKKunqbXd6ijKAlrolU+rON/M/tO+vZuUK3LSY+joNOwpqbM6irKAS4VeRJaIyBERKRKRR7o5/jkR2ev42C4i052OnRCRfSJSICJ57gyvVF+2FFYBcMXERIuTWOviSpY6fOOf+txKUERs2PeBvQYoBXaKyHpjzEGnZseBy40x50TkOmA1MNfp+JXGmCo35lbKJZuPVBA/MoTsJN/dTcoVUWFBZCWOJO+E3iHrj1w5o58DFBljio0xrcAaYJlzA2PMdmPMxVOF94EU98ZUqv/aOzrZerSKyyck+PRuUq7KSY9h16laOjuN1VHUMHOl0CcDJU7PSx2v9eQLwBtOzw2wSUTyRWRFT51EZIWI5IlIXmVlpQuxlOrdntJa6i60ccVE/x6fv2hWegx1F9oormqwOooaZq4U+u5Ohbo9JRCRK7EX+m86vXyZMWYWcB1wv4gs6q6vMWa1MSbXGJObkKC/mGrwNh+pJEBgoZ+tb9OT3IsLnJ3QcXp/40qhLwVSnZ6nAGVdG4nINOBJYJkxpvri68aYMse/FcA67ENBSg25zUcqmZkWQ3RYsNVRPEJmfDgxYUG6wJkfcqXQ7wSyRCRTRIKB5cB65wYikgb8DbjTGFPo9Hq4iERcfAwsBva7K7xSPak838K+03Vc4efTKp2JCDnpMeTrzBu/0+esG2NMu4g8AGwEbMDTxpgDIrLScXwV8B0gDvidYy2RdmNMLjAKWOd4LRD4szHm70PylSjlZEuh/TqPv0+r7ConPZZ/HKqgprGV2HD9S8df9FnoAYwxG4ANXV5b5fT4PuC+bvoVA9O7vq7UUNtcWEn8yGCmjPHvaZVdXdyIZNfJc3wie5TFadRw0Ttjlc/p6DRsPVrJIp1W+THTUqIIDBAdvvEzWuiVzykoqaW2qU2HbboRGmRjSnIU+Trzxq9ooVc+590jFQQIfrdtoKty02PYU1pLa3un1VHUMNFCr3zO5sJKZqRG67TKHuSkx9DS3snB8nqro6hhooVe+ZSqhhb2ltbpsE0vcv5545Sue+MvtNArn/KvaZU6f74noyJDSYkZoStZ+hEt9MqnvHXYvlrlJWOirI7i0XLSY8g/eQ5jdIEzf6CFXvmMlvYONh+u4JrsRJ1W2Yec9BjO1rdQeu6C1VHUMNBCr3zG9qJqGls7WDxltNVRPF5ueiwAO3Wc3i9ooVc+Y9PBM4QH25g/zj83Ae+PSaMjiBoRxPvF1X03Vl5PC73yCR2dhjcPnuWKSYmEBNqsjuPxAgKEOZmxfHBcz+j9gRZ65RN2nzpHVUMr1+qwjcvmZsZysrqJ8jodp/d1WuiVT9h08CxBNtFplf1w6Vj7ENcHxXpW7+u00CuvZ4xh44EzzB8XT2RokNVxvMbkpEgiQgN1nN4PaKFXXm//6XpOVjex5BIdtukPW4AwV8fp/YIWeuX1Xt1bRpBNuE4Lfb9dOjaO41WNnKlrtjqKGkJa6JVX6+w0vLanjEVZCbqI2QDMzXSM0x/X4RtfpoVeebVdp85RVtfM9dPHWB3FK2WPiSQiRMfpfZ1LhV5ElojIEREpEpFHujn+ORHZ6/jYLiLTXe2r1GCs31NGSGCAbos3QLYA4dJxcWw9WqXr3viwPgu9iNiAx4HrgGzgdhHJ7tLsOHC5MWYa8H1gdT/6KjUg7R2dbNhXztWTExkZ4tL2x6obi7LiKT13gRPVTVZHUUPElTP6OUCRMabYGNMKrAGWOTcwxmw3xlxc8/R9IMXVvkoN1I7iaqoaWvm0DtsMysIs+70HW49WWpxEDRVXCn0yUOL0vNTxWk++ALzR374iskJE8kQkr7JSf+BU3/6aX0pkaKBuMjJIGfHhpMWGsaWwyuooaoi4Uui7W++128E8EbkSe6H/Zn/7GmNWG2NyjTG5CQl6d6PqXV1TG2/sP8MNM5MJDdK1bQZrYVY8O45V6T6yPsqVQl8KpDo9TwHKujYSkWnAk8AyY0x1f/oq1V/r95bR2t7JrbmpfTdWfVqYlUBjawe7ddcpn+RKod8JZIlIpogEA8uB9c4NRCQN+BtwpzGmsD99lRqItXklTE6KZMqYSKuj+IR54+KwBQhbj+rwjS/qs9AbY9qBB4CNwCHgRWPMARFZKSIrHc2+A8QBvxORAhHJ663vEHwdyo8cKq9nb2kdt+amIKI7SblD1IggZqRG6wVZH+XSnDRjzAZgQ5fXVjk9vg+4z9W+Sg3G2rxSgmzCshm9zQlQ/bUoK4FfvVVIVUML8SNDrI6j3EjvjFVe5UJrB3/NL2Fx9mhiw3XJA3e6enIixsDbhyqsjqLcTAu98irr95ymvrmdu+alWx3F50wZE0ly9Ag2HTxjdRTlZlroldcwxvDsjpNMHBXBnMxYq+P4HBFh8ZRRbDlaRWNLu9VxlBtpoVdeY9epWg6U1XPnvHS9CDtEFmePprW9Uy/K+hgt9MprPLvjBBEhgdw4Uy/CDpXZGTFEhwWx6cBZq6MoN9JCr7xCxflmNuwr5zM5KYTrAmZDJtAWwNWTRvGPQ2dp69C7ZH2FFnrlFZ7bcZL2TsPn52dYHcXnLZ4yivrmdt003IdooVcer6m1nefeP8ni7FFkxodbHcfnLcpKICzYxuv7dLUSX6GFXnm8v+aXUtvUxhcXjrU6il8YEWxjyZTRvLa3nOa2DqvjKDfQQq88Wken4cmtx5mZFk1OeozVcfzGTbNSON/czj8O6UVZX6CFXnm0TQfOcKqmiRULx+qUymE0b1wcoyNDWbfrtNVRlBtooVce7Q9bi0mLDWPxlNFWR/ErtgDhhpnJbC6spKqhxeo4apC00CuPlX+yhl2narlvYSa2AD2bH243zUqmo9OwvkAvyno7LfTKY63eUkx0WBA356T03Vi53YRREVySHMna/FKM6XZjOOUltNArj3S8qpFNB89yx9x0woL1BimrfHZOOofK69l5Qnee8mZa6JVHempbMUEBAdw1X1eptNKNM5OJDgvi6W3HrY6iBkELvfI4NY2trM0r5caZySRGhFodx6+NCLZx+5w0Nh08Q0lNk9Vx1AC5VOhFZImIHBGRIhF5pJvjk0Rkh4i0iMjDXY6dEJF9zlsMKtWb598/SUt7J/ctzLQ6igLucqwW+uyOE1ZHUQPUZ6EXERvwOHAdkA3cLiLZXZrVAF8Bft7Dp7nSGDPDGJM7mLDK9zW3dfCn7Se4alIiWaMirI6jgKSoESydmsSanSU06Dr1XsmVM/o5QJExptgY0wqsAZY5NzDGVBhjdgJtQ5BR+ZF1u09T3diqZ/Me5gsLMjnf3K5n9V7KlUKfDJQ4PS91vOYqA2wSkXwRWdFTIxFZISJ5IpJXWambHvijzk7DH7YWc0lyJPPGxlkdRzmZkRrNVZMSeeLdYuou6Pmct3Gl0Hd3p0p/JtVeZoyZhX3o534RWdRdI2PMamNMrjEmNyEhoR+fXvmKd45UUFzZyBd1uQOP9PXFE6i70MZTW4utjqL6yZVCXwqkOj1PAVy+Vc4YU+b4twJYh30oSKmPefq94yRFhbJ0apLVUVQ3poyJ4pNTk3hq23GqdVkEr+JKod8JZIlIpogEA8uB9a58chEJF5GIi4+BxcD+gYZVvuvwmXreK6rmznnpBNl01q+n+to1E7jQ1sHj7xyzOorqhz5vOTTGtIvIA8BGwAY8bYw5ICIrHcdXichoIA+IBDpF5CHsM3TigXWOP8MDgT8bY/4+NF+K8mZ/fO8EoUEB3D47zeooqhfjE0dy2+w0/rTjBLfNTmXiaJ0Z5Q1curfcGLMB2NDltVVOj89gH9Lpqh6YPpiAyvfVNLaybvdpbpqVQkx4sNVxVB/+37UTeWN/Od95ZT9rVlyq11O8gP6NrCz3lw9P0dLeyT2XZVgdRbkgNjyYb1w7iQ+O17B+j65s6Q200CtLtXV08uyOEyzMimeC3iDlNW6bncr0lCh+8Poh6pp0uqWn00KvLLVhXzln61v0bN7L2AKEH944lXONrXz31QNWx1F90EKvLPX0eyfIjA/nigmJVkdR/XRJchT3XzmedbtPs/HAGavjqF5ooVeW2XXqHHtKarl7fgYBuoOUV3rgqvFMGRPJt9bto6ax1eo4qgda6JVlnt52nIjQQN1ByosF2QL4xa3TqbvQxrdf3qc7UXkoLfTKEuV1F3hj/xluy00lPER3kPJmk0ZH8tAnJrBh3xle3VtudRzVDS30yhLP7jiJMYbPz8+wOopygy8tGsuM1Gi+88p+Ks43Wx1HdaGFXg27C60d/OXDU1yTPYrU2DCr4yg3CHQM4Vxo7eDRl3QIx9NooVfD7uWC09Q2tXHPZbrmvC8ZlzCSby6ZxFuHK/jLhyV9d1DDRgu9GlbGGJ557zjZSZHMzYy1Oo5ys7vnZ7AwK57vv3aQ4soGq+MoBy30ali9V1RN4dkG7rksQ9dI8UEBAcLPb5lOSFAAD71QQFtHp9WRFFro1TB75r3jxIUHc/30MVZHUUNkVGQoP75pKntL6/j1P45aHUehhV4NoxNVjbx9pILPXZpOaJDN6jhqCC25JIlbc1P43eYidp6osTqO39NCr4bNH7efIDBAuONSXXPeH3zn+imkxobx0JoC6pt14TMraaFXw6K+uY21eSV8atoYEiNCrY6jhsHIkEB+edsMztQ3891XdOEzK2mhV8NibV4pja0d3KtTKv3KrLQYHrhyPH/bfZpXde16y7hU6EVkiYgcEZEiEXmkm+OTRGSHiLSIyMP96at8X3tHJ3/afoLc9BimpkRZHUcNswevGs/MtGi+tW4fZbUXrI7jl/os9CJiAx4HrsO+D+ztIpLdpVkN8BXg5wPoq3zchv1nOFXTxH0Lx1odRVkg0BbAr26bQXun4esv7qGzU++aHW6unNHPAYqMMcXGmFZgDbDMuYExpsIYsxPoesWlz77Ktxlj+P3mY4xLCGdx9iir4yiLpMeF893rp7CjuJontxVbHcfvuFLokwHn+5lLHa+5wuW+IrJCRPJEJK+ystLFT6883buFlRwqr2fl5eN0zXk/d0tuCkumjOZnG49wsKze6jh+xZVC391vp6t/e7nc1xiz2hiTa4zJTUhIcPHTK0/3+83HSIoKZdkMV88NlK8SEX5001RiwoL56prdNLd1WB3Jb7hS6EuBVKfnKYCrl88H01d5ufyTNXxwvIb7Fo4lOFAneCmICQ/mF7dO52hFA9/TvWaHjSu/fTuBLBHJFJFgYDmw3sXPP5i+ysv9YlMh8SODuX1Oat+Nld9YmJXA/VeO4y8flvDiTl3lcjj0ubWPMaZdRB4ANgI24GljzAERWek4vkpERgN5QCTQKSIPAdnGmPru+g7VF6M8x/aiKrYfq+Y7n8omLFh3kFIf9e/XTGRPSR3ffmU/2WMiuSRZp90OJfHEDQJyc3NNXl6e1THUABljuOn32zlT18w7D1+h69qoblU3tHD9b7chIrzywGXEjwyxOpJXE5F8Y0xud8d04FS53duHK9h9qpavXJ2lRV71KG5kCE/cmUt1Ywsrns3Ti7NDSAu9cqu2jk5+/MZh0uPCuDknxeo4ysNNTYnif2+dwa5TtTzy0l7dgnCIaKFXbvXcjpMcrWjgW0snE2TTHy/Vt6VTk/j6NRN4uaCMX2wqtDqOT9KrZMptqhpa+OU/Clk0IYFr9C5Y1Q8PXDWesroLPPZOEbHhwdy7QBe/cyct9Mptfvr3w1xo7eC/rs/WbQJVv4gIP7hhKuca2/jv1w4SGx7MDTP1Jjt30b+tlVu8X1zNi3ml3Lsgk3EJI62Oo7yQLUD41fIZXDo2lofX7mHzkQqrI/kMLfRq0Bpa2nl47R7S48J46BNZVsdRXiw0yMbqu3KZMCqCLz+/i12nzlkdySdooVeD9sPXD3G69gK/uGW63hylBi0yNIg/3TuHxMgQ7v3jTooqzlsdyetpoVeD8s7hCv7y4SlWLBxLbkas1XGUj0iICOG5e+cSZAvgzqc+1A1LBkkLvRqwkpomHnqhgEmjI/jaNROsjqN8TFpcGH+6Zw4Nze3c+dQHnGtstTqS19JCrwakua2DLz2XjzGGJ+7M0Ttg1ZDIHhPJk5/PpeTcBe75404aW9qtjuSVtNCrfjPG8B9/28ehM/X8evlM0uPCrY6kfNjcsXE8dvtM9pbW8uX/20Vre6fVkbyOFnrVb798s5C/7T7N1z4xgSsnJVodR/mBxVNG8+ObprGlsGTFohUAAA56SURBVJJv/HWPLpXQTzpFQvXL8++f5DdvF7F8dioPXjXe6jjKj9w6O5XKhhZ+tvEIYxNG8pWrdSqvq7TQK5e9uqeM77yyn6snJfKDGy7Ru1/VsPu3K8ZxrKKB/32zkHEJI/nktCSrI3kFHbpRLnl9bzkPvVBAbnosv/3sTAJ1wTJlARHhR5+ZSk56DF9fW8De0lqrI3kF/W1VfXpjXzlfWbObmanRPHPPbL0pSlkqJNDGE3fmED8yhPv+lMeZumarI3k8lwq9iCwRkSMiUiQij3RzXETkN47je0VkltOxEyKyT0QKRES3jfIyL+4s4f4/72J6ShR/vHcO4SFa5JX14keG8NTnZ9PY0s59z+6kqVWnXfamz0IvIjbgceA6IBu4XUSyuzS7DshyfKwAft/l+JXGmBk9bXOlPNMT7x7jGy/tZUFWAs/fN5eRWuSVB5k4OoLffnYmB8vqeXitzsTpjStn9HOAImNMsTGmFVgDLOvSZhnwrLF7H4gWEb1K4qWMMfzojUP86I3DfGpaEk/elavDNcojXTVpFI9eN5kN+87w2NtFVsfxWK4U+mSgxOl5qeM1V9sYYJOI5IvIioEGVcOjvaOTR17axxPvFnPHpWn8evlMggP1Uo7yXPctzOTGmcn84s1CNh04Y3Ucj+TKb3B3c+i6/o3UW5vLjDGzsA/v3C8ii7p9E5EVIpInInmVlZUuxFLu1tzWwf1/3sULeSV85eosvr/sEmwBOoVSeTYR4Uc3TWVaShRfe6GAwrO62mVXrhT6UiDV6XkKUOZqG2PMxX8rgHXYh4I+xhiz2hiTa4zJTUhIcC29cpvzzW3c88xONh44y39dn82/XzNB58krrxEaZGP1nbmEhQTyxWfzqG3SBdCcuVLodwJZIpIpIsHAcmB9lzbrgbscs28uBeqMMeUiEi4iEQAiEg4sBva7Mb9yg+qGFj77hw/48EQNv7xtOvdcpvt1Ku8zOiqUVXfkUF7bzIN/2U17h66Jc1Gfhd4Y0w48AGwEDgEvGmMOiMhKEVnpaLYBKAaKgD8A/+Z4fRSwTUT2AB8Crxtj/u7mr0ENwunaC9yyageFZ8/zh7tyuHFmitWRlBqwnPQYfnDDJWw9WsV/v3ZQZ+I4uDSVwhizAXsxd35tldNjA9zfTb9iYPogM6ohUlRxnjuf+pCGlnae+8Jc5mTqxiHK+906O5WiygZWbykmLTaM+xaOtTqS5XTOnJ8qKKnlnmc+xBYQwAsr5pE9JtLqSEq5zSNLJlFS08QPNxwiOXoE103179neOm/OD207WsVn//A+I0MDeenLWuSV7wkIEH552wxmpEbz1TUFbDtaZXUkS2mh9zMb9pVz7x93khYbxksr5+umIcpnhQbZeObu2WTGh/PFZ/PIP1ljdSTLaKH3I3/+4BT3/3kXU1OieGHFPBIjQ62OpNSQig4L5rn75jA6KpS7n9nJ7lPnrI5kCS30fsAYwy/fLOQ/1u3j8gkJPP+FuUSFBVkdS6lhkRgRyvP3zSU2PJjPPfkB24/53zCOFnof19bRyTf+updfv3WUm3NS+MNduYwI1o28lX9Jjh7B2i/NIyVmBHc/s9PvlkrQQu/Dzje3ce8fd7I2v5SvXp3Fz26eRpBuGKL8VGJkKC+smMfkpEi+9Hw+q9495jfz7PW33kedrW/m1ifeZ/uxan76mWl8TZc0UIqY8GBeWHEpn5yaxI/fOMzXX9zjF2vZ6zx6H7Tr1Dn+7fldnG9u4+m7Z3P5BF07SKmLQoNs/Pb2mUwYFcEv/1HIvtN1PPbZWUwcHWF1tCGjZ/Q+xBjDc++f5LYndhAUKKxdOV+LvFLdEBG+cnUWz907l3NNbXz6sW08te04HZ2+OZSjhd5HVJ5v4UvP5fOfL+9nwfh4Xntgod4IpVQfFmTFs+GrC5g/Lo7vv3aQm1dt56gPLnOshd7LdXYaXt59mmt/tYXNhZX8x9JJPPX52Tp9UikXJUaE8vTds/nVbTM4UdXIJ3+zjd+8dZTWdt9Z/VLH6L3Yh8dr+OGGQ+wpqWV6ShQ/v2U6WaN8d5xRqaEiItwwM5kFWfF879WD/O+bhWzYV853Pz2FS8fGWR1v0MQTpxfl5uaavLw8q2N4pJb2DjYeOMsz7x1n96laRkeG8vC1E7lpZjIBuhuUUm7x5sGzfHf9AU7XXmDJlNE8unSSxy8XIiL5xpjcbo9pofd8Z+ub+eB4DW8dOsvbhyo439JOZnw4d8/P4NbcVL0BSqkh0NzWwZNbi/nd5mO0dxjuWZDBA1eOJyLUM4dFtdB7ibqmNk7WNHKiuoljFQ0cPlPPwfJ6SmouABATFsTi7NEsnZbEwvHxegav1DA4W9/MzzYe4a/5pcSEBfHFRWO5a14GI0M8a+RbC72HMMZQ2dDCyeomx0fjv/6taaK2qe2fbUUgMy6cSUkRzEqLYU5mLNlJkQTqna1KWWJfaR2/ePMIm49UEhMWxB2XpnPHpemM8pDFAbXQD6P2jk7K65o5VdPEiepGTlXb/z1Z3cSpmiaaWjv+2TZAIDlmBOmx4aTHhTk+HI9jw3VIRikPVFBSy2NvF/HW4bPYRLh2ymhunJnMogkJBAdadyI26EIvIkuAXwM24EljzI+7HBfH8aVAE3C3MWaXK32746mF3hhDdWMr5bXNnK69QHndBcrrHI9r7Y/P1jfjfM9FsC2A1NgRZMSFkxYX9pF/k6NHWPqDoZQauBNVjTy74yQvF5ymprGVqBFBLMyKZ9GEBOZkxJIeFzasy44MqtCLiA0oBK4BSoGdwO3GmINObZYCD2Iv9HOBXxtj5rrStzvuKPTGGDo6DR3GYAx0dBo6jaGzE9o7O2lq7eBCWwdNrR00tbTb/23roLGlnZrGVmoaW6luaKG6sZXqBsfzxhbaOj76/QoODGBMVChjokeQFDWC5Gj74zTH2fnoyFBsOpaulM9q6+hk69FKXt97hq1HK6k43wJAZGggk5Mi/3lyFz8ymOiwYKJHBBETHkxkaBDBgQEE2YQgWwDBtoBBXXfrrdC7cjVhDlDk2OgbEVkDLAOci/Uy4FnHJuHvi0i0iCQBGS70dZup391IU2uHW25jDgu2ETcymNjwEEZHhTJlTCSxI4MZHRnqKOgjSIoOJS48WBcLU8qPBdkCuGrSKK6aNApjDEfOnqfgVC17T9dxuLyetw6fpaqh1aXPlRARws5vfcLtGV0p9MlAidPzUuxn7X21SXaxLwAisgJY4XjaICJHesgTD3jizgGemgs8N5un5gLPzeapuUCzDcRHcp0E5NsD/lzpPR1wpdB3d7ra9ZS5pzau9LW/aMxqYHWfYUTyevrzxEqemgs8N5un5gLPzeapuUCzDcRw5XKl0JcCqU7PU4AyF9sEu9BXKaXUEHJlysdOIEtEMkUkGFgOrO/SZj1wl9hdCtQZY8pd7KuUUmoI9XlGb4xpF5EHgI3Yp0g+bYw5ICIrHcdXARuwz7gpwj698p7e+g4yc5/DOxbx1Fzgudk8NRd4bjZPzQWabSCGJZdH3jCllFLKffRuHaWU8nFa6JVSysd5RaEXkVtE5ICIdIpIbpdj00Rkh+P4PhEZ1hWGesomIhkickFEChwfqzwhl9PxNBFpEJGHhzNXb9lEZI7T92uPiNzoIbmuEZF8x89XvohcNZy5+sgWJyLvOP5bPjbcuXrL5jj2qIgUicgREbnWinyOHNMddWKfiLwqIh6zz6aIzBCR9x0/93kiMsftb2KM8fgPYDIwEdgM5Dq9HgjsBaY7nscBNg/JlgHs97TvmdPxl4C1wMOekg0IAwIdj5OAiovPLc41ExjjeHwJcNqDvmfhwAJgJfCYJ/2sAdnAHiAEyASODffvp1OWncDljsf3At+3IkcP2TYB1zkeLwU2u/s9PGtB5R4YYw4B3S01sBjYa4zZ42hXPczRestmqd5yicgNQDHQOMyxgJ6zGWOanJ6G0sPNdUOll1y7nZ4eAEJFJMQY0+IB2RqBbSIyfriydNXLz9oyYI3j+3RcRIqwL6myY3gTAvb/EW1xPH4T+0zA/7QgR3cMcPEvjCiG4F4jrxi66cUEwIjIRhHZJSLfsDpQF5kisltE3hWRhVaHARCRcOCbwPesztIdEZkrIgeAfcBKY0y71Zm6+AyweziLvBfraWkUK+wHPu14fAsfvZHTag8BPxOREuDnwKPufgOPOaMXkX8Ao7s59C1jzCs9dAvE/mfrbOzz999yrOD2lgdkKwfSjDHVIpIDvCwiU4wx9Rbn+h7wS2NMw1D+FTLAbBhjPgCmiMhk4E8i8oYxptnqXI6+U4CfYP9L0u0Gk22oDTCby0uguENvGbEP1/xGRL6D/aZN11YZG55sVwNfM8a8JCK3Ak8Bbl3ZzGMKvTFmIF9YKfCuMaYKQEQ2ALMAtxb6gWRznPG1OB7ni8gx7H+BuG2h/QF+z+YCN4vIT4FooFNEmo0xbr2QN8Bszv0PiUgj9jFxq79niEgKsA64yxhzzF15nA32ezaUBvH7OWxLoLiQcTGAiEwAPjlUObrTWzYReRb4quPpWuBJd7+/tw/dbASmiUiYiAQClzNESyD3l4gkiH09fkRkLJCFfVzcUsaYhcaYDGNMBvAr4H/cXeQHSuxLZQQ6HqdjH1c9YWkoQESigdeBR40x71mdx4usB5aLSIiIZGL/HfjQiiAikuj4NwD4NjCss+D6UIa9dgFcBRx1+ztYfcXZxavSN2I/O2gBzgIbnY7dgf0C2X7gp56SDftY7gHssw52Add7Qq4ubb6LNbNuevqe3en4nhU4vmc3eEiub2O/cF3g9JHoCdkcx04ANUCDo022B2X7FvbZNkdwzCyx4gP7GXOh4+PHOFYF8IQP7MPP+Y5a8QGQ4+730CUQlFLKx3n70I1SSqk+aKFXSikfp4VeKaV8nBZ6pZTycVrolVLKx2mhV0opH6eFXimlfNz/B/lVeJVMvu+7AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"kdeplot(trace['center'])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"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.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment