{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Think Bayes\n", "\n", "This notebook presents example code and exercise solutions for Think Bayes.\n", "\n", "Copyright 2016 Allen B. Downey\n", "\n", "MIT License: https://opensource.org/licenses/MIT" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "# Configure Jupyter so figures appear in the notebook\n", "%matplotlib inline\n", "\n", "# Configure Jupyter to display the assigned value after an assignment\n", "%config InteractiveShell.ast_node_interactivity='last_expr_or_assign'\n", "\n", "# import classes from thinkbayes2\n", "from thinkbayes2 import Pmf, Beta, MakeBinomialPmf\n", "\n", "import thinkplot\n", "\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VfWd+P/XOysBQtgCsgfZJOwQ9h1kERXU4hSrVVtbZ2ztTOu037bTab/+nOl82860nS6Ora1trVrXqQqyCYLsIGHfISyBsEjYA2TP+/fHOZzcpNkIuffc3Pt+Ph55eD7nfs7J+x7Jfd/POZ9FVBVjjDEm3MT4HYAxxhhTFUtQxhhjwpIlKGOMMWHJEpQxxpiwZAnKGGNMWLIEZYwxJixZgjLGGBOWLEEZY4wJS5agjDHGhKU4vwMIhbZt22paWprfYRhjjAG2bNlyTlVTa6sXFQkqLS2NzMxMv8MwxhgDiEh2XerZLT5jjDFhyRKUMcaYsBTUBCUiM0XkgIhkich3qng9UUTedF/fJCJp7v5pIrJFRHa5/50ScMzH7jm3uz/tgvkejDHG+CNoz6BEJBZ4HpgG5ACbRWS+qu4NqPYEcFFVe4rIPODHwGeBc8C9qnpKRPoDS4FOAcc9rKr2UMkYYyJYMFtQI4AsVT2iqkXAG8CcSnXmAC+72+8AU0VEVHWbqp5y9+8BmohIYhBjNcYYE2aCmaA6AScCyjlUbAVVqKOqJcBloE2lOp8BtqlqYcC+P7q3974vIlLVLxeRJ0UkU0Qyc3Nzb+V9GGOM8UEwE1RViaPy8r011hGRfji3/f4+4PWHVXUAMN79+XxVv1xVX1TVDFXNSE2ttbu9MWFNVSksKsZWwDbRJJjjoHKALgHlzsCpaurkiEgckAJcABCRzsC7wKOqevjGAap60v1vnoj8BedW4p+D9SaMCbXi4lJ2HTrJ5t3HyD51gQuXrnHhyjVKS8tIiI+jTctmtE5pRu9u7RgxsDu9urWjmhsJxjRqwUxQm4FeItIdOAnMAz5Xqc584DFgAzAXWKGqKiItgYXAd1V13Y3KbhJrqarnRCQeuAdYHsT3YEzIHM05x7sfbWfLnmwKCourrFNUXMLp3Muczr3MnqxTvPvRdlomN2XMkNuZM2UwbVs1D3HUxgSPBPOWgYjMAv4biAX+oKo/FJHngExVnS8iTYBXgCE4Lad5qnpERP4V+C5wKOB004FrwGog3j3ncuAZVS2tKY6MjAy1mSRMuDpx5iJvLNrMxh1HaqwXGxtDaWlZta/HxcUyY2w69985hFYtmjZ0mMY0GBHZoqoZtdaLhnvalqBMOCoqLuG1BZ+wcNXOv3k4e1vbFgzvn8bgvl1o1zqZ1inNSEyI43pBERcuX+fU2Utk7s4mc082V67mVzg2IT6OR+eMYua4fnbrz4QlS1ABLEGZcHP89AV+/vJyjp++UGH/yIHdmTt9KN07t61TcikrK2P7/hzeXJxJ1vGzFV4blt6Nr35uEinJSQ0auzG3yhJUAEtQJpwsW7+Xl/53HcUl5XemB/buzCP3jqRH1/r1OFVVtuw9zqvzN3LizEVvf0pyEt949E4G9K48wsMY/1iCCmAJyoQDVeWNRZt558Ot3r74uFgev28MM8alN8jtuOLiUl5dsIkPVu309sXGxvBPn5/K2CE9bvn8xjSEuiaoqFhuwxi/lZWV8eLba1i2fp+3r2uH1jzz+DS63NaqwX5PfHwsX3hgDIP7duFXr63gcl4+paVl/PxPy7hyNZ+7xvdvsN9lTLDZbObGBFlJSSk//eOyCslpaHpXfvTM/Q2anAIN6duFn/zzZ+jUriXgjH7//TtreXOx3UkwjYclKGOCSFX55Wsr2bjzqLdvQkYvvv3EDBIT4oP6u9u2as4Pv34fvbqVT/j/1pJM5q/cEdTfa0xDsQRlTJCoKn98dz3rtmZ5++6ZOJB/fGQKcXGxIYkhuVkTnv3qvQy+o3xSl5ff28CqzQdD8vuNuRWWoIwJkneXb2fhql1eefrYdB6/f3TIxyY1SYzn21+awR233+bt+/VfPmbr3uMhjcOYm2UJypggWJN5iNc+2OSVRw3szpfnjvNt4GxCfBzf/fJddOnQGnA6bfzkpaUczTnnSzzG1IUlKGMaWPap8zz/+sdeuV/PjvzTo1OJifH3z61500R+8NTdpLZKBqC4pJSfvLSUvGsFvsZlTHUsQRnTgK7lF/KTl5Z6g3A7tWvJt780g4T48BjR0TqlGf/61CyaJDodNM5eyOMXr3xky3iYsGQJypgGoqr88pUVnDl3BYDEhHi+9cQMmiWF12LQndu34msPT/bK2/ad4K0lW3yMyJiqWYIypoH8dfk2Mvdke+Wvfm5S0MY53apRg27n/qmDvfLbSzLZtu9EDUcYE3qWoIxpAIeP5/LGovJBsLMnDwr7qYUeunsE/Xt1BJyBvL9+baU9jzJhxRKUMbeoqLiEX766grIyZ62m3mnteeTekT5HVbvY2BieeWyaN9v5pbzrvPj2Gp+jMqacJShjbtHrCzeT86kzg3hCfBz/+MgUYmMbx59WSnIST82b6JXXbzvM2i1ZNRxhTOg0jr8iY8LUnqxTLAiYOujx+0bTITXFx4hu3vD+aUwZeYdXfvHtNZy/dNXHiIxxWIIypp4KCov51asrvdVwB9/Rhelj032Nqb6+cP8Yb3zUtfxCXnhjlXU9N76zBGVMPb29dAu5F/MAaJaUyFcemthol1hvmpTA1x6ZzI3ot+07wfrtR3yNyRhLUMbUQ/ap88xfWb4o4OP3jaZNy+Y+RnTr+vXsyIxx/bzyH/+6jmv5hT5GZKKdJShjbpKq8tu31ni99tJ7dGDyyD4+R9UwPnfPCFomNwXg4pXrvLFos88RmWhmCcqYm7R8wz4OHD0DOF21n/y7CY321l5lzZIS+cIDY7zy4tW7yco+62NEJppZgjLmJlzOy+eV+eWzlN83ZXDYzhZRX2OH9PDWj1LgN2+t9lqLxoSSJShjbsIbizd7z2Xat2nB3BlDfY6o4YkIX5o7jnh3UcWjOedY+ckBn6My0cgSlDF1lH3qAsvW7fXKX/zM2LCZpbyhdUhN4f47h3jl1z74hPyCIh8jMtHIEpQxdaCqvPzeem/M08DenRmW3tXXmIJtzpRBtE5pBji3Nt9dvt3niEy0sQRlTB1s23eCHQdyABDg8fvHREzHiOo0SYyvMKfg+yt3cPZCno8RmWhjCcqYWpSUlPLyexu88p1j+tKtY2sfIwqdCRm96NElFXCuw6sLNtVyhDENxxKUMbVYvmG/Nxlsk8R45s0a7nNEoSMiPH5/ebfzdVuzOHjsUx8jMtHEEpQxNSgoLOatpeXrPH1m2lBvIGu0SO/RgVGDbvfKr31grSgTGpagjKnBotW7uZyXD0Cbls24Z9IAnyPyxyP3jiTGfea2+9Ap73mcMcFkCcqYalzLL+Td5du88oMzhkVst/LadEhNYero8iU5XluwyWY7N0FnCcqYarz/0Q6uu2N/OqSmMHlEZMy3V18PzhhGnDt49/CJXDbtPOpzRCbSWYIypgqX8q6z4OPy2crn3TXc+3COVm1aNmfW+P5e+fWFm20KJBNUQU1QIjJTRA6ISJaIfKeK1xNF5E339U0ikubunyYiW0Rkl/vfKQHHDHP3Z4nILyXSB6MYX/zvh1spKi4BoFvHNowd2sPniMLD/XcOpkliPAA5n15kdeYhnyMykSxoCUpEYoHngbuAdOAhEam83OgTwEVV7Qn8HPixu/8ccK+qDgAeA14JOOYF4Emgl/szM1jvwUSn85eusjRgSqPP3TMi4gfl1lWL5knMnjzIK7+9dAulpdaKMsERzBbUCCBLVY+oahHwBjCnUp05wMvu9jvAVBERVd2mqqfc/XuAJm5rqwPQQlU3qPOE9s/AfUF8DyYKvb9ih/eh26tbu4if0uhm3TtpIE2bJABw5twV1m7N8jkiE6mCmaA6AScCyjnuvirrqGoJcBloU6nOZ4Btqlro1g/s31rVOQEQkSdFJFNEMnNzc+v9Jkx0uXD5Gh8GtJ4enDHMWk+VNE1K4N7JA73yO0u32LMoExTBTFBV/VVX7pdaYx0R6Ydz2+/vb+Kczk7VF1U1Q1UzUlNT6xCuMTB/xQ6KS0oB6NEllaHWeqrS3RMHkOS2ok7lXmb9tiM+R2QiUTATVA7QJaDcGThVXR0RiQNSgAtuuTPwLvCoqh4OqN+5lnMaUy+X8/JZsnaPV35wprWeqtMsKZG7J5YPWn576RYbF2UaXDAT1Gagl4h0F5EEYB4wv1Kd+TidIADmAitUVUWkJbAQ+K6qrrtRWVVPA3kiMsrtvfco8H4Q34OJIu+v2O61ntI6tSWjXzefIwpv90wcUKFH3/rt1ooyDStoCcp9pvQ0sBTYB7ylqntE5DkRme1WewloIyJZwDPAja7oTwM9ge+LyHb3p5372lPA74Es4DCwOFjvwUSPvGsFLFkb+OxpqLWeapHcrAl3TyhvRb1jrSjTwII6b4uqLgIWVdr3g4DtAuDBKo77d+DfqzlnJtC/qteMqa/Fa3ZTWFQMQNcOrRk5sLvPETUO90wawIKPd1JUXMLx0xfYuvc4w6zlaRqIzSRhol5BYTELV+3yyg/cOcRaT3XUonkS08eUD2+0VXdNQ7IEZaLeRxv3c/V6IQCprZIZM8RmjbgZ904eSEyM81Gy78hp9h8543NEJlJYgjJRraSklPkrd3jlOVMHERtrfxY3o22r5kzI6OWVA2eAN+ZW2F+iiWrrth3m3MWrgHO7asrI6J6xvL7umzrY287ck83x0xd8jMZECktQJmqpaoVv+7Mm9CcxId7HiBqvLre1Ynj/NK/83kf2LMrcOktQJmpt2XucE2cuApCYEM9d461z6K14YNoQb3tN5iGvZWpMfVmCMlFr/oryZ0/Tx/SledNEH6Np/HqntSe9RwcAylQr9Iw0pj4sQZmodORELnuynFmyYkQqTNtj6u/egKU4lm3YR767IrEx9WEJykSl+SvLV8sdPaQHqa2TfYwmcgzv342OqSkA5BcUsXzDfp8jMo2ZJSgTdc5dvMq6bYe98uxJA2uobW6GiHBPwPVcuGqXLWho6s0SlIk6i9fs9tYvSu/RgZ7d2tVyhLkZk0b09p7n5V7MY+POoz5HZBorS1AmquQXFFVYkDDwmYlpGIkJ8cwc188rL1i5wyaRNfViCcpElRWbDnDdfXDfITWF4f1tYtNgmDm+vzcjx6Hssxw4+qnPEZnGyBKUiRqqyqLV5V2f75k40CaFDZJWLZoyMaO3V/7AupyberAEZaLGlr3HOXPuCgBNmyQwaUTvWo4wt+KeSeVd9zftOGIDd81NswRlosbCj8u/xd85uq+3GqwJjm4d29CvZ0fAGbi7ZM1unyMyjY0lKBMVTpy5yM6DOQAIcNcEm9YoFAIHQC/bsI+i4hIfozGNjSUoExUWriofmDtyYHfa2cDckBjev5t3ra9eL2R15iGfIzKNiSUoE/HyrhXw8ScHvfIsm9YoZGJiYiq0Vj9Ytcu6nJs6swRlIt5HG/dTXFIKOM9FbkxoakJj6qg7vGVMTpy+wO5Dp3yOyDQWlqBMRCsrK2PJmj1e+Z6JA6xreYg1S0pk0vDyHpOBXf2NqYklKBPRMvccJ/diHgDNmyYyblhPnyOKTrMmlt/m27zrGLkX8nyMxjQWlqBMRFu8urxr87TRfUmIj/MxmujVuX0rBvbuDIBChemmjKmOJSgTsXI+rdi1fHrA/HAm9AI7S1iXc1MXlqBMxFq6tvzZU0b/NOta7rOMfl1p26o54PSsXB+w5IkxVbEEZSJSfkERKwO6ltvAXP/FxMQwY2x5K3ZxQOcVY6piCcpEpFWbD3nLjXdMTWFg704+R2TA6XJ+Y5bzrONnOZRts5yb6lmCMhFHVVkcMO/bzPH9rWt5mEhJTmLc0PKelNaKMjWxBGUizt7Dp8n59CLgLJ5ns5aHl1njy2+3rtt2mCtX832MxoQzS1Am4gR+K584vBfNkhJ9jMZU1rNbO3p0SQWgpKSUFZsO+ByRCVeWoExEuXD5Gpt2HvXKM61reVgK/P/y4bq9Nj+fqZIlKBNRPtq4n7KyMgDuuP02unVs43NEpipjh/bwWrafnr/Ctn0nfI7IhCNLUCZilJaWsWx9+QwFM8da6ylcJSbEM3lEH68cOGbNmBssQZmIkbknm/OXrgHQonkSowbd7nNEpibTx6V721v2ZHPW5uczlQQ1QYnITBE5ICJZIvKdKl5PFJE33dc3iUiau7+NiKwUkasi8utKx3zsnnO7+9MumO/BNB6B38LvHHUH8fGxPkZjatOpXcsK8/MtX7/P34BM2AlaghKRWOB54C4gHXhIRNIrVXsCuKiqPYGfAz929xcA3we+Wc3pH1bVwe7P2YaP3jQ2p85eYseB8nn3po2t/E/NhKMZAa2oZRv2UVxc6mM0JtwEswU1AshS1SOqWgS8AcypVGcO8LK7/Q4wVUREVa+p6lqcRGVMrZYFfPsemt7N5t1rJIb3T6N1SjMArlzNZ9Ouo7UcYaJJMBNUJyCwa06Ou6/KOqpaAlwG6tLt6o/u7b3vSzVTBIjIkyKSKSKZubm5Nx+9aTSKiktYsWm/Vw78Vm7CW2xsDNPG9PXK1lnCBApmgqoqcVQe7FCXOpU9rKoDgPHuz+erqqSqL6pqhqpmpKam1hqsabw2bD/C1euFAKS2SmZI3y4+R2RuxtRRdxDjfs/ce/g0J85c9DkiEy6CmaBygMBPis7AqerqiEgckAJcqOmkqnrS/W8e8BecW4kmii0NWPzuzjF9iYmxzqmNSZuWzRk+IM0rBw4VMNEtmH/Jm4FeItJdRBKAecD8SnXmA4+523OBFVrDkHIRiRORtu52PHAPsLu6+ibyZZ86z4GjZwBnOYc7R9/hc0SmPqYHdGpZuekAhUXFPkZjwkXQEpT7TOlpYCmwD3hLVfeIyHMiMtut9hLQRkSygGcAryu6iBwDfgY8LiI5bg/ARGCpiOwEtgMngd8F6z2Y8Ld0bfm37ZEDu9MyuamP0Zj6GtSnM7e1bQHA9YIi1m21xQwNxAXz5Kq6CFhUad8PArYLgAerOTatmtMOa6j4TONWUFjMqszyRQlnWueIRktEmDYmnVfmbwSc27ZTRllrONrZzXrTaK3ZcoiCQudWUMfUFPr17OhzROZWTB7Rp8JihkdOWO/baGcJyjRKqsqSgNt708f2s0UJG7mU5CTGDO7hlT+0zhJRzxKUaZQOH8/l2MlzAMTHxTJ5ZJ9ajjCNQWBnidWZWVzPL/IxGuM3S1CmUQrsWj52aE+aN7VFCSNB39tvo3P7VgAUFhWzZsshnyMyfrIEZRqda/mFrN2a5ZVn2Lx7EUNEKrSiltpihlHNEpRpdFZtPkhRcQkAXTu0plc3m9A+kkwa0Zv4OGcm+uxT5zmUbfNBRytLUKZRUVU+DLi9N8M6R0ScZkmJjBvW0ysH3s410cUSlGlUDhz91JurLTEhngkZvXyOyARD4G3bdVuzvLkWTXSxBGUalaXryme7Hj+sJ02TEnyMxgRLz67tSOvUFoDiklJWbjrgc0TGD5agTKORd62A9duPeOXpY6xzRKQSkQqtqGXrrbNENKoxQYnIhwHb3w1+OMZUb+UnBygpcVZc7dEllR5dbRmVSDZ+WE8SE+IBOHn2EnsPn/Y5IhNqtbWgAj8Bqpwzz5hQ+JvOETbvXsRLapLAxOHlzxits0T0qS1BWZvahIXdh05xOvcy4HxwjR3Ss5YjTCQIvM23cccRLufl+xiNCbXaEtTtIjJfRBYEbHs/oQjQGKj47XnS8N40SYz3MRoTKmmd2nrj3EpLy1j5iXWWiCa1LbcxJ2D7v4IZiDHVuXjlOpt2HvXK06xzRFSZMbafN1h32fq9zJkyyMa+RYkaE5SqrrqxLSKp7j6bA9+E1IpN+ykrKwOgT/fb6Naxtc8RmVAaM+R2/vDXdVwvKOLMuSvsPHiSQX06+x2WCYHaevGJiPxfETkH7AcOikiuiPygpuOMaShlZWUsW7fPK9u8e9EnMSG+wmz1H67dU0NtE0lqewb1dWAcMFxV26hqK2AkMFZEvhH06EzU27bvBLkX8wBo3jSR0YNv9zki44fA27qf7DrGhcvXfIzGhEptCepR4CFV9R4AqOoR4BH3NWOCKrBr+ZSRd5AQX9tjUxOJutzWivQeHQAoU2X5hn21HGEiQW0JKl5Vz1Xe6T6Hsm5UJqhyL+SxZU+2V542pq+P0Ri/zRjbz9tevmEfpaVlPkZjQqG2BFXTcpa21KUJquUb9nkD8Qb27kzHdi19jcf4a9Sg7rRongTA+UvX2LL3uM8RmWCrLUENEpErIpLn/ly5UQYGhCJAE51KSkpZvmG/V55unSOiXlxcLFMDO0uss84Ska7GBKWqsaraQlWT3Z8WAWW7xWeC5pPdx7iUdx2AVi2aMrx/N58jMuFg2th0boyA2r7vBGfOXfE1HhNctXUzbyIiXxeRX4vIkyJiT6hNSCwN6Eo8dXRf4twVVk10a9+mBYP7dgGcediWrbf5+SJZbbf4XgYygF3ALOCnQY/IRL2TZy+x+9ApAASYNto6R5hyM8aVd5b4aON+iotLfYzGBFNtCSpdVR9R1d8Cc4HxIYjJRLkP15Z/K87on0bbVs19jMaEm2HpXWnTshngrBG2YcdhnyMywVJbgiq+saGqJUGOxRgKi4pZsam8c0Tgt2VjAGJiYpge0OV8yVq7zRep6tqL70bPvYEBvfrs6aRpcGu3ZnG9wBnBcFvbFgy+w+ZcM39r6qg7iI11Pr4OHD3DsZN/M1zTRIC69uK70XMvLmC7RaiCNNFBVVm8prxzxPSx/WzWalOlVi2aMmpQ+bRXS2x+vohUWwvKmJDJOn6WoznON+H4uFimBIx5MaaywImDV2dmcS2/0MdoTDBYgjJhI/BZwtihPUlu1sTHaEy4S+/RgS63tQKcZ5erNh/0OSLT0CxBmbCQd62AtVuzvPJd1jnC1EJEmDmuv1deunYvqlrDEaaxsQRlwsJHG/dTUuKMZ7m9Syo93WW+janJhIxeJCY4k9rkfHrRGz9nIkNQE5SIzBSRAyKSJSLfqeL1RBF50319k4ikufvbiMhKEbkqIr+udMwwEdnlHvNLsafojV5ZWVmFmSNmjrN590zdNE1KYNLw3l55yZrdPkZjGlrQEpSIxALPA3cB6cBDIlL5k+cJ4KKq9gR+DvzY3V8AfB/4ZhWnfgF4Eujl/sxs+OhNKG3dd4KzF8oXJRw/rJfPEZnGZOb48tvBn+w6xrmLV32MxjSkYLagRgBZqnpEVYuAN4A5lerMwZlOCeAdYKqIiKpeU9W1OInKIyIdgBaqukGdm81/Bu4L4nswIbB4dfm33qmjbFFCc3O6dmhN/14dAWcxw8BFLk3jFswE1Qk4EVDOcfdVWcedqeIy0KaWc+bUck7TiJw8e4nt+51/JoLNHGHqJ7CzxLIN+2x+vggRzARV1bOhyl1s6lKnXvXd2dczRSQzNze3hlMaPwU+exrWrxvt29j4b3PzRgxI8+bnu3I1n/XbbX6+SBDMBJUDdAkodwYqd7Hx6rhLeaQAF2o5Z+DcN1WdEwBVfVFVM1Q1IzU19SZDN6FQUFjMik0HvPJdE/rXUNuY6sXGVpyfb7F1logIwUxQm4FeItJdRBKAecD8SnXmA4+523OBFVrDQAZVPQ3kicgot/feo8D7DR+6CYVVmw+S78671zE1hUF9bN49U3/TRvf15uc7lH2WrOyzPkdkblXQEpT7TOlpYCmwD3hLVfeIyHMiMtut9hLQRkSygGcAryu6iBwDfgY8LiI5AT0AnwJ+D2QBh4HFwXoPJniceffKv+XOHN/f5t0ztyQlOYmxQ3p45UXWimr0gtpdSlUXAYsq7ftBwHYB8GA1x6ZVsz8TsHtBjdzOgyc5ceYiAIkJ8Uwa0buWI4yp3d0TBrA68xDgzIz/6JxRtExu6nNUpr5sJgnji4Uf7/K2p4zsQ7OkRB+jMZGiZ7d29HJnISktLWOprRXVqFmCMiF3OvcyW/dme2XrHGEa0j0TB3rbS9ftsS7njZglKBNyi1bv8sYGDE3vSqd2LX2Nx0SWUYO60zrF6XJ+Oc+6nDdmlqBMSF3PL6rQtfzuiQN8jMZEori42AoDvj9YtctmOW+kLEGZkFqxaT8FhcUAdG7fyrqWm6CYPqYv8XGxABw5kcuBo5/6HJGpD0tQJmTKysoqdC2fNcG6lpvgaNE8qcKkwx+s2lVDbROuLEGZkNm8O5sz564A0LRJAhOHW9dyEzz3TCq/fbxx+2FvxnzTeFiCMiEzf+UOb3vG2HSaJMb7GI2JdN06tmFAb2cuaaXi0AbTOFiCMiFxKPtT9h85Azjzps2yzhEmBGZPHuRtL9uwj2v5hT5GY26WJSgTEvNX7vS2xw7p4XUDNiaYhvTtQuf2rQAoLCpm+Yb9PkdkboYlKBN0Zy/ksXH7Ea88Z8qgGmob03BEhHsnlw/cXbhqJyUlNnC3sbAEZYJu4ce7KHPHoQzo3Ym0Tm19jshEkwkZvWjRPAmA85eusSHgy5IJb5agTFBdyy9k+cZ9XvneSQNrqG1Mw0uIj2NmwMDd91fusIG7jYQlKBNUH67bW2Fg7tD0rj5HZKLRzHH9vIG7R3POsftQleucmjBjCcoETVFxCR8EdO2dPWWgDcw1vkhJTmLyyD5e+d3l23yMxtSVJSgTNKs2H+RS3nUAWrVoyoRhNjDX+Gf25EHc+Hq040AOR07k+hqPqZ0lKBMUZWVlvL+ifGDuvZMHER8f62NEJtp1SE1h1ODyFXff/Wi7j9GYurAEZYJi085jnM69DDjTGk0b3dfniIyBB+4c7G1v2HbY+zdqwpMlKNPgVLXCPf67xvenaVKCjxEZ47i9SyoDezsz6CsVp98y4ccSlGlwuw+d4rB7fz8uLpZZE23FXBM+Hpg2xNtesemA95zUhB9LUKbBvfPhFm97ysg+tExu6mM0xlTUv1dHenRJBaCkpJQFAdNwmfBiCco0qP1HznhjTGJEmDNlcC1HGBNaIlKhFbV4zR7yrhX4GJGpjiWDV3JxAAAX6ElEQVQo06ACW08ThvfmtrYtfIzGmKqNHNidLh1aA84ksh98bK2ocGQJyjSYrOyzbNt3AgCh4r1+Y8KJiDB3+lCvvHD1bluKIwxZgjIN5p0Pt3rbY4b2pFO7lj5GY0zNxgy+nY6pKQDkFxSx0JaFDzuWoEyDOHbyHJt3H/PKgd9OjQlHMTExzJ0xzCsvXLWL/IIiHyMylVmCMg3i7SXlz55GDexOV/f+vjHhbNzQnt5z0qvXC1m0ZrfPEZlAlqDMLTuac46NO4965cBvpcaEs9jYmArPSt//aIc9iwojlqDMLXt94WZve9TA7nTvbAsSmsZjYkZ5b9Nr+YUssB59YcMSlLklB46eYcvebMDpuffZWcP9DciYmxQXF8vfzczwygtW7rRxUWHCEpS5Ja8vKm89jRvWy549mUZp/LCedG7fCoCCwmLet5nOw4IlKFNvuw6eZNfBk4Aza8TfzbRnT6ZxiomJ4bOzyltRH6zaxcUrNkef3yxBmXpR1Qqtp0kj+tDRxj2ZRmz0oNtJ6+Q8Py0uKeWvy7bWcoQJNktQpl42787mwNEzgNMT6kFrPZlGTkR46O7yZ6hL1+3lzLkrPkZkgpqgRGSmiBwQkSwR+U4VryeKyJvu65tEJC3gte+6+w+IyIyA/cdEZJeIbBeRzGDGb6pWUlLKK+9v8MozxqbTrnWyjxEZ0zCGpXelT/fbACgtLeO1Dzb5HFF0C1qCEpFY4HngLiAdeEhE0itVewK4qKo9gZ8DP3aPTQfmAf2AmcD/uOe7YbKqDlbVDEzILd+wn1PuSqRJTRJ40MY9mQghIjw2Z5RXXr/tMIeyP/UxougWzBbUCCBLVY+oahHwBjCnUp05wMvu9jvAVBERd/8bqlqoqkeBLPd8xmf5BUW8uaS84frAnUNo0TzJx4iMaVh9ut/GqIHdvfKf39+IqvoYUfQKZoLqBJwIKOe4+6qso6olwGWgTS3HKvChiGwRkSer++Ui8qSIZIpIZm5u7i29EVPuvRU7uHI1H4A2LZtxz6QBPkdkTMN7+N6RxMQ4H497D58mc0+2zxFFp2AmKKliX+WvIdXVqenYsao6FOfW4VdFZEJVv1xVX1TVDFXNSE1NrWvMpgYXLl9j/oodXvlzd48gIT7Ox4iMCY6O7Voyc1z5E4lX3t9ISUmpjxFFp2AmqBygS0C5M3CqujoiEgekABdqOlZVb/z3LPAudusvZP78/kaKiksA6NaxDRMyevkckTHBM3f6MJKaJABw8uwllqzd43NE0SeYCWoz0EtEuotIAk6nh/mV6swHHnO35wIr1LnZOx+Y5/by6w70Aj4RkWYikgwgIs2A6YBNPxwC+w6fZs2WQ175C/eP8W6BGBOJUpKTKiwb8+biTC7n5fsYUfQJ2ieM+0zpaWApsA94S1X3iMhzIjLbrfYS0EZEsoBngO+4x+4B3gL2AkuAr6pqKdAeWCsiO4BPgIWquiRY78E4ysrK+P3/rvPKowf3YEDvyo8TjYk890wc4C1qeL2giFcXWLfzUJJo6J2SkZGhmZk2ZKq+lqzZw+/eWQNAfFwsv/rePFJt3JOJElv3HueHv13klX/0zP306tbex4gaPxHZUpdhQnaPxtQo71oBry/6xCs/MG2IJScTVYamd2V4/zSv/Lu311q38xCxBGVq9Mr8jVy97izg1q51MvdNHexzRMaE3uP3jyEuzpkr4PCJXD5ct9fniKKDJShTrd2HTvLRxv1e+QsPjLVu5SYq3da2BfdNGeSVX1mwifOXrvoYUXSwBGWqVFRcwm/eXO2VRw7szogBaf4FZIzPPjN9KB3cDhP5BUW8FNBxyASHJShTpXeWbuV0wHx7X5o7zueIjPFXQnwc//DZ8nkBNu08yobtR3yMKPJZgjJ/I/vUed4NWFH00dmjaJ3SzMeIjAkP/Xt1YuqoO7zyS/+7lmv5hT5GFNksQZkKSkpK+dVrKykrKwOg7+0dmDamr89RGRM+Hp0zmpRkZ4Lki1eu88d31/scUeSyBGUqeGvJFo7mnAOchQj/Yd4EnAnmjTEAzZsmVrjlvXLTATbtPOpjRJHLEpTx7D9ypsIy1w/fM5LO7Vv5GJEx4WnM4B6MHdrTK7/wxiouXrnuY0SRyRKUAaCgsJhfvrrCmzK+f6+OzJ480NeYjAlnTz443ns2m3etgBdeX2UDeBuYJSgDwB/+uo5Pz18BnF57T39ust3aM6YGzZsm8vTDk73ylr3ZNoC3gVmCMqzOPFhhQO6X546z6YyMqYNBfTpz98TyRTv/+O567xmuuXWWoKLc8dMXeOGN8gG5Y4f2tHWejLkJj9w7ki4dWgNQXFLKf/7hQ+t63kAsQUWx/IIi/vOlpd4ihB1TU3jqs9Zrz5ibkRAfx7e+OJ0mifEAfHr+Cr96daU9j2oAlqCilKry/OurOOXOFpEQH8e3npjhrSBqjKm7Tu1a8tXPTfLKm3cf472Awe6mfixBRal3PtzKhu2HvfJT8ybQ1b1NYYy5eWMG9+CeieU9X19bsInMPdk+RtT4WYKKQmsyD/HGos1eecbYfkzI6O1jRMZEhs/PHkmf7rcBoMDP/rTcOk3cAktQUWbf4dP86i8rvfKA3p344gNjfIzImMgRFxfLt5+YQWorpxdsYVExP/ztIs5dtKU56sMSVBQ5efYSP/r9EkpLnXn2Ordvxbe+ON1biM0Yc+tSkpP43j/Moqn7PPfilev8x4uLrWdfPViCihJnzl3h2V8v8FbHbdE8iX/5+7tolpToc2TGRJ4utzlf/mJinI/Y7FPn+eFvF5NfUORzZI2LJagokHshj2d/vYALl68BTo+97355Ju3btPA5MmMi18A+nfnKvIle+cDRM/zHi4spLCr2MarGxRJUhDt/6SrPPr+A3It5AMTHxfIvT95F77T2PkdmTOSbPLIPX3xgrFfee/g0P/pd+dhDUzNLUBHs5NlLfO+/3+fMOWeOvdjYGP7PEzMY0LuTz5EZEz3unjiAz88e5ZV3Hszh315YaM+k6sASVIQ6lP0p3/vv97yWU0xMDN/8wnSGpnf1OTJjos99Uwczb9Zwr7z38Gn+9Rfvc/6S9e6riSWoCLR173F+8KsF5F0rAJzbet/+0gxGDEjzNzBjotiDM4ZVaEkdP32Bf/nv9zhx5qKPUYU3S1ARRFV5a0km//HbRd497uZNE3nua7PJ6NfN5+iMMfdNHcw/PjLF69137uJVvv3Tv7Ju2+FajoxOcX4HYBpG3rUCfvHKR2zbd8Lbl9oqme9/5W46tWvpY2TGmEATh/cmuVkT/uuPyygsKqawqJif/WkZB49+yudnj7RxiQGsBRUBtu8/wTf/850KySm9Rwf+3zP3W3IyJgwNTe/Kf3x9Dre1LR/q8cGqnXzvF+/bLb8AEg1TwmdkZGhmZqbfYTS4q9cL+dN761m56UCF/fdPHcxDd48gNta+fxgTzq7lF/KrV1eyefcxb19sbAwPzhjG/VMHR2xrSkS2qGpGrfUsQTU+paVlfLz5AK8v3MzFK9e9/c2bJvKVhyYxcmB3H6MzxtwMVeW9j7bz+qLN3jRkAF07tObROaMZfEfniFujzRJUgEhJUKrK5t3ZvLZgEzmfVrwNMHpwD740dywtk5v6FJ0x5lZkn7rA839ZyeETuRX2D+jdiUfuGUnPbu18iqzhWYIK0NgTVFFxCWu2HOKDj3dx/PSFCq+1TG7Kl+aOY/Tg232KzhjTUMrKyljw8S5eX/gJxSWlFV4b2Lszd08awLD0ro2+RWUJKkBjTFCqyuHjuazbdpiVnxzwxjTdkJgQz31TBzF78iBvqWljTGS4cPkaby3J5KMN+ymr9BndITWFSSP6MHZIDzqkpvgU4a2xBBWgsSSogsJi9h4+za6DJ9m08yifnr/yN3US4uOYNqYvn5k2lJTkJB+iNMaEysmzl3h94WY2bj9MVZ/Ut3dJZXj/bgzs3ZmeXVMbTaeKsEhQIjIT+AUQC/xeVX9U6fVE4M/AMOA88FlVPea+9l3gCaAU+EdVXVqXc1Yl3BKUqnL1eiGfnrvCsVPnOXLiHEdycjl84hxlZWVVHpPaKpm7JvRn6qg7aN7UlsgwJpqcvZDH4tW7Wb5hH9erWbIjMSGe3mnt6N6pLd07t6Fbxza0b9MiLO+w+J6gRCQWOAhMA3KAzcBDqro3oM5XgIGq+g8iMg+4X1U/KyLpwOvACKAjsBy4sSZ5jeesSn0T1Oncy94SFYFuXDNVZ7tMlbIypbSsjJKSMkpLyygsLuZ6fjEFRcVcu17I5av55F0r4MLl65y9kFendWGSmiQwcmB3xg7pwaA+na3buDFRrqCwmMzd2azdmsXWfccr9PqrTovmSaS2ak7L5Ka0SG5CSvMkkpokkJQYT1JiPAnxccTGxhAfH0uMCLGxMQggIohQ7fOuTu1b1rtTVl0TVDBnkhgBZKnqETegN4A5QGAymQM8626/A/xanKsxB3hDVQuBoyKS5Z6POpyzwSxZs4cPVu0Mxqmr1bVDawb27syAPp0Y2LsTCfE22YcxxtEkMZ5xw3oyblhPrl4vZNve4+w8eJJdB096E0NXduVqPleu5jd4LP/0+SlMyOhde8VbEMxPv07AiYByDjCyujqqWiIil4E27v6NlY69sUZEbecEQESeBJ4E6No1/GbwToiPo32bZDrf1prundvQvVNbenRJtedKxpg6ad40kfEZvRif0QtwbgMePp7LsZPnOJpznpNnL5J78WqdWlnhKpgJqqp2YeX7idXVqW5/Vfe4qrxHqaovAi+Cc4uv+jCr175tMn1v7/A3+wNbvDExQozEIAJxsbHEx8UQGxdLYnwcTRLjSEpMIKlJPCnNk2iRnETL5kmktk6mRfMmjb6rqDEmfLRrnUy71skVhpyUlZVx4fJ1zl28ymW3JXX5agEFBUUUFJWQX1hMUXEJpaXO44mS0lJUoUzLKCsr/9is6klQSgjGXAYzQeUAXQLKnYFT1dTJEZE4IAW4UMuxtZ2zwcyaMIBZEwYE6/TGGBNUMTExtG3VnLatmvsdSr0E86n7ZqCXiHQXkQRgHjC/Up35wGPu9lxghTo9EOYD80QkUUS6A72AT+p4TmOMMREgaC0o95nS08BSnC7hf1DVPSLyHJCpqvOBl4BX3E4QF3ASDm69t3A6P5QAX1XVUoCqzhms92CMMcY/NlDXGGNMSNW1m7kNrDHGGBOWLEEZY4wJS5agjDHGhCVLUMYYY8KSJShjjDFhKSp68YlILpB9C6doC5xroHAaK7sGdg3ArgHYNYBbvwbdVDW1tkpRkaBulYhk1qVLZCSza2DXAOwagF0DCN01sFt8xhhjwpIlKGOMMWHJElTdvOh3AGHAroFdA7BrAHYNIETXwJ5BGWOMCUvWgjLGGBOWLEEZY4wJS5agAojITBE5ICJZIvKdKl5PFJE33dc3iUha6KMMrjpcg2dEZK+I7BSRj0Skmx9xBlNt1yCg3lwRURGJuC7HdbkGIvJ37r+FPSLyl1DHGGx1+FvoKiIrRWSb+/cwy484g0VE/iAiZ0VkdzWvi4j80r0+O0VkaIMHoar24zyHiwUOA7cDCcAOIL1Sna8Av3G35wFv+h23D9dgMtDU3X4qGq+BWy8ZWA1sBDL8jtuHfwe9gG1AK7fczu+4fbgGLwJPudvpwDG/427gazABGArsrub1WcBiQIBRwKaGjsFaUOVGAFmqekRVi4A3gDmV6swBXna33wGmioiEMMZgq/UaqOpKVb3uFjcCnUMcY7DV5d8BwL8BPwEKQhlciNTlGnwZeF5VLwKo6tkQxxhsdbkGCrRwt1OAUyGML+hUdTXOQrLVmQP8WR0bgZYi0qEhY7AEVa4TcCKgnOPuq7KOqpYAl4E2IYkuNOpyDQI9gfMNKpLUeg1EZAjQRVU/CGVgIVSXfwe9gd4isk5ENorIzJBFFxp1uQbPAo+ISA6wCPhaaEILGzf7eXHTgrbkeyNUVUuoch/8utRpzOr8/kTkESADmBjUiEKvxmsgIjHAz4HHQxWQD+ry7yAO5zbfJJxW9BoR6a+ql4IcW6jU5Ro8BPxJVX8qIqOBV9xrUBb88MJC0D8PrQVVLgfoElDuzN822b06IhKH06yvqQnc2NTlGiAidwLfA2aramGIYguV2q5BMtAf+FhEjuHce58fYR0l6vq38L6qFqvqUeAATsKKFHW5Bk8AbwGo6gagCc4kqtGiTp8Xt8ISVLnNQC8R6S4iCTidIOZXqjMfeMzdngusUPdpYYSo9Rq4t7d+i5OcIu25A9RyDVT1sqq2VdU0VU3DeQ43W1Uz/Qk3KOryt/AeTocZRKQtzi2/IyGNMrjqcg2OA1MBRKQvToLKDWmU/poPPOr25hsFXFbV0w35C+wWn0tVS0TkaWApTg+eP6jqHhF5DshU1fnASzjN+CycltM8/yJueHW8Bv8JNAfedvuHHFfV2b4F3cDqeA0iWh2vwVJguojsBUqBb6nqef+iblh1vAb/DPxORL6Bc2vr8Uj6wioir+Pcwm3rPmf7v0A8gKr+Bue52ywgC7gOfKHBY4ig62mMMSaC2C0+Y4wxYckSlDHGmLBkCcoYY0xYsgRljDEmLFmCMsYYE5YsQZmIJyKlIrJdRHaLyNsi0vQmj79ahzrPishJ9/fsFZGH6hFnlb9HRL7nzhi+0z3/yJs9dz1iec4dkH0zxwwRkd/XUudpEWnw7sgmMlk3cxPxROSqqjZ3t18DtqjqzwJeF5y/hSqnqAk8vobf8SxwVVX/S0R6AVuANqpaXJ84A/aNBn4GTFLVQndQbIKqht3EpCLyNvDvqrqjhjpNgXWqOiR0kZnGylpQJtqsAXqKSJqI7BOR/wG2Al1E5CER2eW2tH4ceJCI/FREtoqzBlZqTb9AVQ/hDFxs5R7bQ0SWiMgWEVkjIne4+7uLyAYR2Swi/1bN6ToA525MKaWq524kJxGZ6q5FtMtduydRRFLEWcOoj1vndRH5cqX3MkJE/upuzxGRfBFJEJEmInLE3f8nEZnrbh8Tkf/Pff+7bsRf6ZzJwMAbyUmcdYJ+4G7PEJHVIhLjzoR/TERG1HQNjQFLUCaKiDN/4l3ALndXH5zlAoYAxcCPgSnAYGC4iNzn1msGbFXVocAqnBH1Nf2eocChgKmgXgS+pqrDgG8C/+Pu/wXwgqoOB85Uc7oPcZLnQRH5HxGZ6P6OJsCfgM+q6gCcWWGeUtXLwNPAn0RkHs56Tb+rdM6twI0WzHhgNzAcGAlsqiaOc+77f8F9D5VluOe54TvAZ0VkMvBL4AsBLdRM9/caUyNLUCYaJInIdpwPxuM4U1YBZLvr2IDzAf2xqua6S6m8hrNgG0AZ8Ka7/Sowrprf8w0ROYDzIf8sgIg0B8bgTA21HWcewxtr5owFXne3X6nqhKp6FRgGPIkzz9ubIvI4TnI9qqoH3aov34hXVZfhJOHngS9Vcc4SIMudP24Ezi3ECThJY0017+2v7n+3AGlVvN6BgHno3JbSl4FlwK9V9XBA3bNAx2p+jzEem4vPRIN8VR0cuMOdR/Ba4K6bOF91D25/7j6DegD4s4j0wPkSeKny76/DucorqJYCH+PMoL4LZ8Li7dXVF2dJkL5APtAaZ9bpytbgtCaLgeU4rbFYqm4dAdyYtb6Uqj838nEmSw00ADjP3yajJm59Y2pkLShjHJuAiSLSVkRicdb6WeW+FoMzez3A54C1NZ1IVf+K01p7TFWvAEdF5EFwOmSIyCC36jrKJxx+uKpziUgft9PFDYOBbGA/kCYiPd39nw+I9xvAPvc9/EFE4qs49Wrg68AGVc3FWXjzDmBPTe+tBvuAG7EgIt1wJlMdAtxVqedhbyreDjSmSpagjAHcZQK+C6wEduA8c3rfffka0E9EtuA8o3quDqd8DnjGbc08DDwhIjtwEsCNpcP/CfiqiGzGWVusKs2Bl92u6zuBdOBZVS3AmT36bbdVVQb8RkR649zW+2dVXYOTiP61ivNuAtq7rwPsBHbWdzZuVd0PpIhIstsr8iXgm26HjieA37vPzcC5tbm8Pr/HRBfrZm6MaRDiLDuRp6rVjoUSZz2xZ1T186GLzDRW1oIyxjSUFyh/VlWdtsD3QxCLiQDWgjLGGBOWrAVljDEmLFmCMsYYE5YsQRljjAlLlqCMMcaEJUtQxhhjwtL/D4dGK4bg12THAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "beta = Beta(5, 5)\n", "prior = beta.MakePmf()\n", "thinkplot.Pdf(prior)\n", "thinkplot.decorate(xlabel='Prob Red Sox win (x)',\n", " ylabel='PDF')" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [], "source": [ "%psource beta.Update" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd0XPd54P3vMzPoBEAQBEk0EmABAbB3iaREsVklWil2pLXkOJJix46daJO14/escnbjdbQ+Z+PNJk7yypajWG56bUuyXMIksqjCTpESwQo2kCAJgCgkKlGJMpjf+8cdXMyAaCQxmPZ8zoF07517Zx5cDuaZXxdjDEoppVSocQQ7AKWUUmo4mqCUUkqFJE1QSimlQpImKKWUUiFJE5RSSqmQpAlKKaVUSNIEpZRSKiRpglJKKRWSNEEppZQKSa5gBzBRpk+fbvLy8oIdhlJKqTEcPXq00RiTMdZ5EZOg8vLyKCkpCXYYSimlxiAileM5T6v4lFJKhSRNUEoppUKSJiillFIhKWLaoIbT19dHdXU13d3dwQ4lLMXHx5OTk0NMTEywQ1FKRaGITlDV1dUkJyeTl5eHiAQ7nLBijKGpqYnq6mry8/ODHY5SKgpFdBVfd3c36enpmpzugIiQnp6upU+lVNBEdIICNDndBb13St0eYwxud3+ww4gYEV3Fp5RSgdTd08ev3z9O6cVamm900tzWiaffw+89uIqnHl6tX/LukiaoEPH1r3+d+++/n23btgU7FKXUOPT19fOt7+/k1IXqWx57a+dR2ju6+cKTGzVJ3QVNUCGgv7+fF1988bavcTqdAYpIKTWa/n4P//CT94dNTgN2HjzDzZ5env/MZpzOiG9NCYioSVA///nPA/bcTz/99IiPVVRU8NBDD7Fu3TqOHz9OQUEBP/nJTyguLuZzn/sc7777Ls8//zzvvPMOjz76KE888QQffPABX/va13C73axZs4aXX36ZuLg48vLy/K556qmnAvY7KaWGZ4zhu6/v5fCpK/axT21bweZ1C0lNTuCVX+znwNFyAPaVXKTP7eEvntumJak7ENC0LiIPiUiZiJSLyAvDPB4nIm94H/9IRPK8x/NE5KaInPD+fC+QcQZaWVkZX/ziFzl16hQpKSl897vfBaxxRgcOHPBLNN3d3Tz33HO88cYblJaW4na7efnll+3Hh7tGKTV53nznKHs+LrP3H920lM88upasGVNJSojjzz+7he3ri+zHD524xPFzV4MRatgLWIISESfwHeBhoBh4WkSKh5z2eaDFGDMf+DbwLZ/HLhljlnt/vhSoOCdDbm4uGzZsAOCzn/0sBw4cAODTn/70LeeWlZWRn59PQUEBAM8++yz79u2zHx/uGqXU5Gi60cGv3z9u729ZV8hzn7zXr3TkcDj44/98P5vWFNjHfrHzKMaYSY01EgSyim8tUG6MuQwgIq8DjwNnfc55HPiGd/st4CUJUDl4tGq4QBv6Kw3sJyUl3XLuWG/i4a5RSk2OX+w8Sp+3G3l+znS+/NT9w1bdiQif+Z21HDhWTn+/hwsV1zl9sZYlBdmTHXJYC2QVXzbgW66t9h4b9hxjjBtoBdK9j+WLyHER2Ssi9wUwzoCrqqri0KFDgNUWtnHjxhHPLSwspKKigvJyqw77tddeY9OmTZMSp1JqZHUNrXxweLBq7/cfXYfDMfJH6PS0KWxZt9Def+vdowGNLxIFMkENVxIaWjwY6Zw6YLYxZgXwVeBnIpJyywuIfFFESkSkpKGh4a4DDpSioiJ+/OMfs3TpUpqbm/nyl7884rnx8fH88Ic/5Mknn2TJkiU4HA6+9KWwruFUKiK8/tsjeDweAIrnZbK8MGfMaz65bQUObwnr9MVazl++FtAYI00gq/iqgVyf/RygdoRzqkXEBaQCzcaq5+oBMMYcFZFLQAHgtyKhMeYV4BWA1atXh2wFr8Ph4Hvf8+/nUVFR4bf/ox/9yN7eunUrx48fZ6ih1yilJkdlbRMHvT3zAD77n9aNq1fezPQU7l9TYHeqeOvdo/yPL/1OwOKMNIEsQR0BFohIvojEAk8BO4acswN41rv9BLDLGGNEJMPbyQIRmQssAC4HMFallBrRz//jiF39s3rRHBbmzxr3tZ/avsKuKjp+7irllfUTHl+kCliC8rYpPQ/sBM4BbxpjzojIiyLymPe0V4F0ESnHqsob6Ip+P3BKRE5idZ74kjGmOVCxBlJeXh6nT58OdhhKqTtUU3+DI6crAKtN4jOPrr2t67NnTGX9yvn2/t6SCxMYXWQL6EBdY8zbwNtDjn3dZ7sbeHKY634J/DKQsSml1HjsPnze3l69OI85WemjnD28bfcUcvCYVUX40akrfO5TG3Tg7jjo/BtKKTWC/n4Puz8eLPFsvbfwjp6neF4mSQlxADTd6OTy1cYJiS/SaYJSSqkRHDtXxY32LgCmJieysmj2HT2Py+Vk9eI59v5HPtMkqZFpglJKqRHs8qne27y24K4mfV27JM/e1gQ1PpqgwsBvfvMbzp49O/aJQ+zYsYO/+Zu/CUBESkW+G+1dlJypsvc333Nn1XsDVhTlEuOyViCovt5CTf2Nu3q+aKAJKgzcSYJyu9089thjvPDCLXP0jnqNUsqy98hFe2DuwvxZZM+YelfPFxcbw/LCwaGhH2spakxRs9zG7/154CZE/+U/jjzTw0jLbRw6dGjYJTVeeOEFduzYgcvl4hOf+ASf+tSn2LFjB3v37uWb3/wmv/yl1bnxT//0T2loaCAxMZF/+Zd/obCwkOeee45p06Zx/PhxVq5cyZIlSygpKeGll16isrKSz33uczQ0NJCRkcEPf/hDZs+efcs1f/d3fxew+6RUuDDG8MGhc/b+trssPQ1YtzTf7rL+cWkFn9y2YkKeN1JFTYIKprKyMl599VU2bNjA5z73Of7+7/+ef/7nf+aDDz6goKCAZ555hpdffplnnnmGX//615w/fx4R4caNG0ydOpXHHnvMXisKrJkmvve977FgwQI++ugj/uRP/oRdu3YBcOHCBd5//32cTqff7BTPP/88zzzzDM8++yw/+MEP+LM/+zN+85vf3HKNUgouVtbbVXBxsTGsXzFvQp539eI5OETwGMOFius0t3YyLVUngB6JVvFNgqHLbXzwwQfDLqmRkpJCfHw8f/RHf8SvfvUrEhMTb3mujo4OPvzwQ5588kmWL1/OH//xH1NXV2c//uSTTw6baA4dOsRnPvMZAP7gD/7AXvJjtGuUilYHj12yt9evmEt8XMyEPG9yUjzF8zPt/ZLTlRPyvJEqakpQo1XDBdp4B+S5XC4+/vhjPvjgA15//XVeeuklu2Q0wOPxMHXqVE6cODHsc4x3OQ7fmHQJD6UGGWP8etmtXz4xpacBa5fkc/qiNS3px6VX+MSGocvkqQFagpoEQ5fb2LZt27BLanR0dNDa2sojjzzCP/zDP9hJKDk5mfb2dgBSUlLIz8/nF7/4BWD9MZ08eXLMGNavX8/rr78OwE9/+tNRl/xQKppV1DTR0GL9vSXEx7J0gtdwWrVocDzU+SvX7Y4Y6laaoCbB0OU2vvKVrwy7pEZ7ezuPPvooS5cuZdOmTXz7298G4KmnnuJv//ZvWbFiBZcuXeKnP/0pr776KsuWLWPRokX867/+65gx/NM//RM//OEPWbp0Ka+99hr/+I//GOhfW6mwdNin9LSyeDYu18RWf89MTyYtxaq+v9ndS2VtWE4zOikkUpYhXr16tSkp8VuNg3PnzlFUVBSkiCwVFRU8+uijYTthbCjcQ6Um01f+5k2q6qyk8dXntrNhgjpI+Pq/P3yPQyesdq7P/94GHrl/yYS/RigTkaPGmNVjnaclKKWU8qpraLWTk8vlZFXxnU1tNJaiuYPLdZzTRQxHpAkqwHS5DaXCx8elFfb2soKcCeu9N1TxvMGefOcu1REpNVkTLeITlP7D3zm9dyra+PbeW7csL2CvMztzmp38Wtq6aGjpCNhrhbOITlDx8fE0NTXpB+0dMMbQ1NREfHx8sENRalK0tHVx4YpV3SbA6kV5AXstp9PBwryZ9v65S3WjnB29InocVE5ODtXV1TQ0NAQ7lLAUHx9PTk5OsMNQalIcKa2wl3UvmpdJanJCQF+vaF4mJ8uqATh3uY5NawoC+nrhKKITVExMDPn5+cEOQykVBgbmyAO4Z9ncgL+eX0eJS9pRYjgRXcWnlFLj0dvnpvRCjb3vO5g2UBbMmYHDYX0EV19vob2zO+CvGW40QSmlot65y9foc/cDkJmRyqzpKQF/zbjYGOblTrf3z1/RUtRQmqCUUlHvxLmr9vaKotxRzpxYRXMHu5uf1/FQt9AEpZSKesfPDyaoZYWTl6AKdcDuqDRBKaWiWmNLB1e9s0c4nQ4Wz8+atNcuzB9MUOVV9fT26arWvjRBKaWi2smywdJT8bzMgM0eMZzU5AQyM1IB6O/3UKUTx/rRBKWUimrHz1Xb2yuKAjP33mjysgc7SlypaZz01w9lmqCUUlGrv9/DqbLBBLW8cPIHps/N8UlQ1U2T/vqhTBOUUipqlVfV03mzB4C0lERmZ06b9BjyfRLU5Wqd9caXJiilVNTy7b23vCgXEZn0GHxLUJW1zbrCrg9NUEqpqOU//mny25/A6igxsMJub5+bmvrWoMQRijRBKaWiUufNHsor6wFr9vKlBdlBi2VuToa9XVGtHSUGaIJSSkWl0xdr7dnL5+ZmkJwUvKVl8nPS7e3LmqBsmqCUUlHp9MXByWGDWXoC/67m2lFiUEATlIg8JCJlIlIuIi8M83iciLzhffwjEckb8vhsEekQka8FMk6lVPTxnb18cZAT1Nxc/67musiqJWAJSkScwHeAh4Fi4GkRKR5y2ueBFmPMfODbwLeGPP5t4LeBilEpFZ1a2rq4eq0FsKY38l2bKRhmTEsmMT4WsNrGdAl4SyBLUGuBcmPMZWNML/A68PiQcx4HfuzdfgvYKt5+niLyu8Bl4EwAY1RKRaEzF2vt7YV5M4mLnbzpjYYjIn7joa5oOxQQ2ASVDVz12a/2Hhv2HGOMG2gF0kUkCfhvwF+P9gIi8kURKRGREl3WXSk1XqcuDM4esSTI1XsD8v3aoTRBQWAT1HAj3oZWrI50zl8D3zbGjFrONca8YoxZbYxZnZGRMdqpSill821/WrIgNBKUbzuUdjW3uAL43NWA78IqOUDtCOdUi4gLSAWagXXAEyLyf4CpgEdEuo0xLwUwXqVUFLje1EZ9czsAsTEuFsyZEeSILDpp7K0CmaCOAAtEJB+oAZ4CPjPknB3As8Ah4Algl7G6r9w3cIKIfAPo0OSklJoIvt3LF83PxOVyBjGaQTkzpxLjctLn7qfpRidtHTdJmZIQ7LCCKmBVfN42peeBncA54E1jzBkReVFEHvOe9ipWm1M58FXglq7oSik1kU75di8Pkeo9sHoT+k5We6VGZzYPZAkKY8zbwNtDjn3dZ7sbeHKM5/hGQIJTSkUdYwynLwy2NAR7gO5Qc3Onc+mq1eGrsraJZQsnf/mPUKIzSSilokb19RvcaO8CICkhjrzs9DGumFy5swZLUFfrWoIYSWjQBKWUihq+458WL8jC4Qitj8DcWWn29tVruvx7aP3rKKVUAJ0uH0xQi+ZnBTGS4eX4JaiWqJ/ySBOUUioqGGM4U+5fggo1aSmJJCXEAdDd00fTjc4gRxRcmqCUUlGh+voN2jpuAjAlMS4oy7uPRUTIzRwsRVXVRXc1nyYopVRU8G1/WjQ/KyjLu4+HbztU9fXo7iihCUopFRVCvf1pgPbkG6QJSikV8Ya2Py2anxnEaEanPfkGaYJSSkW8mvrB9qekhDjmZIXW+Cdf2pNvkCYopVTE829/ygzZ9ifQnny+NEEppSJeuLQ/wa09+QZW/o1GmqCUUhHt1van0E5QoO1QAzRBKaUiWm1DK63tg+1PoTb/3nC0J59FE5RSKqL5tj8Vzwvt9qcBWoKyaIJSSkW0cGp/GpDjN1j3RtT25NMEpZSKWMYYzob4/HvD8e3Jd7O7N2p78mmCUkpFrLqGVlrarPWfEuNjmZMVevPvDUd78lk0QSmlItbZS3X2dtHczJBb/2k02g6lCUopFcH8upeHSfXeAO3JpwlKKRWhho5/Kp43K4jR3D7fjhI19TeCGEnwaIJSSkWk+uZ2u3NBXGwMc3MyghzR7cmeMdXeronSZTc0QSmlItLZct/2p1k4neH1cZc+NYnYGBcAHV09tHd2BzmiyRde/2JKKTVOvuOfikN4eY2RiAhZfqWo6Kvm0wSllIpIvuOfFs0Lrw4SA7JmpNrbNfXRV82nCUopFXEaWzqob24HIMblZP7s8Gp/GuBbgqqtbw1iJMGhCUopFXF8e+8Vzp2Fy+UMYjR3Lker+JRSKrL4dy8Pv/anAdkzfUtQmqCUUirs+c4gES4TxA7Ht4rvWlMb/f2eIEYz+TRBKaUiSnNrJ3UNVnuNy+VkwZwZQY7ozsXHxTAtNQmA/n4P15vaghzR5NIEpZSKKL7jnxbmzbTHEoUr3558tQ3R1VFCE5RSKqKcLq+xt8O5em9A9gyfKY+irKOEJiilVETxXUE3XNZ/Go1fCSrKOkoENEGJyEMiUiYi5SLywjCPx4nIG97HPxKRPO/xtSJywvtzUkQ+Gcg4lVKRobm1064GC/f2pwHRPJtEwBKUiDiB7wAPA8XA0yJSPOS0zwMtxpj5wLeBb3mPnwZWG2OWAw8B/ywi4V2RrJQKON/2p4I5M8K+/Qn8u5pH26zmgSxBrQXKjTGXjTG9wOvA40POeRz4sXf7LWCriIgxpssY4/YejwdMAONUSkUIv/anCKjeA8hIm0KMd6BxW8dNOm/2BDmiyRPIBJUNXPXZr/YeG/Ycb0JqBdIBRGSdiJwBSoEv+SQsm4h8UURKRKSkoaEhAL+CUiqc+LU/RUAHCQCHw0Fmhs+cfFFUzRfIBCXDHBtaEhrxHGPMR8aYRcAa4C9FJP6WE415xRiz2hizOiMjPOfaUkpNjJa2Lrv9yel0UJA3M8gRTRz/Ofk0QU2EaiDXZz8HqB3pHG8bUyrQ7HuCMeYc0AksDlikSqmw5zu9USSMf/KVHaWTxgYyQR0BFohIvojEAk8BO4acswN41rv9BLDLGGO817gARGQOsBCoCGCsSqkw5zf/XoRU7w3w6ygRRavrBuwrhjHGLSLPAzsBJ/ADY8wZEXkRKDHG7ABeBV4TkXKsktNT3ss3Ai+ISB/gAf7EGNMYqFiVUuEvEtufBvitCxVFs0kEtAxsjHkbeHvIsa/7bHcDTw5z3WvAa4GMTSkVOW60d9ldsK32p/Af/+RraBuUx+PB4Yj8eRYi/zdUSkW80z6lpwVzZhAXGxPEaCZeUkIcqckJgDVpbOONziBHNDk0QSmlwt7pi4PjnxYvGDqaJTJkZURfTz5NUEqpsFd6YTBBLS2IzATlOxZKE5RSSoWBxpYOrjVa6yTFuJwUzImc8U++fBNUXZR0lNAEpZQKa77Ve0VzM4mJcQYxmsDx7cmnCQoQkXd9tv8y8OEopdTtKfXpIBEp8+8NJysKB+uOVYLynT/olu7gSikVTMYYSi9U2/uR2v4EMGt6ij03XH1TG319/UGNZzKMlaB0FnGlVMi61thGk7fLdXxcDPNyI3dOztgYF9PTkgHrg/laU1twA5oEYw3UnSsiO7AmdR3YthljHgtYZEopNQbf9qdF87JwOiO7WT0zI5WGlnbAaofKnZU2xhXhbawE5bt+0/8NZCBKKXW7fNufFhdEbvvTgKwZqZzyVmlGQ0eJUROUMWbvwLaIZHiP6cJLSqmgs9qfBktQSyJ0gK6vaBsLNVYvPhGR/ykijcB54IKINIjI10e7TimlAu3qtRbaOm4C1lRAednpQY4o8KJtXaixKmz/K9bM4muMMenGmDRgHbBBRL4S8OiUUmoE/qWnLESGW/80skTbYN2xEtQzwNPGmCsDB4wxl4HPeh9TSqmg8E1QkTz+ydeMacl2R5CWti66e/qCHFFgjZWgYoZbh8nbDhVZ0wUrFWKMMRijIz2G43b3U+rTg29ZYe4oZ0cOp9PBrPQUez/SS1Fj9eLrvcPHlFJjMMbQ3t5OU1MTbW1ttLe3097eTm9vL263m76+PowxuFwuXC4XMTExJCUlkZycTHJyMlOnTmX69Ok4nZE5tc9oyqsa7NJD+tQksnyqviJdZkaqvfZVTf0N8nOmBzmiwBkrQS0TkTawBzAPfJ0TID5gUSkVoXp6eqiurqauro6Ghga6u7vHvMbtduN2u+nu7qa9vZ1r167ZjzkcDqZNm8aMGTPIzc0lLS0tKtpiTpRdtbeXLcyNit95QNaMqXCmEoj8jhJjdTOPvq9mSk2w/v5+Kisrqaqq4tq1axNabefxeGhsbKSxsZGzZ88yZcoUZs+ezdy5c0lOTp6w1wk1p8p8qvcW5gQxkskXTR0lRk1QIhIPfAmYD5wCfmCMcU9GYEqFu87OTi5evMilS5fo7R2+RjwmJoaMjAymTp1KcnIyKSkpJCQk2FV6ImKXoHp6eujo6KC9vZ22tjYaGxtpa/Of7qajo4OzZ89y9uxZMjMzKSgoIDMzM6JKGJ03e7hYcd3eXxLB8+8NJ5pmNR+riu/HQB+wH3gEWAT8eaCDUiqcdXR0cPr0aSoqKoYtLaWnp5Obm0tmZiapqaljJo+YmBhiYmJISEhg6tSpfo91d3dTX19PTU0N1dXVuN2D3x/r6uqoq6sjJSWFpUuXkpOTExGJ6kx5HR7vfc3Lnm4vhR4tfEtQNddvYIyJiH/X4YyVoIqNMUsARORV4OPAh6RUeOrq6uLMmTNcunTplsSUmJjIggULmDNnDklJSRP2mvHx8cyePZvZs2fT399PXV0dly5dorZ2cAqgtrY2Dhw4wNSpU1m2bFnYl6hOnh9sf1peGF3VewDTUpOIjXHR2+emq7uX9s5uUqZEZpIeK0HZneyNMe5wflMrFSgej4eysjJOnz7tV4IBmDFjBgsXLiQ7OzvgScHpdJKTk0NOTg7t7e129eJATDdu3GDv3r1kZmayatWqsG2jOlXms7xGlLU/AYgIWTOmUlFjjQCqrW+N2gQ10IsPrJ57CT69+owxJmXkS5WKfNevX6ekpOSWtqCMjAyWLl3KjBkzghJXcnIyK1euZNGiRZw7d44LFy7Q32+tH1RXV8fbb79NUVERixYtCqtu6g3N7dR6211iXE6K5s4KckTBkZmRaieouoZWCiP0PmgvPqXugNvt5vjx45SXl/sdT01NZcWKFcyaNSskqtHi4uJYvnw5CxcupLS0lEuXLgFWqe/MmTNUVVVx7733kp4eHvPYnfJZnLBobiaxMWN9x45M2TOiY9LY6PzXVeouNDQ0cOjQITo7O+1jLpeLJUuWUFBQgMMRemsSJSQksHbtWubPn8+RI0dobm4GoL29nffee4/i4mIWL14ckrH7OnF+MEEti8L2pwHRMmmsJiilxskYQ2lpKWfOnPE7npOTw+rVq0lICP12gGnTpvGJT3yCS5cucfz4cdxuN8YYzpw5Q21tLRs2bAjZtimPx+PX/hRt4598+XY1r4ngruaaoJQah5s3b/Lhhx9SX19vH4uJiWHVqlXk5eWFRHXeeIkI8+fPJzMzk0OHDtHQYC3x1tLSws6dO1m3bh25uaE3t115VQMdXT0ApKUkRsXyGiPxLUHVNbTi8XhCvvR7JyLvN1JqgtXX1/POO+/4JaeZM2fyyCOPkJ+fH1bJyVdSUhJbt25lxYoV9odbX18fBw4c4NixY3g8niBH6O/o2Sp7e3lRdE1vNFRSQpzdc8/t7qfxRucYV4QnTVBKjeLixYvs2rXLb868JUuWsHnzZhITE4MY2cQQEQoLC9m2bZvf71NWVsbu3bvp6ekJYnT+jvskqJXFs4MYSWjIioKOEpqglBqGx+OhpKSEkpISe9BtXFwcmzdvZvHixRH37T09PZ2HH36Y7OzBaYPq6+t59913aW0NfhvHjfYuLl21qiIdIlHd/jQgK2Owmq/muiYopaJCb28ve/bs4eLFi/axtLQ0HnroIWbNiszxJgCxsbHcd999LF261D7W0dHBe++9R11dXRAjgxPnBmePKJw7i6SEuCBGExqiYU4+TVBK+ejq6uL999/n+vXByUhzc3NvqQKLVCLCokWL2Lhxoz2At6+vj71793LlypUxrg4c3/anFUVavQf+HSW0BHUHROQhESkTkXIReWGYx+NE5A3v4x+JSJ73+HYROSoipd7/bwlknEoBtLa28t577/lVaS1evJgNGzbgckVXh9fc3Fy2b99uJ2VjDIcPH+bs2bOTvspvf7+Hkz7jn1Yt0gQFt/bki0QBS1Ai4gS+AzwMFANPi0jxkNM+D7QYY+YD3wa+5T3eCPwn70S1zwKvBSpOpQAaGxt5//336erqAqySxD333MOSJUsirr1pvNLS0ti+fTupqYNVSSdPnuTYsWOTmqQuVtbTedPqrDEtNYnZmdMm7bVD2az0FHsl2caWdnr7Im8lpECWoNYC5caYy8aYXuB14PEh5zyOtaQHwFvAVhERY8xxY8zAdMxngHgR0UpnFRDXr19n9+7d9ppNLpeLTZs2kZ+fH+TIgi8xMZFt27b5zSl44cIFPvroo0nrhn7Mr3ovuruX+4qJcTJzujUdqiEyS1GBTFDZwFWf/WrvsWHP8S6E2AoMHX33e8BxY8wt/V1F5IsiUiIiJQODDZW6HTU1NezZs8ee8TsuLo4tW7aQmZkZ5MhCR2xsLA888IDf4N0rV67w4YcfTkqSOqrdy0fkuzZUbb0mqNsx3NecofUCo54jIouwqv3+eLgXMMa8YoxZbYxZnZGRcceBquhUVVXF/v377Q/ZhIQEtm3bFjYTp04mp9PJhg0bmDt3rn3s6tWr7N+/354lPRCaWzvtWbsdDgdLC7R7uS+/jhIROBYqkAmqGvCdLyUHqB3pHBFxAalAs3c/B/g18Iwx5lIA41RRqKqqig8//NBuS0lKSmLbtm2kpOgKMiMREdauXUtBQYF9rLa2NqBJquR0pb1dPG8WiQmxAXmdcOU7Fkqr+G7PEWCBiOSLSCzwFLBjyDk7sDplW0DzAAAgAElEQVRBADwB7DLGGBGZCvwH8JfGmIMBjFFFoaHJKTk5mW3btjFlypQgRxb6RISVK1dSXDzY36muri5gSerj0sGu7WuXaJvgUJE+m0TAEpS3Tel5YCdwDnjTGHNGRF4Ukce8p70KpItIOfBVYKAr+vPAfOCvROSE9yc4K7+piDJcctq6dWtUjHGaKCLC0qVLWbRokX0sEEnqZncvpy7U2PtrluRN2HNHikhfdiOggzuMMW8Dbw859nWf7W7gyWGu+ybwzUDGpqJPdXX1sMkpHJbJCDUiwpIlSwDs5Ufq6uo4ePAgGzdunJCZtY+du0p/v9U+mJc9nRnTQnMZkGBKn5pEbIyL3j43HV09tHd2k5wUH+ywJozOJKGiwrVr1zh48KAmpwk0kKR8q/tqamo4dOjQhIyT8q/ey7vr54tEIhLRM0poglIRr76+nn379tm99aZMmaLJaYIMVPcVFRXZx6qqqvj444/vKkm53f0cOzPYvVwT1Mh8u5pHWkcJTVAqojU3N7N37167bSQxMZEtW7ZocppAIsKyZctYsGCBfezy5ct3NePEmUt1dHVbA6cz0pKjenHCsWT7rq57vSWIkUw8TVAqYrW1tfkNwo2Pj2fLli0kJSUFObLIIyKsWrXKb/aNCxcu2O1Tt+tIaYW9vWbJHJ09YhSRPBZKE5SKSF1dXX4L7sXGxrJ582aSk7WhPVBE5Jbl4ktLS/2WLRkPY4x2L78NOTPT7O3qa1qCUiqk9fT0sGvXLnvi14G59aZOnTrGlepuiQjr16/3WzerpKSEysrKUa7yd6W6kSbvEuZJCXEUzY3cNbgmQvbMwff1tcY23O7Azewx2TRBqYjidrvZu3cv7e3tgDU9zsaNG5k+fXqQI4seDoeD++67z2/KqEOHDnHt2rVxXX/45GDpadWi2bhczgmPMZLEx8WQkWbVDHiMoTaCOkpoglIRw+PxcPDgQZqamuxj99xzj078GgQDpdaBqaOMMezfv5/m5uZRrzPGcPB4ub1/z7K5o5ytBuTMGixFVUdQRwlNUCoiGGM4cuQItbWD0z2uWrWKOXPmBDGq6BYXF8fmzZvtWToGSrcdHR0jXnP5aiPXGtsASIiPZUVR7ojnqkGR2g6lCUpFhNLSUi5fvmzvFxcX+01qqoIjMTGRBx54gJiYGAC6u7vZvXs33d3dw57vW3pas3gOsTHRtZLxncqZ5ZOgImiwriYoFfbKy8v9ujPn5+ezdOnSIEakfKWmprJp0yZ7+qOOjg727dtnd/8fYIzhw+ODXzI2rJw/qXGGMy1BKRWCamtrKSkpsfczMzNZu3atjpsJMRkZGaxfv97eb2pq8psXEayl3RtarM4tSQlxLF+oaz+Nl29Pvpr6G5O22nGgaYJSYau5udlvfr20tLQJm6hUTbzc3FxWrVpl79fU1HD06FH73+/gscFl39Ytzdfee7chOSme1GRrdhS3u5/rTe1Bjmhi6F+yCkudnZ3s3bvXriZKTExk06ZNuFzaZhHKCgoKKCwstPcvXrzI+fPnreq9E4MJasPKecEIL6zl+rRDRcqMEpqgVNjp7e1lz549dkP7wCwROr9eeFi+fLnfbBMnTpxg18HjNLdag3OTk+JZsiA7WOGFrewZgwnqat3o3fnDhSYoFVY8Hg/79++nrc3qijwwKFSXag8fIsK9995LRkaGfezNf99HT481Oew9y/JxOvWj6Xb5j4XSEpRSk8oYw0cffUR9fb197J577mHGDF1sOdw4nU7uu+8+kpOT6fcYyqpu0NzchNvtZsMK7b13J3JnTbO3I2VWc01QKmyUlpZSUVFh7y9btkwH4oaxuLg4Nm3aRHVjD929/Xg8HnpvtrFgti6tcSd8e/JdvdYyIYtGBpsmKBUWLl++7DfWad68eX6L5KnwlJycTIdJsYcFLMhKYv/+/fb6XWr80lISSYyPBaC7p89u0wtnmqBUyLt+/TpHjhyx92fNmsXq1at1rFMEuNHeRVlFI2lpVgN/0ZypNDY2cvjw4YgoAUwmEYm4GSU0QamQ1trayv79++2Bh6mpqTrWKYLsPXIRjzEkJCSwojifqVOsEkBVVRWlpaVBji78RNqMEvpXrkJWd3c3e/fupa+vD4CEhAS/ed1UeDPGsPvjMnv/iUfuZf78wQ4SZ86c8ZtfUY3NvwQV/glKRzWqkNTf38++ffvo7LTq0Z1OJ/fff789M7YKf5eqGuzxOrExLtYvn0dcrIvOzk7q6uoA+Pjjj0lMTPRbAFGNLMeno4SWoJQKAGMMH374od+6Ths2bGDatGmjXKXCjW/paf2KeSTEx+JwONiwYQOpqamA9V44cOAAra2RswhfIPmWoCKhJ58mKBVyTpw4QXV1tb2/cuVKsrN1ZoFI0tvnZv/RwaU1Nq8dXBolJiaGBx54wJ4ZpK+vjz179nDz5s1JjzPczJiWTFysVQXe3tlNS1tXkCO6O5qgVEgZmJttQEFBAQsXLgxiRCoQDh67ROfNHsD6UF00P8vv8aFzK3Z1dQ27RIfyJyLkZQ+OI7tS3RjEaO6eJigVMmpqavyWzsjOzmblypVBjEgFgjGG/9g32EPvwY2Lhh0ykJaWxoYNG+z95ubmW5boULfKyxpMUBW1TaOcGfo0QamQMPDhM2DatGmsX79exzpFoIuV9fY3+xiXk633FI54blZWFqtXr7b3a2pqOHbsWMBjDGf5Ob4lKE1QSt2Vzs5Ov+qbxMRE7r//fl06I0K9ve+0vX3fqgUkJ8WPev6CBQv8lui4cOGCXzWw8udbxVdZo1V8St2xgaUzBhrAhzaQq8hyo73Lb92nh+9bNK7rhi7Rcfz4ca5evTrh8UWC2ZnTGKh3qGtopbunL6jx3A1NUCpo+vv7h106Y6CLsYo87314jv5+a1aQhfmzmJubMcYVloElOqZPn24fO3ToEA0NDQGJM5zFxcaQNcMaD2WAqjBeG0oTlAoKYwyHDx/2Wzpj3bp1zJw5M4hRqUByu/t59+BZe/+R+xbf1vUDg7WnTJkCDA7mHviCowbl5Qwm8oqa8G2HCmiCEpGHRKRMRMpF5IVhHo8TkTe8j38kInne4+kisltEOkTkpUDGqILj5MmTVFVV2fvLli0jLy8veAGpgDt88oo9w3ZqcgL3LMu/7eeIi4vjgQceIC4uDrCqiHfv3q1jpIbw68mnCepWIuIEvgM8DBQDT4tI8ZDTPg+0GGPmA98GvuU93g38FfC1QMWngqesrIxz587Z+/Pnz9elMyKcMYZfvjfY+277+mJcLucdPVdycjKbNm3C6bSu7+rqYs+ePfacjcq/o8SVMO4oEcgS1Fqg3Bhz2RjTC7wOPD7knMeBH3u33wK2iogYYzqNMQewEpWKIJWVlX7dhLOysli1apV2J49wR89W2W0hsTEufuf+26veGyo9Pd1vjNSNGzf8Zr2Pdn49+Wqbw3bsWCATVDbg282m2nts2HOMMW6gFRj3cpoi8kURKRGREm0sDX3Xr1/n8OHD9v7Ah4wunRHZjDH88t3BLyUPbigmZcrd99LMzs5m7dq19v7A+ytcP4wnUlpKon2Pe3r7qGsIz7kMA/nJMNxX4qHvnPGcMyJjzCvGmNXGmNUZGePrDaSCo6Wlxe8bbkpKit9UNipynSmv5ULFdQCcTgf/afPSCXvuefPmsWTJEnu/srKS48ePR32SEpGImFEikAmqGsj12c8Bakc6R0RcQCoQvn0i1bDa29vZvXv3Les6DTR0q8j2y3eP29tb1i0kfeqUCX3+RYsW+a0jVVZWxtmzZ0e5Ijr4zihREaYzSgQyQR0BFohIvojEAk8BO4acswN41rv9BLDLRPtXnwjT1dXFrl276OmxJgYdGIiblJQU5MjUZCivrOfUBWtmegEe37J8wl9DRFi9ejU5OTn2sVOnTlFeXj7KVZHPtx0qXHvyBSxBeduUngd2AueAN40xZ0TkRRF5zHvaq0C6iJQDXwXsrugiUgH8PfCciFQP0wNQhbienh52795NV5c15b/T6WTTpk1MnTp1jCtVpHjjncHJfzesmk9mRmAGYYsI69evZ8aMGfaxI0eO+A1liDZ52T5joWrDsydfQBsAjDFvA28POfZ1n+1u4MkRrs0LZGwqsAbW8BkYRCkibNy4EW0rjB6lF2o4dtZKEAJ8altgZ6YfGMi7a9cumputloJDhw7hcrnIysoa4+rIkz1jKi6XE7e7n6YbnbR3do8572Go0e5TasK53W727Nljf0gA3HvvvVH5IRGtjDG8tmOwx+amtQuZkxX4FZEHqpCTk5MB8Hg87N+/n+vXrwf8tUON0+lgdubgPb90Nfx6OmuCUhNqYPqZxsbBKoU1a9YwZ86cIEalJtvBY5fsD8QYl5OnH1kzaa8dFxfHli1bSExMBKwkNfQ9GS0K5gxWeZ67fC2IkdwZTVBqwng8Hg4ePOj3bXXFihV+PaxU5Ovr6+en//6Rvf/opiVMT5vYnntjSUxMZOvWrfas+AOl+qam8OwscKeK5mba22VXNEGpKOXxeDhw4AA1NTX2sSVLlvit46Oiw86DZ6hvbgdgSmIcn9y+IihxTJkyhS1bttjDGfr6+ti9e7df1XOkK5w7y96+UFFvzyQfLjRBqbs2UHLyTU5FRUUsWjS+tX5U5Ghtv8mb7xy195/4xCqSEoI33i0lJYXNmzcTGxsLDCaplpaWoMU0maanTSF9qjWko6e3L+y6m2uCUndlIDlVV1fbxwoLC1m2bJnOrxeFfvjrD+m8aY15m5mewkMbg/8lJS0tjS1btthJqre3l127dkVNkir0qeY7H2bVfJqg1B0bWHBwaHJavny5JqcodPzcVfYfvWjvf+HJ+4iJubMZyydaWlqaX0mqt7eXDz74ICo6ThTmD66xFm4dJTRBqTvidrvZu3cvtbWDs1ctXLhQk1OU6unt45U399n7G1fNZ0VR7ihXTL5p06axefNmYmJigMHqPt9FMyORb0eJ85frwmqeQk1Q6rYNDML17a1XXFzMihUrNDlFqTffOWp3jEhKiONzn9wwxhXBMW3aNLZu3Wp3nBjo3VdXVxfkyAJnTtY04uOspNzS1kVDS0eQIxo/TVDqtnR3d/P+++/ju7zJ0qVLtc0pipVX1rNj10l7/9nfvYfU5LtfTiNQ0tLS2Lp1K/Hx1qwKA2P3KisrgxxZYDgcDgrmDFbznb8cPslYE5Qat46ODt577z1u3LhhH1uxYoX21otinTd7+LsfvYfHW21UPC+TLetCf2hBamoq27Zt8xvM++GHH3LhwoUgRxYYvt3Nw6kdShOUGpeWlhbee+89OjoGqwfWrl2r45yimDGGl1/fZ1ftJcTH8qef2Rw2Jenk5GS2b99uT4sEcPToUU6dOhVW7TTjUeSToM5rglKRpKamhvfff5/u7m7AmpTzvvvuY968eUGOTAXTB4fPc+jEJXv/y09tYtb0lCBGdPsSExPZvn076emDS1OcOXOGw4cP09/fH8TIJlZB3kwc3i8OV+ua7aEAoU4TlBpVWVkZ+/btw+12A4OTcfquvaOiT2VtM99/64C9v319ERtWhOcXloG5+zIzB3u7VVRUsHv3bnsds3AXHxfDHO/6UAYouxIek+dqglLD8ng8HD16lGPHjtnHBr5t+q65o6JPS1sX//uV39LntkoYubPS+MNPrg9yVHfH5XJx//33+9UKNDQ08O6779pLxoQ732q+M+VDFzcPTZqg1C26u7vZtWuXX4Nxeno6Dz74IKmpgVlwToWH7p4+/vcrv6WhxWp3iouN4avPbScuNibIkd09h8PBmjVrWL58cNXfjo4Odu7c6TcYPVwtXThY6/HRqSth0c6mCUr5aWpqYufOnX7dyHNzc/265aro5PF4+IeffGAvoyHAXzy3zW/NoXAnIhQVFbFx40acTmsWDLfbzf79+8O+88SyhTn2F4m6hlaq6kJ/qidNUAqwemRdvHiR999/316iHawxThs2bLD/WFV0Msbw/bcOcuR0hX3sj564j1WLInOdr9zcXLZv3253Qwer88SePXvszkLhJjbGxcri2fb+4ZOXgxjN+GiCUvT29nLgwAFKSkrweKzp+GNiYti0aROLFi0Km27DKjCMMfzzm/vYefCMfeyxzct46L7IHv+WlpbGQw89xMyZg4Ncr127xjvvvBO2K/Teu3yuva0JSoW8hoYGfvvb3/rVsaempvLggw/qEu0Kj8fDSz/bw3sfnrOPbVg5n2cevyd4QU2iuLg4Nm/eTFFRkX3s5s2b7Nq1i5MnT9pf6MLFyqJcXC6rNqSqrpm6htYgRzQ6TVBRqr+/n+PHj99SpTd//nwefPBBv8GLKjr19fXzD6/tYs/HZfax+1cv4M8/uyWqStUiwvLly9m0aZM9hx/A2bNn2blzp9/MKqEuIT6WFYWDk/iGeilKE1QUamxs5Le//S3nz5+3j8XExLBx40bWrFmj7U2K5tZO/ur//VcOHiu3j21ZV8h/+f3NOJ3R+bGRlZXFww8/7DfM4saNG+zcuZPTp0+HTWlq3dJ8e/vwyStBjGRsrmAHoCZPT08PJ0+e5NKlS37HZ86cybp160hKSgpSZCqUXKy8zre+v5OWtsGS9YMbFvGFJzdGVclpOAkJCWzZsoWysjJOnTpFf38/Ho+H0tJSKisrWbNmTciPE1y9eA4OhwOPx0N5VT2NLR1MT5sS7LCGpQkqChhjuHLlCidOnPAbGe9yuVixYgXz5s2L+g8eZbU3/dueUn72Hx/j9g7CFeDZ313Pow8s0feIl4hQWFhIVlYWhw8fpqnJWka9ra2NDz74gLy8PJYvX05CQmjO6J6cFM/i+VmcumC1O3906gq/s2lJkKManiaoCHft2jWOHz9+Sz15VlYWq1ev1lKTAqC2/gYv/WwPZT5LgiclxPEXf7idZQt1WqvhpKSksH37dsrKyigtLbWnA6uoqKC6uprCwkKKiopwuULvY/be5XPtBLX3yAUeuX9xSH4BCb07pyZEY2MjpaWlXLvmP3NxYmIiq1atIjs7OyTfkGpy9fT28W97Snlr51F76iKAvOzpfO0Pt5OZoTOHjGagNDV79myOHTvG1atXAWtw7+nTpykvL2fRokXMmzcvpNp21y7N4/u/PEB/v4dLVxs4eraK1SE4pk0TVISpr6/n9OnTt4zTcDqdFBYWUlxcHJLf6NTk8ng87D1ykZ/9x8c0t3baxx0OB098YiW/t32F3R1ZjS0xMZGNGzdSV1fH8ePHaW21um93d3dz9OhRzpw5Q2FhIQsWLAiJv7+pyYl8Yn0xv91/GoCf/fvHrCqeHXJfWoN/p9Rd6+/vp6qqiosXL9r14b7mzp3L0qVLQ7ZOXE2e3j43e49c4N92n6Km3r/ad3bmNP7ss1vIz5kepOjCX2ZmJrNmzeLKlSucOnWKmzdvAlaiOnHiBGfPnmXu3LkUFBQEvXr99z6xkg8On6e3z01lbRMHj19i48r5QY1pKE1QYay1tZUrV65w+fLlYZcFyMvLo7i4WCd4VdQ1tLLnyAXePXiWto6bfo+lJifw1MNr2HpPYdR2IZ9IIsLcuXOZM2cO5eXlnDt3zk5Uvb29nD9/nvPnz5OVlcW8efPIysrC4Zj8+56Wksijm5bwq/ePA/D620e4d9nckHoPaIIKM52dnVy9epWKigpaWm6d7NHhcNiJSQfbRrfGlg5KTley50gZFyvrb3k8Pi6GRx9Yyie3Lic+LvxnIw81TqeThQsXMn/+fK5cucLZs2fp7BysTq2traW2tpbY2Fhmz57N7NmzycjImNRk9fjW5bxz4Axd3b3eLzFlbL2naOwLJ4mE8+y8vlavXm1KSkqCHcaE83g8NDU1UVdXR3V1tV23PVRiYiLz589n3rx5Out4lOq82cOFinrOXKzh6Nkqquqahz0vfWoSv7NpKdvuLSQpIW7Yc9TEM8ZQW1vLhQsXbum8NCA2NpasrCyysrKYMWPGpFTLv/XuMX7+Hx8DVqnqm3/+uwFfGVlEjhpjVo95niao0NLd3U1zczNNTU00NDTQ2Ng44tLTDoeDrKws8vPzg1ZNoIKjvbObqrpmKmubqKhp4mJlPVfrmhnpr9nhcLCyKJdNawtYuzhPO0AEWVtbG1euXKGiosJvqrGhUlJSmDFjBunp6UybNo3U1NQJ78jQ3dPHl1/8mV31OzU5kf/xpUcC2hYZEglKRB4C/hFwAt83xvzNkMfjgJ8Aq4Am4NPGmArvY38JfB7oB/7MGLNztNcKpwTV399PZ2cnnZ2ddHR00NbWRmtrK21tbXZd9UgcDgczZ84kJyeH2bNnExsbO0lRq8lgjKHzZi+tHTdp7+impa2L5tZOmm500Hijk2sNrVxrbKPz5thLkTudDhbNy2Lt0jw2rJhHyhTtJBNqjDE0NDRQWVlJbW3tqMkKrGrDlJQUUlNTSUlJITk5maSkJJKSkoiLi7vj5HXsbBXfenWnPUA7IT6W//b5B1lSkH1HzzeWoCcoEXECF4DtQDVwBHjaGHPW55w/AZYaY74kIk8BnzTGfFpEioGfA2uBLOB9oMAYM3xRgrtLUOfLr9LQPPKsvr73aGDb4/GAAYPB4/Hg8XgwxuB29+Px9NPf34+7z01fXx997j76+vro6emhp6fHHtA3XgkJCaSnp5ORkUFGRkZIdFMNJ0Pf48O95Y0xfqUPYwzGDPx/8DGPx9rq7/dg8J7jMfR7PHi8/3e7PfR7rJ8+dz9utwe3u5/evn563W56+/rp6XXT09tHd08fN3v6uNndS0dXD103e/Hc4d+kAHOyp1M0dxZLF+awtCBb25bCiDGGlpYWampquH79Ok1NTbc1v5/D4SA+Pp74+Hji4uKIjY0lJiaG2NhYXC4XLpcLp9OJ0+nE4XDYPyKCiHChsoGXXj/Aze4++zkzM1JZNG8mC2ZnEBfrwul04HI6cDocrF66gNg7XEl5vAkqkJ90a4FyY8xlb0CvA48DZ33OeRz4hnf7LeAlsb4CPA68bozpAa6ISLn3+Q4FItDv/exdDp+qCMRT3zYRISYmhpiYWGJjY4iNjcPlcgPtQEWQo1OhIsblJHtmGnnZ6eRlpZOXnc6COTM0IYUxEWHatGlMmzaNJUuW4Ha7aWxspLGxkebmZpqbm0etYfF4PHR1dY1ZChvN/YUJ/OZAA53d1pfoxsZGSs9dGvbcf3vla2RMC2wP4UAmqGzgqs9+NbBupHOMMW4RaQXSvccPD7n2lrKmiHwR+CLA7Nmzhz4c0qxvMi5cLqf97cblisHlcobcYDk1ueJiY0idEk/KlARSpySQnpbEtNQk0lOTmDk9hcyMVNJSEvV9EuFcLhezZs1i1qxZ9rHu7m6/JoGBZoKuri76+vpGebbxmZ4az39+II89J65RWd9p1xgMG98kdEcPZIIa7q9n6G870jnjuRZjzCvAK2BV8d1ugANmZ6bTeKNjzPPEZ8v6bLD+L9Z/cIgDcVj/dzgEh8OJ0+mw/u9y4vImJGvKE/1wCYTxfmb7frj7XiMiiM+/jcMh9uMOcdjnOxyCQwRxOBABp8Nh/VuL4HQKLqcTp8OBy+XA5XIS4/2JjbF+4mJiiItzER8bQ0JcDAkJsUxJiCMxPpaYGO3AoIY3UIU33Izpbreb7u5u+6evr4/e3l76+vpwu91Ws4P3/wPNEgNNEwM/ABnGMC8vmz53P1XX27lS20pTa7e32tpYP/2G2JjANzUE8hWqgVyf/RygdoRzqkXEBaQCzeO8dsJ89Qu/G6inVkqpSeFyuZgyZQpTpoTm0hl3IpBltCPAAhHJF5FY4Clgx5BzdgDPerefAHYZK43vAJ4SkTgRyQcWAB8HMFallFIhJmAlKG+b0vPATqxu5j8wxpwRkReBEmPMDuBV4DVvJ4hmrCSG97w3sTpUuIE/Ha0Hn1JKqcijA3WVUkpNqvF2M9epB5RSSoUkTVBKKaVCkiYopZRSIUkTlFJKqZCkCUoppVRIiphefCLSAFTe5dNMBxonIJxwp/fBovfBovfBovdh0N3eiznGmIyxToqYBDURRKRkPF0fI53eB4veB4veB4veh0GTdS+0ik8ppVRI0gSllFIqJGmC8vdKsAMIEXofLHofLHofLHofBk3KvdA2KKWUUiFJS1BKKaVCkiYopZRSISkqE5SIPCQiZSJSLiIvDPN4nIi84X38IxHJm/woA28c9+GrInJWRE6JyAciMicYcQbaWPfB57wnRMSISER2NR7PfRCR/+x9T5wRkZ9NdoyTYRx/F7NFZLeIHPf+bTwSjDgDTUR+ICL1InJ6hMdFRP7Je59OicjKCQ/Cd7nfaPjBWpvqEjAXiAVOAsVDzvkT4Hve7aeAN4Idd5Duw2Yg0bv95Wi9D97zkoF9wGFgdbDjDtL7YQFwHEjz7s8IdtxBug+vAF/2bhcDFcGOO0D34n5gJXB6hMcfAX4LCHAP8NFExxCNJai1QLkx5rIxphd4HXh8yDmPAz/2br8FbBURmcQYJ8OY98EYs9sY0+XdPQzkTHKMk2E87weA/wX8H6B7MoObROO5D18AvmOMaQEwxtRPcoyTYTz3wQAp3u1UoHYS45s0xph9WAvJjuRx4CfGchiYKiKZExlDNCaobOCqz36199iw5xhj3EArkD4p0U2e8dwHX5/H+rYUaca8DyKyAsg1xvz7ZAY2ycbzfigACkTkoIgcFpGHJi26yTOe+/AN4LMiUg28DfyXyQkt5NzuZ8htC9iS7yFsuJLQ0L724zkn3I37dxSRzwKrgU0BjSg4Rr0PIuIAvg08N1kBBcl43g8urGq+B7BK0/tFZLEx5kaAY5tM47kPTwM/Msb8nYjcC7zmvQ+ewIcXUgL+ORmNJahqINdnP4dbi+j2OSLiwirGj1bUDUfjuQ+IyDbgvwOPGWN6Jim2yTTWfUgGFgN7RKQCq659RwR2lBjv38W/GmP6jDFXgDKshECFiRAAAAXxSURBVBVJxnMfPg+8CWCMOQTEY02eGm3G9RlyN6IxQR0BFohIvojEYnWC2DHknB3As97tJ4BdxtsqGEHGvA/eqq1/xkpOkdjeAGPcB2NMqzFmujEmzxiTh9UW95gxpiQ44QbMeP4ufoPVcQYRmY5V5Xd5UqMMvPHchypgK4CIFGElqIZJjTI07ACe8fbmuwdoNcbUTeQLRF0VnzHGLSLPAzuxeuz8wBhzRkReBEqMMTuAV7GK7eVYJaenghdxYIzzPvwtMAX4hbePSJUx5rGgBR0A47wPEW+c92En8AkROQv0A/+PMaYpeFFPvHHeh78A/kVEvoJVpfVcBH6BRUR+jlWdO93b3vY/gRgAY8z3sNrfHgHKgS7gDyc8hgi8r0oppSJANFbxKaWUCgOaoJRSSoUkTVBKKaVCkiYopZRSIUkTlFJKqZCkCUpFJBHpF5ETInJaRH4hIom3eX3HOM75hojUeF/nrIg8fQdxDvs6IvLfvTOGn/I+/7rbfe47iOVF78Ds27lmhYh8f4xznheRCe+CrCKfdjNXEUlEOowxU7zbPwWOGmP+3udxwXr/Dzs9je/1o7zGN4AOY8z/FZEFwFEg3RjTdydx+hy7F/h74AFjTI93UGysMSbkJiUVkV8A3zTGnBzlnETgoDFmxeRFpiKBlqBUNNgPzBeRPBE5JyLfBY4BuSLytIiUekta3/K9SET+TkSOibUWVsZoL2CMuYg1WDHNe+08EXlHRI6KyH4RKfQezxeRQyJyRET+1whPlwk0DkwtZYxpHEhOIrLVuw5RqXe9njgRSRVr/aKF3nN+LiJfGPK7rBWRX3m3HxeRmyISKyLxInLZe/xHIvKEd7tCRP7a+/uXDsQ/5DmTgaUDyUmstYG+7t1+UET2iYjDOyN+hYisHe0eKjWUJigV0cSaS/FhoNR7aCHWEgErgD7gW8AWYDmwRkR+13teEnDMGLMS2Is1in6011kJXPSZEuoV4L8YY1YBXwO+6z3+j8DLxpg1wLURnu5drOR5QUS+KyKbvK8RD/wI+LQxZgnWTDBfNsa0As8DPxKRp7DWa/qXIc95DBgowdwHnAbWAOuAj0aIo9H7+7/s/R2GWu19ngEvAJ8Wkc3APwF/6FNCLfG+rlLjpglKRaoEETmB9cFYhTV9FUCld+0asD6g9xhjGrzLqvwUa5E2AA/whnf7/wM2jvA6XxGRMqwP+W8AiMgUYD3WFFEnsOYzHFgnZwPwc+/2a8M9oTGmA1gFfBFrjrc3ROQ5rOR6xRhzwXvqjwfiNca8h5WEvwP80TDP6QbKvXPHrcWqQrwfK2nsH+F3+5X3/0eBvGEez8RnDjpvSekLwHvAS8aYSz7n1gNZI7yOUsOKurn4VNS4aYxZ7nvAO59gp++h23i+kRprv+1tg/oU8BMRmYf1xe/G0Ncfx3MNnmBMP7AHaxb1UqzJi0+MdL5Yy4IUATeBaVgzTQ+1H6s02Qe8j1UaczJ86QhgYPb6fob/rLiJNVGqryVAE7cmo3jv+UqNm5agVDT7CNgkItNFxIm1zs9e72MOrJnsAT4DHBjtiYwxv8IqrT1rjGkDrojIk2B1yBCRZd5TDzI4+fDvD/dcIrLQ2+liwHKgEjgP5InIfO/xP/CJ9yvAOe/v8AMRiRnmqfcB/xU4ZIxpwFqEsxA4M9rvNopzwEAsiMgcrIlUVwAPD+l5WIB/daBSY9IEpaKWd2mAvwR2Ayex2pz+1ftwJ7BIRI5itVG9OI6nfBH4qrc08/vA50XkJFYCGFg2/M+BPxWRI1jrjA1nCvBjb9f1U0Ax8A1jTDfWjNG/8JaqPMD3RKQAq1rvL4wx+7ES0f8Y5nk/AmZ6Hwc4BZy605m4jTHngVQRSfb2inwV+Jq3Q8fnge97283Aqtp8/05eR0Uv7WaulLpjYi050W6MGXEslFjrin3VGPMHkxeZigRaglJK3Y2XGWyrGsl04K8mIRYVYbQEpZRSKiRpCUoppVRI0gSllFIqJGmCUkopFZI0QSmllApJmqCUUkqFpP8f3YkN6zxSaowAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "beta.Update((15, 0))\n", "posterior = beta.MakePmf()\n", "\n", "thinkplot.Pdf(prior, color='gray', label='prior')\n", "thinkplot.Pdf(posterior, label='posterior')\n", "thinkplot.decorate(xlabel='Prob Red Sox win (x)',\n", " ylabel='PDF')" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.8000000040157967" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" } ], "source": [ "posterior.Mean()" ] }, { "cell_type": "code", "execution_count": 29, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.83" ] }, "execution_count": 29, "metadata": {}, "output_type": "execute_result" } ], "source": [ "posterior.MAP()" ] }, { "cell_type": "code", "execution_count": 30, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0.66, 0.91)" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" } ], "source": [ "posterior.CredibleInterval()" ] }, { "cell_type": "code", "execution_count": 31, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.94" ] }, "execution_count": 31, "metadata": {}, "output_type": "execute_result" } ], "source": [ "x = posterior.Random()" ] }, { "cell_type": "code", "execution_count": 32, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "7" ] }, "execution_count": 32, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.sum(np.random.random(7) < x)" ] }, { "cell_type": "code", "execution_count": 33, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "6" ] }, "execution_count": 33, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def simulate(k, dist):\n", " x = dist.Random()\n", " return np.sum(np.random.random(k) <= x)\n", " \n", "simulate(7, posterior)" ] }, { "cell_type": "code", "execution_count": 34, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD8CAYAAACb4nSYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAEpZJREFUeJzt3X+MXedd5/H3BxunS7qFQEaryvbUDpgKl64SdnC0iggrmh+Oiuz80QhbKgqoyHQ3ZouiFaRLlGhNK0GR2P2jZonVeFV+NCZNltUIGUJEWqCCtB4noVk7mE5MqWddNqHOku22JDj97h9zkr2+HeeeGY99b/Z5v6Qrn+c5z3Pu91rWZ46fc+6ZVBWSpDZ8y7gLkCRdOoa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kN6RX6SbYnOZFkPsldS+x/f5KnkzyV5DNJtnb9m5J8vet/Ksmvr/YHkCT1l1HfyE2yBvgr4EZgATgC7K6q4wNj3lJVL3bbO4B/U1Xbk2wCfq+qvr9vQVdeeWVt2rRpmR9Dktp29OjRv6uqqVHj1vY41jZgvqpOAiQ5BOwEXgv9VwO/czmw4mc7bNq0ibm5uZVOl6QmJfmbPuP6LO+sB04NtBe6vuE3vCPJs8BHgH87sGtzkieT/HGSH+pTlCTp4ugT+lmi75vO5Ktqf1V9N/DzwN1d95eB6aq6BrgT+ESSt3zTGyR7kswlmXv++ef7Vy9JWpY+ob8AbBxobwBOv874Q8CtAFX1UlV9pds+CjwLfO/whKo6UFUzVTUzNTVySUqStEJ9Qv8IsCXJ5iTrgF3A7OCAJFsGmu8GvtD1T3UXgklyFbAFOLkahUuSlm/khdyqOptkL/AIsAY4WFXHkuwD5qpqFtib5AbgH4EXgNu76dcD+5KcBV4B3l9VZy7GB5EkjTbyls1LbWZmprx7R5KWJ8nRqpoZNc5v5EpSQwx9SWqIoS9JDenzjVxJumB7P/TAuEt4zUfv3j3uEsbGM31JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ3pFfpJtic5kWQ+yV1L7H9/kqeTPJXkM0m2Duz7YDfvRJKbV7N4SdLyjAz9JGuA/cAtwFZg92Codz5RVe+sqquBjwC/2s3dCuwC3gFsB36tO54kaQz6nOlvA+ar6mRVvQwcAnYODqiqFwealwPVbe8EDlXVS1X118B8dzxJ0his7TFmPXBqoL0AXDs8KMkdwJ3AOuBHBuY+PjR3/RJz9wB7AKanp/vULUlagT5n+lmir76po2p/VX038PPA3cuce6CqZqpqZmpqqkdJkqSV6BP6C8DGgfYG4PTrjD8E3LrCuZKki6hP6B8BtiTZnGQdixdmZwcHJNky0Hw38IVuexbYleSyJJuBLcDnLrxsSdJKjFzTr6qzSfYCjwBrgINVdSzJPmCuqmaBvUluAP4ReAG4vZt7LMmDwHHgLHBHVb1ykT6LJGmEPhdyqarDwOGhvnsGtj/wOnM/DHx4pQVKklaP38iVpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1Jakiv0E+yPcmJJPNJ7lpi/51Jjif5fJI/SvK2gX2vJHmqe82uZvGSpOVZO2pAkjXAfuBGYAE4kmS2qo4PDHsSmKmqryX518BHgB/r9n29qq5e5bolSSvQ50x/GzBfVSer6mXgELBzcEBVfaqqvtY1Hwc2rG6ZkqTV0Cf01wOnBtoLXd/5vA/4/YH2m5LMJXk8ya0rqFGStEpGLu8AWaKvlhyYvBeYAX54oHu6qk4nuQp4LMnTVfXs0Lw9wB6A6enpXoVL0mrY+6EHxl3Caz569+6L/h59zvQXgI0D7Q3A6eFBSW4AfgHYUVUvvdpfVae7P08CnwauGZ5bVQeqaqaqZqamppb1ASRJ/fUJ/SPAliSbk6wDdgHn3IWT5BrgPhYD/7mB/iuSXNZtXwlcBwxeAJYkXUIjl3eq6mySvcAjwBrgYFUdS7IPmKuqWeBXgDcDn0wC8KWq2gF8H3Bfkm+w+APml4bu+pEkXUJ91vSpqsPA4aG+ewa2bzjPvD8D3nkhBUqSVk+v0Jf0xtDaRUktn49hkKSGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWpIr9BPsj3JiSTzSe5aYv+dSY4n+XySP0rytoF9tyf5Qve6fTWLlyQtz8jQT7IG2A/cAmwFdifZOjTsSWCmqv458BDwkW7udwL3AtcC24B7k1yxeuVLkpajz5n+NmC+qk5W1cvAIWDn4ICq+lRVfa1rPg5s6LZvBh6tqjNV9QLwKLB9dUqXJC1Xn9BfD5waaC90fefzPuD3VzhXknQRre0xJkv01ZIDk/cCM8APL2dukj3AHoDp6ekeJUmSVqLPmf4CsHGgvQE4PTwoyQ3ALwA7quql5cytqgNVNVNVM1NTU31rlyQtU5/QPwJsSbI5yTpgFzA7OCDJNcB9LAb+cwO7HgFuSnJFdwH3pq5PkjQGI5d3qupskr0shvUa4GBVHUuyD5irqlngV4A3A59MAvClqtpRVWeS/CKLPzgA9lXVmYvySSRJI/VZ06eqDgOHh/ruGdi+4XXmHgQOrrRASdLq8Ru5ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIb1CP8n2JCeSzCe5a4n91yd5IsnZJO8Z2vdKkqe61+xqFS5JWr61owYkWQPsB24EFoAjSWar6vjAsC8BPwH8uyUO8fWqunoVapUkXaCRoQ9sA+ar6iRAkkPATuC10K+qL3b7vnERapQkrZI+yzvrgVMD7YWur683JZlL8niSW5dVnSRpVfU5088SfbWM95iuqtNJrgIeS/J0VT17zhske4A9ANPT08s4tCRpOfqc6S8AGwfaG4DTfd+gqk53f54EPg1cs8SYA1U1U1UzU1NTfQ8tSVqmPqF/BNiSZHOSdcAuoNddOEmuSHJZt30lcB0D1wIkSZfWyNCvqrPAXuAR4Bngwao6lmRfkh0ASX4wyQJwG3BfkmPd9O8D5pL8BfAp4JeG7vqRJF1Cfdb0qarDwOGhvnsGto+wuOwzPO/PgHdeYI2SpFXiN3IlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkN6PVpZ0rn2fuiBcZcAwEfv3j3uEvQG45m+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kN6RX6SbYnOZFkPsldS+y/PskTSc4mec/QvtuTfKF73b5ahUuSlm9k6CdZA+wHbgG2AruTbB0a9iXgJ4BPDM39TuBe4FpgG3BvkisuvGxJ0kr0OdPfBsxX1cmqehk4BOwcHFBVX6yqzwPfGJp7M/BoVZ2pqheAR4Htq1C3JGkF+oT+euDUQHuh6+uj19wke5LMJZl7/vnnex5akrRcfUI/S/RVz+P3mltVB6pqpqpmpqameh5akrRcfUJ/Adg40N4AnO55/AuZK0laZX1C/wiwJcnmJOuAXcBsz+M/AtyU5IruAu5NXZ8kaQxGhn5VnQX2shjWzwAPVtWxJPuS7ABI8oNJFoDbgPuSHOvmngF+kcUfHEeAfV2fJGkMej1Pv6oOA4eH+u4Z2D7C4tLNUnMPAgcvoEZJ0irxG7mS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDekV+km2JzmRZD7JXUvsvyzJ73T7P5tkU9e/KcnXkzzVvX59dcuXJC3H2lEDkqwB9gM3AgvAkSSzVXV8YNj7gBeq6nuS7AJ+Gfixbt+zVXX1KtctSVqBPmf624D5qjpZVS8Dh4CdQ2N2Ah/vth8C3pUkq1emJGk19An99cCpgfZC17fkmKo6C/w98F3dvs1Jnkzyx0l+6ALrlSRdgJHLO8BSZ+zVc8yXgemq+kqSfwH8tyTvqKoXz5mc7AH2AExPT/coSZK0En3O9BeAjQPtDcDp841Jshb4duBMVb1UVV8BqKqjwLPA9w6/QVUdqKqZqpqZmppa/qeQJPXSJ/SPAFuSbE6yDtgFzA6NmQVu77bfAzxWVZVkqrsQTJKrgC3AydUpXZK0XCOXd6rqbJK9wCPAGuBgVR1Lsg+Yq6pZ4H7gN5PMA2dY/MEAcD2wL8lZ4BXg/VV15mJ8EEnSaH3W9Kmqw8Dhob57Brb/AbhtiXkPAw9fYI2SpFXSK/Slcdj7oQfGXcJrPnr37nGXIK0KH8MgSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTfkauJ+V20/h5a6eLzTF+SGtIr9JNsT3IiyXySu5bYf1mS3+n2fzbJpoF9H+z6TyS5efVKlyQt18jQT7IG2A/cAmwFdifZOjTsfcALVfU9wH8EfrmbuxXYBbwD2A78Wnc8SdIY9DnT3wbMV9XJqnoZOATsHBqzE/h4t/0Q8K4k6foPVdVLVfXXwHx3PEnSGPS5kLseODXQXgCuPd+Yqjqb5O+B7+r6Hx+au37F1b6BTcrFUvCCqdSyVNXrD0huA26uqp/q2j8ObKuqnxkYc6wbs9C1n2XxjH4f8OdV9Vtd//3A4ap6eOg99gB7uubbgROr8NkuxJXA3425hvOZ1NomtS6Y3NomtS6Y3NomtS4Yf21vq6qpUYP6nOkvABsH2huA0+cZs5BkLfDtwJmec6mqA8CBHrVcEknmqmpm3HUsZVJrm9S6YHJrm9S6YHJrm9S6YLJrG9RnTf8IsCXJ5iTrWLwwOzs0Zha4vdt+D/BYLf4XYhbY1d3dsxnYAnxudUqXJC3XyDP9bo1+L/AIsAY4WFXHkuwD5qpqFrgf+M0k8yye4e/q5h5L8iBwHDgL3FFVr1ykzyJJGqHXN3Kr6jBweKjvnoHtfwBuO8/cDwMfvoAax2FilpqWMKm1TWpdMLm1TWpdMLm1TWpdMNm1vWbkhVxJ0v8/fAyDJDXE0B8y6pET45DkYJLnkvz3cdcyLMnGJJ9K8kySY0k+MO6aAJK8KcnnkvxFV9d/GHdNg5KsSfJkkt8bdy2DknwxydNJnkoyN+56BiX5jiQPJfnL7t/bv5yAmt7e/V29+noxyc+Ou67X4/LOgO4REX8F3Mji7aZHgN1VdXzMdV0PfBX4jar6/nHWMizJW4G3VtUTSf4pcBS4dQL+zgJcXlVfTfKtwGeAD1TV4yOmXhJJ7gRmgLdU1Y+Ou55XJfkiMFNVE3cvfJKPA39aVR/r7iT8tqr6X+Ou61VdfvwP4Nqq+ptx13M+numfq88jJy65qvoTFu+KmjhV9eWqeqLb/t/AM0zAt65r0Ve75rd2r4k4w0myAXg38LFx1/JGkeQtwPUs3ilIVb08SYHfeRfw7CQHPhj6w5Z65MTYA+yNonu66jXAZ8dbyaJuCeUp4Dng0aqaiLqA/wT8HPCNcReyhAL+MMnR7pvyk+Iq4Hngv3TLYh9Lcvm4ixqyC5ic562ch6F/rizRNxFnh5MuyZuBh4GfraoXx10PQFW9UlVXs/hN8G1Jxr40luRHgeeq6ui4azmP66rqB1h8qu4d3dLiJFgL/ADwn6vqGuD/ABNxzQ2gW27aAXxy3LWMYuifq9djI3Subs38YeC3q+q/jrueYd0ywKdZfLz3uF0H7OjWzg8BP5Lkt8Zb0v9TVae7P58DfpfJeSruArAw8L+1h1j8ITApbgGeqKr/Oe5CRjH0z9XnkRMa0F0wvR94pqp+ddz1vCrJVJLv6Lb/CXAD8JfjrQqq6oNVtaGqNrH47+uxqnrvmMsCIMnl3cV4uqWTm4CJuGOsqv4WOJXk7V3Xu1j8pv+k2M0bYGkH/B255zjfIyfGXBZJHgD+FXBlkgXg3qq6f7xVveY64MeBp7v1c4B/332Le5zeCny8u6PiW4AHq2qibo+cQP8M+N3Fn+OsBT5RVX8w3pLO8TPAb3cnZCeBnxxzPQAk+TYW7/j76XHX0oe3bEpSQ1zekaSGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXk/wJpDnQVxLfB9QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sample = [simulate(7, posterior) for i in range(100000)];\n", "thinkplot.Hist(Pmf(sample))" ] }, { "cell_type": "code", "execution_count": 35, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.94716" ] }, "execution_count": 35, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.mean(np.array(sample) >= 4)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.5" } }, "nbformat": 4, "nbformat_minor": 1 }