{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Inference: Metropolis Random Walk MCMC\n", "\n", "This example shows you how to perform Bayesian inference on a time series, using [Metropolis Random Walk MCMC](http://pints.readthedocs.io/en/latest/mcmc_samplers/metropolis_mcmc.html).\n", "\n", "It follows on from the [first sampling example](./sampling-first-example.ipynb)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Running...\n", "Done!\n", "R-hat:\n", "[1.0121663209785179, 1.0091279551338153]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1gAAAEYCAYAAABBWFftAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8XNWV6PvfqkFVmiVL8ijZ8ohnwAM2MU0MxIQQghMCwfRNmk4CBhzycumXz21yb9JNE/NI3qMDzQMS3EDs5CZAhyFhbCYb3BADHjDG84QHedQslaRSTev+cUpzSSrbmmyt7+ejj8581inJLq3ae68tqooxxhhjjDHGmDPn6u8AjDHGGGOMMeZcYQmWMcYYY4wxxvQQS7CMMcYYY4wxpodYgmWMMcYYY4wxPcQSLGOMMcYYY4zpIZZgGWOMMcYYY0wPsQTLGGOMMcYYY3qIJVjGGGOMMcYY00MswTLGGGOMMcaYHuLp7wD6Wn5+vhYXF/d3GMYYYxLYuHFjmaoW9Hcc/eVM3qMCoQDBaJD81PyeDcoYYwyQ/HvUoEuwiouL2bBhQ3+HYYwxJgEROdjfMfSnM3mPWluyln1V+/ju9O/2cFTGGGMg+fco6yJojDHGnCMU7e8QjDFm0LMEyxhjzKAnIleJyC4R2SsidyfY7xORZ+P7PxKR4vj2PBFZIyIBEXmk1fGZIrK51VeZiDzUd09kjDGmvwy6LoLGnI7HNj/W5f5lFyzro0iMMT1NRNzAo8AioARYLyIvqer2Vod9H6hU1QkisgT4JXAjEAR+BkyPfwGgqrXABa3usRF4obefxRhjTP+zBMuYU7RuX3mHbcsuSHCgOeuFw2FKSkoIBoP9Hco5x+/3U1hYiNfr7e9QAC4C9qrqfgAReQZYDLROsBYD98SXnwMeERFR1TrgfRGZ0NnFRWQiMBT4r16I3RhjzABjCZYxxnSipKSEzMxMiouLEZH+DuecoaqUl5dTUlLC2LFj+zscgFHA4VbrJcC8zo5R1YiIVAN5QFkS178JeFZVEw6QEpGlwFKA0aNHn1rkxhhjBpxeG4MlIkXxfuk7RGSbiPwovn2IiLwlInvi33Pj20VEHo73b98iIrNaXevm+PF7ROTmVttni8hn8XMeFvsLyBjTg4LBIHl5eZZc9TARIS8vbyC1DCb6AbdPhpI5pjNLgKc726mqK1R1jqrOKSgYtBXqjTHmnNGbRS4iwP+tqlOA+cAPRGQqcDfwjqpOBN6JrwN8BZgY/1oK/BqchAz4Z5xPEy8C/rkpKYsfs7TVeVf14vMYYwYhS656xwB7XUuAolbrhcDRzo4REQ+QDVR0d2EROR/wqOrGngm1a500khljjOlDvZZgqeoxVd0UX64FduB0sVgMrIoftgr4enx5MfA7dXwI5IjICODLwFuqWqGqlcBbwFXxfVmqui7e7eJ3ra5ljDHGJGs9MFFExopICk6L00vtjnkJaOpBcT2wurMuf+3cRBetV8YYY849fTIGK17O9kLgI2CYqh4DJwkTkaHxwxL1gR/VzfaSBNsT3d/6t5seVVjT7sPoNfc73y/7Sd8HY/rMg2/t7tHr3bVoUrfH3Hffffzxj3/E7Xbjcrl4/PHHmTev/fCgnrNw4UIeeOAB5syZc9rX2LlzJ9/97nfZtGkT9913Hz/+8Y97MMKeFx9TdSfwBuAGnlLVbSJyL7BBVV8CngR+LyJ7cVquljSdLyIHgCwgRUS+DlzZqgLht4Cr++5pjDHG9LdeT7BEJAN4HvjvqlrTRbeQzvq3n+r2jhtVVwArAObMmWP9J4wxZ4V169bxyiuvsGnTJnw+H2VlZYRCof4Oq1tDhgzh4Ycf5s9//nN/h5I0VX0NeK3dtn9qtRwEbujk3OIurjuuh0LsliR8WzTGGNPXenWiYRHx4iRXf1DVpvk/TsS79xH/fjK+vbM+8F1tL0yw3RhjzgnHjh0jPz8fn88HQH5+PiNHjgTg3nvvZe7cuUyfPp2lS5c2j71ZuHAhd911F5deeilTpkxh/fr1XHfddUycOJGf/vSnABw4cIDJkydz8803M3PmTK6//nrq6+s73P/NN9/k4osvZtasWdxwww0EAgEA7r77bqZOncrMmTMTtk4NHTqUuXPnDpQS7MYYY0yf6s0qgoLTpWKHqv6q1a7W/dhvBv7SavvfxasJzgeq410J3wCuFJHceHGLK4E34vtqRWR+/F5/1+paxhhz1rvyyis5fPgwkyZNYtmyZbz33nvN++68807Wr1/P1q1baWho4JVXXmnel5KSwtq1a7n99ttZvHgxjz76KFu3bmXlypWUlzvzuO3atYulS5eyZcsWsrKyeOyxtpNpl5WVsXz5ct5++202bdrEnDlz+NWvfkVFRQUvvvgi27ZtY8uWLc1JmxkYNOnChsYYY3pLb7ZgLQC+A1wuIpvjX1cDvwAWicgeYFF8HZyuGfuBvcC/A8sAVLUC+DnOIOT1wL3xbQB3AE/Ez9kHvN6Lz2OMMX0qIyODjRs3smLFCgoKCrjxxhtZuXIlAGvWrGHevHnMmDGD1atXs23btubzrr32WgBmzJjBtGnTGDFiBD6fj3HjxnH4sDOktaioiAULFgDw7W9/m/fff7/NvT/88EO2b9/OggULuOCCC1i1ahUHDx4kKysLv9/PLbfcwgsvvEBaWlofvBLGGGPM2aPXxmCp6vskHicFcEWC4xX4QSfXegp4KsH2DcD0MwjTGGMGNLfbzcKFC1m4cCEzZsxg1apVLFmyhGXLlrFhwwaKioq455572swp1dSl0OVyNS83rUciEaBjmfT266rKokWLePrpjgXwPv74Y9555x2eeeYZHnnkEVavXt1jz2uMMcac7Xp1DJYxxpjTt2vXLvbs2dO8vnnzZsaMGdOcTOXn5xMIBHjuuedO+dqHDh1i3bp1ADz99NNccsklbfbPnz+fDz74gL179wJQX1/P7t27CQQCVFdXc/XVV/PQQw+xefPm0308Y4wx5pzUJ2XajTHmXJBMWfWeFAgE+OEPf0hVVRUej4cJEyawYsUKcnJyuPXWW5kxYwbFxcXMnTv3lK89ZcoUVq1axW233cbEiRO544472uwvKChg5cqV3HTTTTQ2NgKwfPlyMjMzWbx4McFgEFXlwQcf7HDt48ePM2fOHGpqanC5XDz00ENs376drKys03shTFIawlHqGqP9HYYxxgx6MthmfZ8zZ45u2LChv8MwZ5nHNrcUAFi3r7zDPFj3F8fnJbJ5sM4pO3bsYMqUKf0dRo87cOAA11xzDVu3bu3XOBK9viKyUVVPfxKudkTkdVX9Sk9dr7edyXvUox++zKcndrBi8f/o4aiMMcZA8u9R1oJljDHmrCYiszrbBVzQl7H0N6siaIwx/c8SLGOMGWSKi4v7vfWqh60H3iNxYaWcPo7FGGPMIGcJljHGmLPdDuA2Vd3TfoeIHO6HeIwxxgxiVkXQGGPM2e4eOn8/+2EfxtHPOpsZxRhjTF9KKsESEZtryhhjzICkqs+p6q5O9v25r+MxxhgzuCXbgvUbEflYRJaJiPVnN8YYYwYkK3JhjDH9LakxWKp6iYhMBL4HbBCRj4HfqupbvRqdMcYMJGvu79nrJVHW/7777uOPf/wjbrcbl8vF448/zrx583o2jlYWLlzIAw88wJw5p18p/Q9/+AO//OUvAcjIyODXv/41559/fk+FaIwxxgxoSRe5UNU9IvJTYAPwMHChiAjwP1X1hd4K0BhjBqt169bxyiuvsGnTJnw+H2VlZYRCof4Oq1tjx47lvffeIzc3l9dff52lS5fy0Ucf9eo9RcQFzFfVv/bqjYwxxphuJDsGa6aIPIhTqely4GuqOiW+/GAvxmeMMYPWsWPHyM/Px+fzAZCfn8/IkSMBuPfee5k7dy7Tp09n6dKlNE0av3DhQu666y4uvfRSpkyZwvr167nuuuuYOHEiP/3pTwFnouHJkydz8803M3PmTK6//nrq6+s73P/NN9/k4osvZtasWdxwww0EAgEA7r77bqZOncrMmTP58Y9/3OG8L3zhC+Tm5gIwf/58SkpKev7FaUdVY8C/9vqNjDHGmG4kOwbrEWATcL6q/kBVNwGo6lHgp70VnDHGDGZXXnklhw8fZtKkSSxbtoz33nuved+dd97J+vXr2bp1Kw0NDbzyyivN+1JSUli7di233347ixcv5tFHH2Xr1q2sXLmS8vJyAHbt2sXSpUvZsmULWVlZPPbYY23uXVZWxvLly3n77bfZtGkTc+bM4Ve/+hUVFRW8+OKLbNu2jS1btjQnbZ158skn+cpXvtKDr0qX3hSRb8Z7V5wSEblKRHaJyF4RuTvBfp+IPBvf/5GIFMe354nIGhEJiMgj7c5JEZEVIrJbRHaKyDdP98GSegarImiMMQNCsgnW1cAfVbUBnK4YIpIGoKq/763gjDFmMMvIyGDjxo2sWLGCgoICbrzxRlauXAnAmjVrmDdvHjNmzGD16tVs27at+bxrr70WgBkzZjBt2jRGjBiBz+dj3LhxHD7sTAtVVFTEggULAPj2t7/N+++/3+beH374Idu3b2fBggVccMEFrFq1ioMHD5KVlYXf7+eWW27hhRdeIC0trdP416xZw5NPPtk8HqsP/APwJyAkIjUiUisiNd2dJCJu4FHgK8BU4CYRmdrusO8Dlao6AafnRtNDBYGfAR2b8uB/ASdVdVL8uu8lOKbnxCK4Y8FevYUxxpjuJTsG623gS0Agvp4GvAl8oTeCMsYY43C73SxcuJCFCxcyY8YMVq1axZIlS1i2bBkbNmygqKiIe+65h2Cw5Q/rpi6FLpereblpPRKJANC+kaf9uqqyaNEinn766Q4xffzxx7zzzjs888wzPPLII6xevbrDMVu2bOGWW27h9ddfJy8v7/RfgFOgqpmneepFwF5V3Q8gIs8Ai4HtrY5ZjDPfFsBzwCMiIqpaB7wvIhMSXPd7wOR4bDGg7DTjS0re8fcYVnOgN29hjDEmCcm2YPlVtSm5Ir7c+ceWxhhjztiuXbvYs2dP8/rmzZsZM2ZMczKVn59PIBDgueeeO+VrHzp0iHXr1gHw9NNPc8kll7TZP3/+fD744AP27t0LQH19Pbt37yYQCFBdXc3VV1/NQw89xObNmxNe+7rrruP3v/89kyZNOuXYTpc4vi0iP4uvF4nIRUmcOgo43Gq9JL4t4TGqGgGqgU4zx1ZTmvxcRDaJyJ9EZFgnxy4VkQ0isqG0tDSJcBPzhmpP+1xjjDE9J9kWrDoRmdU09kpEZgMNvReWMcYMQEmUVe9JgUCAH/7wh1RVVeHxeJgwYQIrVqwgJyeHW2+9lRkzZlBcXMzcuXNP+dpTpkxh1apV3HbbbUycOJE77rijzf6CggJWrlzJTTfdRGNjIwDLly8nMzOTxYsXEwwGUVUefLBjnaN7772X8vJyli1bBoDH42HDhg2n8QqcsseAGE4Bpp/j9Lp4FOjuBUo0eKn9hFLJHNOaBygEPlDVfxCRfwAeAL7T4SKqK4AVAHPmzLGJrIwx5iyXbIL134E/icjR+PoI4MbeCckYYwzA7Nmz+etfE1cdX758OcuXL++w/d13321ebupa2H7fgQMHcLlc/OY3v+ny/Msvv5z169d3OObjjz/uMu4nnniCJ554ostjesk8VZ0lIp8AqGqliKQkcV4JUNRqvRA42skxJSLiAbKBii6uWQ7UAy/G1/+EM47LGGPMOS6pLoKquh6nH/kdwDJgiqpu7M3AjDHGmFMUjhesUAARKcBp0erOemCiiIyNJ2RLgJfaHfMScHN8+XpgtTbVxk8gvu9lYGF80xW0HdPVe6Jh2P0mBLut72GMMaYXJD3RME4Xi+L4OReKCKr6u16JyhhjTK8pLi5m69at/R1Gb3gYp8VoqIjch5MI/ay7k1Q1IiJ3Am8AbuApVd0mIvcCG1T1JeBJ4Pcishen5WpJ0/kicgDIAlJE5OvAlaq6HfjH+DkPAaXAd3vuUTt5FoDDH8GRjc5XH3drNcYYk2SCJSK/B8YDm4FofLMClmAZY4wZEFT1DyKyEae1SICvq+qOJM99DXit3bZ/arUcBG7o5NziTrYfBC5NKvieVHW447ZDH0HpDpj9930ejjHGDDbJtmDNAaZ21R3CGGOM6U8i8ntV/Q6wM8G2waOmpGU5FgOXC/bFS+mX74Mh46B0F+SNB7e3f2I0xphzWLJl2rcCw3szEGOMMeYMTWu9Eh+PNbufYuk/Hn/L8vEtUHmgZX3Lf0D5Xtj2IhzZ1OehGWPMYJBsC1Y+sF1EPgYamzaq6rW9EpUxxhiTJBH5CfA/gVQRqaGlpHqIePnzwUCaHrsx0LIxEoTN7SaLbowXv2iodApiWCuWMcb0qGQTrHt6MwhjjDkbPLb5sR693rILlnV7zH333ccf//hH3G43LpeLxx9/nHnz5vVoHK0tXLiQBx54gDlz5pz2Nf7yl7/ws5/9DJfLhcfj4aGHHuowkXFPUtX7gftF5H5VHeRVHdr15E/Us79p08ltcPQTmHotDJvW8ThjjDGnJakES1XfE5ExwERVfVtE0nAqLRljjOkl69at45VXXmHTpk34fD7KysoIhUL9HVa3rrjiCq699lpEhC1btvCtb32LnTt3dn/imftfIvJtYKyq/lxEioARqtr1xF3nsv3vdtxWsc/5Hon/LgWr+ywcY4wZDJIagyUitwLPAY/HN40C/txbQRljjIFjx46Rn5+Pz+cDID8/n5EjRwJw7733MnfuXKZPn87SpUtpqkG0cOFC7rrrLi699FKmTJnC+vXrue6665g4cSI//elPAWei4cmTJ3PzzTczc+ZMrr/+eurr6zvc/8033+Tiiy9m1qxZ3HDDDQQCTtezu+++m6lTpzJz5kx+/OMfdzgvIyMDEae7Wl1dXfNyH3gUuBj42/h6IL7NtFa+r+26xuDEdjj2Kax/AmpPOF0H27d+xWLOdmOMMV1KtsjFD4AFQA2Aqu4BhvZWUMYYY+DKK6/k8OHDTJo0iWXLlvHee+8177vzzjtZv349W7dupaGhgVdeeaV5X0pKCmvXruX2229n8eLFPProo2zdupWVK1dSXl4OwK5du1i6dClbtmwhKyuLxx5r2/2xrKyM5cuX8/bbb7Np0ybmzJnDr371KyoqKnjxxRfZtm0bW7ZsaU7a2nvxxReZPHkyX/3qV3nqqad64dVJaJ6q/gAIAqhqJZDSVzc/azUGYPtfYOdrECiFDU/B2gfg3V9AXXnLcVufd/YZY4zpUrIJVqOqNvdLEREPHTp6G2OM6UkZGRls3LiRFStWUFBQwI033sjKlSsBWLNmDfPmzWPGjBmsXr2abdu2NZ937bVO/aEZM2Ywbdo0RowYgc/nY9y4cRw+7MyRVFRUxIIFCwD49re/zfvvv9/m3h9++CHbt29nwYIFXHDBBaxatYqDBw+SlZWF3+/nlltu4YUXXiAtLS1h7N/4xjfYuXMnf/7zn/nZz7qd67enhOOVAxVARAqAWF/dvL/VZ449vRMjwc73bXm2Zbl8L9RXtN1fsd/phmizuBhjTLNki1y8JyJNFZoWAcuAl3svLGOMMQBut5uFCxeycOFCZsyYwapVq1iyZAnLli1jw4YNFBUVcc899xAMtvyR3NSl0OVyNS83rUciEYAO3fbar6sqixYt4umn21WgAz7++GPeeecdnnnmGR555BFWr17dafyXXnop+/bto6ysjPz8/FN/AU7Nw8CLwFARuQ+4HkjcxHYukk4+M80uhOqSxPsAYpHO90XjhYMbqlq2HVwH+RMhPR+2vwThBhg5C/xZzv76CkjNhb7rGmqMMQNKsi1YdwOlwGfAbTiz3Q+eNy1jjOkHu3btYs+ePc3rmzdvZsyYMc3JVH5+PoFAgOeee+6Ur33o0CHWrVsHwNNPP92hyt/8+fP54IMP2Lt3LwD19fXs3r2bQCBAdXU1V199NQ899BCbN2/ucO29e/c2jwnbtGkToVCIvLy8U47xVKnqH4D/AdwPHAO+rqp/6vUbDyAJ25E6S7wuucsp0V62J/F+gHAQTu6Ena+2bNv/Lhz8IL6/wfm+63WnFStQCh89Dgf/ehrRG2PMuSHZKoIx4N/jX8YYMyglU1a9JwUCAX74wx9SVVWFx+NhwoQJrFixgpycHG699VZmzJhBcXExc+fOPeVrT5kyhVWrVnHbbbcxceJE7rjjjjb7CwoKWLlyJTfddBONjU4rxvLly8nMzGTx4sUEg0FUlQcffLDDtZ9//nl+97vf4fV6SU1N5dlnn+3LQhcngP/CeX9LFZFZqjo4Z9TNHA61x8HVqujv8OlwfKuz7PXD2Eth7ztdX2f3f4Ivo+22aNgpetGkYj/seRPyJjjrXbWYGWPMOU40iX7TIvI5CT4YU9VxXZzzFHANcFJVp8e3DQGeBYqBA8C3VLVSnHfefwOuBuqBv296QxSRm2lpLVuuqqvi22cDK4FUnBa1H2kSDzNnzhzdsGFDt89sTGut5z9at6+cwpqNbfbfXxyfl+iyQT4Fzzlmx44dTJkypb/D6HEHDhzgmmuuYevWrf0aR6LXV0Q2quppTcIlIj8H/h7YR8t7lqrq5WcSZ186k/eoP7z2S3ae/JCff/1hqDnqdOPbtxryJ7W0QBVfAgfi4+2a/r/64N8g1LGKZBtpeVBf3vUxrWUXOl9DpziJnjHGnAOSfY9KdgxW6wv5gRuAId2csxJ4BPhdq213A++o6i9E5O74+j8CXwEmxr/mAb8G5sUTsn+O31+BjSLyUrwy1K+BpcCHOAnWVcDrST6PMcaYc8+3gPGtizINJoIQEw/kFDlfADOuh0ij06VPFcZ8oSXBajL3VqfQRfleOP4ZjJrtLLfuOngqyRU4LVjVJU6iN+0bznWzRzlJlzHGnOOSGoOlquWtvo6o6kNAl58IqupaoF25IRYDq+LLq4Cvt9r+O3V8COSIyAjgy8BbqloRT6reAq6K78tS1XXxVqvftbqWMcaYLhQXF/d761Uv2Qrk9HcQA47HB/PvgIuXOd0F8yc6XQObpKRB2hAougjmfh9GXuAkZj2hvhz2r3Fa0nbZZ6DGmMEhqRYsEZnVatWF06KUeRr3G6aqxwBU9ZiINM2lNQo43Oq4kvi2rraXJNjeWfxLcVq7GD169GmEbYwZrFS1L8cPDRrJdE8/DfcDn4jIVqCx1b2u7Y2bDURJvao9kTxd9hOnWuBHj3d9XKgOjm1pWY+GndasjGFOixZALOoU4rB/Z8aYc0SyXQT/tdVyhPj4qR6MI9H/qnoa2xNS1RXACnD6t59OgMaYwcfv91NeXk5eXp4lWT1IVSkvL8fv9/f0pVcBv8SpeDto5r/qNSMvhKOfOAUuGgMd96clGCnQ+tjMYTB0Kuxb46zXlTkTGINTxn3+7U7C9deHIWeMk/iV7nIKc0SCUPw3Tutae5EQeGz+aGPMwJVsFcHLeuh+J0RkRLz1agRwMr69BChqdVwhcDS+fWG77e/GtxcmON4YY3pMYWEhJSUllJaW9nco5xy/309hYY+PxylT1YdP50QRuQqn2JIbeEJVf9Fuvw+nO/psoBy4UVUPiEge8BwwF1ipqne2OuddYAQQr2XOlap6krPFpC/D+Mtgx8vQuAdm/71TLbCpUmBro+fDia1tJxye8CXIGtWSYLUm4rSANdY6CVPZHqcL4dFWZf9TMpzxYKE6mj9DrS+HrS84y63Hc4k4CVnumA63qgvXke5NB5zkfkvZFsZljyMzJfmOOKFIjOPVQUbnJZ5Y+1xnLfkmEVVl8+Eqpo/KxutOduannheMBIlq1Pl3Hos5/x80/b6qQrAKEPBn91lLebJdBP+hq/2q+qsk7/cScDPwi/j3v7TafqeIPINT5KI6noS9Afw/IpIbP+5K4CeqWiEitSIyH/gI+Dvg/08yBmOMSYrX62Xs2LH9HYZJ3kYRuR/nPaV1F8Euy7SLiBt4FFiE8wHe+nhBpe2tDvs+UKmqE0RkCU5L2Y1AEPgZMD3+1d5/U9Wzs3StiDN+a8rXoPYYZI1wvhIZt9AZ1xU4CRtXwvRvQk68S37eeCjf5yyPmu38sVO+r2P3wqbkaszFcOgj+Hyt89VOMBzlZG0jqZEyYr5sANLqj5KZOcJJsA6ug8MfwZzvsb26iic/fYZJmRczIm08DZEAH5x8HTTKVbECPvelE/YPZ4i3iJL6vSwYP4QJeUV8uidI474PEY1QqwHCsSi1vmFcNXkM43Oz2FcNpa4CxhWkk5Pm5UDtHvwePxkykuM1QaobwtQ01lITO8h5w7LJSRlOXTBCut+DW92khOvZ1XCMgtShjMooBISGSIxINEBZ6aeMG/1FYp+/Rl3NIQCGZYykISUNnzeDypwLSY9WUSt1NLg9EI1QuvcdKqoaCeTNIObJINObzpeGpnOk/jAnG6oprW0kw+8hy+8hHIHNgYP8zZivkJmSxZ6qXaiChGqoOXmQ9FQ3aR4vs3PPI9Xj4+MjR3j18GZ+Mvcq/J4UttUeoCAll4LUfI6kplOQMYJgJMjuyt2UVFVyXs50xp74EKkvI6gpRArn48oeSVrgMMFoI6WNFbjFjT+7GPFn4nP7SfWk8trnr1DaUMrXx19HWBvYWb6X0oZSXOLmspGX4dEYZQ1luMSNT7PIdjeQnQJojNLaRupCUVJLt+CtO0qG30uK2+Uk6f4sVldu53iomsuzplLhH42KG4kGSfV6yPB7qGuMkOH3IOLmePlOPOLFlzWPUcWTCEVDHKo5RAxt/vwg3ZNOti8HVfC4vBAT9u3+FFc0RKj8U47UbEXdKWRICtX+PFK9XnLSvYSiyoSia9DMAlxCc1IyxB0iZceLUHsChoyjrPBCDoRrKM4uJj+17QTt0ViUimBFc8LrEheRWIS9lXsZlzOOPP9QItEYm05uAIRsXzYj00eR6kklJVRDSdkWQrEwEu8MJm4PkjXK+Td3chtSNA8yhiORRnC5cfkyGaZu9jXk8NHRMKNyUhmZk0peqIT68qNsPO5j/0cHGR4+zMix04h501GXh0iwnL3hIwzJSME7ZDzFmUU0VlXi8Q4hmjGCaCxMIFyO2xsCFI+4OVwVRBuiTMlzUVl3HCJcmgD5AAAgAElEQVQh0j1+RqcNI6ZKWWOAcFYh26ur2FT6EXnpPk7WNBCtPUymz0VeYxUaqiNcKzRkzSKmMeZ46hmrQSK+GFsby1h08T/iyeu0CHqPOZUqgnNx3rQAvgaspe34qDZE5Gmc1qd8ESnBqQb4C+A/ROT7wCGcaoTgVAG8GtiLU6b9uwDxROrnwPr4cfeqalPhjDtoKdP+OlZB0BhjBrsL49/nt9qmdFOUCbgI2Kuq+wHiH/YtBlonWIuBe+LLzwGPiIioah3wvogkaNbpW732uazHB7nFifcVzXX+iBUBcTsJ2Bf/EVzOH47BcJQdOYvI3LGZvAwfuZOuhIZK6nIm8dd95eSVfkym1vL5iK9wtLKBJZmb8edPgoIpsPt16ssOs6sixuH8Bfgy6qjyHCD7WBonPYWUM745jZ5z7HeMrl9Lxs7/Ii/dB4Cue4TtJ6pwNR4iVnqQXPcIdsT2k6VB3NEgH8cfIejP50AEojHlzUovb6gSqStnrKQQ9EZoQKkNNkIDrPwQREEFgu4sYi5P/LV3Xv2JrjG48RAiTIVWU6HVfObzUhuKoriJurykhSpIidU1v4S56kfFRYMnm8bIMQCGbS8j3PhnUtP9SCyCHHd6vEZjMU7UZeIJV+D1QV32RFyxMFqxj5hC6MQGgp5sMsJV+MjnY+9R6kMRovHEwOsSwjFn5fCBRuq0lvpYLemRSvzhtnXJXgUy3V5qo2EANq+vZLg7jVeDu9scF/O2tOqVBUKsFS9DI0EySaGWELGdr6JADKVBIm3ObfDkICguhJjGEOD/3fIeAC6N4FHnB/yprsKrLsLivA6p6iGTFDJSPBCDYCRGBQ00/SsIpBQwMTuKO9JAOBqjLt6l9PFjr9EQjiDqXMcFpPs81DZGcOEkPB6X0BCJMjm6lrLUQkpix4gQj1tBWmaAaH6OYTE3KdF6gkSoE+f1CnqySInWo9X7m49zaYS39qxGcdHozqDRk4mgpMSCXOr1ENYYBw9upuqTp4lmFpLnyeG8lGJKIieIxaKICEfDJ3G5lFAkRjAcxed1IbEoqY1lZOBlcjifGMo6V+KOXW6BWKuGZrdLcIsQa8oej+7scE5UFXfMRbV/FJsPOtuyg0dafu54OYgLdnzW/KqUidNwn+J24XULwXCUqEK2+nABldL8GVjz/11N574jgtfd9n+0UDTWHHeNbySghGUoWeEKaiLHOSFh0lIyORqtozEaw1P5LpGY8nL8fI/L+Vca3fku1ywYOAlWPjBLVWsBROQe4E+qektnJ6jqTZ3suiLBsQr8oJPrPAU8lWD7BhJ/WmiMMWYQOoPu7IkKKs3r7BhVjYhINZAHlHVz7d+KSBR4Hmcuxw7jgM/qQkwTvkR5oJFX/nqAa2aOIC/D15xcAfz6XaflKif1cor8Mb4EkJrLMb+XrZE0fBlZ5ARLKG3II5bm4ZMxsxkVTuXN7ce5asRsdnwe4pAvh6ORVI6f/JxGCVDoPZ9sz0jG5qfzhQl5HK0KUnnUR0llFYoSJMp2Vxn5vpkENI2wK5XzxozC404hszadTNIBgbpSLvIPp1Ej1CAcOhlkeE0qmziCx6P4c7IZft7leNwpZNX4WbvtZeq99ZCez0WhMFl+L5FojFAkyoFIFZWRBsrYQTSmuERITXExJuRlajiPYGMdojFSvFG2+UOkShoFvgLq04fhDlbiitSjqhwPeQEoq36WqCrZMouq9MmMGuJCQhXI/jcZTgOQiicszC6rxYubPWQy2pvJghFD2d1YxculDbyjhyn3TuK8jDlkSAFZfi97GjcTrd9EdqwGiX1EDuB3uVmSdR6fhP28LcPJS72UmnApIa3Hk+HleOVOVJXfxIAYiPs8AObWHuKENOCPDEEQciWTCt1OqitCalYaR+ui+KJRpuWNIIqbmCeDcN500j2Z1JVvJrexlgy3n0AsRBQlhhKIhQiFhFDQRY6kM3vMODaGjlAarkbdKXhdPrziIdAYIlJ7gnKXhzLPUEKRKEV5fsakpnO8Po2q+hQONx5naGAXDd5cPhNnoux5LmG8+xgjc/x8nncpn5ULKkKaKsOyUgkHyqmOlFIX3MVnDSfIce9lqNdFtsvPAn8x4KKRCJWxIADV2kitBgHB5Uona8wF6NFNzMoaT+Pkq9l0chOFGYVkyDh2lx3h4NE3yJc6RjZUEw4qXreyL1TDiXA9G2Je6r1DyBYvQ6MeKuqOEuAoG9lOUwqSJimkRKL48TBec1G8ZOBBNMyn4qEiJYO1aUK6z0uaFDMncwpuXASiDRypDBIMRwl5s5g1+kJS3G4q64L4Knfi0ghZjcebE/9A5jhqsybh0hANlZ9A3QFEBI+7EbfLRVqKm6CvgEPhIRz25DN2lA9/WhYSOIG7bCexzBEMzRzNrvI0puowqiInadB6QqF9jHGVgbgYnjOBSNhFgzYyxJNNkXcY1cdW81HwOPVpo8jNnEVe6hROhp2aduXB3WQGD/NV8XOyppGUmJtc6lFJJeyawrFxNxBzecmsi+GNHmfCqADjPFC2+3U+bqii0Z+HJ1xHbt7k3vrfsI1kE6zRQOt5RUI4kwUbY4wxA4aIfBWYhjNnIwCqem93pyXY1j4ROqXiSnH/TVWPiEgmToL1HdrODdkUXw8WYurbOk6hSIzfrXM+0j5QXuckWAmkDp/EZ9VBotuO43EJ9aEoAF+dN5UtJaMIBRqprA/z4f6W+ba2BgvYkb2QnXVvMC73GDm12VSFXJQEN3HCtYOvj/5bhmb6GZrpJ3jNd3ji01VsLg0QdaVRkDsKV8Z5HCv3EdQUYkUX8Unlbkh1Jrb+ztTv4Dr2KenVR2H85ZCez8maIPWhKAv8brLT3Hhd3uZYVJXJ479A1FVNRkoGWSlZHZ6xLlxHJOa0dHhdXvweP4K0jF0KN0BdGVd742XxE4wFidWVcfTAu5QEXVQ2RImkXcjxgAuJgoriH5JHhifCVy++kBSvxykSAlzh8TtdL+srmJE/kZTAYWIow9KHMzJjZKvn+ALl9aVEAsdRVeqiDaR70ihIG8aV/iyu9KZ2iKk2VMvW0p24xINqDEGYmHsePnE7rZatniMQriUvLROPy0NFXYhILMbQzETFbK5J+HvSJBSJOS0rLuG8Lo+Mv24xxeVy4ohEYxwor6c+FCEUiZGT5uVbbqWs4STDfXnkxcpJTU1jnDefzAOVRGLK2Px0Jg1rGZMXjoaJ4TxrivsUC6pMv9kZ9+P2cN6QlugvHDUKzr/IWQmcdH52GUPjFTQ7VtFsiDRQ3uD8e3CJi2Fpw3C73ISjznhAgIJMH36vG6IRvhSupzLWSFVjFTGN4XV7Kc4qxiXOBx6qytHqIEMzfe3GS32xZTEadn5P/a1/vzv/WdUGwxwsr2fqiCzn9Vd1xlam5jhTQpyiSOwbTKk9RK4vl1x/bpt9qkpUo3hiMRAXKkJVQxSvx0WGryWdCYajnKwppjA3FRcKw6ZyWUoa0dRsfG4fHleyqc+ZSfYuvwc+FpEXcf73/gYJ3iSMMcaY/iIivwHSgMuAJ4DrobknWFc6K7SU6JgSEfEA2XSc67ENVT0S/14rIn/E6Yp41r93RmNKeV0jlXVhXK3+JszwOQmJqhJojDSneuOHZjAuP53GSIwDZXXNyRXAkPQUvnb+SF7d9xZbtpdR5J/Nl6cN541txzla5XQxqo9WkeEbQqa/nlFkc6y6gfpQLVXR/Wwt87KtbBu14VokN5/x/hz8KS58Hjfp3irEF2Ak2eyubNutLTMlE8Zc0mbb0KzOq1qKCMOy0oH0To9pKqTRKW9qywTQnXCl51M47Xo6K/8SjY3DJSQuOJE5DHA+WTg/Pb/jfpzz8tOHQvrQhPsTyUzJ5OJRc5M6Nt3X8kfxkPTTr/SY4jm1ggmuVr+IHreLCUMzOhxTTFPS4MSYCVwxZVjC63nd3oTbkwsmicQio/vXP9WTSmFmx98Er9tF0ZB2xVbcHjzuLAqAgrSChNcTEUbldEyg217H63wlKdPvZfqo7NY3gfS8pM9vz+PyMC47cfc9EcEjnuYZfAXITe/4Wvu97lbFaATyJ9Dj9WqTkGwVwftE5HXgb+Kbvquqn/ReWMYYY8wp+4KqzhSRLar6LyLyr8ALSZy3HpgoImOBI8AS4G/bHdNUpGkdTuK2OlF3vybxJCxHVctExIvzMfDbp/5I/UNVeWfHScbkpTEx/sn+52V1/PmTI52e89pnx3hv90nqGqNtthflpjJ9VHbzH2JHqxrYUlJNhs/T/MnzgZrd4AvwjQu/htdfzcThHirrXBRk+piemdUmoRiR7fyRuO7oujb3ERGy07ykuFMYnj6c+nA9deGWsU7D04czxD+kQ9GAs4nblagh1Rgz0JxKO1kaUKOqvxWRAhEZq6qf91ZgxhhjzCkKxr/Xi8hInHLq3ZaBjI+puhN4A6dM+1Oquk1E7gU2qOpLwJPA70VkL07L1ZKm80XkAJAFpIjI13Eq3h4E3ognV26c5Orfe+YxO3N6f3yrKnWhaJtuNi99epT9pXVsO1rDj4ZlEo0pH7XqvgcwZUQmhyrqaQjFmgfIt0+ufF4X00Zmt9k2Ml6FrE3kIkwYmonbV8lf9r5Mnj+P78y/EYDHNnf8RH1s9lim503n5f0vd9gXioa4ZpzTrenJz56kMdrIhJwJLBqzyEqNG2P6RLJl2v8Zp5LgecBvAS/wv4EFvReaMcYYc0peFpEc4P8DNuF0aU8qqVHV13Aq2rbe9k+tloO0VL5tf25xJ5edncy9e46e1gisPScDvLrlGIsvGMm4Aqdr1f5Sp+Unpoqqsq80wLHqYJvzrprulGyvDYZ54r+cz1svmzyUaCzGkaog+04GmDg0k+0VW8jwZjA+ZzwxjeHuogtVQ8TpFlgeLOexzY9xx/l3JDyuKLOIwsxCLi28lJJACY2RRo4EOrau5afmcyRwhMLMQkuujDF9JtkWrG/glL/dBKCqR+ODdo0xxph+JyIu4B1VrQKeF5FXAL+qVvdzaH3GJQJot5PCVtSFqKoPMTY/HRHhZI1TLvlwZQNvbT/RYVzK//7oEFNHdP6W33rAfIbPw4ShGRQNCVJVH2LSsAxePfxXAMZUjuFgzUHSvelcPvpyijI7jkd65+A7bWMNJh7mlpXidBucnj+d6flOQeGGSAO/3frbNsctLFrI6kOrGZPVcQJiY4zpLckmWCFVVRFRABHpZiSnMcYY03dUNRYfc3VxfL2RVpMNDwZulxBTp6Ka2915gvX8xhICjRFSPM5g+fx4MYKy2kbqQ1Fe/rRtfY+S6jLWn1xDjqeQbM8oXNK2BcrvdfO180cSDEcpjg8uz01z48n5K7U6s/m4popodeE61paspSiziDnD5rS5Vvs2uG3l2zrEn+JOaVMZr0mqJ5Vrxl2Dz91SyTDbl803Jn6j09fCGGN6Q7IJ1n+IyONAjojcCnyPXu9LbowxxpySN0Xkm8ALXRWgOFc1pVTVDWGGdFIuHSDQ6JQSD0Vi7DsZoGCcU/XrUEV9wuNPNO6gMnyYyvBhsjzDmZK5kPOLctocMybfx6v7X+WzPQ38zai/oTJYSVVjFWtL1jYfUx+pb25t2le1j61lW9latrXLZ2q//4rRVzAxd2Jz6en2RmedZfOIGWPOSclWEXxARBYBNTjjsP5JVd/q1ciMMcaYU/MPOHW0IyLizP7pzGXfcdKic1Caz2lZCjRGOiRYTfnm8Zpgh/O6K4kdJdy8HNJ67rx8YvN6dWM1uyt3MzJjJEcDTsvXltItjMhwxmdlpWRRE6oBwO/2U5hRyLiccfjdfjac2EC6N53Zw2a3ScQ6s7BoYZt5hYwxZqDqNsESETfwhqp+CbCkyhhjzICkqoN6bLC7VavOlpIqth+t4RuzRuHzuPnNe/sRgYZQtMN5TcO1bvviOB5/bz/gJGQnQtv52pQL2P5ZyzA2t6eOD458wPic8QC8ceAN6sJ1fLGwZbLSo3VHyfE7LVxNyRXA5aMvb25hGp01mg0nNjBr6Cym509vTrAWjVlEVKN8Xv05Xy7+MqX1pUQ1yhD/EPye/pjNxhhjTl23CZaqRkWkXkSyB9NgYWOMMWcfEckFJkLL3JKq2n3zyDnAEx8a9dzGkuYiF4+t2UeGz0Mw3DGxatLUuuVxtSRoEQ2SX3CYo5Ew44f62XasFhGYOjKLT0s/5dPST9tcY3vFdsCZlLY2VMvmk5s73EdalZEfnj6cmybfRI6vbVfDiblO69jkIZMBGJaeeCJYY4wZyJIdgxUEPhORt4DmWftU9f/qlaiMMcaYUyQitwA/AgqBzcB8nImBL+/PuPpKSwKjtJ4Tq2nMVWc+2OsUn2g9h+11s0fyZomHimAFw7PSOFrjZXiWv00S1lppfSngdON7/fPXicS6vidArj+322OMMeZslGyC9Wr8yxhjjBmofgTMBT5U1ctEZDLwL/0cU5851Xme6qMVBKJlfHnCHA7Ub2RftTK68BDl9bUMyVgIQExj1EcCTB6eeBib1+UlHGsZo5XiSmH2sNlsKd3SPKdVS4Cdx+L3+JmeN/2U4jfGmIGqywRLREar6iFVXdVXARljjDGnKaiqQRFBRHyqulNEBk1VBLcruQTL53URaKxjX8NactNh0qhZbN97kHcOHYxfCD45+UlS17q08FLeOdQyd1VMY8weNpvZw2bz/O7nSXGnkOfPY3Pp5jbl09v73vTvJXU/Y4w5G3TXgvVnYBaAiDyvqt/s/ZCMGWDW3A9VW5pXC2saujjYGNOPSkQkB+e96y0RqQSOdnPOOSPV6ybV2zJHVUyjCC4q3e8TDQ6nIMUZ3/SlKcP4w65/JxxoJMOfzprDazpca1fFrm7vd0HBBbjbzYnldrWsL56wGEGIapSCtAIKUgtO99GMMeas0l2C1frjsHG9GYgxxhhzJlS1aUbZe0RkDZAN/Gc/htSnohoj0+ehMRbALV4+C7zEnJHTGZHayLp9myhImci4gnSGZ/spyPCR4naRnealurFt/ao0Txr1kcRzYrV28ciL+bzmc8CZ5PeSUZcwNG1o836Py/kTw427uXiFMcYMBt0lWNrJsjHGGDMgiIgfuB2YAHwGPKmq7/VvVH3vSLACv9fNsMKD5KWM4bPPoCJ8gGGpPnLTvZyMfMb3pnyJLJ8Xl0vITU9JeJ1sX3aHBGt6/nSO1x1n7vC5ZKZk4hY3IsKI9BFMy5vGlLwpbZIrY4wZzLpLsM4XkRqclqzU+DIMsskbjTHGDGirgDDwX8BXgKk4BS8GFbe4EIF0n4uiXD8zRmWTHp98eOLQTGJ6jA0nNnDF6CuSup6I4HP7CEaCXDLqElzSsYJgqieVLxZ9McHZxhgzeHWZYKmqu6v9xhhjzAAwVVVnAIjIk8DH/RxPv3DFqwjGNEYkFiHD3/IW73YJboR9VftI86R1eZ0Lh17Isc+P8aXRX6IgtYC6SF3C5MoYY0xiyZZpN8YYYwaq5jrhqho51XLl5w7nuU/UnyDdm57wiEgsknAS4NaKs4u5deateF1eAHLI6fJ4Y4wxbdlHUsYYY85254tITfyrFpjZtNyqa/ugEoqFEm6/edrNXDjswjbbzhtyHlcVXwU446+A5uTKGGPMqbMEyxhjzFlNVd2qmhX/ylRVT6vlpMYKi8hVIrJLRPaKyN0J9vtE5Nn4/o9EpDi+PU9E1ohIQEQe6eTaL4nI1jN5xlNVUlvSYVtxVjHp3nTSPU7rVqonlVtn3srlRZdTkOaUUP/CyC/0ZZjGGHNOsi6CxhhjBjURcQOPAouAEmC9iLykqttbHfZ9oFJVJ4jIEuCXwI1AEPgZMD3+1f7a1wGBXn4E517d7M/yObnmtPxpjMwYSa4/t3lsVWZKJssuWNbLERpjzOBgLVjGGGMGu4uAvaq6X1VDwDPA4nbHLMapVgjwHHCFiIiq1qnq+ziJVhsikgH8A7C890Jv0d1cKifrTwLgEhd5qXlWuMIYY3qJtWAZ0wPW7S8H4MPI7uZtdy2a1F/hGGNOzSjgcKv1EmBeZ8fEC2lUA3lAWRfX/Tnwr0CXs/aKyFJgKcDo0aNPKfDWohrtcn931QONMcb0DPv4yhhjzGCXqHdd+wahZI5pOVjkAmCCqr7Y3c1VdYWqzlHVOQUFBd0d3im/K/HEwU08LvtM1Rhj+oIlWMYYYwa7EqCo1XohcLSzY0TEA2QDFV1c82JgtogcAN4HJonIuz0Ub0Lubrr8XTj0wi73G2OM6Rn2cZYxXXjwrd3MP1TOYVdDl8e95NoLQEnNs83bHtuc17xsg8eNGdDWAxNFZCxwBFgC/G27Y14CbgbWAdcDq1W10xYsVf018GuAeMXBV1R1YU8HnozRWaOZO2wueal53R9sjDHmjFmCZYwxZlCLj6m6E3gDcANPqeo2EbkX2KCqLwFPAr8Xkb04LVdLms6Pt1JlASki8nXgynYVCPvUxSMvZoh/CK/ufxWAa8Zd01+hGGPMoGQJljHGmEFPVV8DXmu37Z9aLQeBGzo5t7ibax8gQQn3nqbxIWHp3nTGZI3p7dsZY4zphI3BMsYYY85BbnH3dwjGGDMoWQuWMT2osGZjy8rnqS3LNgbLGNPHbp52c7el240xxvQ8S7DMgPfY5se63H+mBSS6uv6mmnJOxgtYGGPMQOZzeQFwxTun+D3+/gzHGGMGrbM+wRKRq4B/wxmY/ISq/qKfQzKnqLsE6kzPT5SAPfhWy4TAm2qcSYIvHj9AKmytuT/x9st+0rdxGGPOKhdmjiW15ihFWUXdH2yMMabXnNUJloi4gUeBRThzlKwXkZf6s3rT2ehMW4jOJEFat6+8w7ZEiU6yxyWSKL6mpKq7e/SYzpKmM72GJV3GmLgsTyoXpY4At6+/QzHGmEHtrE6wgIuAvaq6H0BEngEWA5ZgxZ1p69DpXqNXk5U+vMeA1xOJmyVpxpwbqg6Cxvo7CmOMGfSki3kSBzwRuR64SlVvia9/B5inqne2O24psDS+eh6w6wxumw+UncH5fc3i7T1nU6xg8fY2i7dnjFHVgv4Oor+ISClw8AwuMVB/rj3Bnu3sc64+F9izna3O9NmSeo8621uwJMG2Dhmjqq4AVvTIDUU2qOqcnrhWX7B4e8/ZFCtYvL3N4jU94UyTy3P552rPdvY5V58L7NnOVn31bGf7PFglQOvRvIXA0X6KxRhjjDHGGDPIne0J1npgooiMFZEUYAnwUj/HZIwxxhhjjBmkzuougqoaEZE7gTdwyrQ/parbevm2PdLVsA9ZvL3nbIoVLN7eZvGageBc/rnas519ztXnAnu2s1WfPNtZXeTCGGOMMcYYYwaSs72LoDHGGGOMMcYMGJZgGWOMMcYYY0wPGXQJlohcJSK7RGSviNydYL9PRJ6N7/9IRIrj2/NEZI2IBETkkXbn/KeIfCoi20TkNyLijm+/R0SOiMjm+NfVAyHeVue+JCJbW60PEZG3RGRP/HvuAI93QL6+IvJu/JpNcQ3t6loDNNa/F5HSVttvOZVYezHeFBFZISK7RWSniHyzq2sN4HgH3OsrIpmt4tksImUi8lBX1zIDS3e/EwORiBwQkc/iv3Mb4tsSvheJ4+H4820RkVmtrnNz/Pg9InJzPz3LUyJyUpJ4Xz2dZxGR2fHXam/83ERT1fTls3X6HiwiP4nHuUtEvtxqe8LfUXGKlX0Uf+ZnxSlc1hfPVRT//3CHOH/D/Si+/az/uXXxbOfCz80vIh9Ly9/e/9JVPNLFe9ipPnPSVHXQfOEUwtgHjANSgE+Bqe2OWQb8Jr68BHg2vpwOXALcDjzS7pys+HcBngeWxNfvAX480OKN778O+D/s3XmYHGd17/Hv6WX2RZoZSZY1kkaWbcmbbMvyvuAFvAE2BAcImwMEhy0hcHODDQkhEC5OCOGGC8ZAANvEgPEGYrWFsTEGrxJarV2WpdG+zr52n/tH1Ujtcc+oR+qe7p75fZ6nn6l6u6r6nKqarn6r3nrrB8DKlLJ/B24Nh28F/q3A4y3I9Qs8ASxI83lpl1Wgsf5luu1QAPH+C/Cv4XAEaDjWdZuneAty/Q6afzFwWTbWr165f2WyTxTiC9g88H+RUpb2WARcD/yK4Fh7AfBsWF4HbAr/TgyHJ+Yhl8uA+WRwXD2aXIDngAvDeX4FXJfn3D5LmmMwcGq4/5UCs8L9MjrcPgr8mMO/ne4EPjRKeU0F5ofD1cC6MP6i327D5DYWtpsBVeFwHHg23B5p42Ho4+GIc870Nd6uYJ0HbHD3Te7eC/wIuHHQNDcCd4fDDwBXmZm5e4e7PwV0D16ou7eGgzGCDZGtnkNyEq+ZVQGfAP51mGXdDbypwOM9VjmJdxhpl1WgsR6rXMX7PuCLAO6edPeBp7Efy7rNR7zHKqf7g5mdBEwGfj/csrKTimRJJvtEsRjqWHQjcI8HngEmmNlU4Bpgkbvvd/cDwCLg2tEO2t2fBPYPKs5KLuF7Ne7+tAe/DO9h5MfoozZEbkO5EfiRu/e4+0vABoL9M+0+Gn6XXEnw3QJH9/vjqLj7DndfEg63AauBaYyB7TZMbkMppu3m7t4ejsbDlw8Tz1DHsBHlPJIYx1sFaxqwNWW8mVfvbIemcfd+oAWoP9KCzewRYDfQxuGNC/DR8DLyd23kTe5yFe/ngS8DnYPKp7j7jnBZOwh+YBVyvFCY6xfge+Gl939K+SF6tMvKR6wAbwnX7QNmNn3IuUcpXjObEA5+3syWmNn9ZjblaJZVAPFCga3fQf6C4AzfwMmiY12/knuZ7BOFyIFHzWyxmd0Slg11LBoqx0LOPVu5TAuHB5fnW7pj8EhzqwcOht8tqeWjKmw2djbB1ZAxtd0G5QZjYLuZWdTMlhL89l5EcMVpqHiGOobl7DtlvFWw0p1xHXy1KZNpXj2B+zUEl2NLCWrQAN8AZgNnATsIKgkjkfV4zews4ER3f3iEsWRitOMtuPUbeqe7nwFcGr+L/Q4AACAASURBVL7efQzLGkkc2Yz1Z0CTu88DfsPhMz+ZykW8MaAR+IO7zweeBv7jKJc12GjHW4jrN9XbgR9maVkyOop1G10c/n9cB3zEzC4bZtqhcizG3EeaSyHmONQxuOhyC1vKPAj8XUqrpLSTpikrttzGxHZz94S7n0VwnD0POGWYeEY9t/FWwWoGUs8UNwLbh5rGzGJALRleFnf3bmAh4WVEd98V7gBJ4NsEO0C+470QOMfMNgNPASeb2RPhe7vCy9mEf3cXcrwFun5x923h3zaC+8YG4jrqfWu0Y3X3fe7eE072beCcDOPMZbz7CK5iDlS27ye4J+BolpXXeAt0/RJOeyYQc/fFx7osGVWZ7BMFx923h393E/yvnMfQx6Khcizk3LOVS3M4PLg8b4Y5Bo80t70ETe1ig8pHhZnFCSog97r7Q2HxmNhu6XIbK9ttgLsfJLif/IJh4hnqGJaz75TxVsF6HjjJgl5GSgjO0i4cNM1C4OZw+CbgtynNZF7FzKpS/gljBDdArgnHp6ZM+mZg5auXMLrxuvs33P14d28iuNF9nbtfnmZZNwM/LeR4C3H9mlnMzBrC4TjwhpS4RrSsfMY6aN3eQNB2eyRysS84wZWfy8Oiq4AXj2ZZ+Y63ENdvir/glVevjmVZMnoy2ScKiplVmln1wDBwNcF30FDHooXAeyxwAdASNt96BLjazCaGzZ2uDssKQVZyCd9rM7MLwqbc72Hkx+isGuYYvBB4uwU9t80CTiLo6CHtPhp+lzxO8N0CR/f742hzMOA7wGp3/8+Ut4p+uw2V2xjZbpMsbIZvZuXAawmOo0PFM9QxbEQ5jyhIH4XePgrpRVABWkfQVvPTYdnngBvC4TKCM80bwpV8Qsq8mwlqvO0EtdtTgSnhhlgOrAL+H8HZX4DvAyvC9xYCU/Md76BlN/HKHoHqgceA9eHfugKPt+DWL0EPbYtT9of/AqJHWlYBxvrFsGwZwRfW3Hyv27B8JvBkGPNjwIxsrNs8xFuQ6zd8b9PgeLKxfvXK/SvdPlHIL4IeupaFr1Up+3HaYxFBs52vh/mtIKUHVIIOZTaEr/fmKZ8fEjS56gv/r96fzVyABQQ/hjcCXwMsz7kNeQwGPh3GuZaUXvOG2kfDfeG5MOf7gdJRyusSgqZfy4Gl4ev6sbDdhsltLGy3ecCfwhxWAp8ZLh6GPx6OKOdMXxYuRERERERERI7ReGsiKCIiIiIikjOqYImIiIiIiGSJKlgiIiIiIiJZogqWiIiIiIhIlqiCJSIiIiIikiWqYIlkgZklzGypma00s/vNrCLfMQGY2aeysIw/N7NVZpY0swXZiEtERAqPmbWHf5vM7B1ZXvanBo3/MZvLFykkqmCJZEeXu5/l7qcDvcAHM53RzKK5C4sRV7DSxLMS+DOC5zmJiMjY1wSMqIKVwbHsFccjd79ohDGJFA1VsESy7/fAiQBm9hMzWxxeAbplYAIzazezz5nZs8CFZvYZM3s+vAL2rfAJ7JjZE2b2FTN70sxWm9m5ZvaQma03s39NWd67zOy58CraN80sama3A+Vh2b1DTZcuntRk3H21u6/N9UoTEZGCcTtwaXis+Hh4TPlSeJxabmZ/DWBml5vZ42b2A4KH16Y97g1xPBq4Wmbhslea2Qoze1vKsp8wswfMbI2Z3TtwbBQpdLF8ByAylphZDLgO+HVY9D53329m5cDzZvagu+8DKoGV7v6ZcL4X3f1z4fD3gTcAPwuX0evul5nZx4CfAucA+4GNZvYVYDLwNuBid+8zszuAd7r7rWb2UXc/K1zuKemmA+4ZHI+IiIxrtwJ/7+5vAAgrSi3ufq6ZlQJ/MLNHw2nPA05395fC8XTHvVccjwb5M+As4EygIZxnoMXE2cBpwHbgD8DFwFPZT1cku1TBEsmOcjNbGg7/HvhOOPy3ZvbmcHg6cBKwD0gAD6bMf4WZ/QNQAdQBqzhcwVoY/l0BrHL3HQBmtilc5iUEla7nw5N75cDuNDFeNcx0g+MREREZcDUwz8xuCsdrCY5nvcBzKZUrGPq4N5RLgB+6ewLYZWa/A84FWsNlNwOEx9gmVMGSIqAKlkh2dA0+M2dmlwOvBS50904zewIoC9/uDg8mmFkZcAewwN23mtlnU6YD6An/JlOGB8ZjgAF3u/ttR4hxuOkOxSMiIjKIAX/j7o+8ojA4znUMGh/quDfcsoeSesxLoN+tUiR0D5ZI7tQCB8KDzFzggiGmGzj47DWzKuCmIaYbymPATWY2GcDM6sxsZvhen5nFM5hORERkQBtQnTL+CPChgeOJmZ1sZpVp5hvuuJd6PEr1JPC28D6vScBlwHNZyUIkT1TBEsmdXwMxM1sOfB54Jt1E7n4Q+DZBE8CfAM+P5EPc/UXgH4FHw89aBEwN3/4WsNzM7j3CdEMyszebWTNB5xe/MLNHjjSPiIgUteVAv5ktM7OPA/8NvAgsMbOVwDdJfzVpuOPeoePRoHkeDj9vGfBb4B/cfWdWsxEZZebu+Y5BRERERERkTNAVLBERERERkSxRBUtERERERCRLVMESERERERHJElWwREREREREskQVLBERERERkSxRBUtERERERCRLVMESERERERHJElWwREREREREskQVLBERERERkSxRBUtERERERCRLVMESERERERHJElWwREREREREskQVLBERERERkSxRBUtERERERCRLVMESERERERHJkli+AxhtDQ0N3tTUlO8wREQkjcWLF+9190n5jiNfjuUY1dHXQWdf56Hx6pJqymJlWYpMREQyPUaNuwpWU1MTL7zwQr7DEBGRNMzs5XzHkE/ZOEZtad3Czzf9HIAbT7yRaVXTshGaiMi4l+kxSk0ERURExpAZNTN44+w3ArB41+I8RyMiMv7krIJlZtPN7HEzW21mq8zsY2F5nZktMrP14d+JYbmZ2VfNbIOZLTez+SnLujmcfr2Z3ZxSfo6ZrQjn+aqZWa7yERERKRbTq6czvXo6zW3Nr2g2KCIiuZfLJoL9wP9y9yVmVg0sNrNFwF8Cj7n77WZ2K3Ar8EngOuCk8HU+8A3gfDOrA/4ZWAB4uJyF7n4gnOYW4Bngl8C1wK9ymJOMY3csvSNt+YfP+vAoRyIicmQnTjiRrW1beXbHs1wx44p8hyMiMm7k7AqWu+9w9yXhcBuwGpgG3AjcHU52N/CmcPhG4B4PPANMMLOpwDXAInffH1aqFgHXhu/VuPvT7u7APSnLEhERGddOqT+F6pJqVu9fzZbWLfkOR0Rk3BiVe7DMrAk4G3gWmOLuOyCohAGTw8mmAVtTZmsOy4Yrb05TLiIiIsA1TdcAsGzPsjxHIiIyfuS8gmVmVcCDwN+5e+twk6Yp86MoTxfDLWb2gpm9sGfPniOFLCIiMiZMrphMSbSEjr6OfIciIjJu5LSCZWZxgsrVve7+UFi8K2zeR/h3d1jeDExPmb0R2H6E8sY05a/i7t9y9wXuvmDSpHH7eBURERmHpldPZ3/3fnoTvfkORURkXMhlL4IGfAdY7e7/mfLWQmCgJ8CbgZ+mlL8n7E3wAqAlbEL4CHC1mU0Mexy8GngkfK/NzC4IP+s9KcsSERERYGbNTADuW3sfwS3LIiKSS7nsRfBi4N3ACjNbGpZ9Crgd+LGZvR/YAvx5+N4vgeuBDUAn8F4Ad99vZp8Hng+n+5y77w+HPwTcBZQT9B6oHgRFRERSzK2byws7X6C1t5Ufrf0RlfFKJpVP4sLjL8x3aCIiY1LOKlju/hTp75MCuCrN9A58ZIhlfRf4bpryF4DTjyFMERGRMe+tc97KY1seo6u/i+a25kPPx7pq5qsOxyIicoxyeQVLRERECkBJtITrZl0HQGdfJ3etuou1B9YSi8S4rPEygpb2IiKSDaPSTbuIiIgUhop4Be8+9d0ArNq3irtW3cXWtq0kkok8RyYiMjaogiUiIjLOVJdU865T30VtaS1d/V38bOPPeHrH0/kOS0RkTFAFS0REJENmttnMVpjZUjN7ISz7rJltC8uWmtn1KdPfZmYbzGytmV2Tv8hfraakhr+Y+xfcdPJNACzfs1xduYuIZIEqWCIiIiNzhbuf5e4LUsq+Epad5e6/BDCzU4G3A6cB1wJ3mFk0D/EOKWIRJldM5sxJZwLw2JbH8hyRiEjxUwVLREQkN24EfuTuPe7+EsFjSM7L1Yftau3myXV72LC7fcTzXnT8RQC81PIS96y6hyebnyTpyWyHKCIyLqiCJSIi44KZ1ZjZF83s+2b2jkHv3ZHhYhx41MwWm9ktKeUfNbPlZvZdM5sYlk0DtqZM0xyW5cTy5hYWv3yAny3bziOrdo5oXjPjppNv4pS6U2jva2fl3pXcuexO1h9Yn6NoRUTGLlWwRERkvPgewfMZHwTebmYPmllp+N4FGS7jYnefD1wHfMTMLgO+AcwGzgJ2AF8Op03X97kPLjCzW8zsBTN7Yc+ePZlnM8gpU6s5c3otAC9ub+XxNbsJHjGZmckVk7lixhV88MwPcsHUYHUsenkRuzt3H3VMIiLjkSpYIiIyXsx291vd/SfufgOwBPitmdVnugB33x7+3Q08DJzn7rvcPeHuSeDbHG4G2AxMT5m9EdieZpnfcvcF7r5g0qRJR5cZ0DixgivnTuHmi5oAWLr1IKu2t454ORGLMH/KfM6efDYAD6x7gFV7Vx11XCIi440qWCIiMl6Umtmh4567fwH4FvAkcMRKlplVmln1wDBwNbDSzKamTPZmYGU4vJDgSlmpmc0CTgKey0omw6irLOGWy04AYNGLu3h87W66+0b+jKsLj7+QG2bfAMCWti1ZjVFEZCzLqIJlZqfnOhAREZEc+xlwZWqBu98N/C8gk/7JpwBPmdkygorSL9z918C/h123LweuAD4eLnsV8GPgReDXwEfcfVSe5ltZGuPiExsoiUVYuuUg33hiI89v3j/i5TRWNzK5YjIvtbzEwe6DOYhURGTsiWU43Z1mVgLcBfzA3fUtKyIiRcXd/2GI8l8TXF060vybgDPTlL97mHm+AHxhBGFmzXmz6jhvVh1rdrbyqxU7WdHcwrlNdSNeTlNNE7s7d/ODNT/guMrjaChvYF7DPCaUTchB1CIixS+jK1jufgnwToK25C+Y2Q/M7HU5jUxERESO2dzjapjVUElLVx/ff3ozW/d3jmj+BcctYMGU4JFfOzt2snLvSn6w5gfsaN+Rg2hFRIpfxvdguft64B+BTwKvAb5qZmvM7M9yFZyIiIgcu9ecPIlJ1aXsbe/lgcXNdPWOrKXieVPP4wPzPsAHzvgAV0y/AoCHNzzMqn2r6OjryEXIIiJFK9N7sOaZ2VeA1QTt19/o7qeEw1/JYXwiIiJZY2YRM7so33GMtomVJbzrgplcdnIDAA//aduIunAHiEfixKNxTqk/hfOnng/A77b+jrtX3U1LT0vWYxYRKVaZXsH6GkF3tme6+0fcfQkc6q72H3MVnIiISDaFXal/+YgTjlFnTw+egbyrtZvVO9qOejnnTDmHd8x9B2c0nAHAvavv5Zkdz4y40iYiMhZlWsG6nqBziy44dAawAsDdv5+r4ERERHLgUTN7i5mlexDwmBaJGK+fF/Qq/7t1R/9QY4AJZRO4tPFSLmu8DIAlu5bwnZXfYVv7tmOOU0SkmGXai+BvgNcC7eF4BfAoMO6aWYgA8NLvDw8faIErbstfLCIyUp8AKoGEmXUBBri71+Q3rNFx8pRqnq7cx/6OXrYd7GLahPJjWt7pDaczp24OizYvYnPrZn664afUlNQwuWIyVzddnaWoRUSKR6YVrDJ3H6hc4e7tA1ewREREiom7V+c7hny77ozjuPeZLfz4+a1MrS1jR0s3J0yqpDQW4XWnHkc0MrKLe/FInOtPuJ5t7dtYt38d6w6sY8PBDZQ1l3Fq/ak0lDfkKBMRkcKTaQWrw8zmD9x7ZWbnAF25C0tERCQ3wqaB7wRmufvnzWw6MNXdn8tzaKNmcnUZl57UwMv7OtnT3gPApj1Bb4DlJTFec/Kko1rutKppTKuaxoLjFvD9F7/Pyr0rWbl3JQCvmf4aGsoamFQxiYhl3ImxiEjRybSC9XfA/Wa2PRyfCrwtNyGJFK6nN+4DoLG1i+nH2KxGRPLmDiBJ0BPu5wmav38dODefQY22BU11LGiqw91p6+mnqiTGfz22niUvH2BWfSUz6o++oUp1STXvmPsO9nfv59ebfw0EPQ4OuHrm1Zw48cRjzkFEpBBlVMFy9+fNbC4wh6Ct+hp378tpZCIiIrlxvrvPN7M/Abj7ATMryXdQ+WJm1JTFAfiz+dN4aMk2HlzSzF9dOovqsPxoTCibwISyCXz4rA/T3ttOd6Kbbe3b+MO2P/BE8xMs2b2E1818HZXxSgwjHj36zxIRKSSZXsGC4MxeUzjP2WaGu9+Tk6hERERyp8/MooADmNkkgita497M+krOmFbLim0tLH75AJfPmZyV5VaVVFFFFQ3lDcQsxroD69jRsYMfrvnhoWkqYhVcNO0iZtfOJhqJZuVzRUTyIaMKlpl9H5gNLAUGHv/ugCpYIiJSbL4KPAxMNrMvADcB/5TfkArHa0+dwrrdbWzZ35mT5Z/WcBqn1p/KugPr6OoPbuf+4/Y/0tnfyW9e/g2P8RhvnfNW6svrc/L5IiK5lukVrAXAqa4nCIqISJFz93vNbDFwFUGz9ze5++o8h1VQKkti7GvvZcu+TuqqSqgqHUmDlyMzM+bUzTk0XlNSQ1+yj2V7lrG3ay/3rb2PuXVzuXLGlVn9XBGR0ZBpNz4rgeNyGYiIiMhoMLPvu/sad/+6u3/N3VeHLTUkdPaMCQA8uKSZbz+5iZau3N52fcKEE5hTN4e3znkr1826DoA1+9eQSCaOMKeISOHJtILVALxoZo+Y2cKBVy4DExERyZHTUkfC+7HOyVMsBemMabVcc9pxnDi5CoBnN+1jtBqxzKqdxSXTLgHggXUPjNrniohkS6bX/D+byyBERERyzcxuAz4FlJtZK0HzQIBe4Ft5C6wAmRmnHl/DCZMq2bC7nVXbW5laW84ZjbWj8vmnN5zOszueZV/3Pr65/JucWn8qlzVeNiqfLSJyrDK6guXuvwM2A/Fw+HlgSQ7jEhERySp3/6K7VwNfcvcad68OX/Xuflu+4ytEZfEoN1/UBMBvVu/i5X0do/K5EYvwzlPeyZmTziTpSVbuXcnOjp2j8tkiIscqowqWmX0AeAD4Zlg0DfhJroISERHJoU+b2bvM7J8AzGy6mZ2X76AKVV1lCTeedTwAK7a1jNrnVsQruHjaxVw/63oAHlr/EJ19uenZUEQkmzK9B+sjwMVAK4C7rwey83AMERGR0fV14ELgHeF4e1gmQzhhUhVl8Sjrd7WzblfbqH52U20T8ybNA+CuVXexvX07Sddjy0SkcGVawepx996BETOLET6gUUREpMic7+4fAboB3P0AUJLfkArfxScGz6X6xfIdrBzFK1kAl0y7hBk1MwD4yYafHLqa1ZfIbe+GIiJHI9NOLn5nZgM3Br8O+DDws9yFJSIikjN9Yc+BDmBmkwBdEjmCeY0TOGFSFd9+chOLXtzF42t2c/aMicyoq2B6XTlmduSFHINrZl7Djo4d/HzTz9nduZu7Vt0FwOwJs7mm6ZqcfraIyEhkegXrVmAPsAL4a+CXwD/mKigREZEc+irwMDDZzL4APAX8n/yGVByqSmO8+expnDSliv6k8/zm/Ty4pJlfrtiZ8+7U49E4M2pmcMu8W7hqxlVcOu1SADYe3Mi3l3+bh9c/zJr9a3Iag4hIJjK6guXuSeDb4SsjZvZd4A3Abnc/PSyrA+4Dmgh6JXyrux+w4LTXfwHXA53AX7r7knCemzlcmftXd787LD8HuAsoJ6jwfcz1sAwRETkCd7/XzBYDVxF01f4md1+d57CKRlNDJU0NlfQnkuxo6eaBxc2s29XG5n0dlMYiJN25aHYDJ0+ppiSW6XnczMUiMebUzQGCq1dLdy+lpbeFl1peYkfHDpbvWc5NJ99ExLL/2SIimciogmVmL5Hmnit3P2GY2e4Cvgbck1J2K/CYu99uZreG458ErgNOCl/nA98Azg8rZP8MLAg/f7GZLQzby38DuAV4hqCCdS3wq0zyERGRcW8X8HuC42C5mc0fOLEnmYlFI0yvq+BDl89m0Yu76O1P0p9Msv1gN4te3MWiF3exoGkiZ8+YSFVppnckjExFvIKLpl0EwK6OXTy4/kH2du3lzmV3UhItwTDqy+u5cfaNOW/CKCIyINNvvAUpw2XAnwN1w83g7k+aWdOg4huBy8Phu4EnCCpYNwL3hFegnjGzCWY2NZx2kbvvBzCzRcC1ZvYEUOPuT4fl9wBvQhUsERE5AjP7PPCXwEYOnzx04Mp8xVTMyuJR3njm8YfGO3v7eXbTftbsbOOFzQd4cXsrrz11Cic0VOa0kjOlcgrvPvXdLN29lIM9B6mKV7HuwDq2t2/nia1PcPG0iymJqi8TEcm9TJsI7htU9H/N7CngMyP8vCnuviNc5g4zG+jqfRqwNWW65rBsuPLmNOUiIiJH8lZgdmrvuJI9FSUxrpg7mSvmTmbhsu1s3N3OwqXbqSqN0TixnCm1ZdSUxZk2oZzykmhWP7u6pJpLGy89ND5v0jzuW3sfq/evZvX+1fz5yX/OpIpJWf1MEZHBMm0iOD9lNEJwRas6i3GkO6XlR1GefuFmtxA0J2TGjBlHE58IvPR7Glu78h2FiBy7lcAEYHe+Axnr3jhvKnvae/j9ur1s2d/Jmp1trNl5+Dla9VUlxCIRmhoqOL62nJn1FVm9ylVfXs8t827hly/9kua2Zu5fdz+XTruUMyadkbXPEBEZLNMmgl9OGe4n7KDiKD5vl5lNDa9eTeXwwa0ZmJ4yXSOwPSy/fFD5E2F5Y5rp03L3bwHfAliwYIE6wpCs2HowqGw9vX8fz/Sv4+OvOznPEYlIhr4I/MnMVgI9A4XufkP+QhqbzIzJ1WW85ZxGkkmnN5Gko6efzfs6aD7QRVdvgt1tPexq7T40TzRixKLGacfX0tufZHpdOXOPqznqGGKRGDfMvoG1+9fy2JbH+P2233Naw2nqBENEcibTJoJXZOnzFgI3A7eHf3+aUv5RM/sRQScXLWEl7BHg/5jZxHC6q4Hb3H2/mbWZ2QXAs8B7gP+XpRhFRGRsuxv4N4JHj4zo+VdmthloAxJAv7svOJoecsejSMQoi0Qpi0epryrlnJmH3+vuS7B2Zxtt3f0k3Vn88gGWvHwAgJXbWnhi7R5ikeDK1jkzJxKNGKdMrSEeDSpJnb39JD1o3lIejxKJvPoq2Jy6Oezr2sfSPUt5ZvszhzrHEBHJtkybCH5iuPfd/T/TzPNDgqtPDWbWTNAb4O3Aj83s/cAWgs4yIOgF8HpgA8FB6L3hcveHNyM/H073uYEOL4APcbib9l+hDi5ERCQze939q8cw/xXuvjdlfEQ95B7D545ZZfEoZ06fcGj83KY62nv62dnSzZ724OrW3vZeth3o4om1ewB4bPVu6qtK6OxN0NWbeMXy6iqDzizqq0p47SlTKIsH93rNnzKfpXuWsnTPUja3bmZu3Vxm1sykvrwegKQH9e18X93q6U/Q3t2PAxUlUSpKctMLo4jkxkh6ETyX4EoTwBuBJ3llBxSv4O5/McRbV6WZ1oGPDLGc7wLfTVP+AnD6sFGLiIi82mIz+yLBMS21ieDRXl0aUQ+5A509ydDKS6KUl0SZVF0K1B4q7+jpZ+uBTp7euC98D+oqg/dm1FWwv6OXjp4E8d6DvLxzDzv2w7M9u3jNyUHHFmXu3DDhNJ7avZj9+xbzzPbFPFNSQXlJDRGL0JnowYEpFZM4Y/JZ7En2YmZELELEoaRtJ7VeyjSLUIIDTl8iqJQZTiwSoaO3j86eBHWVcaJm9CX78b7OIMhkP5GDW4lZBE57E9Qc7n1x054OtmxYzL7di2jr6mcbbfRGynCLcVL8eGqtnM39+ymNwNSKBBESTI1UEnfo8h6Or25kenU5eDJ4AUyYCWa4O119CdxhVdtGtrRsYmJNcE/6CeVTmFIygb6qqSQmnUhdZQ1723vYHjaD70/2EY/1MrEqysGeFiJmzKieScmmx+HgVvrNiViUHpJs7W2hqqSMhtKJROe8nnhpFVGLsmF3O+1dwfIcSCZ66e3ej+N0JTpJegIwqmJVVEQrcYx+T9BqUBKvoKzf6e7Zy46ebbjDcZVTmdG/FevvpjJWRW+/cVxtJdFYnB0l0ymtPY7qshgtXX3sbOmmL9nL9s6XMWDgFr+IQWVpnPrKwz1L9iWSJB3KbAI7D8Cq1ido6z9ISaSEeOTwdNFIlOMrUu9sMeJRO7RPRizCSRNOojRWyu7O3azc9TJtPd2097USMcOIEDGjvjJOWSxOY9XxzKk7hb7uFjZveYGd+3cxMVqFhV0OGFAWKSGJk/AEbcl2/tS+jp7+4MSCeT/9iV68bCKdyR5i0QgRi3D1pMuIWZyDfa0kHQ509BKtnUZ5ZSUlkSh7O3qJJMooi5RwXLyLaRMriba8DJ7AejuItLxMpydo9QQ++RSqGqYHK3DgBETXAVq6eimP1BApq2N7335qW7cy2UrwiU10RyqwWCnEKzGDiBkV1kv0wCbwJGZBk2CvqIeyWhLubO/dR1W8mpJoCV19XbywYw04HOw5wO7WDVwSq6PGuim3OBOo4eUDvUyffw3ltQ0AtPf0s7+jl8a6KFUlwT2dXX09dPX1UhotpTQWIxbN/QmUTCtYDcB8d28DMLPPAve7+1/lKjAREZEcOTv8e0FKWabdtDvwqJk58M3wHt+R9pCrCtZRqiyNMfe4GuZOqYZk/+E3dq2Ctp1wcEtQueg6QH95kiVbDpLY5axYG6M8HmFyTRnVHlxa7PZK9iS62NK6h9bEdva19VARibPZWmiJR9i4rxL3RAAAHH9JREFU+lE8VkbUYiSSCdo720gkHQdqvZTKaCVGhKD+DDFiTLQy+h166Ocg3bRHuummj343nCjxKMQTQUUjuuUFGqz8UK9dvf1J9tBJNGp0l9VTEYkxMea09nXwEutxoKuvHXNnU0+U7mgN2AFwpzTZjh/cTHmsFieCJXuJJTuJE6PEoyTDGHvop4OO4AN3biUZ1jamJOPssTAuN6JEmEQl5cRZb0FTzSBQC3sacyLJJG7QHasF9/AdJ5bspTaeoP1P9wNQRw1RnIgn6KCPPpL0ksCH6MtksldQToyXrRUAtyiR5NDTO5FwHSYp8xjgtJdMBoKYEhajK5y3rL+VaNh5qAfVLSKRCGAkLE4fUXpi1cQS3cSSPRhO1CLUU0t3OFcH3XR4F+tSYrBB/awlidIbq8SJkIjED5XHE13UUEppfwsl/a1sAXqi/bjDJMrZTSfJEfQUUOZRJlBGH0kO0k0ZMUqIUEGcbdbOD19Of86oNxqclTB34slO4j50haM9VkFpoh02/pZYxIhFjETSSXpwQqEncbiVdcJKDq3fAXVeRjJaSYwYXfQSTXQces+BNgumr/I47dZH0mIkLE468WQXC4GkRUlYCdMSRowIbF9ER+l0aiMTINHONmvlgPVSbn1U9rbQ5ZC0YN+48vT3cv0512Swdo9NphWsGUDqGuslaGsuIiJSVI7xvuKL3X17WIlaZGZrhpk2ox5v1dPtCO1eAy/+9PBVmnSmnUOs7gSOb+ph7c42Vu8LryD1gVuE3mgFJclOemMVRKL9mDslZUksAk3JfnZ1v0h51xaOS0yiIlJOkiheasRrnaXxGM0dbRA1qstiwZn5aA+9iSSb3ImYsaeth7J4lNJYjFk1J9PRVUJLVx8NtUmaO9ZQ1t9NhcdoJXlohzDgZNvLlJrJnH/Oh/CIMbF0Its7trOzYycRizCxpI5Sm0wsEqEzuY+u/qBStGrvWjbvbTt0ZcTMSHg/PYlO+syIYMRiESZEI8yI1XJieSM11bN4qXUdm9vWYX2dNHXspL8/SXmklBfbthIvj5MsizGht4ZkT5wFtbOoipTQ3LuXSFjZ2lM1lcryKcQicdydpzavoyY6haq+7UzpeZ4W66asfjIejVNTVU482U2CJLMrjmdK5VTi0TJikQhVsXIO9rWz9OAa3CHhCU5I9FCZdBrL6olblGhZLXUVk7GksXTbPsyMiuNm0pHsZm97Dwd6t1PW9jKVLdupjsWC6l6yn+pYH+XxGBOtlBmxWXgkSmLCbLr7+tjb2gWeJJLoIZroIN7XQX80wUGPEY3EKamayAnxOsoi8Vf84x4eDqt2buzr6MEx4j0HaO7fz8ZkC7uSnVjSuDQyjROraiiN1oS9ZE7kYGcfB6INvFw+gY0dKyjt2sFMqimbcDp1jecSq4DKeCWl0TK6E90kkgmiFg0+00uYUFbLzNrjiUfDRx10t0B/Dx29/SSTsOzAMlr6WolYhEll9dSX1tG+dRX1kXaSScdxWhJddFkPyd5OuloP4NE4PXVz6a2dFaYXBYtQnXS697XQ3Rf8rxhJDnb0YvESOuMJGuPbqWjdDEBFtJyOxsvoa29mf1878b4WurwPx6mkigo7Ho9X0lbWSHtPgsrkNqodamJV9LXspMf7OLF8MlGCK8eVkVLm1DaQTCZpTXazNTmZl0qPY/XBxezwBLRspaJnPWXRjbQZeDRKvKeXqRGYSx29sRoMKC2J0et9NJaMztM5Mq1gfR94zsweJtiv3gzck7OoREREcsjMXg+cBpQNlLn75440n7tvD//uDo+J5zHyHnIHL1M93aaTTEL3weAVXoHBk7Dq4WB40slQM/AITIOGk6C3A6omQyxoqtXYAI2zg3uadhzsPtRMCYJONyIGRvA3Ho0woSKOmdHddxk9/clDzckMKI1FKYlFuC5NqIlkgva+9kPjJdESymPl6dPy5Iju8ZpWNY1pVeke9Xm4eeHc+tkwJ+NFHnL6tHMJ7gDJjnfMC85dmBmJRDJYv5HDubo7/d5PPPLqKxQzgTPJrCPP2SPptDfRd3j/AYjGD7UTnDV42o590NsGkRhUToZY5g+mPnQHoTsz+7q4GIfuVvDU+wMt2D+jcSoJLmdn7V6XkuA+wopw9Orj05ysaVqQrU97tWQS+jogWjqi9TYSDcAJwGuA1P22q7+LLa1bqCqpoipeRU28+pWPe8jhA86Hkmkvgl8ws18BA0/ve6+7/yl3YYmIiOSGmd1J8DvkCuC/gZuA5zKYrxKIuHtbOHw18DlG2ENultMpfl0HYO966G2H9j2Q6IX+bujYO/Q8c18PU+e9uryiLu3kpbEoTQ2VGYdUFo8e6hgjE9FIlNrS2iNPSP470Mil1B+10TT3uZgZ8SGaf+VMdASfV1kfvI6FGZSE1ZySzPe5oheJQGk2H5GbufJYOXPqjuIMQw6NpFuaCqDV3b9nZpPMbJa7v5SrwERERHLkInefZ2bL3f1fzOzLwEMZzDcFeDj8ERkDfuDuvzaz5xlBD7njXqIfusIOgVu2wrpHD78XjUFFffDDdOA16ZTDP1gBIvHgKoCISIHKtJv2fyboSXAO8D0gDvwPcHHuQhMREcmJgafadprZ8cA+0rQWGszdNwFnpinfxwh7yB2XEv2wbz2s+smr3zvxtTBtPkQyv2okIlKoMr2C9WaCXpeWQNAG3czycx1QRETk2PzMzCYAXyI4rjnw7fyGNMb198Dz3wluxAeY2BRUqCC4GlU+MW+hiYhkW6YVrF5397Bb2oF26CIiIkXFzCIEDwU+CDxoZj8Hyty9Jc+hjV3blsC6R4LhWAkseB+UTcjLjeciIqMh0wrWj83sm8AEM/sA8D50tk9ERIqMuyfDe64uDMd7SHnYsGRRf09wv9RA5er4s+Ckq9UMUETGvEx7EfwPM3sd0EpwH9Zn3H1RTiMTERHJjUfN7C3AQ+6ubtFzYd0jwZWrASddDY3n5C8eEZFRdMQKlplFgUfc/bWAKlUiIlLsPgFUAv1m1k3wmCN395r8hjUGuMPaX8KO5cF4/eyg6+bjz85vXCIio+iIFSx3T5hZp5nVqo26iIgUO3dXJ0254A47lh2uXF30N1Bald+YRETyINN7sLqBFWa2COgYKHT3v81JVCIiIjlkZhOBk4CygTJ3fzJ/ERW5g1vgT/ceHj/3/apcici4lWkF6xfhS0REpKiZ2V8BHwMagaXABcDTwJX5jKso7V0Pq38WdGgBMP28oDlgRV1+4xIRyaNhK1hmNsPdt7j73aMVkIiISI59DDgXeMbdrzCzucC/5Dmm4tLbCXvXwdpfBeMTm4L7raafl9ewREQKwZGuYP0EmA9gZg+6+1tyH5KIiEhOdbt7t5lhZqXuvsbM5uQ7qKKQ6INNT0DzC4fLmi6BWZfmLSQRkUJzpApW6lMAT8hlICIiIqOk2cwmEJxEXGRmB4DteY6pcLlDsh+2L4UNvzlc3ngunPAaiMbzF5uISAE6UgXLhxgWEREpSu7+5nDws2b2OFAL/DqPIRWmZAJ2LIX1vwFPHi6vOR7mvQ3iZUPPKyIyjh2pgnWmmbUSXMkqD4dBzwwREZEiY2ZlwAeBE4EVwHfc/Xf5jaqAPXsndLceHm+6BI47Hcon5i8mEZEiMGwFy92joxWIiIhIjt0N9AG/B64DTiXo8EIGuMP+TbD8x4fLXvNJiETyF5OISJHJtJt2ERGRYnequ58BYGbfAZ7LczyFI5mAnSsO9woIMOlkOOUGVa5EREZIFSwRERkv+gYG3L3fzIabduxLJuGlJ6B9N+x/6XD55FNg2jkwYXreQhMRKWaqYImIyHhx5qB7ictT7jMef/cVr/sV7FgeDNdMhfoT4bgzoKw2v3GJiBQ5VbBERGRc0H3FIXdY9+vDlavL/l5drYuIZJEaVouIiIwnO1cEz7QCOPudqlyJiGSZrmCJiIiMF8kkrP1lMHzRR6G0Or/xiIiMQapgiRyjhZENNLfexx1L6w+VffisD+cxIhGRIexeFTQRnDpPlSsRkRxRBUtERGSs62mH1T+DnrZg/MTX5jceEZExTBUsERGRsSyZhMV3Ha5cNZ4LsdK8hiQiMpapgiUiIjKWrXooqFyV1cKFar4sIpJrqmCJiIiMRb2d8PIfYe/6YPzc9+c3HhGRcULdtIuIiIxFW/4Izc8Hw3OuVbNAEZFRoitYIiIiY82WZ2BrWLma/x6onprfeERExhFVsERERMaS/h7Y+HgwfPpboHZafuMRERlnVMESyZKnN+47NPzhs/IYiIiMXy89CZv/EAxPOwcmnZzfeERExqGivwfLzK41s7VmtsHMbs13PFKEHv/i4ZeIyDDMLGpmfzKzn4fjd5nZS2a2NHydFZabmX01PDYtN7P5OQ+ur+tw5aq2EU56Xc4/UkREXq2or2CZWRT4OvA6oBl43swWuvuL+Y2sSKVWMK64LX9x5NpwFanB7x3levjKonVHNd/HX5f+bPMdS+94xRWyC2fX8+GzDne3nPp5pZN+86r5U6fNpTuW3pHRZw+e7umN+5hf87ZD40Oth2zINMZCU6xxj0EfA1YDNSll/9vdHxg03XXASeHrfOAb4d/c2bcx+Nu4QJUrEZE8KuoKFnAesMHdNwGY2Y+AG4H8VLAyraAM9wM/F/MNt4yh5svn1ZzBMWWady5iPsplLmm976jme/eDmU339MZ9PL3x82nfu3BS/ZDTXji7Pt0sh6ZLZ37N215V4Xn3g4c/e7hlQvqKQbrPS11ndywdfplDLSM1lqHyGTwdvLKCOpBrurhHssxMKpajKdPPPta4x3Klz8wagdcDXwA+cYTJbwTucXcHnjGzCWY21d135CzA9Y8Ef2delLOPEBGRI7Pgu784mdlNwLXu/lfh+LuB8939o4OmuwW4JRydA6wd1UCPTgOwN99BZNlYy0n5FDblU9iGymemu08a7WAyYWYPAF8EqoG/d/c3mNldwIVAD/AYcKu794RNCG9396fCeR8DPunuL6RZbjaPUWNtP0ml3IrPWM0LlFuxOtbcMjpGFfsVLEtT9qoao7t/C/hW7sPJHjN7wd0X5DuObBprOSmfwqZ8Clux5WNmbwB2u/tiM7s85a3bgJ1ACcFx5pPA58jw+ATZPUYV23odCeVWfMZqXqDcitVo5VbsnVw0A9NTxhuB7XmKRURExq6LgRvMbDPwI+BKM/sfd9/hgR7gewRN10HHJxGRcavYK1jPAyeZ2SwzKwHeDizMc0wiIjLGuPtt7t7o7k0Ex5rfuvu7zGwqBL0GAm8CVoazLATeE/YmeAHQktP7r0REpGAUdRNBd+83s48CjwBR4LvuvirPYWVLUTVpzNBYy0n5FDblU9jGSj73mtkkgiaBS4EPhuW/BK4HNgCdwHtHKZ6xsl7TUW7FZ6zmBcqtWI1KbkXdyYWIiIiIiEghKfYmgiIiIiIiIgVDFSwREREREZEsUQUrj8wsamZ/Cp+XQngz9BfMbJ2ZrTazv00p/6qZbTCz5WY2P7+RDy1NTleZ2RIzW2pmT5nZiWF5qZndF+b0rJk15TPudMxss5mtCGN/ISyrM7NFZrY+/DsxLC/4bTREPl8yszVhzA+b2YSU6W8L81lrZtfkL/L00uWT8t7fm5mbWUM4XpTbJyz/m3AbrDKzf08pL+jtA0Puc2eZ2TMDZWZ2Xlhe8Nuo0JnZteH+sMHMbs13PJnI1vesmd0cTr/ezG7OUy7fNbPdZrYypSxruZjZOeG62hDOm+5RAKOZ22fNbFu47Zaa2fUp76X9fhpqH7Wgs7Jnw5zvs6DjstHIa7qZPW7Bb65VZvaxsLzot9swuY2F7VZmZs+Z2bIwt38ZLh4b5jfnSHPOmLvrlacX8AngB8DPw/H3AvcAkXB8cvj3euBXBDdRXwA8m+/YR5DTOuCUcPjDwF0pw3eGw28H7st37Gly2Qw0DCr7d4IHiQLcCvxbsWyjIfK5GoiFw/+Wks+pwDKgFJgFbASi+c7hSPmE5dMJOr55eeD9It4+VwC/AUrD8YHvhILfPsPk9ChwXcp2eaJYtlEhvwg6etoInEDwTK5lwKn5juso95ERfc8CdcCm8O/EcHhiHnK5DJgPrMxFLsBzBA+1tnDe6/Kc22cJHrg9eNq030/D7aPAj4G3h8N3Ah8apbymAvPD4WqC3yynjoXtNkxuY2G7GVAVDseBZ8PtkTYehvjNeTQ5Z/rSFaw8MbNG4PXAf6cUfwj4nLsnAdx9d1h+I3CPB54BJljYNXAhGSInB2rC4VoOPwfmRuDucPgB4KrRPBt3DFLjvpugW+aB8oLfRoO5+6Pu3h+OPkPwrB4I8vmRu/e4+0sEPaGdl24ZBegrwD/wyoe6FuX2IfhOuN2DZywN/k4o1u0z3HdCMW6jQnEesMHdN7l7L8Gzum7Mc0xHa6Tfs9cAi9x9v7sfABYB14520O7+JLB/UHFWcgnfq3H3pz34ZXhPyrJybojchjLU91PafTQ89l9J8FsAXrmecsqD59gtCYfbgNXANMbAdhsmt6EU03Zzd28PR+Phy4eJZ6jfnCPKeSQxqoKVP/+X4EdgMqVsNvA2C5rN/MrMTgrLpwFbU6ZrZvh/knxJl9NfAb80s2bg3cDtYfmhnMIf+C1A/eiFmhEHHjWzxWZ2S1g2xcNn2YR/J4flxbCN0uWT6n0EZ9egSPMxsxuAbe6+bNC0RZkPcDJwadik4Xdmdm5YXgz5QPqc/g74kpltBf4DuC0sL5acClWxrr9sfM8Wcu7ZymVaODy4PN8+GjaV++5AMzpGnls9cDDlZF9ecgubjZ1NcDVkTG23QbnBGNhuFtySshTYTVCh3ThMPEP95szZd4oqWHlgZm8Adrv74kFvlQLd7r4A+Dbw3YFZ0iymoPrXHyanjwPXu3sj8D3gPwdmSbOYgsoJuNjd5wPXAR8xs8uGmbao8zGzTwP9wL0DRWnmL4Z8Pg18Js20xZpPjKC5yQXA/wZ+HJ51K4Z8IH1OHwI+7u7TCb4fvhNOWyw5FapiXX/Z+J4txtxHmksh5vgNghPDZwE7gC+H5UWXm5lVAQ8Cf+furcNNmqas2HIbE9vN3RPufhZBy5vzgFOGiWfUc1MFKz8uBm4ws80Elx2vNLP/IaghPxhO8zAwLxxuJrivZEAjh5vVFIp0Of0CONPdB86Y3AdcFA4fysnMYgRNhTJtfjAq3H17+Hc3wfY4D9g10Gwp/DvQZKvgt9EQ+RDejPsG4J1hEwYoznxeQ9CGelm4HzYCS8zsOIozn/MI4n4obA7xHMHV4QaKIB8YMqebgYfCSe7ncNPGosipgBXl+svS92wh556tXJo53IQ7tTxv3H1X+CM3SXBS+Ej/y0OV7yVoahcbVD4qzCxO8NvrXncf+G4aE9stXW5jZbsNcPeDwBMEJyKHimeo35w5+05RBSsP3P02d2909yaCm+1+6+7vAn5C0H4Ugh+L68LhhcB7LHAB0DJw6bpQpMuJoL1qrZmdHE72OoI2wBDkNNDLzk0E6yDfZ+MOMbNKM6seGCboDGIlr4z7ZuCn4XBBb6Oh8jGza4FPAje4e2fKLAuBt4c978wCTiK4UbcgDJHP8+4+2d2bwv2wmeAG350U6fYh5Tsh/D8qITioFfT2gWFz2k7w/QZBbuvD4YLeRkXgeeAkC3rRKiH4Hl6Y55iGlcXv2UeAq81sYtjc6eqwrBBkJZfwvTYzuyC8iv2elGXlhb3yHsk3E2w7GPr7Ke0+Gh77Hyf4LQCvXE+5zsEIrqKvdvf/THmr6LfbULmNke02ycJej82sHHgtwe/LoeIZ6jfniHIeUZA+Cr196DVsTyiXc7jHvQnAL4AVwNMEV38guFT5dYL2pSuABfmOewQ5vTmMeRnBGYYTwvIygrPXG8Kd+YR8xz0ohxPCmJcBq4BPh+X1wGMEPwofA+qKYRsNk88GgnbGS8PXnSnzfDrMZy2j2FvVseQzaJrNHO5FsFi3TwnwPwQHwCXAlcWwfY6Q0yXA4rD8WeCcYthGxfAi6OFsXbgOX/U/UWivbH7PEtxDuiF8vTdP+fyQoMlVH8EJnvdnMxdgQfhdsBH4GmB5zu37YezLCX58Tk2ZPu3301D7aLgvPBfmfD9hz6mjkNclBE2/lnP4OHj9WNhuw+Q2FrbbPOBPYQ4rgc8MFw/D/OYcac6ZvixciIiIiIiIiBwjNREUERERERHJElWwREREREREskQVLBERERERkSxRBUtERERERCRLVMESERERERHJElWwRLLAzBJmttTMVprZ/WZWke+YAMzsU1lYxpfMbI2ZLTezhweePSEiImOLmbWHf5vM7B1ZXvanBo3/MZvLFykkqmCJZEeXu5/l7qcDvcAHM53RzKK5C4sRV7DSxLMION3d5xE8E+K2bAQmIiIFqwkYUQUrg2PZK45H7n7RCGMSKRqqYIlk3++BEwHM7CdmttjMVpnZLQMTmFm7mX3OzJ4FLjSzz5jZ8+EVsG+FT2DHzJ4ws6+Y2ZNmttrMzjWzh8xsvZn9a8ry3mVmz4VX0b5pZlEzux0oD8vuHWq6dPGkJuPuj7p7fzj6DNCYu1UnIiIF4Hbg0vBY8fHwmPKl8Di13Mz+GsDMLjezx83sBwQPr0173BvieDRwtczCZa80sxVm9raUZT9hZg+ErSjuHTg2ihQ6VbBEssjMYsB1hAca4H3ufg7Bk9z/1szqw/JKYKW7n+/uTwFfc/dzwytg5cAbUhbb6+6XAXf+/3buJtSKMo7j+PcXRtoLglHgInqBoiTKKKNQojCiVVRccOGiKIIijFYm2KJ9iwikRRHUIlpFVIuyCCo3oSB5w6BNLxRUWvRG9qK3f4t5rh3vvXPUmIudy/ezmXlm/jPzzNn8z/+ZZwZ4DXgYuBK4N8m5Sa4ANgHrq2otMANsrqpt/PtkbXNfXE9/+twHvPmffyBJ0iTYBuxq+eMp4H7g56paB6wDHkhycYu9HtheVWtae17em5uP5lzrbmAtcDVwK/BkktVt3zXAo8Aa4BJg/aLcrTSwZae6A9ISsSLJR219F/B8W38kyV1t/QLgUuAHuuLmlZHjb0myFTgTWAXsB95o+15vy4+B/VX1DUCSz9o5NwDXAnva4N4K4MACfdw4Jm5uf+ZJsh04Arw0Lk6StOTcBlyVZKq1V9Lls7+A3VX1+UhsX97rswF4uapmgO+SvE9XxP3Szv01QMuxFwHjBgGl/wULLGkYv7enQkcluZluNO7GqjqU5D1gedv9R0smJFkOPANcV1VfJXliJA7gz7b8e2R9tr0MCPBiVR3v3ahxcUf7s+CByT10T9U2VlUd5zqSpKUlwJaq2nnMxi7P/Tan3Zf3xp27z2jOm8H/rZoQThGUFs9K4MeWZC4HbuiJm00+3yc5G5jqievzLjCV5HyAJKuSXNj2HU5y+gnE9UpyO/AYcEdVHTrJvkmSJs+vwDkj7Z3AQ7P5JMllSc5a4LhxeW80H436ANjU3vM6D7gJ2D3IXUiniCMB0uJ5C3gwyTTwKd0HIuapqp+SPEc3BfALYM/JXKSqPknyOPB2ktOAw3TvaX0JPAtMJ9nb3sPqixtnB3AG8E6bWvhhVZ3wVxIlSRNnGjiSZB/wAvA03fS8ve1DEweBOxc4blzeOyYfjWx/le7jSvuAArZW1betQJMmUpztI0mSJEnDcIqgJEmSJA3EAkuSJEmSBmKBJUmSJEkDscCSJEmSpIFYYEmSJEnSQCywJEmSJGkgFliSJEmSNJB/AIQN1WKcVIbwAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAGRCAYAAACtyi4tAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYnFWZ/vHvTcImeyAshoQwEEFQCNAsjozDImuURVFAhAwDZFRUHBw16CgqRuCniDIsCsMuGgGFBIJAQCLDDFsDgbAEEzBKEiDBsAqE7fn9cU5B0anurup+q6qr6/5cV11V76l3OaeXeuos7zmKCMzMzIqwXLMzYGZmg4eDipmZFcZBxczMCuOgYmZmhXFQMTOzwjiomJlZYRxUzMysMA4qZmZWGAcVMzMrzNBmZ6DR1llnnRg9enSzs2GD1D333PNMRAxvdj7MmqXtgsro0aPp7OxsdjZskJL0l2bnwayZ3PxlZmaFcVAxM7PCOKiYmVlhHFTMzKwwDipmZlYYBxUzMytM2w0ptmKMnjitYvq8U8Y1OCdmNpC4pmJmZoVxUDEzs8I4qJiZWWEcVMzMrDAOKmZmVhgHFTMzK0zdgoqklSTdJel+SQ9J+l5O31jSnZLmSPqNpBVy+op5e25+f3TZuU7I6Y9K2qssfe+cNlfSxHqVxczMqlPP+1SWArtFxEuSlgduk/R74Hjg9IiYLOnnwFHAOfn52YjYVNIhwKnAwZK2AA4BtgTeC9wk6X35GmcBewDzgbslTY2Ih+tYprbT3f0oZmaV1K2mEslLeXP5/AhgN+DKnH4xcEB+vX/eJr+/uyTl9MkRsTQi/gzMBXbIj7kR8XhEvAZMzvuamVmT1LVPRdIQSTOBRcB04DHguYh4I+8yHxiRX48AngDI7z8PrF2e3uWY7tLNzKxJ6hpUIuLNiBgLbEiqWby/0m75Wd28V2v6MiRNkNQpqXPx4sW9Z9zMzPqkIaO/IuI5YAawE7CmpFJfzobAwvx6PjASIL+/BrCkPL3LMd2lV7r+uRHREREdw4d7+XAzs3qp5+iv4ZLWzK9XBj4KPALcAhyUdxsPTMmvp+Zt8vt/iIjI6Yfk0WEbA2OAu4C7gTF5NNkKpM78qfUqj5mZ9a6eo782AC6WNIQUvC6PiGslPQxMlvQD4D7g/Lz/+cClkuaSaiiHAETEQ5IuBx4G3gCOjYg3ASR9EbgBGAJcEBEP1bE8ZmbWi7oFlYh4ANimQvrjpP6VrumvAp/q5lyTgEkV0q8Drut3Zs3MrBC+o97MzArjoGJmZoVxUDEzs8I4qJiZWWEcVMzMrDAOKmZmVhgHFTMzK4yDipmZFcZBxczMCuOgYmZmhXFQMTOzwjiomJlZYRxUzMysMPVcT2WkpFskPSLpIUnH5fTvSlogaWZ+7Ft2zAmS5kp6VNJeZel757S5kiaWpW8s6U5JcyT9Jq+rYmZmTVLPmsobwFcj4v2kFR+PlbRFfu/0iBibH9cB5PcOAbYE9gbOzmvcDwHOAvYBtgAOLTvPqflcY4BngaPqWB4zM+tF3YJKRDwZEffm1y+SVn0c0cMh+wOTI2JpRPwZmEtad2UHYG5EPB4RrwGTgf0lCdgNuDIffzFwQH1KY2Zm1WhIn4qk0aQFu+7MSV+U9ICkCyStldNGAE+UHTY/p3WXvjbwXES80SW90vUnSOqU1Ll48eICSmRmZpXUPahIWhX4LfCViHgBOAfYBBgLPAmcVtq1wuHRh/RlEyPOjYiOiOgYPnx4jSUwM7Nq1XONeiQtTwool0XE7wAi4umy988Drs2b84GRZYdvCCzMryulPwOsKWlorq2U729mZk1Qz9FfAs4HHomIn5Slb1C224HAg/n1VOAQSStK2hgYA9wF3A2MySO9ViB15k+NiABuAQ7Kx48HptSrPGZm1rt61lQ+DBwOzJI0M6d9kzR6ayypqWoe8G8AEfGQpMuBh0kjx46NiDcBJH0RuAEYAlwQEQ/l830DmCzpB8B9pCBmZmZNUregEhG3Ubnf47oejpkETKqQfl2l4yLicdLoMDMzGwB8R72ZmRXGQcXMzArjoGJmZoWp65BiG1hGT5xWMX3eKeManBMzG6yqqqlI+kC9M2JmZq2v2uavn0u6S9IXJK1Z1xyZmVnLqiqoRMTOwGGkO9s7Jf1K0h51zZmZmbWcqjvqI2IO8J+kGw7/GThD0mxJn6hX5szMrLVU26eylaTTSdPX7wZ8PK+Tshtweh3zZ2ZmLaTa0V9nAucB34yIV0qJEbFQ0n/WJWdmZtZyqg0q+wKvlM3FtRywUkS8HBGX1i131ifdDR02M6u3avtUbgJWLtt+T04zMzN7W7VBZaWIeKm0kV+/p6cDJI2UdIukRyQ9JOm4nD5M0nRJc/LzWjldks6QNDevCrlt2bnG5/3nSBpflr6dpFn5mDPydPtmZtYk1QaVv3f5kN8OeKWH/SFNX//V3KG/E3CspC2AicDNETEGuDlvA+xDWkNlDDCBtEIkkoYBJwI7kmYkPrFsCeJz8r6l4/ausjxmZlYH1fapfAW4QlJpZcUNgIN7OiAiniQtF0xEvCjpEdIa8vsDu+TdLgZmkIYp7w9ckhffukPSmnlBr12A6RGxBEDSdGBvSTOA1SPi9px+CXAA8Psqy2RmZgWrKqhExN2SNgc2I62RMjsiXq/2IpJGA9sAdwLr5YBDRDwpad282wjgibLD5ue0ntLnV0g3M7MmqWVCye2B0fmYbSQREZf0dpCkVUnr1H8lIl7ooduj0hvRh/RKeZhAaiZj1KhRvWXZzMz6qNqbHy8FfgzsTAou2wMdVRy3PCmgXBYRv8vJT5fWqc/Pi3L6fNI0MCUbAgt7Sd+wQvoyIuLciOiIiI7hw4f3lm0zM+ujamsqHcAWub+jKnkk1vnAIxHxk7K3pgLjgVPy85Sy9C9KmkzqlH8+N4/dAPywrHN+T+CEiFgi6UVJO5Ga1Y4A/qva/JmZWfGqDSoPAuuTO96r9GHgcGCWpJk57ZukYHK5pKOAvwKfyu9dR7rJci7wMnAkQA4eJwF35/2+X+q0Bz4PXES6h+b3uJPezKypqg0q6wAPS7oLWFpKjIj9ujsgIm6jcr8HwO4V9g/g2G7OdQFwQYX0TsBrvZiZDRDVBpXv1jMTZmY2OFQ7pPiPkjYCxkTETZLeAwypb9bMzKzVVDv66xjgSuAXOWkEcHW9MmVmZq2p2mlajiV1vL8Aby/YtW6PR5iZWdupNqgsjYjXShuShtLNjYZmZta+qg0qf5T0TWDlvDb9FcA19cuWmZm1omqDykRgMTAL+DfSPSVe8dHMzN6l2tFfb5GWEz6vvtkxM7NWVlVQkfRnKvShRMQ/FJ4jMzNrWbXM/VWyEmlqlWHFZ8fMzFpZVX0qEfG3sseCiPgpsFud82ZmZi2m2uavbcs2lyPVXFarS47MzKxlVdv8dVrZ6zeAecCnC8+NmZm1tGqbv3Yte+wREcdExKM9HSPpAkmLJD1YlvZdSQskzcyPfcveO0HSXEmPStqrLH3vnDZX0sSy9I0l3SlpjqTfSFqhtqKbmVnRqm3+Or6n97sswlVyEXAm0HXJ4dMj4sddzr8FcAiwJfBe4CZJ78tvnwXsQVrp8W5JUyPiYeDUfK7Jkn4OHAWcU015zMysPqq9+bGDtCDWiPz4HLAFqV+lYt9KRNwKLKn0XgX7A5MjYmlE/Jm0UNcO+TE3Ih7P08RMBvbPq0ruRprkEuBi4IAqr2VmZnVSyyJd20bEi5CasYArIuLoPlzzi5KOADqBr0bEs6RAdUfZPvNzGsATXdJ3BNYGnouINyrsb2ZmTVJtTWUU8FrZ9mvA6D5c7xxgE2AsaWni0gCASitERh/SK5I0QVKnpM7FixfXlmMzM6tatTWVS4G7JF1F+vA+kGX7SnoVEU+XXks6D7g2b84HRpbtuiGwML+ulP4MsKakobm2Ur5/peueC5wL0NHR4dmVzczqpNrRX5OAI4FngeeAIyPih7VeTNIGZZsHAqWRYVOBQyStKGljYAxwF3A3MCaP9FqB1Jk/Na9nfwtwUD5+PDCl1vyYmVmxqq2pALwHeCEiLpQ0XNLGuVO9Ikm/BnYB1pE0HzgR2EXSWFJtZx5pxmMi4iFJlwMPk+6DOTYi3szn+SJwA2n54gsi4qF8iW8AkyX9ALgPOL+GspiZWR1UO6T4RNIIsM2AC4HlgV+SVoOsKCIOrZDc7Qd/rg1NqpB+HWmq/a7pj5NGh9kAMnritG7fm3fKuAbmxMyaodqO+gOB/YC/A0TEQjxNi5mZdVFtUHkt92MEgKRV6pclMzNrVdUGlcsl/YI04uoY4Ca8YJeZmXVR7cqPP85r079A6lf5TkRMr2vOzMys5fQaVCQNAW6IiI8CDiRmZtatXpu/8tDelyWt0YD8mJlZC6v2PpVXgVmSppNHgAFExJfrkiszM2tJ1QaVaflhZmbWrR6DiqRREfHXiLi4URkyM7PW1VufytWlF5J+W+e8mJlZi+stqJRPMf8P9cyImZm1vt6CSnTz2szMbBm9ddRvLekFUo1l5fyavB0RsXpdc2dmZi2lx5pKRAyJiNUjYrWIGJpfl7Z7DCiSLpC0SNKDZWnDJE2XNCc/r5XTJekMSXMlPSBp27Jjxuf950gaX5a+naRZ+Zgz8rr1ZmbWRNXO/dUXFwF7d0mbCNwcEWOAm/M2wD6khbnGABNIyw4jaRhpHZYdSdPcn1gKRHmfCWXHdb2WmZk1WN2CSkTcCizpkrw/UBqefDFwQFn6JZHcQZq4cgNgL2B6RCyJiGdJ08Tsnd9bPSJuz7MnX1J2LjMza5JaVn4swnoR8SRARDwpad2cPgJ4omy/+Tmtp/T5FdIrkjSBVKth1KhR/SxCY3W36JUXvDKzgaiezV+1qNQfEn1Irygizo2IjojoGD58eB+zaGZmvWl0TeVpSRvkWsoGwKKcPh8YWbbfhsDCnL5Ll/QZOX3DCvu3jZ6W7R2oXOsyG/waXVOZCpRGcI0HppSlH5FHge0EPJ+byW4A9pS0Vu6g35M0Df+TwIuSdsqjvo4oO5eZmTVJ3Woqkn5NqmWsI2k+aRTXKaRVJI8C/gp8Ku9+HbAvMBd4GTgSICKWSDoJuDvv9/2IKHX+f540wmxl4Pf5YX3QirUeMxuY6hZUIuLQbt7avcK+ARzbzXkuAC6okN4JfKA/eTQzs2I1uk/FbBnuazEbPAbK6C8zMxsEHFTMzKwwDipmZlYYBxUzMyuMg4qZmRXGQcXMzArjoGJmZoVxUDEzs8I4qJiZWWF8R/0A4fm3zGwwcE3FzMwK05SgImmepFmSZkrqzGnDJE2XNCc/r5XTJekMSXMlPSBp27LzjM/7z5E0vrvrmZlZYzSzprJrRIyNiI68PRG4OSLGADfnbYB9gDH5MQE4B1IQIk2nvyOwA3BiKRCZmVlzDKTmr/2Bi/Pri4EDytIvieQOYM28auRewPSIWBIRzwLTgb0bnWkzM3tHs4JKADdKukfShJy2Xl7Rkfy8bk4fATxRduz8nNZdupmZNUmzRn99OCIWSloXmC5pdg/7qkJa9JC+7AlS4JoAMGrUqFrzamZmVWpKTSUiFubnRcBVpD6Rp3OzFvl5Ud59PjCy7PANgYU9pFe63rkR0RERHcOHDy+yKGZmVqbhQUXSKpJWK70G9gQeBKYCpRFc44Ep+fVU4Ig8Cmwn4PncPHYDsKektXIH/Z45zczMmqQZzV/rAVdJKl3/VxFxvaS7gcslHQX8FfhU3v86YF9gLvAycCRARCyRdBJwd97v+xGxpHHFMDOzrhoeVCLicWDrCul/A3avkB7Asd2c6wLggqLzaGZmfTOQhhSbmVmLc1AxM7PCeELJBvKkkWY22LmmYmZmhXFNJeuuFjHvlHENzomZWetyTcXMzArjoGJmZoVxUDEzs8I4qJiZWWHcUV8HHjpsZu3KQaUfHDyapy8/e4/kM6s/N3+ZmVlhHFTMzKwwLR9UJO0t6VFJcyVNbHZ+zMzaWUsHFUlDgLOAfYAtgEMlbdHcXJmZta+WDiqkZYjnRsTjEfEaMBnYv8l5MjNrW60++msE8ETZ9nxgxyblxdqIR/6ZVdbqQUUV0mKZnaQJwIS8+ZKkR6s8/zo6lWf6mrkWtA4MnPLq1Lqfrx7l3ajg85m1lFYPKvOBkWXbGwILu+4UEecC59Z6ckmdEdHR9+y1FpfXzPqr1ftU7gbGSNpY0grAIcDUJufJzKxttXRNJSLekPRF4AZgCHBBRDzU5GyZmbWtlg4qABFxHXBdnU5fc5NZi3N5zaxfFLFMv7aZmVmftHqfipmZDSAOKmZmVhgHFTMzK4yDipmZFcZBxczMCuOgYmZmhXFQMTOzwjiomJlZYRxUzMysMA4qZmZWGAcVMzMrjIOKmZkVxkHFzMwK46BiZmaFcVAxM7PCOKiYmVlhHFTMzKwwDipmZlYYBxUzMyuMg4qZmRXGQcXMzArjoGJmZoVxUDEzs8IMbXYGGm2dddaJ0aNHN/SaS15ayoLnX317+4Mj1mjo9a1x7rnnnmciYnijr9uMv2trH/fcc88LwO0RsXdv+7ZdUBk9ejSdnZ0Nu96v/u8xvjl1NhfutznfnDqbeaeMa9i1rfEk/aUZ123037W1F0lzqgko0IZBpdE+84+bvP1cem1mNli5T6UBHEzMrF04qJiZWWEcVMzMrDAOKmZmVhh31JtZSxg9cVq373lU5cDhmoqZmRXGQcXMzArjoGJmZoVxn4qZWZnu+m7cb1Md11TMzKwwDipmZlYYBxUzMyuMg4qZmRXGQcXMzArjoGJmZoVxUDEzs8I4qJiZWWEcVMxalKQJkjoldS5evLjZ2TEDHFR61NOsqGbNFhHnRkRHRHQMHz682dkxAxxUulUKKA4sZmbVc1DpRmmeH8/3Y2ZWPQeVHjigmJnVxkHFzMwK07CgImmepFmSZkrqzGnDJE2XNCc/r5XTN5d0u6Slkv6jh3NuLOnOfPxvJK3QqPKYmdmyGl1T2TUixkZER96eCNwcEWOAm/M2wBLgy8CPeznfqcDp+fhngaPqkGczM6tSs5u/9gcuzq8vBg4AiIhFEXE38Hp3B0oSsBtwZdfjzcysORoZVAK4UdI9kibktPUi4kmA/LxuDedbG3guIt7I2/OBEYXl1szMatbI5YQ/HBELJa0LTJc0u5/nU4W0qLhjCmITAEaNGtXPy5qZWXcaVlOJiIX5eRFwFbAD8LSkDQDy86IaTvkMsKakUmDcEFjYzbV957GZWQM0JKhIWkXSaqXXwJ7Ag8BUYHzebTwwpdpzRkQAtwAH9eV4MzMrXqOav9YDrkp96wwFfhUR10u6G7hc0lHAX4FPAUhaH+gEVgfekvQVYIuIeEHSdcDRuebzDWCypB8A9wHnN6g8ZmZWQUOCSkQ8DmxdIf1vwO4V0p8iNWdVOte+Xc67Q3E5NTOz/mj2kGIzMxtEHFTMzKwwDipmZlYYBxUzMyuMg4qZmRXGQcWa4rbZFe9TNbMW56BifdbXpZZvm72Qz192nwNLP0maIKlTUufixYubnR0zwEHF+qgUUPoSWHbe/L2cc9g27Lz5e4vOVlvx9EM2EDmoWJ+Ullru65LLDihmg1MjZym2QaavAcWs2fradGu9c02lzfifyczqqaqgIml1SSdLulTSZ7q8d3Z9smZF608/iJlZNaqtqVxIWhTrt8Ahkn4racX83k51yZkVrr/9IGZmvak2qGwSERMj4uqI2A+4F/iDpLXrmDerAwcUM6unajvqV5S0XES8BRARkyTNB24FVq1b7szMrKVUW1O5BtitPCEiLga+CrxWdKbMzKw1VRVUIuLrEXFThfTrI2JM8dmywWrGwwuanQUzqyMPKbaGmfHwAr5w2UwHFrNBzDc/WsPsssUIzj4sPZvZ4FR1UJG0HLBTRPxfHfMz6Hx00jRu+pZHXJU4oFhvfB9Va6u6+SuP/DqtjnkZdD46aRpzX0zPA9U2/gc2swLV2qdyo6RPSlJdcjPI3PStcWy6GgO2prLNxGk8S/8Ci79Vmlm5WoPK8cAVwGuSXpD0oqQXqjlQ0jxJsyTNlNSZ04ZJmi5pTn5eK6dL0hmS5kp6QNK23ZxzhqRH8zlnSlq3xvLU3UANKAD3nTKOtfJzX3jaFzPrqqagEhGrRcRyEbF8RKyet1ev4RS7RsTYiOjI2xOBm/Ow5JvzNsA+wJj8mACc08M5D8vnHBsRi2opj/U9oICnfTGzZdUUVHIN4rOSvp23R0raoR/X3x+4OL++GDigLP2SSO4A1pS0QT+uY3XigGJm5Wpt/job+BBQmqn4JeCsKo8NUp/MPZIm5LT1IuJJgPxcar4aATxRduz8nFbJhbnp69vd9fW047KrbpIys2aoNajsGBHHAq8CRMSzwApVHvvhiNiW1LR1rKSP9LBvpeAQFdIOi4gPAv+UH4dXOln5sqvD1l6nyuy2rmr6OvoTdBywzKw7tQaV1yUNIX/ASxoOvFXNgRGxMD8vAq4CdgCeLjVr5edSn8h8YGTZ4RsCCyucc0F+fhH4VT5nj557+TWef2VwT1fWW19HfzrY3TlvZj2pNaicQQoI60qaBNwGnNzbQZJWkbRa6TWwJ/AgMBUYn3cbD0zJr6cCR+Q+nJ2A50vNZGXnHCppnfx6eeBj+Zw9WvM9K7DGytVWrlpXT30d/ekHcef8wNGOzbo28NU6+usy4OukQPIkcEBEXF7FoesBt0m6H7gLmBYR1wOnAHtImgPskbcBrgMeB+YC5wFfKJ1I0sz8ckXgBkkPADOBBXnfHg1ZbuDeYnPM+be+/bpRNYG+XMcBZWAob9YdPnx4s7NjBtQ++uvSiJgdEWdFxJkR8YikS3s7LiIej4it82PLiJiU0/8WEbtHxJj8vCSnR0QcGxGbRMQHI6Kz7Fxj8/PfI2K7iNgqn/O4iHiztuIPHMecfyvT57zIMeffWnUT09iy9y+9bW6v1yg/n2scZlYPtTZ/bVm+kftXtisuO+3rvKM+wh5jVuO8oz5S1Qf+2InTeC4/X3rbXL597aM9BpZKgcoBxcyKVlVQkXSCpBeBrcrupH+R1LE+pZfDrUrnHfXOgLjePvBnnjKONYE9Pjicw3felJM+thmH77xpt/v3pWZSTe2nO2fe+EifjzWz1lXtIl0nR8RqwI/K7qRfLSLWjogT6pxH68YeHxzOFbMW87XL7uoxoJTUGlB6q/0A7PK9ZZvozrzxEX78h8fZ/QceIWbWbmpt/vpWwXfUWz/86LAd+NQHh/Ojw4r/Fezxwfdy0sc249vXPgpU7t/Z5XvTmPfKsoHli3u+n01Whcdegi9d7JUSzNpJrUHlLPp+R73VQdEBZdLVD/DU8y9zxvQ5ywSUroFl3iswemWYceKyNaDHXoKPv38t/mv8PxaaPzMb2Bp5R70NYFPv/QuTrn6A8+54gvNvmcuUzvnL9MOUN5+VAsy8V1gmrfR8zSPPNiLrZjaANOyOehu4pt77F467/EE+OGoNjtlpJFfc8QR/J40s666Dv2t6eUDxcGWz9lXEHfU/LDxXVjeV+kb223YjfvbpD7DfthvxrQO2entk2cxegkJ50OipVmNm7aOIO+qvqEfGrHg93VS537YbvWu7u4DS0w2ZDiRmVmtNBeBp4H+A/wNW7m5VRht4+luLaKXJJFshj2aDUa3TtJwEPEBqBjstP35ch3y1tXp+ILbDZJKtFPzMBpuhNe7/aWCTiBjcc8c3UaUO7/6ca/dNV+X8o/+5iKwVkqdGmHfKuJbJqxWjuy8Q/htovFqbvx4E1qxHRiwpqjZQ+ie7ee5LHPXff+x3vlrt278/TMyao9agcjJwn6QbJE0tPeqRsXZWxAdi6RxF1VQqBbtDz7mp3+c1s8Gl1uavi4FTgVn4/pQBr+hv610Dyu1/Wcqh59zErz//0W6P+dpld9VlGhkzG5hqrak8ExFnRMQtEfHH0qMuObMB7def/ygf2mjFXgNKacJLM2sPtdZU7pF0Mmm536WlxIi4t9BcWUvoKaBAnpfMNRWztlJrTWUbYCfSXfQeUtyiKnW216sD3gHFrL3UVFOJiF3rlRFrjEpDloscxlzN9T0yqxiSJgATAEaNGtXk3NSuESMJW2W04mBS8x31ksZJ+rqk75Qe9ciYFWOZ6eorjOJq1E2NrTYseaCLiHMjoiMiOoYPH97s7JgBtd9R/3PgYOBLgIBPARv1eJA1TbfroFQIHo2oPbTKHflm1ne11lT+MSKOAJ6NiO+RFuwaWc2BkuZJmiVppqTOnDZM0nRJc/LzWjldks6QNFfSA93NLyZpu3zOuXl/1VieAedbk+8p7FzVfIhv3uBagwOK2eBWa1B5NT+/LOm9wOvAxjUcv2tEjI2Ijrw9Ebg5IsYAN+dtgH2AMfkxATinm/Odk98v7bt3DXkZcL41+R4um/lUXQJLpSanzSdO41X6F1ga2elvZgNfrUHlGklrAj8C7gXmAb/ux/X3J91QSX4+oCz9kkjuANaUtEH5gXl79Yi4PSICuKTs+JY06ZDtOGzs+kw6ZLt+n6v8g727ZrDZp4xjpfzcn2tUcy0zaw9VBxVJy5FqFc9FxG9JfSmbR0S1HfUB3CjpnjxqBWC9iHgSID+vm9NHAE+UHTs/p5UbkdN72qflFBlQ3l7yt4dmsNl58sW+aGanv5kNTFUHlYh4i3RfSml7aUQ8X8O1PhwR25Kato6V9JEe9q3UNxJ92CftKE2Q1Cmpc/HixdXldgDrLQj09GHf3blKz7MXLKkpL83q9DezganW5q8bJX2yLx3iEbEwPy8iLUm8A/B0qVkrPy/Ku8/n3QMANgQWdjnl/Jze0z6law/ooZe11BSqbV6q9oO9PADNXrCE439zP7MXLOH5V95Z3cBNWWZWrVqDyvHAFcBSSS9IelHSC70dJGkVSauVXgN7kqbRnwqMz7uNB6bk11OBI/IosJ2A50vNZCV5+0VJO+Ugd0TZ8S2j1j6IejQvlc61+Yhh/OTgrdlg2KpMuW8Bz7/yWrf5+8l1D/V4zu9dObOw/JlZ66h1jfrVImK5iFghIlbP26tXceh6wG2S7gfuAqZFxPXAKcAekuYAe+RtgOuAx4G5wHm/+oGqAAAZL0lEQVTAF0onklT+afV54L/zfo8Bv6+lPANBX4JEfwJKb8Fr8xHDWGPlFdh/mxGssfIKFfP3k+se4oxb53UbWL535Uwu7FzgwGLWhmqdUJJ8L8kYYKVSWkTc2tMxEfE4sHWF9L8Bu1dID+DYbs41tux1J/CBavM+UFUTJIpaCbLac62x8grd5u/4fbd813PX85140Nh3PZtZ+6j1jvqjgVuBG4Dv5efvFp8tK1dtE9noidN63KfIprPygFIpb10Dyo5V5N3MWl+tfSrHAdsDf8mTS24DtP5wqgGummBQ6V6Rns5VlO1HLv+u50p2nDiNp+k+sPjeFrPBo+Y76iPiVQBJK0bEbGCz4rPVno696H+7fa+3YFDN8OF6uOLYPdl+5PJcceyeb6d1dAkOT5M61e7sJl++t8Vs8Kg1qMzPd9RfDUyXNIVuhvFabY696H+ZNvu5HgNLb+adMq4pH8xdA8ozvBNYSrWPp3s5hwOK2eBQ63oqB+aX35V0C7AGcH3huWpDZ/3Lh+Gi/03P/dTMNUs6TxlHx8RpdJbVPryGig0G3TXP+m/73aqqqUhaSdJXJJ0p6d8kDc3r00+NiNd6P4NVo6iAUv7cDJ1d/slK/3Q/u/7hZfYtcvJMM2u+apu/LgY6gFmkaVZO63l3a5aB2j/xs+sf5vQZf35XYKnHrMxm1lzVBpUtIuKzEfEL4CDgn+qYp5Y2EEYwVRtQqs3rt39zb3+yA8Bxe2/Bv++yMcftvcXbaeWzMg+En5uZ9V+1QeX10ouIeKNOeWl5vTU9bTaAPjirbSb79m/u5dL7niwssAAced6Mt9PKA0pfA4sDktnAUW1Q2TrP9fWCpBeBrWqZ+6td9NT0tNnEaSxl4ASWapvJTjp4Ww7fZgNOOrji4pvL2OuUnst35HkzuOWxv78rsPSnyW4g9CGZ2TuqCioRMSTP9VWa72tojXN/tY3uPhgfPWUcK+bnZpj71HPLpFX7IV5LQHn0uWUDy3GX3P726wuP2YVdN1mFC4/ZpU956Wqg9iE1wmBb0sEGh1rvU7F+aFZA2euUaXzjygcqBpYi3TBxHJutmZ5LjrvkdqY8vGSZwFKkdgwoMPCXdLD25KAyyO3/01R7eOmlF9l0/TXrco3ypqfygALwsyM+xP5bDONnR3yoLtc2s4HFQWWQm/KVcWy9/rIf9kXpqU9jv9PTBJcOKGbtw0GlDUz5Sv2ahyr1aTz/ymvsd/o0Hshzs7gT3ax9OKhYn+xYVkPpGlCm3LeASz+3B1utl9J66/P42GkOOmaDhYOK1aw0lX2lpq/yVSOn/nvvE1x+7LRpPLi4/oFlxsML6nr+tw1Zvvs1AMzagIOK1ezOU8axHt0P5y1fNbI31351HB8Ynp7rZcbDC/jCZTPrHljmLX6Boauvu1FdL2I2wNW8nLAZvLM2ShHDeesZUAB22WIEZx+Wnutp9PDVeeOFRX+p60XMBjjXVKwt1DugvO3N11/vfSezwctBxRrKI8HMBjcHFWuYvs7T9fqbb9UjO2ZWBw0NKpKGSLpP0rV5ezdJ90p6UNLFkobm9LUkXSXpAUl3SfpAN+e7SNKfJc3Mj7GNLE9ffO/Kmc3OQtN07dj/2mV39XrM62++xbxn/u7AYtYiGl1TOQ54BEDScqTFvw6JiA8AfwHG5/2+CcyMiK2AI4Cf9XDOr0XE2PwY0J/Y37tyJhd2LnBgIQWUK2Yt7jWwLD9kOUavswrLD3Gl2qwVNOw/VdKGwDjgv3PS2sDSiPhT3p4OfDK/3gK4GSAiZgOjJa3XqLzWy4kHjeXIjhGceFBzKlQn/LqzKdet5IEnFr/ruScOKGato5H/rT8Fvg6U2jGeAZaX1JG3DwJG5tf3A58AkLQDsBGwYTfnnZSbyU6XtGJdcl6gZgaUX9//9LsCy+wFS/p0rqeef7nf+SnNaDzpoO36fS4zGzgaElQkfQxYFBFvL0YeEQEcApwu6S7gRaC0quQpwFqSZgJfAu4re6/cCcDmwPbAMOAb3Vy/7dedOPnQDg7dej1OPjTF8NkLlnD8b+6vObA89fzLnDF9TiGBZdJB2/GFX82kc+5T/T6XmQ0MSp/tdb6IdDJwOCkwrASsDvwuIj5bts+ewNER8ekuxwr4M7BVRHS7yqSkXYD/iIiP9ZSXjo6O6OwcOM1AzTR7wRI2HzGs5uOeev5l1l/jPYXkoXPuU3Rsun4h5xoIJN0TER2971msVvy7HizDy9thPZ9a/q4bckd9RJxAqlWUf/h/VtK6EbEoN1t9A5iU91kTeDkiXgOOBm6tFFAkbRART+bAcwDwYCPKM1j0JaAAhQUUYFAFFFvWYAkcVr1m94B+TdIjwAPANRHxh5z+fuAhSbOBfUijxgCQdJ2k9+bNyyTNAmYB6wA/aFzWzcysq4bP/RURM4AZ+fXXgK9V2Od2YEw3x+9b9nq3umTSzMz6xBNKWsMMtv6TZpM0AZgAMGrUqKbmxc1cVtLs5i9rE51zn/JIr4JFxLkR0RERHcOHD292dswABxVrkI5N1+fsz4x1TcVskHNQsZr1tamjY9P1OX/Gn3rf0cxalvtUrCaju1mbvhrnz/gTJ10/B4Cjdnlf4Xmz2nT35aAd7ruw+nFQsYYpBRIHFLPBy81fVpPu1qWvlgOK2eDmoGI1c/OImXXHzV9mZi1soPWNuaZiZmaFcU3FzKzBBlrtokiuqZiZWWFcUzGzqnmOr8GtiN+vaypmZlYY11TMzPrBtbd3c03F+qXWNe7NbHBzULE+m71gCcf/5n4HFjN7m4OK9dnmI4bxk4O37vNa92Y2+DioWL84oJhZOQcVMzMrjIOKmZkVxkOKrc9ef/Mtlh/i7yUDwawFzxc29YeHyFp/KCKanYeGkrQY+EuVu68DPFPH7Aw0tZRXGrriivHG0qVAq/4R1eP3u1FEDC/4nBVJmgBMyJubAY824rpZM/83mv1/2Y5lHwPcHhF797Zj2wWVWkjqjIiOZuejUVxeq1Yzf3bN/r21c9mr4bYLMzMrjIOKmZkVxkGlZ+c2OwMN5vJatZr5s2v2762dy94r96mYmVlhXFMxM7PCOKiYmVlhHFSsTySp2XmwgUPSkGbnwQbG/6WDSj9J+mdJBzc7H40g6b2SRgJEG3TGtdPvtj8k7Q0c2azAIul9kjaXtGqTrr+PpAOadO2xkraXtDGk/0tJTf1cd1DpB0l7AacDTzQ7L/Um6RPADOB8SVdK2kbSik3OVt200++2P3JA+TEwNyLe7PJe3b81SzoQuBw4C5gk6Uv1vmaX6+8J/AR4rpHXzdf+OPAr4BvAdySdChARbzU1sESEH314ALsALwAfzNurAu9pdr7qVNa1gRuB7fL2/wN+DuwDrNDs/Pl327Sf07bAY8DBeXtNYBQwqmwf1fH67wH+APwTsEr+e/wlcGKDyv+PwFxgj7y9GmkaldUacO0VgGuA/fP2+4GpwH834mff08M1lb57E3gWWC9/Y58M/FLSVZLGNDdrhXsVWBEYDhARXwfmAPsD74OB0ZZboHb63fbHe4AHgRck7USqMfwX8GNJ34G6N5O+CSwB/h4RfyfVpH8GbCLpc3W8bskw0pePhfnvYjJwMXCapCPreeGIeI30P/hq3n4E+FdgHUk/yGlNaaJ2UOmjiPgf4LOkf6Kngd8DRwFPkqrDg0b+h/0VsJ2kf8hppwGvAN/N24OmjyX/bg+nDX63/RERt5Ganf4FuBC4kvRz+ymwlaQt63z9pcDtwNmSRkTEK8BDwFXAZpJWqPP1rwV+RLoh8SrgauBzwB+BnSWtX8/rAw8D50kanfPzDOn/cQNJG9T52t3yzY81kLQb8BFSlff+iJiV07aOiNPL9rse+JeIeKpJWe03SeOAccBTpGr2IuA7wP3ADRHxWN5vGqmsi5uV1yLk8n6MFDhuiIg7Je1K+t3+tGy/lv/dFkGSSl8kJO0CjIyIS8ve/y3wvYh4oAHX/y6wBfDvEbFA0urAFOBfI+LPDbj+J4D3RsSZeXs1Uq3tC/W6flk+vgUcAHwqIubltGuBEyJiVj2v3R2vp1Kl3CF3GqkvYXXgt5KOiog/SLq1bL/PkNqWX2lOTvsvN2X8BDgZWAu4ATgIuIj0rXR9SfcAKwMbAa81JaMF6VLeNYFrJR0dEVMG2++2SKUP1oiYIWloWfoBwEjgb/W8PO8sufAL4N+AKbmjfnNSP9jf63h9JC0XEW9FxO+6jHzbHVijntcv+9lPkhTA/0j6Mun/cRRNGDhQ4qBSvZ2AH0XEJQD5W+zVkvaLiP/N/1SHAhOBT0fE803Ma3+NAP4vIi4CkPRn4BzSP+63gX1J1fzXgcNbvKxQubwnS3ojIqZJWh44hMHxu+0TSf9M+oKxfERcUfYtvfTh9kbe/iJwNPDZiFhQ4PV3AdYFhkbEryLirdJ7EfEk8F1JS0jry6wLHBMRi+p4/QCirPxv5v2+TGoqPayo60vaEHix/O+uvLk5In4oaT6wNanD/rMR0bxRi80YHdBqD9K3oh8BvyhLm0hqO74T2JBUezka2LzZ+S2gvNsC/w1sUJZ2IKkJbGzeXglYpdl5rWN598/l3Y400uZo4P3NzmuTfj4fJS0AdjRptNX5wDrd7LtP0f8DwK6kZsn/ADqBs0nNTZX2XYEU+Jp1/cOALQu89gHAbOD4Sj9zYLku20Oa/ffiPpUeSNqZ9Ad6i6RhwF358RZpJNS+wBnAuRFxv6Shkb+xtRpJ+wArRsTVecTTJcAC0j9SRERI+jqwNCJ+1sy8FqHK8n4NeC0iftbKv9v+yLW0S4BpEfHLfIPhn4DrgK9HxJK8317AzIh4uuDrCzgVeDIiTpe0EimoPQP8sHQ9SR+PiGuKvHaN1/9YpI77Iq89nDSi7K/AfNKgkcmROuTL96vLz76vPPqrAiXrAucBl0raN//zbEca4fE7YFykKu+qpDZMWvVDp+sNXJFG1UwAtiHdALhR3nVFcllbWQ3lXan0ulV/t31VGiIeEa+Thq4iaeWIeInUCb0l8P2yQzpIQ4wLFelb772k0VzrRcSrwDHAesCJZbt+VNImTbz+nnW4/vPAV0hNzTNJS/oemj+bKLvBcUdS/+bA0Oyq0kB+AN8Cfkga8XRQhfePJDULjGp03gosY6UbuDbIr1cmNQtdBPyWNFzzA83Os8vbkJ/TumWvDySNAPwhaejwRfnndhO5ObQO1x9J+hJTGgxyGbAHsHLZ7+oe8s1/g+n6pC9uK9Dlhlvgk6SWkS/l7W2a/XdS6eGO+gryN4AgfVP9E/C/wEn5m8ibEfFjSTsAe5E6bv/avNz2W9cbuH4KIGkR6Q7dY0jfkDYFHok6D5FsgHYrb81y0+B/5BF+L5F+RgtJo6qWAGdFxCt5QMPLdbj+OFKT0/+RgtfxwK9J39olaVZEPCnpZqDwGmQzr9/l2mtIOjEiZgNExG9zBXJrSVcDu0vaPAocEFGIZke1gfjgnft3tgEm5tffJ929+v28vTywRrPzWlB5DyUFzgdJH6ojSR2OFwHDm50/l7ehP5stSTW5fwLGkqbk+QPpPhTIHcOku7fvBjYs8NrKv4tZpKly1gO+TupTGEG6j+iS/DiZ1M/wvsFw/W6u/VVSMN+yy76/BOaRpxEaaA/XVDJJW5GCyf2Rf3PZppL2Bw4mTQFxuKQ7IuI6UptnyykvK0BE/FrSUtKIlvPyPlNJMwasCrT6jY1tVd5+GgrcFmlWASSNAj4MXCLp8IiYn2vpXyLdBDq/qAvn/7snJN1OaiFYFBH/T9LrpG/uO5H6N7YnDZ/dPSL+NBiu3821T8vXvlHSrhHxJ6U75bcFDogm3dzYq2ZHtYHwIH0D+RvpG8CHu7z3s/zegXl7f+Afmp3nOpV1SNnrA0j/SOs2Mn8ub9N+TjuTanCrAI+T7k6H1I/yZdL0H5/OaasDaxV8/Y8D/05qAZgMfLPL+ycAF5BG7NWj/E27fhXX/jqpFl3qz1m12X8vPZan2Rlo9oPUIXY2MIl009Ivyj98gN14Z3bepsz62aiylu33ZdLghJbupG638vbxZ7QcqXb2EOkb8kHAZqQBKJOB6aTay+eBk+uUhz1Jo5v2ytujSU1O3yjbZzRpjq3C/webef1ml70ej7Zv/oqI1yT9J7AUWJ901/AR+b6EP0bEH5qbw+L0UNblIjd3ZE8Dn4mIh5qRz6K0W3n7ItKd6S9Jupg06+8BpHuzNpO0RuS7uCW9BQwt+n4dSf8IXAp8PCLukrQOqa/iAGBabv65ljRqb1vSNDnPDobrN7vs9dK2QUXSWNKHDZGmjQZ4TNIUUhPXZyQ9AfwD8EREPBr5a0OrqaKsn5W0gFTWv0bEb5qT02K0W3kL8gZpKOuFwARJ2wNLlSYsPAb4Gqkdv+jRVn8jTfezgaS1gStyXh4iDe/ejjQarwM4MiKK/lBt5vWbXfb6aHZVqRkP0lQSTwBnkkYAHdnl/U1Jo1s6SUMqxzQ7zw0q699buaztWN4Cf26b8M5Ix6+Shgqfmbe/DGxWx2tvTerHmU8KYMuRbkY9i3dGnRXahzNQrt/sstelTM3OQEMLm4btrUqaYmK/nLYTaQjl57rs+0PSsL3C5vFxWV3egfoA3kuqpRxDuoP+O6Sml8826PpbAMd2SbsB2Lb0+x2s12922Yt+tFXzV6Tf0EuSOoHVJS0fEXdIOgS4QtKrEXGR0jTWmwOfiBZtZ2+nskL7lbdoEbEwN/d+m/QBd43STNxzG3T9h0mLTgEg6ZOkpXkX5Pfr2vTczOs3u+yFa3ZUa8YD+ALpW9nqZWk7k5pENm12/lxWl7dJP7uR5JGOeXu5JuRBpObJh2lCTbKZ12922Yt6tNWEkmWT5J1Nmvzu55LWyN9qbwMeoA7TPjRDO5UV2q+89RART0TEPWU/y7d6O6ZOHqe5NclmXr/ZZe+3QT/1vaTNSPM9dQJvRV5MJ783mbSK3x2kkXDHA/8cBd4l3EjtVFZov/KatYJBHVSU1o7+IaltcgHpw+eiiHihbJ9/JXVSbg18t1W/IbRTWaH9ymvWKgZtUFFaXOiXwBmRlvv9JGk00FLSssDPd9l/xUjrarScdiortF95zVrJYO9TWZ108xDAVaQhkiuQ5jhC0g6Sts3vv9b47BWqncoK7Vdes5YwaINKpBXrfgJ8QtI/5U7H20jz7HxE0sqk2VcX5v1btsrWTmWF9iuvWSsZtM1fAErrSR8NbAX8MiJuzekzgKMi4rEmZq9Q7VRWaL/yDmaS3iStIzIUeAQYHxGFL/5VK0nfjIgf9vMcPyLNQvwa8BhphofnisjfQDWob36MiFclXUZaxfEESZuT2t2Hk6ZfGTTaqazQfuUd5F6JiLEA+Xf6OVJNtFeShpSP+ivYN0mDQapWIT/TgRMi4g1Jp5Km0P9GgXkccAZt81dJpEnYziOtYLcbsCtp6omnm5qxOminskL7lbdN/A9pfjYkXS3pHkkPSZpQ2kHSS5K+L+lO4EOSviPpbkkPSjq3dJ+NpBmSTpd0q6RHJG0v6XeS5kj6Qdn5PivpLkkzJf1C0hBJpwAr57TLutuvUn7KCxMRN8Y7k3DeAWxYvx/dANHsuy8b+QCG0IS7hF1Wl9ePHn93L+XnocAU4PN5e1h+Xpk0OejaeTvIC4aV75dfl6aSB5gBnJpfH0fqY9sAWJE0gePawPuBa0jT/UNaf+eI8nzl1z3t96789FDOa2jQXGrNfAzq5q+uon7V5AGnncoK7VfeQWZlSTPz6/8Bzs+vvyzpwPx6JGm0399I6778tuz4XSV9nTSTwjDS1PHX5Pem5udZwEMR8SSApMfzOXcmTTF/d67grAwsqpDH3XvYr2t+lqG0hMAbwGU97TcYtFVQMbMB6e0+lRJJuwAfBT4UES/nARgr5bdfLX2JyAM2zgY6IuIJSd8t2w/yujrAW2WvS9tDSfNtXRwRJ/SSx572e7WnLzWSxpOWtd49cpVlMBv0fSpm1pLWAJ7NAWVz0s2tlZQCyDOSViUth1yLm4GDJK0LIGmYpI3ye6/nG217269bkvYmdczvFwNgRFsjuKZiZgPR9cDnJD0APErq5F5GRDwn6TxS89Y84O5aLhIRDystOX2jpOVIKzEeC/yFtC78A5LujYjDetivJ2eS+nCm52azOyLic7XksdUM6vtUzMyssdz81SIkvZmHMj4o6QpJ72l2niDdIFbAOT6Vh42+JamjiHyZWXM4qLSOVyJibER8gHR3btVV6NJ4+jqpOahUyM+DwCeAWwvJkZk1jYNKaxpsN4g9EhGP1vuHZmb156DSYiQNBfYhdUwC/GtEbAd0kMb1r53TVwEejIgdI618eGZEbJ9rOiuThjiWvBYRHwF+Trr57FjgA8C/SFpb0vuBg4EP56GfbwKHRcRE3qlBHdbdft3kx8wGIY/+ah2D/gYxM2t9DiqtY1DfIGZmg4Obv1rboLhBzMwGDweV1nY9MDTfIHYSPdwgRprNdxZwNX24QQwo3fj1AGk67w3y26UbxC7rZb9uSTpQ0nxSB/40STfUkj8zGzh886OZmRXGNRUzMyuMg4qZmRXGQcXMzArjoGJmZoVxUDEzs8I4qJiZWWEcVMzMrDAOKmZmVpj/D2waT1Sn+IV6AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from __future__ import print_function\n", "import pints\n", "import pints.toy as toy\n", "import pints.plot\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# Load a forward model\n", "model = toy.LogisticModel()\n", "\n", "# Create some toy data\n", "real_parameters = np.array([0.015, 500])\n", "times = np.linspace(0, 1000, 1000)\n", "org_values = model.simulate(real_parameters, times)\n", "\n", "# Add noise\n", "noise = 10\n", "values = org_values + np.random.normal(0, noise, org_values.shape)\n", "\n", "# Create an object with links to the model and time series\n", "problem = pints.SingleOutputProblem(model, times, values)\n", "\n", "# Create a log-likelihood function\n", "log_likelihood = pints.KnownNoiseLogLikelihood(problem, noise)\n", "\n", "# Create a uniform prior over the parameters\n", "log_prior = pints.UniformLogPrior(\n", " [0.01, 400],\n", " [0.02, 600]\n", ")\n", "\n", "# Create a posterior log-likelihood (log(likelihood * prior))\n", "log_posterior = pints.LogPosterior(log_likelihood, log_prior)\n", "\n", "# Choose starting points for 3 mcmc chains\n", "xs = [\n", " real_parameters * 1.1,\n", " real_parameters * 0.9,\n", " real_parameters * 1.15,\n", "]\n", "\n", "# Choose a covariance matrix for the proposal step\n", "sigma0 = np.abs(real_parameters) * 5e-4\n", "\n", "# Create mcmc routine\n", "mcmc = pints.MCMCSampling(log_posterior, 3, xs, sigma0=sigma0, method=pints.MetropolisRandomWalkMCMC)\n", "\n", "# Add stopping criterion\n", "mcmc.set_max_iterations(30000)\n", "\n", "# Disable logging mode\n", "mcmc.set_log_to_screen(False)\n", "\n", "# Run!\n", "print('Running...')\n", "chains = mcmc.run()\n", "print('Done!')\n", "\n", "# Show traces and histograms\n", "pints.plot.trace(chains)\n", "\n", "# Discard warm up\n", "chains = chains[:, 10000:, :]\n", "\n", "# Check convergence using rhat criterion\n", "print('R-hat:')\n", "print(pints.rhat_all_params(chains))\n", "\n", "# Look at distribution in chain 0\n", "pints.plot.pairwise(chains[0])\n", "\n", "# Show graphs\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }