Skip to content

Instantly share code, notes, and snippets.

@Sayam753
Created July 7, 2020 08:33
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save Sayam753/1a014bbc1afcf4dea0bb5e946e2e103f to your computer and use it in GitHub Desktop.
Save Sayam753/1a014bbc1afcf4dea0bb5e946e2e103f to your computer and use it in GitHub Desktop.
Testing transformations in PyMC4
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import arviz as az\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import pymc3 as pm3\n",
"import pymc4 as pm4\n",
"import seaborn as sns\n",
"import tensorflow as tf\n",
"import tensorflow_probability as tfp\n",
"\n",
"tfd = tfp.distributions\n",
"tfb = tfp.bijectors"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"experiments = np.array([1., 1., 1., 0., 0., 0., 0., 0., 0., 0., 0.])"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"@pm4.model\n",
"def model():\n",
" prob = yield pm4.Uniform('p', 0., 1.)\n",
" ll = yield pm4.Bernoulli('ll', prob, observed=experiments)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"WARNING:tensorflow:From /usr/local/lib/python3.7/site-packages/tensorflow_probability/python/math/minimize.py:77: calling <lambda> (from tensorflow_probability.python.vi.optimization) with loss is deprecated and will be removed after 2020-07-01.\n",
"Instructions for updating:\n",
"The signature for `trace_fn`s passed to `minimize` has changed. Trace functions now take a single `traceable_quantities` argument, which is a `tfp.math.MinimizeTraceableQuantities` namedtuple containing `traceable_quantities.loss`, `traceable_quantities.gradients`, etc. Please update your `trace_fn` definition.\n"
]
}
],
"source": [
"trace = pm4.fit(model())"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3dd3gU1foH8O9JI6EFQgmdJXQQgRCBiEgVUFDvtaPXdr0i9nJFg/qzcVW8eu2ocEWsqIidIFx6b4kIBAIkhAQSQghJCCGQfn5/7Oxme0l2d3Znv5/n4WFn5uzMmUzy7tkz75wjpJQgIiJtCVG7AkRE5HkM7kREGsTgTkSkQQzuREQaxOBORKRBYWpXAADatm0rdTqd2tUgIgooqampp6WU7Wxt84vgrtPpkJKSonY1iIgCihAix942dssQEWkQgzsRkQYxuBMRaRCDOxGRBjG4ExFpEIM7EZEGMbgTEWmQpoJ7ak4x0vPPql0NIiLVaSq4X//RNlz57iYcKTyHZXtPqF0dIiLV+MUTqp424T8bAADTLu6kck2IiNQR0C33qpo6PLj4Dxw9XW63TNG5SnC2KSIKNgEd3HdlFyN5bz6e/Wmfze2Zp8ow7F+r8eV2u8MvEBFpUkAHd2eyCvUt+o2HT6tcEyIi39J0cCciClZBHdwzCspwsrRC7WoQEXmcJrNlXHXF2xsBANlzp6pcEyIiz9J0y31F2km1q0BEpApVg7sQ4mohxILS0lKv7P/H3XlW62rrJNLyvHM8IiJ/oWpwl1L+JqWcER0d7bNjfrA2E9Pe34y9uWd8dkwiIl/TRLeMlHD5QaV9SqudN1KJSMsCOrgLk9c9Zi936718ZpWItCygg7urAXp1eoHxtRAOChIRaURAB3cDBmwiInOaCO4N8faqw2pXgYjIa4I2uB88WaZ2FYiIvCbogntGAYM6EWlfQAf3hsy2lF103gs1ISLyLwEd3L/ZeRyAPs+9Mf48zgeaiEhbAjq4e8pf5m1RuwpERB6lieCeX3rBY/vKPHXOY/siIlKLJoK7p/rRf9+Xj4lvbcCKtHyP7I+ISC2aCO6ekq6kRzJNkogCHYO7A0tTc6FLSkZlTa3aVSEickvQzMS06kABKqpdC9IXqmuRnn8Wb6w8CAAoKa9Gh+hQb1aPiMijgia43/tFistl52/IwvwNWV6sDRGRd7FbhohIgxjcXSAdDC78x7ES5BSV+7A2RETOBU23jLdc9+FWAED23Kkq14SIqB5b7qbsjGOwN7fUow9KERF5G1vuLrjvy1SECCDrNbbOiSgwsOXuojpOukpEAYTBnYhIgxjc3aBLSsbZimoAwMnSCjy0+A+Va0REZBuDu5uOKYOU/Sv5AJbt5QBjROSfGNwVu7KLsSr9lMvl2QVPRP6M2TKKGz/e5lZ52djpn4iIvIgtdzcJof+fsZ2I/BmDewNZBvfKmlq8s/qwyyNPEhF5E4N7A1mON/P51my8szoDn2ziaJJEpD4G9wYoLq9CZU2d2bqK6jqz/w10Scl4JfmAz+pGRAQwuDdI/JxVWH+o0Gyd0hVvcwTJ/2466oNaERHVY3D3EFdutFbX1qGW4xgQkQ+oGtyFEFcLIRaUlpaqWQ23CGMb3WK9Et1/3p1nN02y97O/49p5m71WNyIiA1WDu5TyNynljOjoaDWr4RZHE3cAwInSCqxIO2l3e1reWU9XiYjICrtlPESYNOjPXKhWryJERGBwd5vdbhk764mI1MDg7iFpJ+rvGxjCfGWN7QeadEnJ+HJ7jg9qRUTBisHdQ5ItRoisrZPo+9wK43J5ZY3Z9s+2MD2SiLyHwd1LaurMH2b6cH2mSjUhomDE4O6mpam5DXpfda15lk1Ds90PF5ThZGlFA99NRMGCQ/666VOVu1Mmvb0RAJA9l5N1E5F9bLl7QZ0EdmQVq10NIgpibLl7wTM/7XNeiKMQEJEXseVORKRBDO4+UmUxRDARkTcxuPvIZ1uzzZZNe2XS8koxfcF2uw892XLbJ+6VJ6LgwuCusiOF53DT/G3YllWEAydcH1RsS2YRDuaXebFmRBTIGNxVNuE/G3C+qr4FXlFdi1nf78Hpc5WorKnlnKxE1CDMllHJ0dPl2Hi40Gr9L3/m4XvlQamNGYUoOFvJnHYichtb7iq649OddrcJARScrfRhbYhISxjciYg0iMHdzxhm6FuS0rAxbIiIAAZ3IiJNYnD3IxeqalHJh52IyAOYLeNHbv1kh1vlBWf2IyI72HLXkDPnq1Bdy5Y/ETG4a8qQl1fhocV/qF0NIvIDDO4BTNoYNnjl/gLfV4SI/A773DUgp6gcURGhaleDiPwIW+4BQJeUjM8tRpU0NeaN9Rj+yhrfVYiI/B6De4B44df9aleBiAIIg3sA0SUlmy0LAeSWnFepNkTkzxjcA9z4NzeoXQUi8kMM7gGuinntRGQDgzsRkQYxuAeJF3/db9VnT0TaxeCuQVJKvLs6A4Vl9ZN9WE7QTUTaxoeYAlh5pe35VXvMXg4ASMkpxpf3jPBllYjIT7DlHsC+TznucDsn1yYKXgzuAezH3XkOt9sae4aIggODOxGRBjG4a1hZRY3aVSAilTC4a1hlDfvciYIVgzsRkQYxuBMRaRCDu4ZlF3HESKJgxeCucWfOV6ldBSJSAYO7xi3eeUztKhCRChjcNY4PMhEFJwZ3IiINYnDXuDdWHuJUfERBiME9CFz2+jq1q0BEPsbgHmTWHTqldhWIyAcY3IPM3Yt2oaTcOj3yp925OFxQpkKNiMgbGNyDkK1JtR//bg8mvb3RpfdvzyrC26sOu3XMlOxi1NYxdYfIVxjcg5BpkD1VVoFPNmW59f5bFmzHu2syXC6/K7sYN3y8DR+szXTrOETUcJxmLwjtOFqEnUeL8UNqns1WPAAcLihDbMtIREeFAwDWHizARZ2ikX6yvutm/H/WY/kjoxEZHurweCdLK/T7PMVuHyJfYXAPQo9/t8dpmUlvb0Sf2Ob43+NjAAB//ywF3ds0RfsWTYxlsgrLkVFwDoO6RDvcFztjiHzP48FdCBEH4FkA0VLKGzy9f/KdwwXnzJZzis7jWHHDc+ZFYytERC5zqc9dCPGpEOKUECLNYv0UIcQhIUSmECIJAKSUWVLKe7xRWfK+onOVKDpXabZOmoxhwOEMiAKDqzdUPwMwxXSFECIUwDwAVwIYAGC6EGKAR2tHPjfsX6sx7F+rjctHT5fjQrXjGZ22HSnCle9uMpv56b01GThSeM7Bu4jIm1wK7lLKjQCKLVYPB5CptNSrAHwL4FoP149UNu7N9Rjw/EqHZZ77eR/S88/iuNJlc+Z8Fd5adRjTF2z3RRWJyIbGpEJ2BnDcZDkXQGchRBshxMcAhgohZtt7sxBihhAiRQiRUlhY2IhqkJqEAI4UlgMA/jxeCqC+68aQiSOd9OVMX7AdX2zL9lYViYKSx/PcpZRFUsqZUsqeUsrXHJRbIKVMkFImtGvXztPVIBU8+b3tLJxHv/0TAFBhp3tnW1YRnv9lv9fqRRSMGhPc8wB0NVnuoqyjIGbaFWPZYF+dznFtiHylMcF9F4DeQogeQogIALcA+NUz1aJA8dTSvWbL27KKUF1n+8EoIvIdV1MhvwGwDUBfIUSuEOIeKWUNgIcArASQDmCJlJLfrYPMgfyzVuseWrwbgL4/nojU4dJDTFLK6XbWLwew3KM1ooC386g+scpZTvz0BduxM9syCYuIPIHDD5DX1NVJfL412+72bVlFvqsMUZDhqJDkNWWVNXjhV/bUEamBwZ1UMentDVbrdEnJyC05j4rqWqshEIjIPap2ywghrgZwda9evdSsBvmYlNJqUDKDrUeK8N2u40jNKcH/TRuA3JLzuGpQR1yii/FxLYkCm6otdynlb1LKGdHRjoeMJW1xdKNVAEjNKQEAzFl2AIu2ZOPGj7eZldmUUWhzqkB/sDnjNJbvy/fqMVak5UOXlIyzFdVePQ4FNnbLkE+tOlCAjFP2BxQTTvInK2tqcfvCnbj90x1OyzlzoaoWx4oaPoSxLX9buAMPfP2HR/dp6YN1+hmtsk+Xe/U4FNgY3Mmn7v0iBZPfsT9Xq7PUeEOr//BJ+x8Qu4+VoO9zK7D+kOMnYu/5fBcuf2Od2brPt2YjLa/USS3UcaGqFrqkZKTlWT9b4GtSStTYmcWL/AODO/kVew13XVIyBr24Evd/lQoAqHHwFOwuJXd+c8Zps/Ul5VW4ef42FJzVT/u39Yh1KuYLv+7HtPc3N6TqXpdfesFs2bJ762xFNYa+/D/s8EGK6Uu/HUCvZ39HHSc991sM7uRXLAOyqbKKGqw7pB9B1DSmnCytwGkb2TWWHxRLUo5jx9FiLNx81Gz9p5uP4ufd7g2LtCmjEJkmc8LW1UnM+CLFuLz2YAF+2p0LXVIyzlXWuLVve0JDzE/IMqzuyy1FyflqtyYvb6jPt2XbrIOprZmncaqswut1IdsY3Mmv/OhmkAWAka+tQYLJBCOuKDO5GfnysgN47Ls/7ZZ95qd9+O/GLONySXkVbl+4ExPfqu9eKquowf8OFBiX//5ZinGuWk/164c4uR9haMm7OuxDXZ2ELinZ4XDLxeVVGDV3LUrPu3/z9tZPduDaD7YYl7MKz+FClfN7IeQZDO4U0G78eKvZ8vHi85j7+0Grck9+vwcfbTgCQP/tYNCL/3P5GIt3HMMry9ONyy8vO2BdyEFAnb/xiM31UkqsSMtHtZ2+679+uAW6pGTj8hNLzD+ALMfJl0o7Wrg4W61hvP3nf9lv9u1i+b585BTpb9bGz1mFvDMX8Oh3u83eazjCrKV7kJZXavcc8kv1Lfea2jqM/88GPPB1qlWZX/ecwDYbXWTUOAzuFNB2ZZcYX9fWSYz+9zpjl41p5s3S1FycUVqftgY7s7Rw81Gz1r3BpoxCZJh0xxg4ai3/8ucJs2VDGufq9FOY+dUfmKdkv1jafeyM2bLpuQLmXSKFZZVm2TPrD53C4QLrepoyTaW86IWVxhTUB77+A1Pe2WRWtqzCvGvJcOwf/8jDtPc3m32gPv9LmlWwNlyTTSbdbj/vzsPx4vN45JvdmP7f7aiprcPK/SedTu5CruHYMhSwxr+53mw5xWIQsuLyKrOWrzvmLDuAORYt9NScYty+cKfN8s7aynlnLqBzqyjsPlaCv364Fe9PH2oMrgVnre8XuBLgfttzAgM6tkRkeCgueaW+W0oI4K5FuwAA2XOn2n2/ZcbR9R9tNZa3nDfXWX325eozjLIKz+GLbTn4YluOzXKGvUgp8dh3f6Jt8ybGbR+sy8Q7qzOw8M4ETOgfa7/eBWUIDw1Bj7bNzNYvTc1FVHgopl7cEQBw6mwF2reMdFhvWzZlFOJY8XncNqK707Lnq2pw9kINOkS7fxxvU7XlLoS4WgixoLTUP1PPyL9lWeR5z/zK/Cv/0tRct/Z30/xtDrdf/5H19tyS89AlJWPtQcdpl6PmrgUA7FPSLHccLUJNrT7UfbPzmLE1vyLtJGb/uM/shvG187Ygy8Zk44u2ZFt9ANmjS0rGw9/Ud62UlFfZnPi8qsZ294qhOtPe34QRr662ytQxfHMZ/x/rYSUAoM9zvwPQf7sqLq8yvt/0RnhuiT4bqMjJA2qT3t6IcRYf7IC+6+3BxfpnDNYdPIXhr67BOifpsAafbTmKeesyUVsncfvCnXj2pzTklljfK6mrk/htzwljltD0/+7AyNfWON3/7B/3QpeUjGd+2of9J0ptXk9P4xOqpBklDbjpZ8owVLE7DN0M36c4/yCprKnF+2v1XTChQph1+zyvDLA286tUfLPzmFmq557jZ/DOatsZMCfOXEBhmf1xeGb/uNeY+vnbnhO4UFWLL7fnYOicVbjXJLvHwFAWgFkeu5T6wJyWd9bmNw1H3VK7j5l3J8XPWeUwy8YT9uSeUY5d37WVU1SONekFmLcuE1e+W9/t9EryAbz42wG8sfIQej5TP4L5zfO3I/PUOSxNzUXmqXP48/gZfLThCB7+Zje+3qH/ZrLnuHnXmamtR05Dl5SMVQcK8M1O/XTTi3ccw9T3Ntv9EPQk4Q/9WwkJCTIlxfoXzZmGfuUmUkNkeAgqqm23jMf2bYdFd12CHrP1weXaIZ2s+uobqmVkGM5WuJaOufqJMZj4lvcDz1+GdMLPHjo/Z/71l4tw24huiHtmudk3ju2zJ+CvH24x3vR118E5U9Dv/1YAABK6t0ZKTglmjumJv43shj+PnzFOWmPP6ifGoFf75g06toEQIlVKmWBzG4M7EZF6Nj01Dl1jmjbovY6CO7NliIhUdNxG374nMLgTEakoMa6NV/bL4E7kAkcphbbKGv7Z22ZveUDHljj62lXG5fSXp1i9f8HtwwAAE/u3NzvWLw+Owru3DLE6bvbcqdj9f1e4fH6bnx6Hg3OmIO2lyVj9xBhsfnqccdt3M0Yie+5UxDSLAAD8+/qLkT13KmZN7mu1z4n9Y5E9dypWPX6502Oa+ui2eHx4WzwA4MlJfRAWYvtubbsWTWyud9aP/eiE3sbXCd1b2y13c0JXm+tjmkWgm41ulOenDXB4XADooKRmThqgT/WMCAtxOhJqQzHPnYLShlljMeaN9QCAIV1bYXy/9rh+WBdc+8FmnD5XhTsSu9vN1XbHwjsT0KppBDYcLkTLSPt/bhFhIaiqqUPyI5dBCIH1T45Fm+YRiIoItSrbprk+sPa0CGKDu7bC4K6tbO6/VdNwjIyLQeapcpw+V4mkK/vh7lE69H1uhVXZLq3rA5chUG6cNQ7J+/IxQmllNo0IRXE5cIUSpB4c1wsPjuuFd1dnIKeoHD/uzsP704fq69muOe66VIfPbMyn+8kdCejWpil2HC1GdFQ4rhncybjNEPwfGt8b13ywGXtzS/Hzg6NQcr4KPds2R7c2+noa7r1lvHIlKqpr0SIyHOn5Z5FTVI6ZX1kPv2waS01fZ8+dij+OleC6D/VPPY/q3RZ3jdJh3rpMLNubj1mT+6JfhxbGHHzTe36L7r7E7CGyN28cjKd/2IsHx/XCe8pYP+P7tceHt8VDSv3N9ddXHMI0JSffGxjcSRNGxsVg3q3xGGZnjJk7E7tDCIG/DO2Mf684iK4mAeznB0dZlR/Vq63d4L7g9mEY3iMGb686jM+VMv06tMB/70jAwZPmT4UaAsEwixbi3aN0Zrnxyx6+DJsyThtbcTqTB3S+nTESST/sxaAurZR9xeDLe4ZjpIOv88sevgzT3t+MtsoHgRAC385IRFlFNf694hDuTNShSZj1B4c93do0xf1jexqXv7pnBFbsP4nWSgve4NGJ+lbxWzcPMa4LCRF48ZqBxuD+91E98LeR3QAAE5UPhz6xLRwe/65LdXhiyR7o2jTFEIsPsKem9MWBE2cRHhqC8FB9Z0T/ji3Rv2NLm/syHaPHcqiG+G6tMfXijkjem2/cj6F8l9ZRdh+u6hQdhaOF5ca63jCsC24Y1gUAjMH9igGxiAyv/5knXdnP4Tk3FoM7BZR3bxmCR7/Vj7Hy+vWDkJ5fZgwabZo3wZp/jsEEGznEL117kfH14ntH2t2/o+SxuLbNkHW6HJMGdtCXVdbfdakO/5zUBy0iw13Oenjh6oF44eqBxuU+sS3sBriRcW2wftY4s3Wje7dzuP/2xi4L8+DVIjIcc/5S/7NY8dhot4K8ga5tM8wc09N5QRObnhqH0BCBTq2i3D7edfFdcF18F5vbHhhrf5rOT+5IwD8s8vnt9PI0inQhc987nS/2sc+dVPP0lPqWy1QbX09bWHRjHH3tKlw7pLOx5XbzJd0waaB5S6pnu+Z4bmr/BtfJ0Z/oiscux8E59X3ghg+CHm2boUVkeIOP6Q2G83AWyPp1aGn1GL+3dI1p2qDA3hgTB8Qivpt5S79jtEkdHPx8PB2MvdS1bheHHyCfsHdzyiDSVuvR6hF3/V/H4ntHYMczE+zua3y/9m7Xz3hIB033iLAQs6/VNynn1JjjeYu7w/9q2aK7hpt1KV0X3xnTh+u7hYZ2s32PwhZH3+ps/v4qDN0z3rpxao9mhh+Y4Id/YKQX0ywCr103yGEZW7/3t47oZrNs04gwxCpZB4ZW2Ige3kknA4BfbPTJA8CgLtHInju1wQ+geJO7w/9qWXTTcNw+Uj8I2HvTh0IIgdeuG4TfHx2NWZOss3wsGxWX99F3gfWONb+BHR0VjhABzLs1Hrq2zdCvg75bzfKegFo00+c+oX8s1jgZvIm8477L4zDfZDILSx9MH4qQEIGnp/TD6ysOIq5dM2QVmg/6JaC/sZmWV4ohXVshJbsYd16qwz8n9TUOOmVLj7bNsHHWOHRpXf9V27SFdP/YnsaWkyuGdW+N1emn0FTJUhndu63dDBR/1r5FJCYPjMWMy93rF9eqTq2irFIv7d1wNTD8Gt0wrAsm9GtvdfM49bmJAIAw5Sbupb3aNuppU0/TTHAn9TwxqQ/atWiCc5U1Nge4urRXWwDAzDFxuPPS7mgaof+1+2h9/SQWQuhbPIZWz0Wd9d/mIsKctzwNKXG23Dq8m1t/bO9NH4qswnIM7NQS/7yiD24Zbvvbg78LDRGYf7vNp9LJCVs3Ry0DO1Af1E35S2AHeEOV7LhmcCdkvXoVnrmq/qbnJbr6dL4nruhjfN0kLBT/GB2Hxyb2gSNCCGNgB7z39dXwoMjdo3Ru/7E1jQjDRZ2jIYTAwxN6231QhrQv0Lu02HInm94zeQDF4NW/DsJdi3Yh78wFXBffGS0iw5DY035f99zrBjnMjkjs2cbYVePJP6SoiFCXnii9YVgXTOzPezXB7rmp/bHbZOhePxhL0SMY3MllQgADOrVE3pkLiAwPxd2jetgt2zQi1KUujVZN1UshfPPGwaodm/zHP0bH2Vwf6JlGmumWCfQLoTZXx0555+YhWDoz0Wx6NFOzlafuXB0M6apBHTFcF4OHxtt/EOWNGy5uVO46kTu00nLXTHDXygXxhsFdou0O3uSM+c9VoFmTMCToYuyWNzwQ42pOb3RUOJbMTHTYN35jQle7rSsiTzM8GNe3g+MhEVxlSOkdpSQW+Ipmgjs5IAR6WzzavvfFSS69lZ+ZFGyui++Cg3OmmN1vaoz4bq2RPXcqOvv46Vz2uWvcsO6tbQ5F2tLFx+XdnamLHwakBaZPIgcqttz91B2J3V0ue1Fn2w9jNAkLwQ/3X+qxh3BiW7qeFsh7IETq0kxw11owucnJWCym3B3Vz3RUQGcMHwyL7r7EpcGxeO8jcCVd2Q//uMx+BhQFFs0Ed60RAnjEZMYYR561k0li7wPPMM6GK2JbRiJ77lSM6+tePrjGPmuDwswxPfGcC7MJUWDgqJB+qlXTCDxxRR+XUhTju9meKizQn7AjoobTzKiQWhLXtplP76z/cH8iABhn7Wkc9ssQ+QN2y/ghZ6PVGdyR2B27ntWPTLfkvkTj+tG99fm0zu5DGIL5sO76vPXNT493t6pWOI44kX/QTHDXVCxxcjKGiQeio8KNA1sN71H/YNECZTTAuHaOZ9hZ88+x2JJUH9Ajw0ONg241VKgy9U9Dpm4jIs9hnnsjLbkvETfN3+aRfU0Z2AEr9p+0uz08VGDjU+NQUyvx9fYcu+OUR0WEYtFdl+DiLo67u6KjwhEdZZ4Bs+Gpsairc7/uBhP6x+KBsT1xL58oJVIVg3sjmU7TteOZCRjx6poG72tcv3YoPl+FJ23NDgMgLCTEOPPQ3hcnO9lXw0Y7bGyLOzRE4Kkp3p3VnYic00y3jKe50z2R0F2frRLbyC6NZk3CsOS+RJ9NWExE2hXQwX1L0ngMUwKrKzfwHncymYSpGBszr9giAHx97wjsecG1sVoawtFk0O5YeGcCpgzs4JF9EZF/C+jg3rlVFHo6uWloqlkT210O22ZbZ4kYPiwiw53/iJqEhVr1XVuKa0RrvEWkZ3rPJvSPxce3D/PIvojIvwV0cPcUQz+2qYY+Rv/69YOs1k3s3x7LHrnMuGzIdnFVuDJXozsTPRNRcGNwt8NWbO/XoQU2zhqH5Y+Mtvu+mxK64vuZiWbrYppFmM0dam9QMHsfKOGhITjw8mS8dM1Ap/UmIgI0kC0TpQzNGW5jJvLGMHThmwbcqwZ1RLc25pNKWE5KIYTAJQ4mswBsf1NwxvTDgYjImYBvuc+a0g+PjO+FawZ3UrsqRER+I+CDe/MmYXhiUl+EebrlrjTIGzpSygtXqzO6nqtZPkSkbQEf3L3FECQfdXHYXUvOZnL58LZ4JJvcZAUaPx5L+stTsDWp8ePDEFHg45C/dlwxIBbZc6fi5ktcnzTDmZ8euBTvTx8KQN9/P7BTtEcnR4iKCNXE9GBE1Hgc8tcOw4QWzlIi7TW2p13c0Wrd0G6tcbXFvYHnpg3A1EHWZYmIGiOoumXcyV03ZME0tKukRWS4zZx3RzhFHRF5SlAFd0uuTIjRxhc3KDU1XjER+YOgTp5e9+RYSCf5MEIIPDy+F95fm+mjWhERNV5QB/eIsMZ/cXHUbTO4q3444An9Yxt9HCIidwR1cPe2fh1aIuOVKz3+9CwRkTNBFXWcdcF4AwM7EakhaCNPiBs3MY2ZM16qCxGRpwVlcL+oc0tsmDXO5fIzLo/DrSO64e8efODI1GRlAo3+HVt6Zf9EFHyCss89Ma4NusY0dV5Q0bxJGF79q+2cdctRIRvimsGdMHlgbKPnLyUiMgiqlrs/PyTEwE5EnhRUwd3AE61tIiJ/FpTBnYhI6xjcbWDDnogCXVAF9zF92wGoz06x5dsZI7HpKdczaYiI/FFQZcv069AS2XOnOiwzMq6Nj2pDROQ9mm+5u/Owkj/sl4jIEzTVcl987wjENIvAlHc2AQCWPzIabZpHYMSrazx+rN8evgxr0095fL9ERJ6gqeB+ac+2ZssDOnnvic+BnaIxsJP/zSBFRARwDlUiIk3iHKpERBqk+RuqRETBiMGdiEiDNBncrxncSe0qEBGpSpPB/e2bhyD95SlqV4OISDWaDO6hIQJRERxCl4iCl6by3Bujb2wLr+bFExH5EoO7YuXjl6tdBSIij9FktwwRUbBjcCci0iAGdyIiDWJwJyLSIAZ3IiINYnAnItKgoEiFXHJfInKKyivPtM8AAAU4SURBVNWuBhGRzwRFcB/eIwbDe8SoXQ0iIp9htwwRkQYxuBMRaRCDOxGRBjG4ExFpEIM7EZEGMbgTEWkQgzsRkQapGtyFEFcLIRaUlpaqWQ0iIs0RUkq16wAhRCGAnAa+vS2A0x6sTiDgOQcHnrP2NfZ8u0sp29na4BfBvTGEEClSygS16+FLPOfgwHPWPm+eL/vciYg0iMGdiEiDtBDcF6hdARXwnIMDz1n7vHa+Ad/nTkRE1rTQciciIgsM7kREGhTQwV0IMUUIcUgIkSmESFK7Pg0lhOgqhFgnhDgghNgvhHhUWR8jhFglhMhQ/m+trBdCiPeU894rhIg32dedSvkMIcSdap2Tq4QQoUKI3UKIZcpyDyHEDuXcvhNCRCjrmyjLmcp2nck+ZivrDwkhJqtzJq4RQrQSQiwVQhwUQqQLIRK1fp2FEI8rv9dpQohvhBCRWrvOQohPhRCnhBBpJus8dl2FEMOEEPuU97wnhBBOKyWlDMh/AEIBHAEQByACwB4AA9SuVwPPpSOAeOV1CwCHAQwA8G8AScr6JACvK6+vAvA7AAFgJIAdyvoYAFnK/62V163VPj8n5/4EgMUAlinLSwDcorz+GMD9yusHAHysvL4FwHfK6wHKtW8CoIfyOxGq9nk5ON/PAfxDeR0BoJWWrzOAzgCOAogyub53ae06A7gcQDyANJN1HruuAHYqZYXy3iud1kntH0ojfpiJAFaaLM8GMFvtenno3H4BcAWAQwA6Kus6AjikvJ4PYLpJ+UPK9ukA5pusNyvnb/8AdAGwBsB4AMuUX9zTAMIsrzGAlQASlddhSjlhed1Ny/nbPwDRSqATFus1e52V4H5cCVhhynWerMXrDEBnEdw9cl2VbQdN1puVs/cvkLtlDL80BrnKuoCmfA0dCmAHgFgpZb6y6SSAWOW1vXMPtJ/JOwCeAlCnLLcBcEZKWaMsm9bfeG7K9lKlfCCdcw8AhQAWKV1RnwghmkHD11lKmQfgTQDHAORDf91Soe3rbOCp69pZeW253qFADu6aI4RoDuAHAI9JKc+abpP6j2zN5K0KIaYBOCWlTFW7Lj4UBv1X94+klEMBlEP/dd1Ig9e5NYBrof9g6wSgGYApqlZKBWpc10AO7nkAuposd1HWBSQhRDj0gf1rKeWPyuoCIURHZXtHAKeU9fbOPZB+JqMAXCOEyAbwLfRdM+8CaCWECFPKmNbfeG7K9mgARQisc84FkCul3KEsL4U+2Gv5Ok8EcFRKWSilrAbwI/TXXsvX2cBT1zVPeW253qFADu67APRW7rpHQH/z5VeV69Qgyp3vhQDSpZRvmWz6FYDhjvmd0PfFG9bfodx1HwmgVPn6txLAJCFEa6XFNElZ53eklLOllF2klDror91aKeVtANYBuEEpZnnOhp/FDUp5qay/Rcmy6AGgN/Q3n/yOlPIkgONCiL7KqgkADkDD1xn67piRQoimyu+54Zw1e51NeOS6KtvOCiFGKj/DO0z2ZZ/aNyEaeQPjKugzS44AeFbt+jTiPC6D/ivbXgB/Kv+ugr6vcQ2ADACrAcQo5QWAecp57wOQYLKvvwPIVP7drfa5uXj+Y1GfLRMH/R9tJoDvATRR1kcqy5nK9jiT9z+r/CwOwYUsApXPdQiAFOVa/wx9VoSmrzOAlwAcBJAG4EvoM140dZ0BfAP9PYVq6L+h3ePJ6wogQfn5HQHwASxuytv6x+EHiIg0KJC7ZYiIyA4GdyIiDWJwJyLSIAZ3IiINYnAnItIgBnciIg1icCci0qD/B/l0XelGnbkVAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(trace.losses)\n",
"plt.yscale('log')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([<matplotlib.axes._subplots.AxesSubplot object at 0x147831650>,\n",
" <matplotlib.axes._subplots.AxesSubplot object at 0x147a76c50>],\n",
" dtype=object)"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA+kAAAFTCAYAAABF801FAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdd3hUVf4G8PfMZNJ7g/QKKXQJvVcpIkVBFFFEmm0Lll1RERXb7lrW/eGKooiKiCCoGEDpCSC9Q0IgvZHee2bu74+ELCA1meRMeT/PM0+YO7e8N0DO/eaee45QFAVEREREREREJJ9KdgAiIiIiIiIiasAinYiIiIiIiMhAsEgnIiIiIiIiMhAs0omIiIiIiIgMBIt0IiIiIiIiIgPBIp2IiIiIiIjIQLBIJ2oGIcQsIYQihAhsxraKEGLJdZZ7CSF0QoiBeogIIUSKEOLLZmwXKIRYIoQI1keO1tTccyQiIjJWxnANQkQtYyE7ABE1mQggD8B+Pe1vMoDSZmwXCOBVAHsBJOkpCxERERkufV+DEFELsEgnMhyTAGxSFEWnj50pinJcH/vRByGEGoBQFKVedhYiIiL6A71egxBRy7C7Oxm1xm7ZihAiXAjxqxCiQgiRJoR4rPHzmUKIeCFEuRBilxAi5IptNUKIpY1dpmsbvy4VQmiuOUawECJaCFEphMgTQvwbgNUN8swTQpwUQlQLIfKFEJ8LIVxv4zwcAQwD8GOLviFX7/OqruBXdI/rK4RYLYQoFUJkCSE+EkJYN64zFMCuxk22Na6vNC6/7XNs3OZNIcTfhRDJAGoB9BJC1Ash/nSdrC8IIeqEEB6N70cLITYLIbIbv+9nhBDPNhb7LfmeDG3Mdp8Q4kshRFHj92G1EMKtJfsmIiLzYqrXII2PvSlCiCeFEO8LIXIbj/9Lc7rYE9GdY5FOpmIdgGg0/Cb4KIAvhBBvAXgCwN8BPAYgDMC3V2yzqvGzrwDcA+BLAH9rXA4AEEJYAtgGoAeApwDMAhAE4OVrAwgh3gGwDMB2APcCeB7AGABbbqO4HIeGQnb7bZ9x830NIBHAFAD/RcN5vdj42bHG9wDwJwD9Gl/HgDs+x1kAxgN4rvFrUuN2D18n00wAWxVFyWt8HwxgB4DZjduuArAEwJvNOuM/+hCAAuBBAC81nst6Pe2biIjMi6leg7wIoENj/qcA9ATw27W/SCCiVqAoCl98Ge0LDYWbAuCRK5a5AKgHUADA8Yrlf2pcNwBA58Y/L7lmfy83Lu/a+H5u4/u+V6yjAnC2cXlg47JAAFoAi6/Z34DG9SZdsex6x/0OwHo9f29SAHx5xftZjcd+7Zr1fgGQcMX7oY3rjbxmvTs9xywANtesO6Pxs7ArlnVvXDbtBuch0PBozksAigCobnSOt/E9uXxuW2+Qa4Tsf9N88cUXX3wZx8tUr0Ea96cAOHdNm3t5f4/L/t7zxZepv3gnnUzFlst/UBSlCEAugAOKolw5cFp841c/AIMb//zNNfu5/H5I49d+ANIVRTlwxf51AL6/ZrtRaGg4VwshLC6/ABwEUHbF8f6g8TflY6HHru63EH3N+9MA/G9juzs9x62KolRds2wjgHI03Dm/bCaAEgA/X14gGkaZXS6ESEXDb/frACwF4AzA8zay3sq1f3/rAOjQ8PdNRER0J0z1GmS9csUz6oqi7AOQAbaVRK2OA8eRqSi65n3tDZYBgDWAy89oZV+zzqXGr5c/9wKQc53jXbvscuF48Qb5bva883AAtmi4o90WCq95X4MbPN92jTs9x2u/t1AUpVII8QOAGUKIV9BwUfEggHWKolQDgBBChYaC3RsNdyniAVShoRvhS2j4+2upq/7+FEWpFUIUAfDRw76JiMi8mOo1yI2OzbaSqJWxSCdzdblQbY+G57NxxfsrP88G0Ok627e75n1B49fR+GPDfOXn1zMJwB5FUYpvso4huNNzVG6wn68BPApgIAAbNFyEfH3F5yEAogDMVBSl6S6DEGJCMzLfyFV/f413ElwAZOrxGERERNdjLNcg1x7n8rITN9kfEekBu7uTuYpp/Dr9muUzGr/ubvz6OwA/IUTfyys03umdds1229DQXdpfUZQj13klXy+EEEKgYYCXturqfjtqGr/aXLO8Wed4HbvQ0F1uZuMrBUDsFZ/bNn6tu7ygcZCaGdCfa//+pqLh5+HvejwGERHR9RjLNcj9jce7vP4AAL5gW0nU6ngnncySoihnhBBrACxpfG5rPxqesXoFwBpFUU43rnp59NUNQohFaHjObAEAx2v2lyiEeBfA/wkhwgDsAVCNhmfPRgFYoSjKLvxRHzTcSTakIj0BDYPezBZCFKKhaD/fgnO8iqIoOiHEagDzAWgAfKAoypV33eMApAJ4UwihRUOx/lf9nR4AoJMQYiUaBsvpiIZR43crirJDz8chIiK6ihFdgzgA+FEIsRyAB4C3AVxAw4j0RNSKWKSTOZuFhmnBZqNhRNUsAO8CeO3yCo3PKo8C8H8APgZQgYYpVKIBfHLlzhRFWSSEiEPDNCVPoaG7dzoaphK7cIMMkwAcVRQlQ29n1UKKohQIIZ5Gw1QwewCo0TB/6u5mnuP1fN24/8t/vvL4tUKISWj4nn+Fhm5/XwBIA/BZc8/rGn9Gw92DtWg4v01oGHmXiIioLcyC4V+DvA0gFA3Tw9mhoSfc04qi1N1gfSLSE3H1DSwiaktCiHgA3yiKslR2FnMghBiKhouMUYqitMWc9ERERAbpRtcgQohAAMkA5iqKskJCNCKzxzvpRBIpihIuOwMRERGZH16DEBkuFulEBkYIoQYgbrKKDjcf9FFRFEWr31SGr/G5vpsxu+8JERERERkfdncnMjBCiN0AhtxklVVomMLsRvYoijJUn5mMgRDiVj/MHlMU5cu2yEJERERE1Fws0okMTOPIrA43WSUfgPtNPi9TFOW8flMZPiFE1C1WSVYU5WZzxRIRERERSccinYiIiIiIiMhA3OoZTlbwREREhuNm41XcDrbrREREhuO67frNBp8iIiIiIiIiojbEIp2IiIiIiIjIQLBIJyIiIiIiIjIQLNKJiIiIiIiIDASLdCIiIiIiIiIDwSKdiIiIiIiIyECwSCciIiIiIiIyECzSiYiIiIiIiAwEi3QiIiIiIiIiA8EinYiIiIiIiMhAWMgOQEStR1EUVNZqUVmrhUYt4GxrKTsSERGZCZ1OQWFlLeq1CnSKAmdbDWwteelJRHQr/ElJZGJ0OgXb4nLw29kc7EnIQ355TdNn3k7W6O7vjOm9/DGogzuEEBKTEhGRqdDpFJzMKMaBpEIcTytC/KUyZJdUoU6rXLWeg5UFQjzt0d3PGb2DXDEszBM2lmpJqYmIDJNQFOVmn9/0QyIyHIqi4NezOfhgWwLO55TB2VaDQR08EOnlCHsrNSprtTibVYr9iQXIL69BhJcjXrknAv1D3GVHJ6Lb19LfrLFdJ706k1mCtYfTsfXsJeSVNfxSOMjdDp28HeHnaov2jtbQqFUQAiiqrEVOSTXiL5XhdGYJKmu1sLVUY0yn9pg7OBgRXo6Sz4aIqM1dt11nkU5kAoora/HihtPYcuYSgt3t8OeRHXBPV2+oVX/8f19Tr8VPJ7KwbNdFpBZUYs7AIDx3dxisNbyTQWQEWKSTdIqi4LdzOfh4dyJOphfDykKFERGeuLtTewzq4AFXu1s/WlWv1eFQSiE2nczGppNZKK+px6jIdlg0LgJB7nZtcBZERAaBRTqRKTqeVoQnvjmGgooaLBwVhrmDgmChvvWYkJW19Xh7czy+PpCK7n7OWDmrF1xu48KKiKRikU5SxV7Iw9ub43EuuxSBbrZ4tH8gpvTwhZOtptn7LKmsw5f7U7BibxJq63X466iOmDPw9toyIiIjxyKdyNTsPp+LJ745Bg8HKyx76C508XW6431sPZONP313Av6utvhqdm94O9u0QlIi0hMW6SRFRlEl3vjlHH49mwN/V1v8ZWQH3NvNW6+FdE5pNV758Qx+O5eDPkGuWDbjLrjbW+lt/0REBohFOpEpiT6VjT9/dxwd2zlg1eze8HBo/oXMgaQCzF11BE62Gmx4sj88Haz1mJSI9IhFOrUpRVGw9nA63vjlHHQK8PTwUMwZFAQri9Z5REpRFGw8nokXN5yGq50lls/sia6+zq1yLCIiA8AinchUxCTk4fFVh9Hdzxmfz+oFR+vmdzO87ER6MR789ABCPO2wdl4/2Flx8gciA8QindpMcWUtnlt3EtvjctEv2A3/nNoVvi62bXLsM5klmP/1URRX1uLzWb3QN9itTY5LRNTGWKQTmYIT6cV46LMDCHCzw9r5ffVSoF+2Mz4Hc786ioGh7vhiVq/rDjxHRFKxSKc2cTqjBE+sPoqc0mq8ODYCs/oHQtXGbUJOaTVmrDiIjKJKLJ8ZhSEdPdr0+EREbeC6P1g5IgeREckuqcKcVYfhZm+JVY/p5w76lYaHt8MbEztjT0IePtyeoNd9ExGRcfjlVBbu+2Q/dDoF6xb0x+yBQW1eoANAO0drfDevL4Lc7THvqyM4mlrY5hmIiGRgkU5kJGrqtXjim2OoqtVi5axe8HRsnefGH+rjj2lRvvjPzovYfT63VY5BRESGR1EULNt1EU9/exxdfZyw6ZmB6O4n93lwd3srrJ7TB15O1piz6giS8sql5iEiagss0omMxJKfz+JEejHem9YNoZ4OrXqs1+7tjPD2Dvjr2hPILqlq1WMREZF8Op2CxT+dxT9/PY+J3b2xem4fuBnIyOqudpZYNbs3VEJg1srDKKyolR2JiKhVsUgnMgI/ncjEmkPpeGJoCMZ09mr149lYqvHxjLtQXafDC+tP4RZjVxARkRGr0+rwl7Un8PWBVMwfHIwPH+jeaqO3N1eAmx1WPBqFSyXV+PN3x6HVsV0iItPFIp3IwGUUVeLlH8/gLn9nPDuqY5sdN9jDHovGRyD2Qj5WH0xrs+MSEVHbqa3X4cnVx/DzySz8bUw4XhwXASEMc9DQHv4ueH1iJ8ReyMcH2zhuChGZLhbpRAZMq1Pw17UnoCjAhw/0gIW6bf/LPtzHH4M6uOOtzXFILaho02MTEVHrqq3X4elvj2HbuRy8dm8nPDE0RHakW5re2x8PRPnh/3ZdxK54jptCRKaJRTqRAVsRm4TDKUV4fWIn+Lu1zdy0VxJC4N37ukItBBZtPM1u70REJkKrU/CXtcfx27kcLJkQiUf7B8qOdNtem9gJ4e0d8Pz6Uygor5Edh4hI71ikExmopLxyvL8tAaMi22FyDx9pObydbfDC2HDsu1iAH09kSstBRET6oSgKXv7xDDafvoSXx0dg1oAg2ZHuiLVGjQ+nd0dpVR3+voG/QCYi08MincgA6XQK/v7DaVhaqLB0UmfpzwfO6O2P7n7OWPpLHIorOaouEZExe39bAtYcSsOTQ0MwZ1Cw7DjNEt7eES+MCcO2czn4/ki67DhERHrFIp3IAH17KA2HUgrxyvhItGul+dDvhEol8NbkLiiuqsO7W+NlxyEiomb6/kg6/rPzIqb38sPzd4fJjtMiswcEoW+wK5ZGxyGntFp2HCIivWGRTmRg8str8I+t8egf4oapUb6y4zSJ9HbEY/0D8d3hdJzJLJEdh4iI7tDviQVYtOE0BnVwxxsG0EurpVQqgbendEVtvQ6v/nRWdhwiIr1hkU5kYN7eHI+qOi1en2h4F1DPjOgAF1tLvL7pHJ8BJCIyIsn5FXhi9VEEuNni/x66C5o2ni2ktQS52+EvIzti69lL2HomW3YcIiK9MI2f0EQm4lByIX44loE5g4IR6mkvO84fONlo8NzoMBxKKcTm05dkxyEiottQUlmHx788DAHgi1m94GSjkR1Jr+YMCkKklyNe/fksKmrqZcchImoxFulEBkKrU7D4pzPwdrLGM8NDZce5oQd6+SHCyxFvbY5DdZ1WdhwiIrqJOq0OT6w+ioyiKnz6SBQC3OxkR9I7jVqFNyZ1Rk5pDT7Zkyg7DhFRi7FIJzIQ3x9JR/ylMiwaHwFbSwvZcW5IrRJYfE8kMoursCI2SXYcIiK6iTej47A/sQBvT+mCXoGusuO0mp4BLpjU3RvLY5KQXlgpOw4RUYuwSCcyAGXVdXjvt/OICnDB+C5esuPcUr8QN4zt3B7LdiXiUglH1CUiMkSbTmbhy/0pmD0gCPf1NJyBSFvL38aGQy0E3tocJzsKEVGLsEgnMgDLdiUiv7wWiydEGtxgcTeyaFwEtIqCf3BKNiIig3Mxtxx//+EUega44MVx4bLjtAkvJxs8OTQEW85cwv7EfNlxiIiajUU6kWRZxVX4Yl8yptzlg66+zrLj3DY/V1vMGRiEDcczcSqjWHYcIiJqVFFTjye+OQprjRrLTGgk99sxd3AwfJxt8Pqmc6jX6mTHISJqFvP5qU1koP69/QKgAM+ODpMd5Y49MTQErnaWeJd304mIDIKiKFi08TQS88rx0YM90N7JWnakNmWtUeOl8RGIv1SG7w6ny45DRNQsLNKJJErMK8e6o+mY0dcfPs42suPcMQdrDZ4eFop9FwsQeyFPdhwiIrP3zcE0/HQiCwtHdcSAUHfZcaQY27k9+gS54r3fzqOksk52HCKiO8YinUii939LgLVGjaeGGe6Ua7cyo68/fF1s8O7WeOh0iuw4RERmK/5SKd745RyGhnngyaHG2660lBACiydEoriqDv/llGxEZIRYpBNJcjqjBNGnszFnYBDc7a1kx2k2Kws1Fo7qiDOZpYg+nS07DhGRWaqu0+LPa07A0VqDf03tBpXKOAYhbS2dvJ0wsZs3Vu5L5iwkRGR0WKQTSfKPX+PhbKvBnMHBsqO02MTuPghv74D3fjuPOg7UQ0TU5t7dGo/zOWX419SuRv2LX316dnQYdIqCf++4IDsKEdEdYZFOJMHviQWIvZCPJ4eGwNFaIzvODe3btw99+vSBtbU1goKC8NFHH113PbVK4IUxYUgpqMTaw+nYtm0bBgwYACcnJ7Rr1w6TJ0/G+fPn/7BdamoqHnzwQbi6usLW1hbdunXD1q1bW/u0iIhMyu7zuVi5LwWz+gdiaJin7DgGw8/VFg/c5YX/vvsq3Dw8YGdnh/HjxyMlJeWW2y5YsADh4eGwt7eHi4sLBg8ejO3bt1+1Tl5eHv70pz+hd+/esLS0RGBgYOucCBGZHRbpRG1MURT849d4tHe0xiP9AmXHuaGLFy/i7rvvRlBQEDZv3oz58+dj4cKFWLFixXXXHxbmid6Brnjrq2iMHz8ePj4+WLduHT7++GMkJSVh5MiRKC0tbVo/PT0d/fr1Q3FxMVauXImff/4ZM2fORFVVVVudIhGR0Ssor8Fz604hrJ0D/j7WPOZDvxNp0ctQfno7ut/3NNavX4/8/HyMGjUK1dU37wJfVVWFp59+Ghs3bsQ333wDd3d3jB07FgcOHGhaJzMzE2vXrkX79u3RvXv31j4VIjIjQlFuOtATR4Ei0rPt53Iw56sjeGtyFzzUx192nBuaP38+du3ahXPnzsHCwgIA8OSTT2LTpk1IS0uDEH983vFoaiGGPzAfSsJuFOZmN2136tQpdOvWDZs3b8bYsWMBANOnT0dmZib27NkDlYq/LyS6TS190JjtuglRFAVzvzqCmIR8/PT0AER4OcqOZFAyMjIQGBiIaQvfwn5VJ/z01AC4qyoQFBSEjz/+GHPmzLntfWm1WgQFBWHSpElNvcp0Ol1T+/Xcc89h/fr1t3WXnojoCtdt13llTHSbZs2ahaioKERHRyMyMhK2trYYP348CgsLcfHiRQwbNgx2dnaIiorCqVOnmrbT6XR45513EBoaCisrK0waGgXL5BhMjfJtWic6OhqjRo2Cp6cnHB0d0bdvX/z2229XHX/JkiVwd3fH8ePH0bdvX9ja2qJHjx6IjY1tlfPdsmULpkyZ0lRoAw2FdUZGBs6cOXPdbXoGuMLf2RLVigXKav73bLqzszOAhgtKACgpKcGGDRvw5JNPskAnImqm74+kY3tcLl4YE9aiAl1f7VvHjh2xatWqq/Yts327fJz3XpgLVztL/OPXePj4+GDgwIHYsmXLHe1LrVbD2dkZtbW1TcvYfhFRa+FPF6I7kJaWhsWLF2Pp0qX49NNPsX//fsybNw/Tp0/H9OnTsX79etTX12P69OlNBekzzzyDpUuXYt68eXjrk2+gDu6Di+v/hV+3bG7ab3JyMiZMmICvv/4aP/zwA/r374+xY8di3759Vx2/srISjz76KObPn48ffvgBVlZWmDJlCiorK5vW0el0qK+vv+lLq9Xe9DwrKiqQnp6O8PCru05GREQAAOLj42+47RvPPYH6sgJMf/rvKCoqQnp6OhYuXIjw8HCMGDECAHDs2DHU1dVBCIEBAwZAo9HA19cXb7/9Nm7Ru4eIiABkFVdh6S9x6BvsitkDglq8v5a2b9HR0Zg8eTJmz56NX375pWm/Mtu3+Ph4+Pr6wsvdBU8NC8W+iwWIvZCHiIiIm7ZjlymKgvr6ehQUFOCDDz7AhQsXMHv27JZ+q4mIbk1RlJu9iKjRo48+qqjVauXixYtNy55//nkFgLJq1aqmZdHR0QoA5dy5c8qFCxcUIYTy5ZdfKjqdTpm8bK/S763tyowZDytRUVHXPY5Wq1Xq6uqU0aNHK4899ljT8ldffVUBoOzYsaNp2fHjxxUAypYtW67KiYYurTd8DRky5KbnmpGRoQBQNm7ceNXyuro6BYCyfPnym25//8ufKCpr+6bjhYeHK6mpqU2ff/vttwoAxdHRUfnb3/6m7Ny5U3nllVcUlUqlLFu27Kb7JjJzt2q32a6bAZ1Op8z8/KAS/vIWJTW/osX7a2n7dqWZM2caTPs2Z84cpVu3boqiKEp1Xb3S/+0dyj0fxSqLFi1SvLy8bvl9WbNmTdN+7ezslJ9++umG6z777LNKQEDALfdJRHSN67bX/+vHSkS3FBgYiJCQkKb3oaGhAIDhw4f/YVlmZiYSExOhUqkwefJkxJzPwdGUArw2IRLCeyTWrv0OWq0WarUaGRkZeOmll7B9+3ZkZ2c33aUYMGDAVce3tLTE0KFDm95HRkYCaHju7rIlS5bg6aefvul5ODg4NP1Zq9Vedff6yu7tzXH27Fns+uQV2If1x9333oeHo9rjnXfewbhx47B//344Ojo2HW/s2LF45513AADDhg1DRkYG3n77bTz55JMtykBEZMrWHclATEIeXru3E/zdbPWyz5a0b/X19U3rjBgxAmvWrDGI9u1KVhZq/HVURzy37iRsC8tvuo/L7r77bhw+fBj5+flYvXo1pk+fjs2bN1+Vk4ioNbBIJ7oDl5+tvszS0vIPyy8vq66uRn5+PrRaLZycnJo+f+yf/9s+Ozsb3t7euPfee1FWVobXX38doaGhsLOzw+LFi5Gbm3vV8RwcHK56Bu7KY13m7+8PX19f3MyVg76FhIQgNTW16X1ycjI8PDwANDw7fqWioiIAgIuLyw33/corryCsYweMefEDrNyXjLf7DUF09CD4+vpixYoVWLhwYdP2w4YNu2rb4cOHY+XKlSgtLYWjIwdAIiK6VlZxFd745Rz6BLliZt8Ave1XH+3blQyhfXNxcbmqHZvcwwefxSRhb0wqnG/Sjl25fVRUFABgzJgxyMrKwuLFixETE3PLbYmIWoJFOlErcnV1hYWFBZZ/vxkv/XgWcwYFYWJ3n6bPPT09cfHiRRw/fhxbtmzBmDFjmj5r7lRks2fP/sPAPdcaMmQIdu/eDQDYtGkTampqmj7z9vaGpaUl/Pz8/vDM3uX31z6rfu06w4YNwxNDQ7DmUBo+3H4B/3mwBwICApCYmAjgf8+2X3kH/8r3HIyHiOiPFEXBixtOo16n4B/3d4VK1dLB/pvvcvu2b9++6/7MNoT2LTw8HOnp6aioqICdnR3UKoHn7w7DhH8nwyfozmdX6dGjB7777rvmRCciuiMs0ola0fDhw6HVarE6Ng4+HTph0SPDYWOpvmqdyxcrVlZWTctSU1Oxb98+dO3a9Y6PeafdAbt06XLddcaOHYuNGzdi6dKlUKsbMq9duxZ+fn7o3LnzDfcdEBCA48ePw93eCrP6B+K/exLxUDcXpKSkIDAwEEBDt8pOnTph586dWLBgQdO2O3bsQEhICOzt7W/3dImIzMa6oxnYk5CHJRMiEeBmJzXL5fatpKQEo0aNuu46stu30aNHAwA2btyIhx9+GAAQ4VSPmoxzyOsxAlW12j+0yTeiKAp+//13BAW1fJA+IqJbYZFO1IrCwsJw34zHsOG/L2PiIwuwP1aF6upqnD17FgkJCVixYgXCw8Ph6+uLZ599Fm+88QbKysrw6quvwsfH59YHuI7AwMCmYrglnn/+eaxevRozZ87E3LlzcfjwYSxfvhz//e9/r+pOaGFhgcWLF2Px4sUAgAULFmDSpEmYNWsWJky+H9rEg5g44W+wtLTEjBkzmrZ74403cN999+H555/H6NGjsXv3bnz99df46quvWpydiMjU5JfX4M3oOPQOdMUj/QJlx0FYWBgWLFiA6dOn44UXXkBUVJTBtW++vr54/PHH8Ze//AWKosDDwwNLliyBr58/dCGD8MW+ZDw1LBSvv/46Xn/99aZn62NjY/H+++9j8uTJ8Pf3R0FBAVatWoUDBw5g06ZNVx1j/fr1AICEhARUVlY2vR8yZEjTo2NERHeKRTpRK7MfPh/t8zWIi43GuC//DUdHR0RGRuLxxx8H0HCHYcOGDXjqqadw//33w9fXFy+99BJ27959w/nI20JoaCi2bt2KhQsXYuzYsWjfvj3ee+89zJkz56r1tFotdLr/zYk+ceJErF27Fv/85z+xYcMGKGoN6l2D8NXaTfD29m5ab/Lkyfjqq6/w5ptv4t///jf8/f2xbNmyqwp5IiJq8PbmeFTW1uOtKZ2ldnO/0rJly9CxY0d89tlnWLx4sUG2bx999BHs7OywcOFCVFZWYsiQIVizZg3ejMnHJ3sS8VBvf+h0uqumbvPz84NarcaiRYuQl5cHDw8PdO/eHXv37kW/fv2u2v/UqVOv+37Xrl0cYI6Imk1c+0zoNThhMVELJOSUYfQHMRAPM60AACAASURBVPjryI7488gOsuNIUVpdh8H/2IUefs5Y+Vhv2XGIjF1LqzO260boQFIBpn96AE8ODcELY248JgjdvvOXyjDm3zGYMzAIL42PlB2HiMzXddt1js5E1IpWxCbBWqPCzH76G4HX2DhaazBvcDB2nc/D0dQi2XGIiIxKbb0OL/94Br4uNnhmuHn+src1hLV3wH13+WLV76nILG7eQHZERK2FRTpRK8ktrcaPx7MwtacfXO0sZceRalb/QLjbW+L9bedlRyEiMiqfxSbhYm45Xp/Y6bYHOaPb89dRHQEAH2xLkJyEiOhqLNKJWsmq31NQp9Ph8YEcCdbW0gILhoRg38UCHEgqkB2HiMgopBdW4j87L+DuTu0wPLyd7Dgmx8fZBo/0DcCGYxlIyCmTHYeIqAmLdKJWUFFTj28OpOHuyPYIdJc7TY6heLhvADwcrPDRjguyoxARGTxFUfDqz2ehEgKvTugkO47JempYKOwsLfCPrezpRUSGg0U6UStYdyQdJVV1mDs4WHYUg2GtUWP+4GDsTyzA4ZRC2XGIiAzarvO52Bmfi7+O7AhvZxvZcUyWi50lFgwNwfa4HBxKZttERIaBRTqRntVrdfh8XzJ6BrigZ4CL7DgGZUafALjbW/JuOhHRTdRpdXgzOg7B7naYNSBQdhyTN3tAENo7WmNp9DnodJwAgYjkY5FOpGe/ns1BemEV5g7iXfRr2ViqMW9wMGIv5HOkdyKiG/juUBoS8yrw4rgIaNS8VGttNpZqvDAmDKcySrDxeKbsOERELNKJ9ElRFHwak4hAN1uMiuQgP9czo08AXO14N52I6HpKqurwwfYL6BfshpERnrLjmI1J3X3QzdcJ//g1HpW19bLjEJGZY5FOpEeHU4pwMqMEjw8KhlolZMcxSHZWFpgzKAh7EvJwIr1YdhwiIoPy8a6LKKqsxUvjIyAE25G2olIJLJ4QiZzSGnyyJ0l2HCIycyzSifTo05gkuNhqcP9dvrKjGLRH+gXC2VbDu+lERFdIK6jEyn0puO8uX3T2cZIdx+z0DHDFhG7eWL4nEZnFVbLjEJEZY5FOpCeJeeXYHpeDmX0DYGOplh3HoNlbWWDOwCDsjM/F6YwS2XGIiAzCu1vjoVYJPDc6THYUs/W3MQ3f+39sjZechIjMGYt0Ij35fG8yLC1UmNkvUHYUo/BI/0A4Wlvgo528m05EdCSlENGnszF/SDDaO1nLjmO2fF1sMXdQMH46kYVjaRzglIjkYJFOpAf55TX44WgG7rvLBx4OVrLjGAVHaw1m9Q/EtnM5uJhbJjsOEZE0iqLg3a3x8HSwwrzBnBlEtieGhsDDwQqvbeKUbEQkB4t0Ij34+vdU1NTr8PhAXlzdiUf7B8Jao8JyDtJDRGZs78V8HE4pwjPDQ2FraSE7jtmzs7LAonHhOJlejDWH02THISIzxCKdqIWqarX4+kAqRkZ4ItTTXnYco+Jmb4VpUX748UQmsks4SA8RmR9FUfD+tgR4O1ljWi8/2XGo0aTuPugX7IZ3t8Qjv7xGdhwiMjMs0ola6IdjGSisqMXcQbyL3hxzBwVDpwBf7E2WHYWIqM3tPp+H42nFeHp4B1hZcNBRQyGEwBuTOqOqTou3NsfJjkNEZoZFOlELaHUKPt+bjG6+Tugd5Co7jlHyc7XF+C5e+PZgGkoq62THISJqM5fvovu52mBqFKfuNDShnvaYNzgYG45l4kBSgew4RGRGWKQTtcD2uBwk51dg7uBgCCFkxzFa84cEo6JWi28OpsqOQkTUZrbH5eJ0ZgmeGd4BGjUvyQzR08M6wNfFBi//eAa19TrZcYjITLBFIGqBz2KS4OtigzGd2suOYtQ6eTthcEcPrNyXjOo6rew4REStTqdruIse6GaLKT18ZMehG7CxVOO1ezvhYm45PudjWUTURlikEzXT0dQiHEktwuMDg2DBOyAttmBIMPLLa7H+aIbsKEREre7Xs5cQl12KP4/swDbEwI2IaIfRke3w7x0JSC+slB2HiMwAWwWiZloRmwRHawtMi+JovPrQL9gN3fyc8VlsErScl5aITJhOp+DD7RcQ7GGHe7vxLroxWHJvJ1ioVPjbD6egKGyjiKh1sUgnaobUggpsPXsJD/cNgJ0V57TVByEEFgwORmpBJbady5Edh4io1eyIz8X5nDL8aXgHqFUcz8QYeDvbYNG4COxPLMC3hzh3OhG1LhbpRM3wxd5kWKgEZvUPlB3FpIzu1B6+Ljb4Yh+f+yMi0/VpTCJ8nG0wvquX7Ch0Bx7s7YeBoe54KzoOGUXs9k5ErYdFOtEdKqqoxfdHMjCpuw88Ha1lxzEpapXAo/0CcSi5EGcyS2THISLSu2NpRTicUoTZA4M4oruREULg7SldAAB//+E0u70TUath60B0h1YfTEVVnRZzBwfLjmKSpvXyg62lGiv3pciOQkSkd5/uSYKTjQbTe3E8E2Pk52qLF8dFYO/FfHx3OF12HCIyUSzSie5AdZ0WX+5PxdAwD3Rs5yA7jklystHg/p6+2HQyC3llNbLjEBHpTXJ+BX49dwkP9/XneCZG7KHe/ugf4oY3o+OQWVwlOw4RmSAW6UR34KcTmcgvr8G8QbyL3poe7R+IWq0O3x7k4DxEZDpWxCZBo1LhUY5nYtRUKoF37+sKnaLghfUnoeOMJESkZyzSiW6TTqfgs9hkRHo5ol+Im+w4Ji3Ewx7Dwjzw9YFU1NRrZcchImqx/PIarD+agSl3+cDTgeOZGDs/V1u8PD4S+y4WYOX+FNlxiMjEsEgnuk27E3JxMbcc8wYHQwhOmdPaHhsQhPzyGkSfypYdhYioxb76PRU19TrMYU8sk/Fgbz+MjGiHd7fGI/5Sqew4RGRCWKQT3aZPY5Lg5WTNKXPayKAO7gj1tMcX+5I5gi4RGbWqWi2+/j0FIyPaIdTTXnYc0hMhBN69rwscrTX485oTqK5jzy8i0g8W6US34VRGMQ4kFWL2AE6Z01aEEHhsQCDOZJbiSGqR7DhERM22/mg6iirrMH8I76KbGjd7K/zz/q44n1OGf/56XnYcIjIRrDaIbsNnsclwsLLA9N6cMqctTenhCycbDb7kdGxEZKQURcGq31PR1dcJUQEusuNQKxgW7olH+gXg873JiL2QJzsOEZkAFulEt5BeWInNp7PxYB9/OFhrZMcxKzaWakzt6Ytfz15Cbmm17DhERHfs96QCXMwtx8y+ARzPxIS9ODYCoZ72eG7dSRRV1MqOQ0RGjkU60S2s3JcCAWAWp8yR4qE+/qjXKVh7OF12FCKiO/bV/lS42GowoZu37CjUimws1fjwge4orKjFoo2nOZYKEbUIi3SimyipqsPaw2mY0M0b3s42suOYpWAPewwIdcOaQ2nQci5aIjIi2SVV2BaXg2m9/GCtUcuOQ62ss48Tnh0dhi1nLmH90QzZcYjIiLFIJ7qJNYfSUFGrxZxBQbKjmLWH+wQgq6QaO+NzZUchIrpt3x5Mg05R8HCfANlRqI3MHRSMPkGuWPLzWaQWVMiOQ0RGikU60Q3U1uuwcl8yBoa6o5O3k+w4Zm1kZDt4OljhmwOpsqMQEd2Weq0Oaw+nY2hHD/i52sqOQ21ErRJ4/4HuUKkEFn5/EvVanexIRGSEWKQT3cCmk1nIKa3B3MGcMkc2jVqF6b39EXMhD2kFlbLjEBHd0p6EPOSW1eCBXv6yo1Ab83G2wRsTO+NoahE+2ZMoOw4RGSEW6UTXoSgKPotNQlg7Bwzu4C47DgF4sLcfVEJg9SHeTSciw7f2cDrc7S0xIsJTdhSSYGJ3b9zT1Qsfbr+AUxnFsuMQkZFhkU50HbEX8hF/qQxzBwdzyhwD4eVkgxHhnlh3JAM19VrZcYiIbiivrAY743Nx312+0Kh5qWWOhBB4c1IXuNtb4S9rT6Cqlu0WEd0+thxE1/FZbBLaOVrhXk6ZY1Ae7huAwopabDl9SXYUIqIb2nAsA/U6BVOj/GRHIYmcbDV4b1o3JOVV4J0tcbLjEJERYZFOdI1zWaWIvZCPWf2DYGnB/yKGZGCoOwLcbDmAHBEZLEVRsPZIOqICXBDqaS87Dkk2INQdswcEYdXvqdiTkCc7DhEZCVYgRNdYEZsEO0s1HurDwX4MjUolMKOPP46kFiEhp0x2HCKiPziSWoSkvApM68W76NTghTFh6OBpj+fXnURRRa3sOERkBFikE10hu6QKP5/MwrRefnCy0ciOQ9cx5S5faNQCaw+ny45CRPQHaw+nw85SjfFdvGRHIQNhrVHjw+ndUVRZi0UbT0NRFNmRiMjAsUgnusKX+1KgUxTMHhAkOwrdgLu9FUZFtsOGYxxAjogMS1l1HaJPZWNCN2/YWVnIjkMGpJO3ExaOCsOWM5fww7FM2XGIyMCxSCdqVFZdh28PpmFcFy/4udrKjkM38UAvfxRV1mHbuRzZUYiImvxyKhtVdVo8wK7udB3zBgejd6Arlvx8FlnFVbLjEJEBY5FO1Gjt4XSU1dRj3uBg2VHoFgaGusPH2YZd3onIoKw9nI6O7ezR3c9ZdhQyQGqVwL+mdoNWp7DbOxHdFIt0IgB1Wh1W7ktBnyBXdPXlxZWhU6sEpkX5IfZCPtILK2XHISLCxdwynEgvxrQoPwghZMchA+XvZosXxoRh9/k8bGC3dyK6ARbpRAB+OZWFzOIq3kU3IlOjfCEE8P0R3k0nIvl+PpEFlQDu7e4tOwoZuEf7BSIqwAWvbTqL3NJq2XGIyACxSCezpygKlu9JQsd29hgW5ik7Dt0mb2cbDOnogXVHMlCv1cmOQ0RmTFEU/HQyC/1D3OHpYC07Dhk4lUrg3fu7orpeh1d+OsNu70T0ByzSyeztTshD/KUyzB8cApWKXRSNyfRefrhUWo2YC3myoxCRGTuVUYLUgkrc24130en2hHjYY+Gojvj1bA42n74kOw4RGRgW6WT2PtmdCG8na3ZRNEIjItrB3d4S3x1il3cikufnk1mwVKtwd+f2sqOQEZkzMAhdfJyw+KczKKyolR2HiAwIi3Qya8fTinAwuRCzBwZBo+Z/B2OjUatwX09f7IjP5XN9RCSFVqdg08ksDA3zgJONRnYcMiIWahX+ObUrSqrqsDT6nOw4RGRAWJWQWVu+JwlONho82NtfdhRqpgei/KDVKVh/LEN2FCIyQweTC5BbVsPeWNQs4e0dMX9IMDYcy8T+xHzZcYjIQLBIJ7OVmFeOX89dwsy+AbCzspAdh5op2MMevQNdsf5IBgffIaI29/OJLNhZqjEivJ3sKGSknhneAf6utnh54xnU1GtlxyEiA8AinczWitgkWKpVmDUgUHYUaqGpUb5Iyq/AkdQi2VGIyIzU1Gux5cwljO7UHjaWatlxyEhZa9R4Y1JnJOVX4JPdSbLjEJEBYJFOZim3tBo/HM3E1ChfuNtbyY5DLTSuixfsLNX4/jAHkCOithOTkI+SqjqO6k4tNqSjB+7p6oVluy4iKa9cdhwikoxFOpmlL/aloF6nw9xBwbKjkB7YWVngnq7eiD6djfKaetlxiMhM/HwyCy62Ggzs4C47CpmAxfdEwspChZd/5NzpROaORTqZnbLqOqw+kIqxXbwQ4GYnOw7pybRevqis1WLzqWzZUYjIDFTU1GP7uRyM6+LF2UFILzwdrfHCmDDsTyzAjycyZcchIonYqpDZ+fZgGspq6rFgcIjsKKRHd/m7INjDDt8fYZd3Imp92+NyUFWnxcTuPrKjkAl5qE8Auvs5Y+kvcSipqpMdh4gkYZFOZqWmXovP9yZjQKgbuvg6yY5DeiSEwLQoPxxJLUIin+cjolb2y6lstHe0RlSAi+woZELUKoGlkzqjsLIWH25PkB2HiCRhkU5m5afjWcgtq8GCIbyLboqm3OUDtUpg3RHOmU5Eraeiph4xCXkY07k9VCohOw6ZmM4+Tniwtz+++j0V5y+VyY5DRBKwSCezodMp+CQmEZ28HTEwlIP8mCJPB2sMC/PED8cyUK/VyY5DRCZq9/k81NTrMKZze9lRyEQ9PzoMDtYWePVnDiJHZI5YpJPZ2BaXg6S8CswfEgIheOfDVE2L8kVeWQ32JOTJjkJEJmrLmWy42VmiV6Cr7ChkolzsLPHs6DAcSCrELxwQlcjssEgns6AoCj7Zkwg/VxuM450PkzYs3BPu9pYcQI6IWkV1nRa74nMxulM7qNnVnVrRQ739EenliLc2x6GC04sSmRUW6WQWDqcU4XhaMeYNCoYFp8oxaRq1ClPu8sWOuFzkl9fIjkNEJmbvhXxU1GoxprOX7Chk4tQqgdcndkJ2STWW7booOw4RtSFWK2QWPtmTCFc7S9zf0092FGoDU3v6ol6n4MfjnGeWiPRr69lLcLC2QL9gN9lRyAxEBbpiSg8frIhNRnJ+hew4RNRGWKSTyTt/qQw743Mxq38gbCzVsuNQG+jQzgE9/J2x9nA6B9whIr2p0+qw7VwORkW0g6UFL6Gobfx9bDgsLVR4bdNZtmlEZoItDJm85TGJsNGoMbNvgOwo1IamRfnhQm45TmaUyI5CRCbiQFIBSqrqcDfHNqE25Olojb+M7IDd5/OwIy5XdhwiagMs0smkZRVX4ecTWZje2w8udpay41AbuqerF6w1Kg4gR0R6s/XMJdho1BjS0UN2FDIzj/YPRKinPV7/5Ryq67Sy4xBRK2ORTibt873JUADMGRQsOwq1MQdrDcZ18cKmE1moquUFDRG1jFan4NezORgW7gFrDR+doralUauwZEInpBVW4rOYJNlxiKiVsUgnk1VcWYs1h9Jwbzdv+DjbyI5DEkyL8kNZTT22nuUcs0TUMsfSipBfXsNR3UmagR3cMbZze3y8OxFZxVWy4xBRK2KRTibrmwOpqKzVYv4Q3kU3V32CXBHgZovvD2fIjkJERm7L6UuwVKswPNxTdhQyY4vGRUCnKHh7S7zsKETUilikk0mqrtNi5b4UDAvzQHh7R9lxSBIhBKb29MXvSQVIK6iUHYeIjJSiKPj17CUM6uAOeysL2XHIjPm52mLBkBBsOpmFg0kFsuMQUSthkU4mad3RDBRU1GL+kBDZUUiy+3r6Qghg/VEOIEdEzROXXYbM4iqM7tROdhQiLBgSAh9nG7z681nUa3Wy4xBRK2CRTiZHq1PwWUwSuvs5o0+Qq+w4JJmXkw0Gd/DA+qMZ0Oo4vywR3bkdcTkAgGHs6k4GwMZSjUXjIhB/qQxrDvMX0ESmiEU6mZwtZ7KRVliJBUNCIISQHYcMwLQoP2SVVGPfxXzZUYjICO2Iz0U3P2d4OljLjkIEABjXpT36Brvivd/Oo7iyVnYcItIzFulkUhRFwSd7EhHsbodRkeyWSA1GRnrC2VbDOdOJ6I7lldXgZEYxRvAuOhkQIQSW3NsJpVV1eH9bguw4RKRnLNLJpOxPLMCZzFLMGxwMtYp30amBlYUak7r74LezObzjQER3ZNf5XCgKMCKCRToZlvD2jpjZNwDfHEhFXHap7DhEpEcs0smkfLInER4OVph8l4/sKGRgpkX5oVarw08nsmRHISIjsiMuB15O1oj04kwhZHj+OqojnGw0WPLzWSgKx10hMhUs0slknMksQeyFfMweEAQrC7XsOGRgIr0d0dnHkV3eiei21dRrEXshH8PDPTnGCRkkZ1tLPHd3GA4mFyL6dLbsOESkJyzSyWQsj0mCg5UFZvT1lx2FDNS0KD+czSrFmcwS2VGIyAgcTCpEZa2WXd3JoE3v5Y9IL0e8FR2Hytp62XGISA9YpJNJSCuoRPSpLDzU1x+O1hrZcchA3dvNG5YWKqzj3XQiug074nJgrVGhf4i77ChEN6RWNQwil1VSjU92J8qOQ0R6wCKdTMKKvUmwUKkwe0CQ7ChkwJxtLXF3p/b48UQWquu0suMQkQFTFAU74nMxMNQd1ho+QkWGrXeQK+7t5o1PYpKQXlgpOw4RtRCLdDJ6BeU1+P5IOib38EE7R85hSzf3QJQfSqrqsO1cjuwoRGTAEnLKkVFUheHhnM6TjMOL48KhFgJvRsfJjkJELcQinYzeqv0pqKnXYe7gYNlRyAj0D3GDj7MNB5AjopvaEd/wi7zhnB+djISXkw2eHh6KrWcvYe+FfNlxiKgFWKSTUausrcdXB1IxMqIdQj3tZcchI6BSCdzf0xd7L+Yjs7hKdhwiMlA74nLR2ccR7Z3YQ4uMx+MDg+DvaovXNp1FnVYnOw4RNROLdDJq3x9OR3FlHRYMCZEdhYzI/T19oSjAD0czZEchIgNUWFGLY2lFGMGu7mRkrDVqvHJPJC7kluPr31NlxyGiZmKRTkarXqvDZ7HJ6BXogp4BLrLjkBHxc7XFgFA3rDuaDp1OkR2HiAzM7vO5UBRw6jUySiMjPDG4owc+2J6A/PIa2XGIqBlYpJPRij6djcziKswbzLvodOemRfkhvbAKB5ILZEchIgOzIy4Xng5W6OztJDsK0R0TQmDxPZGoqtXiX7+elx2HiJqBRToZJUVR8GlMEkI87DCCg/pQM9zdqT0crC2w7gi7vBPR/9TW6xCTkIfh4Z5QqYTsOETNEuppj1n9A7H2SDpOZ5TIjkNEd4hFOhmlfRcLcDarFPMHh/AiiprFWqPGxO7e2Hw6G6XVdbLjEJGBOJxSiLKaeoyI4PPoZNz+NLID3Ows8erPZ6AofLSLyJiwSCejtDwmEZ4OVpjYw1t2FDJi06L8UFOvw6aTWbKjEJGB2BGXC0sLFQaEusmOQtQijtYavDAmHMfSivHjiUzZcYjoDrBIJ6NzJrMEsRfy8diAIFhZqGXHISPWxccJ4e0d8D27vBMRGh6l2hGfgwEhbrC1tJAdh6jF7r/LF918nfDW5nj2GiMyIizSyeh8FpsEeysLPNTHX3YUMnJCCEyN8sPJ9GKcv1QmOw4RSZaYV4HUgkoMZ1d3MhEqlcDSSV1QUF6D9ziIHJHRYJFORiW9sBK/nMrGg7394GSjkR2HTMCk7t7QqAXWHUmXHYWIJNsRlwMAGM4BScmEdPF1wiP9AvHVgVScyiiWHYeIbgOLdDIqn+9NhgAwe2CQ7ChkItzsrTAyoh02Hs9Ebb1OdhwikmhHfC4ivBzh42wjOwqRXi0c3REe9lZYtPE06rVs64gMHYt0MhpFFbVYezgdE7v7wMuJF1CkP9Oi/FBQUYud8bmyoxCRJMWVtTiaWsRpPckkOVprsHhCJM5kluLrA6my4xDRLbBIJ6PxzYFUVNVpMW9wsOwoZGIGdXBHO0crdnknMmN7EvKg1SkYEcEinUzT+C5eGNzRA+/9loBLJdWy4xDRTbBIJ6NQXafFl/tTMCzMA2HtHWTHIRNjoVbhvrt8set8LnJKeeFCZI52xOXC3d4S3XydZUchahVCCLwxsRNqtTq88cs52XGI6CZYpJNR+OFYBgoqajF/SIjsKGSipkX5QacA3x/m3XQic1On1WH3+VwMC/OESiVkxyFqNQFudnhmWCiiT2c3DZRIRIaHRToZPK1OwWcxSejm54w+Qa6y45CJCnS3w6AO7lhzKA1anSI7DhG1oaOpRSitrmdXdzIL84YEo2M7e7y08QzKOHc6kUFikU4G77ezl5BSUIn5g4MhBO9wUOuZ0ccfWSXV2MUB5IjMyo64HFiqVRjYwUN2FKJWZ2Whxrv3dUVuWTXe3hIvOw4RXQeLdDJoiqLgk5gkBLjZ4u5O7WXHIRM3IqIdPB2ssPogR74lMic74nPRJ9gV9lYWsqMQtYke/i6YPSAI3x5Mw++JBbLjENE1WKSTQTuUXIiT6cWYOygYaj4nSK1Mo1Zhem9/7E7IQ3phpew4RNQGkvMrkJRXgZER7WRHIWpTz44OQ4CbLf6+4RSqarWy4xDRFVikk0FbHpMENztL3N/TV3YUMhPTe/lBAFhzKE12FCJqA5cHzxrO+dHJzNhYqvHOlK5ILajEe7+dlx2HiK7AIp0MVkJOGXbG5+LR/oGw1qhlxyEz4e1sg+Hh7fD9kXTU1utkxyGiVrYzPhcd29nDz9VWdhSiNtcvxA0P9fHHF/uScTytSHYcImrEIp0M1qcxSbDRqDGzb4DsKGRmZvT1R355LX47d0l2FCJqRaXVdTiUXIgR7OpOZuzFseFo52iNF9afQk09u70TGQIW6WSQskuq8NOJTDzQyw8udpay45CZGdzBA74uNlh9gF3eiUxZTEIe6nUKRrCrO5kxB2sN3prcBRdyy7Fs50XZcYgILNLJQK3clwKdAjw+MEh2FDJDapXAQ3388XtSAS7mlsuOQ0StZEdcLlxsNejh7yI7CpFUw8I9MbmHDz7enYi47FLZcYjMHot0Mjil1XX49mAaxnfx4jOCJM3Unn7QqAWnYyMyUVqdgl3nczEszJOzhxABWHxPJJxsNHhh/SnUaTkmC5FMLNLJ4Hx7MA3lNfWYNzhYdhQyYx4OVhjb2Qvrj2SgvKZedhwi0rPjaUUorqzj8+hEjVzsLPHGpM44nVmCj3clyo5DZNZYpJNBqanX4ou9yRgY6o7OPk6y45CZe2xAIMpq6rHhWIbsKESkZ9vjcmGhEhjU0V12FCKDMa6LF+7t5o3/7LyAM5klsuMQmS0W6WRQfjqehdyyGswfwrvoJF8Pfxd083PGl/tSoNMpsuMQkR7tjM9Bn2BXOFprZEchMiivT+wEVztLLPz+BKrrONo7kQws0slg6HQKPo1NQqSXIwaG8s4GGYbZAwKRlF+BmAt5sqMQkZ6kF1YiIaccw8PZ1Z3oWs62lnj3/q5IyCnHB9sTZMchMkss0slg7IzPxcXccswfEgwhOIgPGYaxnb3g6WCFlftSZEchIj3ZEZcDABgZwanXiK5nWJgnHuzth09jknAkpVB2HCKzwyKdDMbymET4ONtgXBcv2VGImlhaqPBw3wDsSchDYh6nYyMyBTvicxHiYYcANzvZUYgM1kvjI+HrYoNn151Eu0F7hQAAIABJREFUZS0HUCVqSyzSySAcTS3C4ZQizBkUBI2a/yzJsDzUxx+WahVW7U+RHYWIWqi8ph4Hkgo4qjvRLdhbWeCf93dDWmEl3t4cLzsOkVlhNUQG4dOYRDjZaDAtyk92FKI/cLe3woRu3lh/NAOl1XWy4xBRC8Qm5KFOq2B4OLu6E91K32A3zB4QhK8PpCKWY7MQtRkW6SRdUl45fjuXg0f6BcDOykJ2HKLremxAICprtfj+cLrsKETUAtvjcuFsq0FUgIvsKERG4fm7wxDqaY8X1p9CSRV/UU3UFlikk3SfxiTBUq3C/7d33/FNlfsfwD9PutM96KJ00NKWvcGy93JeFUVBGVdRVEARxZ/CtaDXgV7RKwoqCIoyFBQQ2bL3LrRAGWV00L13kzy/P1J6qcymSU/afN6v13m1TU+Sz5Pk5Mk35znPeTYqWOkoRLfVqrErOge744f9l6Hl6diI6iWtTmLb2TT0i/CGNQ+tIron9jZW+M/wtkgvKEP02jil4xBZBPZQpKi0/FL8diwZwzsFoJGzndJxiO5oXPcQJGaXYHNcqtJRiMgAx67mIKe4AgNa8Hh0oppo28QNE/uF4ffjyVgbk6J0HKIGj0U6Ker7PZeg0ekwvmeo0lGI7mpQS1+EeDli/s6LkJJ704nqmy2n02BrpUKv8EZKRyGqd17pG4YOgW6Y/vsppOSWKB2HqEFjkU6KySuuwE8HruCBNv4I9FQrHYforqxUAs/3bIqYpDzsT8hSOg4R1dDW02m4L9QTTpz/hKjGrK1UmPNkO2h1ElN+OcFDv4hMiEU6Keang1dQVK7FC72bKh3FpFavXo02bdrAzs4OISEh+Oyzz+64/muvvQYhBKZOnVrt8rNnz6Jr165wdXXFiBEjUFhY/Zzdu3btQuPGjW+6/FYWL14MIcQt142OjoaXl1fV35cvX4YQompxdHREaGgoRo4cid27d990/TFjxqBTp053zVBfPdqhMbyc7PDNzgSloxBRDVzMKERCZhEGNues7sbGfs5yBHk64t2HWuJAQja+281+kMhUWKSTIkortFi09xJ6hzdCS39XpeOYzN69e/Hoo4+iS5cu+OOPPzBu3DhMmzYNn3/++S3XP336NBYuXAgXF5eb/jdmzBiEhYXhl19+wenTp/HBBx9U/U+n02Hy5Mn48MMP4eTkZJK2fPrpp9i/fz/Wr1+PGTNmICsrC7169cLMmTNNcn/myt7GCmO7B2PnuQycTslXOg4R3aOtp9MAgOdHNzL2c5ZneMcADG3li/9sjkdscp7ScYgaJBbppIhfjyYhs7AcL/Zu2Meiz5o1C927d8eCBQswaNAgzJgxA5MmTcKsWbNQXl5+0/oTJ07E5MmT4e5e/dRAhYWFOHjwID7//HMMHjwY77zzDrZs2VL1/++//x42NjZ45plnTNaWiIgI3HfffejduzfGjBmDjRs3YsaMGYiOjsaOHTtMdr/maFTXIDjaWuHbXReVjkJE92jrmTS09HeBv5uD0lEaFPZzlkcIgQ/+0RoejraYvPw4Ssq1SkcianBYpFOd02h1+HbXRbRr4ob7mnooHcekTpw4gYEDB1a7bNCgQcjJycH+/furXb5y5UqcPXsWb7311k23c/2DjoOD/sOlWq2uuiw/Px/Tp0/HF198ASGEKZpxW++++y78/f0xf/78Or1fpbmqbfB010D8cfIaErOLlY5DRHeRVViGo1dyMIB70Y2O/Zxlcne0xX+Gt8PFjCJ8sP6M0nGIGhwW6VTn1semIjG7BBP6hNZ5Z1vXSktLYWtrW+2y63+fOfO/Tq2kpASvv/46PvroIzg6Ot50Ox4eHggODsaXX36J7OxsfPvtt1XHw7333nsYMGAAoqKiapxPq9VCo9FUW3Q63T1f38rKCv369cOBAwdqfN/13bgeIVAJYOGeS0pHIaK72B6fAZ0EBvLUa0bHfs5y9Wjmhed6hGDJgSvYdjZN6ThEDQqnN6U6JaXEvB0XEdrIEQMtYI9GWFgYDh8+XO2yQ4cOAQCys7OrLvvwww/h5+eHUaNG3fa2vv76awwfPhxvv/02mjVrhq+++goXLlzAggULcOrUKYPyubm53fJyT0/Pe76NgIAApKVZXufs5+qAh9s1xvLDVzGpfzN4ONre/UpEpIitp9Pg62KPlv43HwdNtcN+zrK9MSQCey5k4s2VJ7Fhci80crZTOhJRg8A96VSndp3PxJlr+XihdyhUqoa9Fx0AXnzxRaxevRrfffcdcnJysGnTpqpZb1Uq/eZ36dIlfPrpp3cdxjd06FCkp6cjPj4eZ86cQWBgIKZMmYLXXnsNAQEB+OqrrxAYGIjAwEB8/fXX95Rv165dOHz4cLXl+eefr1EbLfl84S/2borSCh0W77usdBQiuo3SCi12nc/AgBbeDX70lhLYz1k2O2sr/Pep9igo1eD1X2Og42nZiIyCe9KpTs3bcQG+LvZ4pF1jpaPUiXHjxiEmJgYTJkzA+PHjoVar8fHHH2PixInw9fUFALz11lsYOnQoIiIikJubC0A/i21ZWRlyc3Ph6upa9aFGrVYjPDwcALBlyxbExMRgxYoViImJwYwZM7Bv3z4AQFRUFHr06IE2bdrcMV/79u1vmiV33bp1NWpjcnIyfHwa/qiIWwnzdsaQlr5YtOcS/tk9BK5qG6UjEdHf7E/IQnG5lsejmwj7OQr3ccb0B1pgxupYLNxzCc/3atin1iWqC9yTTnXm+NUcHEjIxnM9Q2BrbRkvPSsrK8ydOxcZGRk4efIk0tLScN999wFA1c/4+Hj89ttvcHd3r1oSExMxd+5cuLu7Izk5+abb1Wq1eO211zB79mw4ODhgx44d6NevHyIjIxEZGYn+/ftj586dJm+fRqPBtm3bDDpOsKGYPKAZCso0WLiH54slMkdbT6fB0dYKUaH3PryZ7h37OQKAUV0DMbilD2ZvOouTSblKxyGq97gnnerM/J0X4WJvjRFdApWOUueufygB9MfcdevWDZGRkQCABQsWoLCwsNr6I0aMQO/evTFhwgQ0atToptubN28e3N3d8eSTT1ZdVlz8v1nGi4qK6mR43qxZs5CSkoIXX3zR5Pdlrpr7uWBYa198v/cyxvUIgZuax6YTmQudTmLL6TT0Cm8EO2srpeM0aOznLJsQAh8/1gbDvtiNScuOY92knnCyY5lBZChuPVQn4lMLsCkuDRP7hVnUm/aBAwewZ88etGvXDvn5+Vi2bBk2bdqEPXv2VK1zffbaG9nb26NJkybo06fPTf/Lzs7GzJkzsWnTpqrLevXqhTfffBPff/89pJTYtm0bPvroI6O2JT4+Hl5eXigvL8elS5ewfPlybNy4EdHR0ejdu7dR76u+mdw/HBtiU7Fg9yVMHRyhdBwiqnQ8MQfpBWUY0spX6SgNFvs5us5NbYvPR7THiG/341+rY/HZk+2UjkRUb1lOtUSKmrv9AtS2VhjXPUTpKHXKxsYGK1asQHR0NFQqFXr27Im9e/eidevWBt9mdHQ0HnroIXTo0KHqsvbt22P27Nl45513AACffvop2rZtW+v8N5o6dSoA/QcrPz8/REVFYdeuXejZs6dR76c+ivB1xrDWfli09xL+2SME7pzpncgsbIxNha2VCv0ivZWO0mCxn6MbdQnxwOT+4Ziz9Rx6NPPCox0ClI5EVC+JuwwV4hSNVGsXMwox4LOdeKFXKN4aGql0HCKTOJdWgMGf78KE3qF4cwhf52QytZ2e3GL6dSkles7ejmbeTlg0tovScYgshlYn8dR3BxCXnId1k3oixMtR6UhE5uyW/bplzN5Fivpq+wXYWavwXE/L2otOliXcxxn3t/bDD/suI7uoXOk4RBYvLiUfSTklGNrKT+koRBbFSiXwxYh2sLFWYdKy4yjX6JSORFTvsEgnk7qSVYQ1J1IwqmsQvJzslI5DZFKT+zdDcYUW3+7iTO9EStsYmwqVAAa04KmziOqan6sDZj/WBqeS8zB741ml4xDVOyzSyaS+3n4RViqB8TxnJlmAZj7OeLCNP37Ydxnp+aVKxyGyaBvjUtE1xBMenCOCSBGDWvri2aggLNhzCZvjUpWOQ1SvsEgnk0nKKcaqY0l4uksgvF3slY5DVCemDAyHRqfDnK3nlI5CZLEupBfgQnohhrbmrO5ESnrn/uZoE+CK13+NwZWsIqXjENUbLNLJZObtuAiVEHihN/eik+UI9nLEyK5BWHE4EefSCpSOQ2SRNsbq99oNasEinUhJdtZW+OrpDlAJgQk/HUNphVbpSET1Aot0MonUvFL8eiQJwzsFwM/VQek4RHVqUv9mcLS1xkcbeBwekRI2xqWiQ6AbfF05iotIaU081JjzZFucvpaP6LVxSschqhdYpJNJzN1+Hjop8WLvUKWjENU5D0dbvNQ3DNvOpmPfxUyl4xBZlMTsYsQm52NIK+5FJzIX/SJ98HLfUCw/nIhfjyQqHYfI7LFIJ6NLzC7G8kOJGNGlCZp4qJWOQ6SIsd2D0djNAR+sPwOdzmJOTU2kuE2VE1QNaclTrxGZkykDI9At1BPTV8ciLiVP6ThEZo1FOhnd51vPw0ol8ErfZkpHIVKMvY0Vpg4OR2xyPtbGpCgdh8hibIhNRQs/FwR68ktiInOiP396e7irbTH+x6PILCxTOhKR2WKRTkZ1Ib0Avx9PwjP3BfFYQLJ4D7dtjFaNXfDJpnhOlkNUB9LyS3H0Sg6HuhOZqUbOdvju2U7IKirDi0uOokzDvpHoVlikk1HN2XIeDjZWmNCHx6ITqVQCbw9tjuTcEizcc0npOEQN3rqT1wAAw1pzqDuRuWod4Ir/DG+HI1dyMP33WEjJQ8KI/o5FOhlNbHIe/jx1DeN6hMDTyU7pOERmoVuYFwa39MHcbReQnFuidByiBm1tTApa+rsgzNtJ6ShEdAf3t/HD5P7N8OvRJH6JTXQLLNLJaD7bcg4u9tZ4rifPi050oxkPtICExPvrTisdhajBupxZhJjEXDzczl/pKER0Dyb3b4ZhrX3xwfoz2H42Xek4RGaFRToZxdErOdh2Nh0v9A6Fq4ON0nGIzEqAuxqv9A3DhthU7DqXoXQcogZpbUwKhAAebMsinag+UKkEPh3eFs39XPDy0mM4kZirdCQis8EinWpNSolPNp2Fl5MtxnYPVjoOkVl6vldTBHuqEb02jhPlEBmZlBKrTySjS7AH/FwdlI5DRPdIbWuNRWM6w8vJDmMWHcL5tAKlIxGZBRbpVGvb49NxICEbr/QNg9rWWuk4RGbJztoK0Q+1REJmERbs5vF3RMYUl5KPhIwiPNyusdJRiKiGvF3s8dM/u8LGSoVnFh5CUk6x0pGIFMcinWpFo9Xhg/VnEeLliKe7Bikdh8is9Ynw5iRyRCawNiYFNlYCQ3nqNaJ6KdBTjR/HdUFxuQbPLDzEc6iTxWORTrWy4kgiLqQXYtqQSNha8+VEdDfXJ5Gb9Uec0lGIGgSdTmLtiRT0Dm8Ed0dbpeMQkYGa+7ng+zGdcS2vBKO/P4S8kgqlIxEphlUVGaywTIM5W86hc7A7Brf0UToOUb0Q4K7GpP7NsCkuDRtjrykdh6jeO3Q5G6n5pZwwjqgB6BTsgXmjOuJcWgGeXXiQhTpZLBbpZLBvdl5EZmE53h7WHEIIpeMQ1RvP92yKlv4umL46DrnF5UrHIarX1pxIgYONFQa24JfFRA1B3whvzBvZEWeuFWDUgoPsJ8kisUgng1zLK8F3uxPwYFt/tA90VzoOUb1iY6XC7MfbILe4HO+tO6N0HKJ6q1yjw/pT1zCopQ8nLiVqQAa08ME3z3REfGoBHp+/n5PJkcVhkU4G+c/mc9DpgDcHRygdhaheaunvihd7h2LVsSTsiE9XOg5RvbTrXAbySirwcDsOdSdqaPpGeuOHcV2Qll+KR7/eh7iUPKUjEdUZFulUY7HJeVh1LAljugejiYda6ThE9dYr/cIQ2sgR7/wei8IyjdJxiOqd1SeS4aa2Qc9mjZSOQkQmEBXqiZUvdoOVSuDJbw5g9/kMpSMR1QkW6VQjOp3Ev9bEwkNti5f7hCkdh6hes7exwuzH2yIlrwSzN55VOg5RvZJbXI7NcWl4pF1j2Fjx4wxRQxXh64zfXuqGAHcHjF10GKuOJikdicjk2KtRjaw8loRjV3Px1tBIuKptlI5DVO91DHLHmG7B+HH/FRxIyFI6DlG9sTYmBeVaHYZ3ClA6ChGZmJ+rA355MQqdgz3w+q8xeG/daVRodUrHIjIZFul0z3KLy/HRhrPoGOSOxzrwQxGRsbwxOALBnmq8/ksM8kt5uhmie/HLkUS09HdBS39XpaMQUR1wsbfBD+O6YHRUEBbuuYRRCw4is7BM6VhEJsEine7Zp5vj9bNRP9wKKhVPuUZkLGpba8x5sh1S80sRvSZO6ThEZu90Sj5ik/MxvCO/MCayJLbWKsx8uBX+M7wtTiTm4sEv9+BEYq7SsYiMjkU63ZNTSXn4+eBVPBsVjBb+LkrHIWpw2ge64+W+YfjteDL+PHlN6ThEZu3Xo4mwtVLh4XaNlY5CRAp4rGMAVk3oBpUQeGL+fvx04AqklErHIjIaFul0VzqdxPQ1sfB0tMOUQeFKxyFqsCb2C0PbAFe8s/oU0vJLlY5DZJZKK7T4/XgyBrbwgbujrdJxiEghrRq7Yt3EHogK9cT01bGYuOw4CnjIGDUQLNLprlYcSURMYi7euT8SLvacLI7IVGysVPjsyXYordBi6q8x3CtAdAvrT11DbnEFRnYNVDoKESnM3dEWi8Z0xptDIrAhNhUPfrkHsck8nzrVfyzS6Y7S8kvx4foz6BrigUc4rJDI5EIbOeGd+1tg9/lM/LDvstJxiMzOzwevoqmXI6JCPZWOQkRmQKUSeKlPGJY9fx9KK3R4dN4+LOHwd6rnWKTTbUkp8c7vsSjT6PDRY20gBCeLI6oLo7oGom9EI3yw4SziUrhHgOi6M9fycfRKDp7uGsg+iYiq6RLigT8n9UBUU0/MWB2LVzj8neoxFul0W3+cvIatZ9IwdVAEQrwclY5DZDGEEPh0eFu4q23wytLjKCzTKB2JyCwsPXgVttYqPM5Z3YnoFjyd7KqGv2/k8Heqx1ik0y1lFpbh3TWxaNvEDeN6hCgdh8jieDrZ4b8j2uNKVhHe/u0Uh+2RxSss0+D348l4oI0f3NScMI6Ibu368Pfl4yuHv3/N4e9U/7BIp1t6d20cisq0+OTxNrDiOdGJFNG1qSdeGxCOtTEpWHE4Uek4RIr69UgiCss0eOa+IKWjEFE90DnYA+sn90RUKIe/U/3DIp1usjH2Gv48eQ2T+och3MdZ6ThEFu2lvmHoEeaFd9fG4WxqvtJxiBSh1Uks3ncZHQLd0D7QXek4RFRPeFTO/j5tSCQ2xqbiAQ5/p3qCRTpVk11Ujumr49DS3wUv9A5VOg6RxbNSCcx5sh2c7W3w8s/HUMTj08kCbTubjitZxTz8iohqTKUSmNAnFMvH34cyDn+neoJFOlWRUuKNX2OQX1KBTx5vCxsrvjyIzEEjZzv8d0Q7XMoswhsref50sjzf77kEf1d7DGnpq3QUIqqnrg9/7xbG4e9k/liFUZUf9l3GX2fT8dbQSLTwd1E6DhHdoFuYF6YNicT6U6mYt/Oi0nGI6kxcSh72J2RhdLdgWPPLYyKqBQ9HW3w/msPfyfyxtyMAwOmUfHyw/iz6RXpjbPdgpeMQ0S2M79UUD7Txwyeb4rEjPl3pOER1Yv7OBDjaWmFE50CloxBRA3B9+PuKG4e/77/MUWpkVlikE4rLNZi47Bjc1Db45PE2EIKzuROZIyEEZj/eBhE+zpi07DiuZBUpHYnIpC5lFuHPkykYFRUEV7WN0nGIqAHpVDn8vXuYJ2asicMrSzn8ncwHi3TCrD9OIyGzCHOebAdPJzul4xDRHahtrfHtM50ghMALS46iuJwTyVHD9c3Oi7C2UuGfnDCOiEzAw9EWC68Pf49LxcNz9yI+tUDpWEQs0i3dHzEpWH44ERN6h6J7mJfScYjoHgR6qjH36fY4l1aAKStioNNxiB41PKl5pVh1LAlPdAqAt7O90nGIqIG6Pvx96XNdUVCmwSNf7cXq48lKxyILxyLdgsWl5OHNlSfRMcgdrw0MVzoOEdVAz2aN8M79LbAxLhUfbTyrdBwio5u/8yJ0EnihF08HSkSm17WpJ/6c2AOtA1zx6ooTmL76FMo0WqVjkYVikW6hsgrLMP7Ho3B1sMG8UR14ujWiemhc92CM6RaMb3clYMn+y0rHITKapJxiLD14FcM7BqCJh1rpOERkIbxd7LH0ua54oVdT/HTgKp6Yvx9JOcVKxyILxMrMAlVodXjp52PILCzDt8925DBConpKCIEZD7TAgObeeHdtHLadTVM6EpFR/Pev8wCASf2bKZyEiCyNtZUK/zesOeaP6oiEjCI88OUebOcZVaiOsUi3QDP/iMPBS9n4+LE2aBPgpnQcIqoFK5XAf59qj5b+rnhl6XGe75XqvYsZhVh5NAmj7guCv5uD0nGIyEINaeWLtRN7wNfFHmMXHcb7605z+DvVGRbpFubng1fw04GreKFXUzzSvrHScYjICNS21lg4uhPc1bYYt/gwrmZxaB7VX//ZHA97Gyu81JfHohORskK8HLH65e54NioIC/ZcwmPz9iEho1DpWGQBWKRbkJ3nMvDumjj0iWiEN4dEKh2HiIzI28Uei8Z2RrlWh5ELDyA1r1TpSEQ1tv9iFtafSsX4Xk3hxVOCEpEZsLexwqyHW+HbZzoiKacED3y5B78cSYSUPLMKmY64ywuMr74G4mRSLkZ8ewBBno5Y8cJ9cLG3UToSEZlATGIuRi44CF9Xe6wYfx88Weg0NKKW1zfbfl2j1eGBL/egoFSDrVN6w8HWSulIRETVpOaV4tUVx3EgIRsDmvvgg3+0grcL53aiWrllv8496RbgYkYhxi46DHe1LX4Y2/meC/Rr167hjTfeQNu2beHk5IQmTZpg9OjRSElJMXFiovrv9OnT6N+/P9RqNfz9/fGvf/0LWu2dj2U7fPgwxo4di7CwMKjVakRERGDmzJkoLa2+Vzw4OBhCiFsu3tYlWDi6ExKzi/Hs94eQX1phymYS1djtto1lhxNxNrUAbw9rfssCPS8vD2PHjoW7uztcXV0xcuRIZGVlVVtHSol///vfCAwMhL29PTp06IBNmzZVWycuLg5DhgyBv78/7OzsEBgYiOeeew7Xrl0zabuJ6NYM6S8vX758yz5wxIgRt73OmjVrIIRAp06dDM7q62qPn5+7D9Pvb47d5zMw4LOd+O1YEveqk9FZKx2ATOtiRiGe+vYAhAB+/GeXGn3bd/ToUfz+++947rnn0LVrV6SlpSE6OhrdunVDbGwsnJycTJicqP7KycnBgAED0KJFC6xZswYXL17E66+/Dp1Oh/fff/+211uxYgUuXryIadOmoVmzZjh58iRmzJiBkydPYtWqVVXr/f777ygrK6t23eeffx5WVlbw8/ODH4D5z3TE+B+PYNyiw/jxn12gtuXbPSnvdttGYUk5tjn3Q9cQDwxr7XvL6z7xxBM4d+4cFixYAJVKhWnTpuGRRx7B7t27q9b56KOPMGvWLMyaNQvt2rXDTz/9hAcffBB79+5F586dAeiL/ZCQEDz77LPw9/fHpUuXMHPmTBw9ehSHDx+GtTW3FaK6Ymh/ed2nn36K7t27V/3t5eV1y/VKS0vx2muvwcfHp9aZrVQCz/Vsin6R3nhz5UlM+SUGf568hn//ozV8XblXnYxESnmnheqxC+kFsvP7W2TH9zbL+NT8Gl8/JydHVlRUVLssPj5eApCLFy82VkyiBueDDz6Qbm5uMi8vr+qyjz/+WDo4OFS77O8yMjJuuuybb76RAOTly5dve71r165JKysr+dFHH1W7fF1Migx5a50cPm+fzC8pN6AlZIbu1m+bdb9+u23D2tZONn19pTx3m75q3759EoDcuXNn1WUHDx6UAOSWLVuklFKWlZVJZ2dnOX369GrX7dChg7z//vvvmGvz5s0SgDx69KihTSMiAxjaX166dEkCkH/88cc93c+sWbNkjx495OjRo2XHjh1rnfs6jVYnF+xOkBHT18sWMzbI73ZdlBUardFunyzCLftrDndvoBIq96BrdRJLn78P4T7ONb4NNze3m/YohIeHQ61Wc8g70R1s2LABgwcPhouLS9VlI0aMQElJCXbu3Hnb691qD0D79u0B4I7b3C+//AKdTnfTML/72/jh8xHtcfRqDkYtPITc4vKaNoXIqG61bXi36wtNeRkGumei2W36qg0bNsDHxwe9evWquqxLly4ICQnBhg0bAAAXL15EQUEBBg4cWO26gwYNwpYtW1BefvvXv6enJwDccR0iMj5D+8uauHr1KmbPno0vvvjCKLd3IyuVwD97hGDTq73QJcQD7/95Bg98uQeHL2cb/b7IsrBIb4ASMgoxorJAXzbesAL9dk6ePIni4mKEh4cb7TaJGpqzZ88iMrL6GRQCAwOhVqtx9uzZGt3W/v37oVKpEBp6+9NRLV++HFFRUQgKCrrpfw+19ce8kR1wJiUfI749gMzCslvcAlHd+Pu2kVlYhrmHcmFla48Q67x7vt51zZs3r9qmrs/dYGtrW20dW1tblJeXIyEhodrlOp0O5eXliI+Px1tvvYXOnTujS5cuBreNiGqutv3l2LFjqw71mjJlCkpKSm5a5/XXX8cTTzyBDh06GC333wV5OuL7MZ0xf1RH5JdUYPj8/Xj9lxik5fNMK2QYFukNzLGrOXh8/n6TFOg6nQ6TJ09Gs2bN8NBDDxntdokampycHLi5ud10ubu7O3Jycu75dlJTU/H+++/jmWeegbe39y3XuXLlCg4cOHDHyXIGtfTFwjGdcDmrCE98s5+nZyPF3LhtaHUSry4/gfzYoKXgAAAcXUlEQVRSDbw8PZCXl3tP17vRjdtU06ZNIYTA4cOHq61z6NAhAEB2dvU9W8OGDYOdnR0iIyORnZ2NdevWQaXixyKiumRof2lnZ4eXX34ZCxcuxF9//YUXXngB8+bNu6kv3LZtGzZv3owPPvjA6Nn/TgiBIa18sfX13pjQJxRrY5LR55MdmLPlHIrKNCa/f2pY2Bs1IBtOXcNT3x6As701fn0xqkYFupQSGo2marnVrJr/93//h/3792PJkiWwseEp3IhMqby8HE888QScnJwwZ86c2663fPlyqFQqPPHEE3e8vZ7NGuHHcV2Rnl+Gx+btw7m0AmNHJqqRudsuYM+FTMx6qCWsVbU9sxzg6uqKp556Cv/+97+xfft2ZGdn48svv8TWrVsB4KYC/Msvv8SBAwewZMkSFBYWYujQoTedSYGIzJOfnx/mzp2Lhx56CH369EF0dDQ+++wzrF27FjExMQAAjUaDSZMm4Z133jHKhHH3Sm1rjWlDIvHXlD7o19wbX/x1Hn0+3YHlh65Cq+Ms8HRvWKQ3AFJKfLcrAS8tPYaW/i74bUI3NG1Us5nXd+7cCRsbm6qlf//+1f7/9ddf45NPPsEPP/yArl27GjM+UYPj7u6OvLybh+7m5OTA3d39rteXUuLZZ59FXFwc1q9ff8frLF++HH379r2nDyBdQjyw7Pn7UK7V4bGv92HXuYy7XofImK5vG9vj0/H5X+fwaPvGeLJzk7tuG/e6TX3++edo0aIF+vXrB09PT3zyySeYPn06AMDXt/qs8c2aNUPXrl0xatQobNq0CcePH8fSpUuN1FIiuhe17S9v9PjjjwPQn50IAL777jvk5eVhzJgxyM3NRW5uLsrLy6HVapGbm4uKCtOeojTQU42vnu6AVRO6IdBDjbd+O4VhX+zGTva9dA9YpNdzFVod/rUmDv9efwZDW/li6fP3wdPJrsa307FjRxw+fLhq+eabb6r+t2rVKkycOBGzZ8/Gk08+acz4RA1SZGTkTcfSJSYmori4+JbH1f7dq6++ijVr1mDNmjV3XD8+Ph4nTpzAU089dc/ZWge4YvXL3dHY3QFjFx/GkgNX7vm6RLUVGRmJA8dO4eWfjyHS1wXv/6MVkpKS7rpt3GqbAm4+nrVRo0bYtm0bEhMTERsbi4SEBDg6OsLX1xfBwcG3vf2goCB4eHjcdNw6EZlWbfvLGwkhqv2Mj49HUlISfHx84O7uDnd3dyxbtgwnTpyAu7s7VqxYYZxG3EXHIHesfDEK80Z2QKlGi9HfH8LIBQdw7Oq9H/5GlodFej2Wll+Kp787gCUHruCFXk0x96kOsLexMui2nJ2d0alTp6olIiICALBjxw6MHDkSEydOxNSpU40Zn6jBGjp0KDZt2oSCgv8NKV+xYgUcHBzQu3fvO173ww8/xNy5c/HTTz+hR48ed1x32bJlsLW1xaOPPlqjfI3dHLByQjf0Dm+EGatjMeuP0xyCR3Uiqnd/bNmyGS5WGiwe2xlqW+t72jaGDh2K1NRU7Nmzp+qyI0eOICEhAUOHDr1p/YCAALRs2RIajQbff/89xo0bd8dc8fHxyMrKQkhIiOGNI6Iaq01/+XcrV64EoN/xBACvvPIKtm/fXm0ZPHgwwsPDsX379pvOBGFKQggMbe2HLa/1xr8eaIH41AI8+vU+jFt8GLHJt580kyyXkPKOH8z4qc1MbY9Pxxu/xqC4XIsPH22Nh9s1Nvp9nDlzBlFRUQgODsa8efOqvpkE9Hsr7jTbNJEly8nJQYsWLdCqVStMmzYNCQkJmDJlCl599VW8//77VeuFhYWhd+/eWLhwIQBg6dKlGDlyJMaMGYMXXnih2m2GhoaiUaNG1S5r3rw5IiIisHr1aoNyanUS7/95Gov2XkZUU098MaIdvF3sDbotqjO1PXhbsX79QnoBnpr7F07M+Se6dGiLmTPevudtAwAGDx6M8+fP49NPP4VKpcK0adPg7e2N3bt3V62zZMkSVFRUoGnTprh69SrmzJkDjUaD/fv3w8lJfxjY1KlTYW1tja5du8LNzQ1nzpzB7NmzYW1tjZiYGDg6Otbdg0Jk4QztL6Ojo1FQUIDu3bvDxcUFu3btwieffIJhw4Zh1apVt72/MWPGIDY2FkeOHDF52+6kuFyDH/ZdwfydF5FXUoHBLX3w2sBwRPq63P3K1NDcul+/3QnUKxcyMwWlFXLayhgZNG2dHPTZTnkuNd9k97Vo0SIJ/Qe6m5bRo0eb7H6JGoK4uDjZt29faW9vL319feX06dOlRqOptk5QUFC1bWn06NG33eYWLVpU7brHjx+XAOSyZctqnfWXw1dlxPT1suN7m+Wuc+m1vj0yqbv122bZrx+/miPbzdwkO763Ra7ZfrDG24aUUubk5MgxY8ZIV1dX6ezsLJ966imZkZFRbZ3FixfL8PBwaWdnJ729veX48eNlZmZmtXWWLVsmu3XrJt3d3aWDg4OMiIiQU6ZMuem2iKhuGNJfLlu2THbs2FG6uLhIGxsbGRoaKmfMmCFLS0vveF+jR4+WHTt2NEUzDJJXUi7nbImXrf61UQa/tU6+svSYvJBeoHQsqlu37K+5J70e2X8xC2+sjEFKbgnG9wrFawObwc7asOHtREQ3Op9WgJeXHsP59EK83CcMrw5oBmsrHhFlhurVnnSdTmLxvsv4aONZ+LjY4ad/dkWQJ/dUExHdKLe4HN/uSsDifZdRWqHFP9oHYFL/ML5fWoZb9uss0uuB3OJyzNlyDj/sv4JgTzX+80RbdAzyUDoWETUwJeVaRK+Nw4ojiegQ6IaPH2uDZjU4lSPViXpTpJ9Nzce//zyD3eczMaC5Nz5+rI1BE5sSEVmKzMIyzN9xEUsOXEGFVof72/jjxd5N0dLfVeloZDos0usbjVaHnw9exZyt55BfUoFno4Lx5pAIqG2tlY5GRA3YmhPJeHdtHIrKNHipTxhe6hvKUTvmw6yL9DKNFvsuZmHV0ST8eeoanGyt8ebQSIzqGlhtXhMiIrq99PxSLNx7CT8fuIrCMg36RDTCS33C0CWEO+kaIBbp9cmucxl4b91pnE8vRLdQT8x4oAWa+3EyCSKqG5mFZXhv3WmsOZGCMG8nfPRoa3QK5ocDM2B2RfqC3Qk4eCkbafmluJheiKJyLZzsrDGmWzCe79kUrmobY98lEZFFyCuuwJIDl/H93svILipHpyB3vNQ3FH3CvaFS8YvPBoJFurmTUmL3+UzM3XYBhy5nI8hTjbeHNcegFj7cA0FEitgen47pv8ciObcE97fxw+sDw9G0kZPSsSyZ2RXp01aexInEXPi42iPYU42+kd6Iaupp8ClBiYioupJyLX45kohvdyUgObcEwZ5qPBMVjMc7BsDVgV+E1nMs0s2VlBJbz6Rj7rbziEnKg5+rPcb3aoqnuwZyiCkRKa6oTIP5Oy9i4Z5LKNPoMLxjACb1bwZ/Nwelo1kisyvSiYioblRodVh/6hp+3H8FR6/kwMHGCo+098fIrkFo6e/CnXr1E4t0c5NXXIFVx5Lw88EruJhRhCYeDnipTxge7dCYxTkRmZ2MgjJ8tf0Clh68CgjgyU5NMLpbEMK8OblcHWKRTkREiE3Ow5L9V7D6RDLKNDqE+zjh0Q4BeKRdY/i62isdj+4di3RzIKXE8cRcLD14FX/EpKBMo0P7QDc8GxWEB9v485RHRGT2knKK8d+/zmP18RSUa3XoHuaJZ6OCMaC5D6x4jJypsUgnIqIqucXlWHfyGn4/noyjV3IgBBDV1BODWvigf3MfNPFQKx2R7oxFupLOpRVg7YkU/HEyBVeyiqG2tcIj7Rvj6S6BaNWYp1Ugovons7AMKw4n4ucDV5CSVwp/V3vc38YP97fxR9sAVw67Mw0W6UREdEuXM4vw+/Fk/HEyBQkZRQCACB9n9G/ujR7NvNC+iTscbDla18ywSK9LOp1EXEo+/jqbho2xqTibWgCVALqFeuGhtv4Y2toXzvac6IGI6j+NVoetZ9Lxy5FE7D6fgQqtRGM3Bwxr7Yv+zX3QPtCNh/AYD4t0IiK6q0uZRfjrTBq2nknD4cs50OokrFUCrRq7okuIBzoEuqOlvwsC3B34pbqyWKSbWl5JBQ4mZGF7fDr+OpOO9IIyCAF0CHTHg238MKyNH7ydeYwIETVceSUV2HI6DetPXasq2O1tVOgS4onuoZ6ICvVEcz8X2PDQHkOxSCciohrJK6nAsSs5OHQ5G4cvZeNkUh7KtToAgLOdNSL9nNHczwVNvRwR5OmIQE81mrirYWvNvroOsEg3toLSChy/mot9F7Ow/2ImTiXnQScBJztr9Ar3Qv9IH/SJaARPJzuloxIR1bn80gocTMjG3guZ2HMhExfSCwEAdtYqtPR3QdsmbmjXxA1tAtwQ7KnmN/n3hkU6ERHVSmmFFqev5eNM1VKAs9fyUVSurVpHJQA/VwcEeaoR6KFGYzcHNHZ3gL+bAxq7OcDX1Z5fuBsHi/Ta0OokzqcX4PjVXJy4movjiTk4n14IKQEbK4H2TdwRVbmXqEOge4P65ik6OhozZ85UOgaRxXj33XcRHR2tdAyjS8svxcFL2TiZmIuTSXk4lZyHkgr9BwJXBxu09HdBuI8zInz1S7iPM5zsrBVObXZYpN+A/ROR5WqofaVSpJTIKCzD1axiXMkqxpXsYlzNKsKV7GIkZhcjs7C82vpCAD7O9lWFu7+bPQLcrv+uL+hdeGjvvbhlv85PP7eQV1KB82kFiE8rwLnUApxNLUBscl7Vt0tuahu0b+KGB9r4o32gGzoGuUNty4eSiOhOfFzs8VBbfzzU1h+A/lj28+mFiEnMRUxSLk5fK8AvRxJRfMM3+QHuDojwcUa4rzMifZ0R2sgJwV6OLN6JiIiMSAgBb2d7eDvbo1Owx03/L63Q4lpeKVJyS5CcU4LkXP2SkluCU0m52BRbWjWE/jpnO+tqRXxjN3XlT30R7+1sz7PC3IbFfsrJL61AUnYJEnOKkZRTgsTsYiRkFuFcagFS80ur1nO0tUIzH2c83jEA7QLd0L6JO4I4LJOIqNasrVRo7ueC5n4uGNElEIB+0s3k3BKcTS3AuTT9l6TnUguw81wGNLr/7QT2crJDiJcaIV6OCPZyRIin/mewpyNnriUiIjIyexsrhHg5IsTL8Zb/1+kkMovKkJxTgpTcymL++pJTgmNXc5BbXFHtOtYqAV9Xe/i7OVTthW/srh9Of31YvaX26Q1muLuUEkXlWhSUViC/RIOswjJkFJYhs7AcmYVlyCwoQ2ZhGdILypCYXYz8Uk2166ttrRDs6YhIX/0em3AfJ4T7OKOxG2c8JCJSWrlGh0uZRbiUWYiEzCJczizC5cxiXMoqQkZBWbV1vZzs0NjNHn6uN357rz9+zsvJDp5OtvV59BOHuxMRUb1UVKa5qXhPydUX9cm5JUjNL4VWV72b8nC0rSza7eHjou/HGznb3fDTFl5OdrC3qbfFvHkck34+rQCnr+VDq5P/W6SErvJ3jU5CJyUqtBJlFVqUanQoKdeitPL30gr972UVOhSVa5BfWZQXlmluelKvs1YJeFY+gd7OdghwVyPA3QFNPCp/uqvhprZhMU5EVA8VlFbgSlZxZRFfpO/wK4fkpeSWVBs+f529jQqejvqC3dPRFq4ONnC0s4aTvTWcbK31v9vpf9paq2BtJWBrpYKN1f9+t7YSsLFSwUalghD64/NUQkAIwN7aCu6OtqZoLot0IiJqkDRaHdIKyqqK9+TcEiTd8HtGQRnySipueV1baxVc7G3g4mANVwebyt9t4GRnBTtrKzjYWsHe2gr2Nqqq361UAtZWAlYqAStR+bNysVapoFIBAvp+vUnl5HkmYB5F+tc7LmD2xvh7Xt/eRgV7G/0D6WBrBTvryr9tVFDbWsPF3houlU+Ec+XvzvbW8HC0RSMn/bcsrg42UPF4ByIiiyOlRF5Jhf4b+rxSZBWVI7uoHFmFZcgqLEdWUTmyivSdflGZFoVlGpRrdHe/4bvoH+mNhWM6G6EFN2GRTkREFqtMo0VW5UjpjMqR0pmF5VU7bvNLKip/r9D37eValJZrUarRokJreBc4bUgkJvQJNWJLqphHkZ5TVI7s4vKbvq1Q3fTthYCdtYp7t4mIqE5VaHUoKtOP0Coq06Jco0OFTocKjQ4anUS5VgeNVqJCq6tcJKSUkBKQkNBJwN/NAb3DG5kiHot0IiIiA2h1EqUVWpRUjsyuGsVd+VP7t5HdGq2ErOw2Az3UCHBXmyKWeRTpREREZDAW6URERA3HLfv1hnMybyIiIiIiIqJ6jkU6ERERERERkZlgkU5ERERERERkJlikExEREREREZkJFulEREREREREZoJFOhEREREREZGZYJFOREREREREZCZYpBMRERERERGZCRbpRERERERERGaCRToRERERERGRmWCRTkRERERERGQmhJTytv+cOXPmRgBedRfnnvkDSFE6hAXj468cPvbK4WOvLD7+epnvvvvuEEOvbMb9urFYwuvEEtoIWEY72caGwRLaCFhGO5Vo4637dSllvVuio6Ol0hkseeHjz8feEhc+9nz8uZj/YgmvE0too6W0k21sGIsltNFS2mlObeRwdyIiIiIiIiIzUV+L9JlKB7BwfPyVw8deOXzslcXHn+6FJbxOLKGNgGW0k21sGCyhjYBltNNs2njHY9KJiIiIiIiIqO7U1z3pRERERERERA0Oi3QiIiIiIiIiM8EinYiIiIiIiMhM1OsiXQjxnRDiohCiRAiRIYRYI4RornQuSyCE8BBCfCmEOFv5+CcKIeYJITyVzmYJhBDjhRDbhRC5QggphAhWOlNDJoR4SQhxSQhRKoQ4KoToqXQmSyCE6CWEWCuESK58nY9ROhMpqybbohDCTwixtLKf0gohFtdhVIPVsI2PCiE2V34GKhBCHBRCPFSXeQ1Rwzb2FkLsE0JkVX7eOCuEmFqXeQ1laN8hhOghhNAIIWJNnbG2avhc9ql8L//7ElmXmWuqps+jEMJWCDGr8jplQoirQohJdZXXUDV8Lhff5rksqsvMNWXAc/m0EOKEEKJYCJEqhPhJCOFbF1nrdZEO4AiAMQCaAxgMQADYKoSwUTKUhfAH0BjAmwBaAxgFoBeAZUqGsiBqAJsBRCuco8ETQjwJ4AsAHwBoD2AfgA1CiEBFg1kGJwCxACYDKFE4CynMgG3RDkAmgI8AHKyTkLVkQBt7A9gG4P7K9dcD+N2cv0g0oI2FAP4L/WeMFgDeBzBTCPFSHcQ1mKF9hxDCHcCPAP4yechaqkX/2BKA3w3LeVPmrA0D27gcwBAA4wFEABgO4KSJo9aKAe2cjOrPoR+ABAC/mD6tYWraRiFEdwBLAPwA/Wv2Eejfg36uk7wNaXZ3IUQbADEAIqWU8UrnsTRCiGEA1gFwk1LmK53HEgghOgE4DCBESnlZ4TgNkhDiIICTUsrnb7jsPICVUsr/Uy6ZZRFCFAJ4RUq5WOkspIzabItCiHUAMqWUY0ybsnaM8X4jhDgEYLeU8nUTxawVI7XxNwBlUsqnTBSz1gxtZ2XbYqDf8fS4lLKVycMaqKZtFEL0AbAdQCMpZWadBa0FA9o4CMCvAELrSxuB2m+XlQXtHgDdpZT7TJfUcAY8l1MBTJRSBt1w2VgAX0opnUydt77vSa8ihHAEMBbAVQCXlU1jsVwAlAEoVjoIkTEIIWwBdIR+1MKNNgPoVveJiCyTJWyLRmyjM4AcY+UyJmO0UQjRvnLdncZNZzyGtrNydIAP9KMFzFotn8sjQohrQoi/hBB9TRLQCAxs4yPQ7zyZIoRIEkKcF0L8Vwhh8qLOUEZ673keQJwZF+iGtHEvAD8hxINCzwvACOhHLJlcvS/SK48tKIR+ONRQAP2llGUKx7I4Qgg3AO8B+E5KqVE6D5GReAGwApD2t8vTANTJMUlEBMAytsVat1EI8TKAAOiHaJojg9tYWfCUQX+o49dSyvmmiWgUNW6nEKI1gHcBjJJSak0bzygMeS6vAZgA4DEAjwKIB/CXGR+eYUgbmwLoAaAt9O18Bfqh74tNE9EoavXeI4RwBfAEgO+MH81oatxGKeV+6IvynwGUA8iAfoTLaNPF/B+zK9KFEO/fZiKCG5c+N1zlZ+iPK+gN4ByAX4UQaiWyNwQGPP6o/HbwDwDJ0B+jTgYw5LEnIiICACHEYwA+AfC0lPKK0nlMoCeATgBeBPCqEOIZhfMYjRDCDsAKAFOllJeUzmMqUsp4KeV8KeVRKeV+KeVLADYCeEPpbEakAiCh3w4PSik3QV+oPyaE8FE2msmMgr7d5vrloEGEEC0AfAn9TsiO0H/Z4gvgm7q4f+u6uJMa+hzAT3dZ5+r1X6SUeQDyAJwXQhyAfojXY2hgL5Q6VKPHv7JAvz7s4wEpZampglmAGj32VCcyAWihH354Ix8AqXUfh8hiWcK2aHAbhRCPQz/Z2LNSyj9ME88oDG7jDcXrqcpiJxrm+1mvpu30g34S5EVCiEWVl6kACCGEBsAwKeXfh+kqzVjb5EHo91aaI0PaeA1AcmV9ct2Zyp+BuHlPrjmo7XP5PIBVUspsYwczIkPa+H8ADkkpP6n8+6TQz16/WwjxtpQyyTRR9cyuSK+cZMHQiRZE5WJnvESWpSaPvxDCGcAG6B/zIVLKQlNma+hq+donE5BSlgshjgIYCP1EMNcNBLBKmVRElscStkVD2yiEeAL62YdHSylXmjZl7RjxeVTBjD/rGdDOZOjPlHOjlyrX/wfMcK4lIz6X7aAvbM2OgW3cC2C4EMLphs/F4ZU/zXKES22eSyFEF+iH9r9quoS1Z2Ab1dAX9je6/rfJR6ObXZF+r4QQYdDvMd8K/TECAQDegn7isnUKRrMIlQX6Zugni3sEgGPl5H0AkC2lLFcsnAUQ+nM0+uJ/b/wtKucFuGrm32TWR58BWFI5Y/Je6Ida+gMw5+MhG4TKkTphlX+qAAQKIdpB/x7DUSWW547bohDiRwCQUj57/QqVrxdA31fpKv8ul1KersvgNVCjNgohRkC/N3kqgF3if+fvLTfjvqCmbZwI4BL0xy8D+lOxTQXwdd3GrrF7bqeUsgL6001WEUKkQz+DvTmfK72mz+Wr0H/hEAfAFvph0o9A/3neXNX0fWcpgBnQj4qIBuAG/Wm/Vkop0+s2eo3U+P210ngA56WUO+ouqsFq2sY/AHwnhJgAYBP0I14+B3CsLj6D1NsiHfpivA+A16HfANIA7AIQJaVsKEPfzFlHAPdV/n7ub//rC2BHnaaxPC9CP8HMdX9W/hwL856cpN6RUq4QQngCmA79G3Qs9EMPzfIb8QamE/Sn67luZuXyA4AxSgQi5dzDtnirc90e/9vfD0K/NyvYVDlrw4A2vgj9Z7nPK5frdkL/GcnsGNBGKwAfQ/+caQBchH6njFl/UWrg67VeMaCNttDPmxAAoAT6Yv1+KWWdzJZtiJq2UUpZKIQYAP2xzIehPwx3NfSvWbNlyOu1cofdCACz6ixoLRjwXC6ubOMrAP4D/eHV2wBMq4u8Deo86URERERERET1mdnN7k5ERERERERkqVikExEREREREZkJFulEREREREREZoJFOhEREREREZGZYJFOREREREREZCZYpBMRERERERGZCRbpRERERERERGaCRToRERERERGRmWCRTkRERERERGQm/h/pLScv8I0jxgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 993.6x331.2 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"az.plot_posterior(trace.approximation.sample(10000))"
]
}
],
"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.7"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment