{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline\n", "import sympy as sym\n", "from sympy import init_printing ; init_printing()\n", "from IPython.display import display, Math, Latex\n", "maths = lambda s: display(Math(s))\n", "latex = lambda s: display(Latex(s))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Define the Beta distribution (symbolically using sympy) which will act as our prior distribution" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\operatorname{Beta}(\\mu|a,b) = $" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATkAAAAcCAYAAAAa98jAAAAABHNCSVQICAgIfAhkiAAAB+VJREFUeJztnFmMFEUYx39cAoJZFBFRuTxAosgCcgXBRdEXogghKsSYFR40KJgoiUc8QCNG5YGoaEKErICCQURDCILhCMQDXVBBRUV0Ua7AKocoKiI+fNXZ3t7qc6qne6B+yaRnurqq/v19VdV19YDFYrFYLCXGM8AHWYuwWPJA46wFRGQosAzYA5wExmQrJ/eUA1+mnIf1iaUkKJVGrhVSae/LWkiJUA58kXIe1ieW04YngC3AUeAAMBdokWJ+UXsN84D9SGU8leiL2GCCT3g7FT4OWA38iTRGA1PUlEZPrtT8NxGxw/iQ68L8lwYmbJmW7iBtUW0KKdq1ETANGAx0BoYDu4GHTWfkIkqFuhr4D3jAc34M8BKwATii0lpgWmARWArsBVprwm5E7mszcC1wObAC2Ak0TUlPmE9ORvi0cV3v57888xpyH70iXBvkPzftgRPAi65zWdoyTLdpbXFs6quv0EJ/EnjS9XsnsBypWFkyHWnEXvWcfwwx2FFgF9nrTMqzwEZgMnKvbsqBf4BRiD8AHkJ6c12AH4ojUcu0gLC/XN/9/Jdn+iL38HWEa4P852YkMqW0VBOWhS2j6jalLY5N4+iLRUfkKbMV+A1pPI4DT/tcP5Xwlr4iJM+wXkM35OkwWxM2DLgM6YFWULo9OYBtSCPWxHN+IfCG59wlyL1206QzlfR94qQThSD/pUUl0e7Tj+bIg2VjjDh+/nOzAqj1XJO1LYN0m9SWxKZafd6Fh8n4j4HLlKg16ndb4DPgfGAKMAToB/yN/8rey0CPkM+n8e6pAeORRuwtTdhaYDvRHVEs4tjdYRHQCZkicFMOfO451wc4DPyoSb8YPolDkP8cxiEN+ffA78BBpCzelbo6PT2BZsAmxP5LgF+BP4BVwJWaOH7+cygDrkNWsE8k1JWGLcN0m9KWxKZafd7hah913KSJ3EeJcsJGIAsMt1HXaFQiE4jeSuZQqz5pMhwpFJ+knI9J4tjd4UN1vAFYqb6fifRUm7muawTcD8wH/tWkXwyfxCHMf2cBrwPVyNzqfuA8ZGg3V31/Ln2Z9eirjp0Rv6wA5gCDEP+sQXouh1xxdP5zMwI4A3inAF1p2DJMtyltSWwaSd8W4Bj6ubopSGM2Vv0egVSaUcClwCRk0u8QUrFM0hppzcuVhinqeyfPda2Upq0R0qwgP8PVOHZ3KFPn3b2sgch0wTfAAMQvC4Aa4ByjiqP7BKIPY6L4rzUyevDSAemJfBshHy+VFDZcna3i1wK9PWHzVNijnvM6/7lZjEz/eHcqZG3LIN0mtSWxaZg+WiAVxG8MvJD68zqNgFnIxOF+ZG5uBrAuQHhSKtDPFVV5ruumzq+KkWbWjVxcu7s5Buxz/b4HKZijgZ+R7SNL0BfkQqkgmk8geuGP4z8d3yFlMS6VFNbIVeM/3TBAhS3UhHn959ACaWQWa8LyYEs/3Sa1JbVpkD76q4iv+ETcjjRopntpJhmE3EPQHIRDBckauRqiLZU7n7D0C7H7bvRD0LwRtfBH8d/ZwOPIMOcgMuRx23tzSB41xPNfVUh6zZB56Br0PfHOKh3dCqmf/25WccZpwvJgSz/dprQVYtMG+twJBM0LlSErdBvI36S9m2PqmOZm5B3UXwYPY09IeCF2b0ndPZ8KhPnvKuTp3x4ZjixCVvWPA12BOwl/nW0m9fdqgQyzRyLzUzWesLA3R3oic2fL0Ff8zuq4UxPm579RyMri8pC8g0jTloWWuzBthdgUPPqiNnKDkZ5E2FMya5zuddsU87jecHpJ7d4Yqaw/GdaTJWH+m4/c8zAaTos8pY7VIXnM1JyrRBq5Kk26YTj+q/EJH6WO3qGZn/+aADchE+uHY2pxk5YtTZS7MG1JbQoafe4tJE7CuzURb1dHXUXME3uRV8u6Zy0kBknt3h1pANN+R7WYBPmvI9L7WEfDStmGui0PxS6jziqgbmGnA3A3MuXgrZB+/huKVH6/oVhU0rKliXIXVk+T2lSrz2nkmlG37+QaT6RbgTvU9yirlllyElgPnIusLOadQuzuvIu6NgVdWRHkP2eK4GLqb5Fpi8ztXIQMbYrd6DsVciz1379sjawCNkf2QXqHXX7+G43si3yvQF1p2dJEuQurp0ltGqivt8p4H2KAN5F3PNcjqzz7VPgS0n3R2wRjEa33asJuQYYkVcD76rodrnMziqDPTSF2X4g4uWOxxBZAnJ3wQf5bTd2E+PPIkKsWWYU8QfK/l6ok2epqU2TuZzOyqr0deAHZabALaawm+8T1898vyBysH1nbMqjcmdBWiE0D9U1QGU5EKnotsrN4DfIS7YPqdzVwQcSbyIozkMZBtyVjKsEraTVFUVhHUruXIQXh3WKKLYA4hT/If+2QJ/kB5CHwETJB7jws5iTUV0myRq6XijcbuBDZuHtEfVYGpOfnv34qvaCX6bO0ZVi5M6EtqU1D9c1SCfePKDDvPILcj3cTYd5IavdJKt4Q44ryQan4Lyl+/puuznc1mJdJW5oud6b9HKjvY6SL19JQZlnTAlleXpa1kBCS2L0lsi3l7VQU5YNS8V8Sgvy3DfNziqZsmUa5M+nnQH1NkCHRVwYyyhNDkb+ByuufLia1ew9k2N3FsJ68kXf/JSUL/5mwZVq6Tfk5UN8V5OP1ptMNa3eLxWKxWCwWi8VisVgsFovFYrFYLBaLpZT4H0YXK1vLOFD5AAAAAElFTkSuQmCC\n", "text/latex": [ "$\\displaystyle \\frac{\\mu^{a - 1} \\left(1 - \\mu\\right)^{b - 1} \\Gamma\\left(a + b\\right)}{\\Gamma\\left(a\\right) \\Gamma\\left(b\\right)}$" ], "text/plain": [ " a - 1 b - 1 \n", "μ ⋅(1 - μ) ⋅Γ(a + b)\n", "────────────────────────────\n", " Γ(a)⋅Γ(b) " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "a,b,mu = sym.symbols('a b mu')\n", "G = sym.gamma\n", "\n", "# The normalisation factor\n", "Beta_norm = G(a + b)/(G(a)*G(b))\n", "# The functional form, note the similarity\n", "# to the Binomial likelihood\n", "Beta_f = mu**(a-1) * (1-mu)**(b-1)\n", "\n", "# Turn Beta into a function\n", "Beta = sym.Lambda((mu,a,b), Beta_norm * Beta_f)\n", "\n", "maths(r\"\\operatorname{Beta}(\\mu|a,b) = \")\n", "display(Beta.expr)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, '$\\\\mu$')" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEQCAYAAACgBo8fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3iTVfvHPyfdLW2hpWUVKKvsUjZlyVRQhggK8iqCKIp7i6/i1lffn+ALLkARcAGCiIiAirIpQkF2KXuU1QHdTZo05/fH06JC6UzyJOn5XFeups2Tc75dd+7c5x5CSolCoVAoXB+D3gIUCoVCYRuUQVcoFAo3QRl0hUKhcBOUQVcoFAo3QRl0hUKhcBM89dq4Zs2aMjIyUq/tFQqFwiXZuXNnqpQyrLjHdDPokZGRxMfH67W9QqFQuCRCiFPXe0yFXBQKhcJNUAZdoVAo3ARl0BUKhcJN0C2GXhxms5mkpCSMRqPeUqo8vr6+RERE4OXlpbcUhUJRRpzKoCclJREYGEhkZCRCCL3lVFmklKSlpZGUlESjRo30lqNQKMqIU4VcjEYjoaGhypjrjBCC0NBQ9U5JoXAxnMqgA8qYOwnq96BQuB5OZ9DLg8lSwIUMI/mWAr2lKBQKRZmYsfYIm46k2GVtlzboZoskOctIvkX/nu4mk4nRo0fTtGlTunbtysmTJ4u97t577yU8PJw2bdpcd63x48ezdOnSCunYuXMnbdu2pWnTpjz22GMU1+/+0KFDxMbG4uPjw3vvvVehfRQKRcX44PcjxB1Ls8vaLm3Qi6ICEv0N+ty5c6lRowZHjx7lySef5Pnnny/2uvHjx7NmzRq76Zg8eTJz5szhyJEjHDlypNi9QkJCmDlzJs8884zddCgUimuxFFixWCU+nh52Wd+lDbqh0KJbbWjPb731Vjp27Ejr1q2ZM2dOmZ/3ww8/cM899wAwatQofvvtt2K94969exMSElLqemvXrqVXr15ERUWxcuXKMmk4f/48mZmZxMbGIoRg3LhxLF++/JrrwsPD6dy5s0pJVCgcjMliBcDXyz6m16nSFv/Oaz8e4OC5zBKvkVKSm1+Aj5cHnobSD/Fa1Q3ilaGtS7zm888/JyQkhLy8PDp37szIkSMJDQ1l9OjRJCYmXnP9U089xbhx4zh79iz169cHwNPTk+DgYNLS0qhZs2apuorj5MmTbNiwgWPHjtG3b1+OHj3KqVOnGD16dLHXr1+/nrNnzxIREXHlaxEREZw9e7ZC+ysUCttjNGvnfb5e9vHQndagl40iIy7/dr9yzJw5k++//x6AM2fOcOTIEUJDQ1m8eHGJzyvOG69Mpsgdd9yBwWCgWbNmNG7cmEOHDhETE8Pu3bsdpkGhUNiWIg/dx7OKeeiledKgxaMOns+kbnU/albzqfSe69evZ+3atcTFxeHv70+fPn2u5GKX5qFHRERw5swZIiIisFgsZGRklCm0cj2uNsRCCBITE0v00CMiIkhKSrrytaSkJOrWrVthDQqFwrYoD70EioxecZ5pRcjIyKBGjRr4+/tz6NAhtm3bduWx0jz0YcOGsWDBAmJjY1m6dCn9+vWrlHe8ZMkS7rnnHk6cOMHx48dp3rw5vr6+JXro1atXJzAwkG3bttG1a1e++OILHn300QprUCgUtsVotm8M3cUPRbWPtjoUHTRoEBaLhejoaKZOnUq3bt3K/NyJEyeSlpZG06ZNmT59Ou+88w4A586d4+abb75y3Z133klsbCyJiYlEREQwd+7cYtdr3rw5N9xwA4MHD2bWrFn4+vqWSccnn3zCfffdR9OmTWnSpAmDBw8GYNasWcyaNQuACxcuEBERwfTp03nzzTeJiIggM7Pk8wqFQlF5jIU1M/bKchG28m7LS6dOneTVAy4SEhJo2bJludbZl5RBWKA3tYP9bClPQcV+HwqF4vrEHUvjzk+38c39XenepGIJE0KInVLKTsU95tIeOmheui3TFhUKhcJeFHno9oqhu7xBF0LYLIauUCgU9sRkLgq56BRDF0L4CiG2CyH2CCEOCCFeK+YaHyHEYiHEUSHEH0KISHuILQ7loSsUClfhr8Ii/Tx0E9BPStkOiAEGCSGuPi2cCFyWUjYF3gfeta3M66M8dIVC4SrYO22xVIMuNbILP/UqvF1tQYcDCwrvLwX6CwdVtAjloSsUChehKG1Rt5ALgBDCQwixG0gGfpVS/nHVJfWAMwBSSguQAYQWs84kIUS8ECI+JcU27SMNQmBVHrpCoXABTM5wKCqlLJBSxgARQBchxNW9X4vzxq+xslLKOVLKTlLKTmFhYeVXWwwGAc5gz0+fPk3fvn1p37490dHRrFq16ppr1q9fz5AhQyq1z9KlSxFCcHXKp0KhcH6uFBbp6aEXIaVMB9YDg656KAmoDyCE8ASCgUs20FcqQgisTtA+98033+SOO+7gzz//ZNGiRTz00EM23yMrK4uZM2fStWtXm6+tUCjsj9FcgIdB4OmhX5ZLmBCieuF9P2AAcOiqy1YA9xTeHwX8Lh10UmlrD72i7XOFEFeqLTMyMq7bQyUzM5MRI0bQqlUrHnzwQaxWa5n3mDp1Ks8991yZq0YVijJhzICzu2D/d7BlBuz+Bk5vg+xkvZW5HSaL1W7eOZStl0sdYIEQwgPtBeBbKeVKIcTrQLyUcgUwF/hSCHEUzTMfU2llq6fAhX2lXlbLUoDVKsG7DN9K7bYw+J0SL6lo+9xXX32VG2+8kQ8++ICcnBzWrl1b7Prbt2/n4MGDNGzYkEGDBrFs2TJGjRpV6vp//vknZ86cYciQIWrKkKLyWAvgyC8QPw+O/gryOo5FRGfo8gC0Gg6e3o7V6IYYzQV2i59DGQy6lHIv0L6Yr7/8t/tG4HbbSisbgmKC9ZWgou1zFy5cyPjx43n66aeJi4vj7rvvZv/+/RgM/3w17tKlC40bNwa0vi6bN29m1KhRJa5vtVp58sknmT9/fuW+OYVCSji4HH6ZChlnoFot6PE41OsIIY0hOAJyUuHScbi4H3Z9Ccvug19ehH5Tof1df40KU5Qbo9lqtwwXcOZui6V40kWkpueRnptP67rBld6yMu1z586de2XcW2xsLEajkdTUVMLDw/9xfXFtcUtbf/jw4ezfv58+ffoAWnOtYcOGsWLFCjp1Kralg0JxLVkXYdXTkPAj1ImBm96G5oPB46rJVb7BENoEmg2E7o/Dsd9g0zRY8QgcXQtD/wd+NfT5Hlwck0VnD93ZsWUMvTLtcxs0aMBvv/3G+PHjSUhIwGg0Ulwmz/bt2zlx4gQNGzZk8eLFTJo0qUzrp6amXrnfp08f3nvvPWXMFWXnyFr4biKY82DAaxD7CHiU4d/fYNAMe5P+sHUG/P4mnN0Jty+AiI721+1mGM1WfOxo0F2/lwtaHrotzmAr0z532rRpfPrpp7Rr144777yT+fPnF9sPPTY2lilTptCmTRsaNWrEiBEjKq1boSiRP7+Gb+6A4PoweQv0fKJsxvzvGAzQ80m49xcQBvhiOJzZYR+9bozJUlBFQy5lpKgnui2G0Pn4+LB69eoKPbdVq1Zs2bKlxGv69OlzJWxSGdavX1/pNRRVAClh03uaV924D9zxJfgGVW7NiI5w7xqYdzN8dRuMW67F3xVlwmS22m24BbiDh17oBatqUYXiKn5/U7tFj4axSypvzIsIqgvjV2px9C9HwLnrT9FS/BOjnWPoLm/Qr3joyp4rFH/xxxzNO+8wDm6dZfuUw+AIzaj7BMHCO7XMGEWpGM0F+NppWhG4gUG39VxRhcLlObAcVj8HzW+GW97X4t/2oHoDGPM15KbBsvu13HZFiZgsVnxUyOX62HquqELh0pzaqhnX+l1g1OflP/wsL3XaweB34djvsFEVvJWG8tBLwaBi6AqFRtYF+PYezXO+cxF4OWjObsfxWpx+/X/g2DrH7OmiGNWhaMkIFUNXKKDAAksngilLy2bxD3Hc3kLAkPehZhQsf0jToCgWk6VA5aGXhLN46NOnT6dVq1ZER0fTv39/Tp06dc01J0+epE2bqzsPl40ffviB6OhoYmJi6NSpE5s3b66sZIU7se4tOLVZq+Ks1crx+3sHwPCPIOs8/P6W4/d3AaSUmoduxzx0lzfozuKht2/fnvj4ePbu3cuoUaN47rnnbLp+//792bNnD7t37+bzzz/nvvvus+n6Chfm8C+weTp0uAfaVb4vXoWp3xk63QvbZ8O5P/XT4aQUzRNVHnoJ2NpDr2j73L59++Lv7w9At27dSEpKKvY6i8XCPffcQ3R0NKNGjSI3N7dM61erVu1KRk9OTk6xVaiKKkjuJfjhYajVBgb/V2810P9lCAiDH59QWS9XccWgV8VK0Xe3v8uhS1e3Xb8Wq5Tk5Rfg4+mBp0fJRq5FSAue7/J8iddUtH3u35k7dy6DBw8udv3ExETmzp1Ljx49uPfee/n444955plnePLJJ1m37toDpTFjxjBlyhQAvv/+e1544QWSk5P56aefSvw+FFWE1c9B3iW4exl4OUGffL/qMOg/sPRe2P4pdHtQb0VOg8nOA6LBiQ16WRGFBf/SJsX/FW+fW8RXX31FfHw8GzZsKPbx+vXr06NHDwDuuusuZs6cyTPPPMP7779f6tojRoxgxIgRbNy4kalTp16357qiipDwI+xbAn3+rfX6dxZa36b1j1n3NkTf4dgDWifmyvi5qmjQS/Oki7AUWDl4PpM6wX6EBfpUas/KtM8FWLt2LW+99RYbNmzAx6d4Lddrn1sWD72I3r17c+zYMVJTU6lZs2b5v1GF65OTBiufhNrR0OspvdX8EyHgxjfgkx7aBKSBr+mtyCkwFg6IrpIhl7JSFEOXNhhzUZn2uX/++ScPPPAAa9asuaYH+t85ffo0cXFxxMbGsnDhQnr27AlQqod+9OhRmjRpghCCXbt2kZ+fT2hoaDm+O4VbsWYK5KXDuB+u7WfuDNRqDW1vhz9mQ9cHIaiO3op0x+QAD93lD0WFDStFK9M+99lnnyU7O5vbb7+dmJgYhg0bVux1LVu2ZMGCBURHR3Pp0iUmT55cpvW/++472rRpQ0xMDA8//DCLFy9WB6NVlRMbYd+3WjvbWq31VnN9+v4brGbY6ASHtU5AkYduz8Iil/fQhRAYhLBJL5fKtM8tSzw7MjKSgwcPVmj9559/nuefL1sYSuHGWPLhp2egekPnC7VcTUgjrYp053zo/qg24q4KYzQXhVyUh14iwoZTixQKp2bbx5CaqKUoOqq0vzL0fhYMXtoBaRXnr5CLKiwqEYMQuleKKhR2JyMJNryrdVFsPkhvNWUjsDZ0nQT7lkLaMb3V6MpfIRfloZeIEKrboqIK8POL2lvRQWUboO40dHtYO7jd+oHeSnSlKG3Rnlkupa4shKgvhFgnhEgQQhwQQjxezDV9hBAZQojdhbeX7SO3eAzYJoauUDgtp7fBweXaPNAaDfVWUz4Ca0G7O2H3N5CdrLca3TA5iYduAZ6WUrYEugEPCyGK6/6zSUoZU3h73aYqS0F56Aq3RkrNOw+sox0uuiLdH4OCfC2NsYpypbBIz0NRKeV5KeWuwvtZQAJQz26KKoCtslwUCqdk/3dwNh76TdW6GroiNZtCyyGw49Mq2173SpaLsxyKCiEigfbAH8U8HCuE2COEWC2EKDY5VggxSQgRL4SIT0lJKbfY6+vS30PfuHEjHTp0wNPTk6VLlxZ7TWXa5x46dIjY2Fh8fHx47z01GabKYDbC2te00v52d+qtpnL0eAKMGbDrC72V6IIjmnOVeWUhRDXgO+AJKWXmVQ/vAhpKKdsBHwDLi1tDSjlHStlJStkpLCysopqvwRk89AYNGjB//nzGjh1rl/VDQkKu9H1RVCH+mAUZp+HGt+w3G9RRRHSChj0h7iMoMOutxuGYzAX4eBrsWhBYpr8QIYQXmjH/Wkq57OrHpZSZUsrswvurAC8hhMOajNjSQ69o+9zIyEiio6MxlPJPV9H2ueHh4XTu3BkvLycs81bYh7zLsGk6NLsJGt+gtxrbEPswZJ6FxIoV8LkyxkKDbk9KrRQV2svJXCBBSjn9OtfUBi5KKaUQogvaC0VaZYRdePttTAmlt88FMFqsWK2SU94lHzb4tGxB7X//u8RrbNE+tyQq0z5XUcXYMhNMmVqPcXch6iYIrg87PoNWxbfHcFdMFqtdM1ygbKX/PYC7gX1CiN2FX/s30ABASjkLGAVMFkJYgDxgjHRgDESADVpzaVS2fW5pVKZ9rqIKkXVRC7e0HQW1K3bu4pQYPLR2AL+/ASmHISxKb0UOw2gu0N+gSyk3U0qjcSnlh8CHthIFlOpJ/51z6XlczsmnYb3gSu1Z2fa5ZcEW7XMVVYBN08Bigj4v6K3E9nQYB+vfgfjPYbCLFUlVAqPZqn/IxRUQAqw2WKcy7XPLSkXb5yqqEOmnNWPX4W4IbaK3GttTLRxaDdcKjfq7cCpmOTFZ7O+hu/ixuUZRlktlozyVaZ+7Y8cOIiIiWLJkCQ888ACtWxff1rSi7XMvXLhAREQE06dP58033yQiIoLMzKuTjRRuwfp3QRigt20HjTsVne8DU4bW46WKYDRb7dqYC9zEQzf8rSd6KWNFS6Qy7XM7d+583cHQRVSmfW7t2rVLXV/hBqQdgz0LoesDEOxU9Xu2pUE3CG+tFRp1GPfXYAM3xmgpIMDbvibXLTz0oji03rnoCkWl2fie1siqxxN6K7EvQkDne+HCPji/u/Tr3QCTAzx0tzDoBhtOLVIodCPtGOxdDJ3u1RpauTttRoGHjzZQugpgtBTgU9Vi6BXxspWHbnvUz1IHNk0r9M6vaWjqnvhV1/q77FuiZfS4OSYHZLk4lUH39fUlLS2t3Mak6JtQHrptkFKSlpaGr6+v3lKqDpeOw55Fhd55bb3VOI6Yf4ExHRJX6a3E7jgiy8WpDkUjIiJISkqivI27jOYCUrPzkZd98LbzK2BVwdfXl4iICL1lVB02VjHvvIjGfSConhZ2aT1CbzV2xWi22rV1LjiZQffy8qJRo0blft7WY6nc/80fLLy/G+2ahNpBmUJhRy6f0jJbutxftbxz0CpH242Bze9D5jkIqqu3IrthNBfYtXUuOFnIpaIUTdEumgiiULgUW2ZoeefdH9NbiT7E/AukVTsQdlMsBVYsVml3D90tDHpRKlDRRBCFwmXIugB/fgUxY90777wkQptAg1gt7OKmh/FFvdBV2mIZKDpoUB66wuWI+xCsZm1WaFUmZiykHYGzu/RWYheuTCuqSlkuFaXoh1T0Q1MoXILcS7Djcy0fO6Sx3mr0peUw8PCG/e7ZCuAvD12FXErlLw9dhVwULsQfs8CcAz2f1FuJ/vhVh2Y3avNTre7nmBU5m8qgl4GiH5Ly0BUugylLM+gthkCtVnqrcQ7a3g7ZF+HkJr2V2Jyi8z0VcikDf4VclIeucBF2ztcGJvd8Sm8lzkPUTeAd6JYdGI0W5aGXGS8PAx4GoQ5FFa6BxaQNSo7sBREd9VbjPHj5aa0ADq5wu1YApiIPXWW5lA1fT4Py0BWuwd5vIeu8ymwpjrajtD7pR37VW4lNKfLQfVQeetnw8fJQMXSF82O1aoVEtdtCk/56q3E+GvUB/5pul+1iunIoqjz0MuHraVBZLgrnJ/EnLd+6xxNVYqhDufHw1Hq6JK7WDo7dBJW2WE58lYeucHak1HqW1IiEVrfqrcZ5aTsKLEZIXKO3EpvhNIVFQoj6Qoh1QogEIcQBIcQ17eCExkwhxFEhxF4hRAf7yL0+3iqGrnB2Tm2Bszuh+6OaJ6oonoguEFgHDi7XW4nNKLJNzuChW4CnpZQtgW7Aw0KIqxNnBwPNCm+TgE9sqrIM+Hp5qCwXhXOzZaYWH475l95KnBuDQascPboWTNl6q7EJJmdJW5RSnpdS7iq8nwUkAFd3ERoOfCE1tgHVhRB1bK62BHy9DFdSgxQKpyP5EBz5GbpM0tLzFCXT+lYt7HLkZ72V2IQrHrreIZe/I4SIBNoDf1z1UD3gzN8+T+Jao29XfDw9rqQGKRROx9YPwNMPOt+ntxLXoH5XqFYLDrhH2MVoLsDDIPD0cBKDLoSoBnwHPCGlzLz64WKeck0fTCHEJCFEvBAivrxTiUpDeegKpyXznNbru/1dEKAGsJQJg4cWdjnyK+Tn6K2m0pgsVrt751BGgy6E8EIz5l9LKZcVc0kSUP9vn0cA566+SEo5R0rZSUrZKSwsrCJ6r4uvl/LQFU7KH7NAFkDsw3orcS1aDQdLHhz5RW8llcZotv88UShblosA5gIJUsrp17lsBTCuMNulG5AhpTxvQ52l4uNpUGmLCufDmAnx8zTjFFL+8YpVmobdISAMDv6gt5JKYzRb7Z6yCGWbKdoDuBvYJ4TYXfi1fwMNAKSUs4BVwM3AUSAXmGB7qSWjZbmokIvCydj1BZgytVRFRfkweEDLobBnEeTngre/3ooqjMniGA+9VIMupdxM8THyv18jAV3fT6rCIoXTUWCGbZ9Aw55QTzXhqhCtboX4z+HYb5pxd1GMZis+zhBycRV8CguLpJvOJFS4IAe+h8wk6FFFhz/bgoY9wK8GJKzUW0mlMFkKHBJycRuDXvR2Jr9AhV0UToCUsHUm1GwOTQfqrcZ18fCEqMFweLX2jsdFMZmtdm/MBW5k0NWQC4VTcWIDXNgH3R/RKh8VFafFLdowkJOb9VZSYYwOiqG7zV9aUXzKpOLoCmdg6wcQEA5t79BbievTpJ9WlHXoJ72VVBijWYVcykVR0r7KdFHozsWDWh+SrpPAy1dvNa6Ptz807a8ZdKtr/n+bLFbloZcHNSha4TTEfQhe/tBpot5K3IcWQyDrHJz7U28lFcJoLsDXztOKwI0MuoqhK5yCzPPaiLn2d4F/iN5q3Ieom0B4wKEf9VZSIbS0RecoLHIJijx01ULXPkgpyTBlcDrrNCm5KaQZ00gzppGTn0OeJQ9jgRGrtGIQBgQCX09fArwCCPAKoLpPdcL8wgj3D6dWQC1CfUMR7jqtZ/tsrcy/22S9lbgX/iEQ2VMLuwx4VW815cZpCotchb9CLspDryxmq5nES4kcTDtIwqUEEi8lcjLzJFn5144E8/P0w8/TD18PXwzCgERSIAvIL8gnOz+bfGv+Nc8J8AqgQWADIoMjaV6jOS1CWtAipAWhfi7euMqUpRXBtBwKIY31VuN+tBgCq5+FlMMQFqW3mjIjpcRodkxzLrcx6H+FXJSHXl6klBy+fJhNZzex48IO/kz+kzxLHgCB3oG0DGnJzY1upkFgA+oH1r/iZYf4huDl4VXi2vkF+aSb0knJTSE5N5lzOec4nXmaU1mn2J28m9UnVl+5tn5gfWLCYmhfqz3danejflD9ElZ2Qv78Skuv664KiexCi1s0g574k0sZ9KJEDUdUirqNQQ/y0wxLRp7rFh84Eiklu1M0g7r+zHrO52i91JpWb8rwJsPpWLsjbWu2pW5A3UqFR7w9vAn3DyfcP5zWtL7m8QxTBocvH2Z/6n72pOxhy7kt/Hhci5PWq1aP2Lqx9InoQ9c6XfH1dOKMkQILbPsYGsRCRCe91bgnwfWgTjttgHTPJ/VWU2YyC21SkK/9za3bGPTwQB8AkrNMOitxbs5mn+X7I9+z8vhKzmafxdfDl251u/FA9AP0juhNmL9t2xqXRrBPMJ1rd6Zz7c6A9kJzMvMk285vI+5cHKtPrGbp4aX4evgSWzeWGyNvpG/9vgR4BThUZ6kk/ADpp+Gm/+itxL1pfgus/w9kJ0O1cL3VlIkimxQWaH+HxG0MeoCPJ/7eHqQog34NVmkl7lwcCw8tZGPSRoQQdK3dlYdjHqZfg35OZRyFEDQKbkSj4Ebc2eJOzAVmdlzcwfoz6/n99O+sO7MOHw8fekf0ZmjjofSM6ImXoeSwj92RUpsXGtIEmt+srxZ3p/lgWP82HP4ZOtytt5oykXLFoPvYfS+3MeigeenJWUa9ZTgNFquFNSfXMHffXI6mHyXUN5T7o+/n9qjbqR1QW295ZcLLw4vudbvTvW53pnSZwu7k3aw5uYafT/7Mr6d+JdQ3lCGNhzAyaiSNgnXqN35qC5zfDUPeV2X+9qZ2Wwiur4VdXMSgF9mkcGXQy0dYoI/y0IECawE/Hv+R2Xtmk5SdRNPqTXm759sMihxU6iGmM2MQBjrU6kCHWh14tvOzbE7azPKjy/k64WsWHFxA1zpdGd18NH3r98XT4MA/7S0zwb8mtLvTcXtWVYTQvPRdX4I5zyUGbisPvYKEB/qScP7qcadVByklG5M28r9d/+No+lFah7bm2c7P0qd+HwzCvTxHL4MXfRv0pW+DvqTmpfL9ke9ZcngJT61/ijoBdRjbYiy3Rd1GkHeQfYUkH9Im0/f5t0sYF7eg+WDYPgeOr9fuOznJWSaCfD1V6X95qcoe+rH0Y9z/y/088vsjWKwWpveZzsJbFtKvQT+3M+ZXU9OvJvdH38/q21Yzo+8MIgIjmLZzGgOXDOT/dvwfF3Mu2m/zuA+1xlGd77PfHop/0rAn+ARB4iq9lZSJlCwT4UGOydByKw89LNCHLJOFvPwC/Lzt/2roDOSac5m1ZxZfHvwSfy9/Xuz6IiOjRup/UKgDHgYP+jXoR78G/TiYdpAvDn7B1wlfs/DQQoY1GcbEthOpH2jD3PasC7B3MXQYBwEuXhTlSnh6Q9MBkLhGa9bl5OcWyVkmwqrZP9wCbuah/5W6WDUORuPOxXHrD7cy78A8hjUdxsoRKxnTYkyVNOZX0yq0Fe/0eoeVI1ZyW7Pb+PHYjwz9fiivbH2Fs9lnbbPJH7O1oQvdHrLNeoqy0/xmyEmGc7v0VlIqyVlGwoOUQS83RYcO7h52yc7P5rW415j06yR8PHz4cvCXvNb9NWr41tBbmtMRERjBS91eYvXI1YxpMYaVx1YyZNkQ3oh7g+Tc5IovbMqCHXOh1TAIbWI7wYqy0WyA1qzLycMuUkpSlIdeMcILE/fdubhoT8oeRv04imVHljG+9XiWDF1CTHiM3rKcnnD/cKZ0mcKq21YxMmoky44u45Zlt/C/nf8jw5RR/gV3LgBTBvR43PZiFaXjVwMadtfCLhbsD7AAACAASURBVE5MlsmC0WxVHnpFKPqhuaOHbpVWPtv3GeNXj0dKyYJBC3i609POXQ7vhNQKqMVL3V5ixa0r6N+wP5/v/5xbvr+FrxO+xmwtY9uIArNW5h/ZC+p1tK9gxfWJGgTJB7QKXSelyBaFO6BKFMpg0IUQnwshkoUQ+6/zeB8hRIYQYnfh7WXbyywbIf7eeBiE28XQ043pTF47mRm7ZtCvQT+WDFNeeWWpH1ifd3q9w7dDv6VlSEve2f4OI34YwbrT65BSlvzk/d9B5lnlnetN1CDtoxN76cmZjstBh7J56POBQaVcs0lKGVN4e73ysiqGwSCoWc3brTz0xEuJjPlpDDsu7ODl2Jd574b37J9bXYVoEdKCOQPn8FH/j/AQHjy27jEeXPsgx9OPF/8EKWHLDAhvpWVaKPSjZlMIbQqHndegp2QXeehOYtCllBuBSw7QYhPCA33dJoa++sRq7lp1F2armQWDFnB71O3uOxhCR4QQ9I7ozdJhS5nSZQr7UvYxcsVI/m/H/5FjzvnnxUd+heSDmneufhf6EzUITm7SDqmdkOTMorJ/Jwm5lJFYIcQeIcRqIcS1PVILEUJMEkLECyHiU1JSbLT1PwkL9LnyNsdVkVIya88sntv4HK1CW7F4yGLahrXVW5bb42Xw4l8t/8XK21YyvOlwvjz4JcO+H8aak2v+CsNs+R8E1YM2I/UVq9BoPhgK8uHYOr2VFEtKlglvTwNBfo4p+bGFQd8FNJRStgM+AJZf70Ip5RwpZScpZaewMPu0aQ0P9LnyNscVMReYmbplKh/t/oihjYfy2Y2fUdOvpt6yqhQhviG82v1Vvrr5K0L9Qnl2w7M8uPZBziT+qDXi6v4ouHBPHLeifjfwre60YZeilEVHvbOutEGXUmZKKbML768CvIQQulmgsEAf0rJNFFhLOdhyQnLMOUz+bTI/HPuBye0m81bPt1y6mZarEx0WzcJbFvJClxfYk7KHEdte5LOatTDHjNVbmqIID09oNlBrp2t1vmllyVkmhx2Igg0MuhCitih8+RFCdClcM62y61aU8EAfrBLSclzLS083pnPfz/cRfyGeN3q8wUMxD6l4uRPgYfBgbMux/NBzGr1ycpgR6MPoXyZwIPWA3tIURUQNgtxUOLtTbyXXkJJlctiBKJQtbXEhEAc0F0IkCSEmCiEeFEI8WHjJKGC/EGIPMBMYI0vN+7IfRVNBXCmOnpybzISfJ3D48mHe7/M+tza9VW9JiquotfNL3r+cy/9i3yDdmM7YVWOZvnM6Rot7pci6JE2LqkZXl36tg3Fk2T+UoTmXlLLEJs9Syg+BD22mqJJcKf93kTj6+ezz3PvzvaQZ0/h4wMd0rdNVb0mKq0k/A/uWQOf76R91K50j+zEtfhrz9s9j3el1vNHjDVUXoCd+1bWq0cM/w4BX9FZzhXyLlcu5ZsKqOa74z60qReGvfM8UF/DQL+Rc4N6f7yXDlMFnN36mjLmzElfor8Q+DECQdxCvdX+NOQPnkF+Qz7jV45gWPw1TgfP/zbktTlg1mlqUg+5AD93tDHqYi3RcLDLm6aZ0Zg+cTXRYtN6SFMWRnaL1bYkeDdX/2Xo3tm4sy4YvY1TUKOYfmM/tP97O/tRiC6oV9qaoavTwz/rq+BtXhkM7qDEXuKFB9/XyIMjX06mrRVNyU5j480QuGS8xe+BslWPuzGz7CCxG6PlksQ8HeAXwcuzLzB44m1xzLnetuotPdn9S9r4wCtvghFWjV/q4KA+9coQHOW+1aIYpgwfWPkBKXgqzBsxSnrkzk3cZtn8GrW+Fms1KvLR73e4sG76MwY0G8/Gejxm3ahwnMk44SKgC0Lz0ExvBlK23EuDvw6FVDL1ShFVzzlF0ueZcHvntEU5mnGRG3xnqIM3Z2f4p5GdBr6fLdHmQdxD/6fUfpt0wjTPZZxi9cjTfJn5berMvhW2IGqRVjR53jqrRlCwTQkBoNW+H7emWBj08yMfpPHSz1cxTG55ib+pe3u39LrF1Y/WWpCgJU7bWIjdqENQuX0jsxsgbWTZsGTFhMbyx7Q0eW/cYl4wu0w7JdWnQDXyDnSbskpxlIsTfGy8Px5lZ9zTohcOincUzklLyetzrbDm7hZe7vczAhgP1lqQojZ3ztJBLr2cq9PRw/3BmDZzFc52fY8vZLYxaMYq4c3E2Fqn4Bx5eWk764V+0WaM6k+LgKlFwU4MeFuhDnrmAbJNFbykAzNk7h+VHlzO53WRGRqmmTk5Pfi5smQmNboD6nSu8jEEYuLvV3Sy8ZSGB3oFM+nUS0+OnYy5QB6Z2I2pw4azRP/VW4vCyf3BTg+5Mo+h+PPYjH+7+kKGNhzK53WS95SjKws75mlHoM8UmyzUPac6iIYsYFTWKeQfmMW71OM5knbHJ2oqraNpfqxo9rH/VaEqm0aEHouCmBt1ZhkXvuriLl7e+TOfanXmt+2uqN4srYM7TWuRG9tKqD22En6cfr8S+wrQbpnEq8xR3/HgHa046R6zXrfAP0WLpOk8xklKSkq08dJsQfqW4SD+Dfj77PE+uf5J61erxfp/3VddEV2HXF5B9EW543i7L3xh5I0uGLaFx9cY8u+FZXo97XfWDsTVRg+DiPl2rRtNzzZgLpEMbc4GbGvR6NfwQAk6k5JR+sR3Is+Tx+LrHyS/IZ2a/mQT7BOuiQ1FOzEbY/D407AGNetltm3rV6jF/0HwmtJnAksNLuGvVXZzMOGm3/aoczW/WPupYNXo8VbM99UP8HbqvWxp0f29PGoUGcPB8hsP3llLyypZXOHTpEO/2fpfGwY0drkFRQf78ErLOww3P2X0rL4MXT3V8io/6f8TF3IuMXjma1Sf0j/u6BUVVozp2Xzx4PhOAVnUdO//XLQ06QMu6QSScd/ycwS8OfsHqk6t5rMNj9I7o7fD9FRXEbIRN07UJOI1ucNi2vSN6s2ToEpqHNOe5jc/x5rY3VZMvW9B8sK6zRhPOZxLs50XdYHUoahNa1Qni9KVcsoyOSxHbdXEX7+98n4ENBzKxzUSH7auwATvnQdY56Peiw4c/1w6ozdyb5jKh9QQWJy7m7lV3qyyYyhJVNGv0d122P3guk5Z1Ah2eCOHWBh3g0AXHvEKn5qXyzIZniAiM4PXur6uMFlciPwc2TdMyWxrp867Ky+DFU52eYmbfmSRlJzF65WjWn1mvixa3oH5XbdaoDmGXAqsk8UIWreo4/uzMbQ16y0KDfvBcpt33KrAWMGXjFDLzM5l2wzSqeVez+54KG7L9U8hJgX4v6a2Evg368u2Qb4moFsGjvz/K/3b+D4vVOQrkXAoPT4i6SZdZoyfTcsgzF9CyTqBD9wU3Nui1gnwICfB2iEGfvXc2f1z4gxe7vkjzkOZ2309hQ4yZsGWGVjLeoJveagCICIzgy5u/ZFTUKObun8ukXyeRmpeqtyzXI2oQ5F2CM9sdum2RzXH0gSi4sUEXQtCyTiAJF+xr0Hde3MnsvbMZ2ngoI5qNsOteCjvwxyztn77vi3or+Qc+Hj68EvsKb/Z4k70pexn942h2J+/WW5Zr0bQ/GDwhcZVDt004n4mnQdA03PHv1N3WoIMWRz90IQtLgX0a9WSYMpiyaQr1qtXjxW7OZRAUZSD3Emz9AJrfAvU66K2mWIY3Hc7XN3+Nj6cPE9ZM4OuEr52m6ZzT4xsMkT0dHkc/eD6TpuHV8PH0cOi+4OYGvWWdIPItVk6k2r7ASErJa3GvkZqbyn97/5cArwCb76GwM5umQX429J+qt5ISKeoF07NeT97Z/g5TNk0h15yrtyzXoMUQSDsCKYcdtmXC+cwrSRmOxq0NelEMqyjJ35Z8f/R7fj31K4+0f4Q2NdvYfH2FnUk/ox2GtrsTwlvqraZUgryDmNFvBo+2f5TVJ1Zz1+q7OJ3pPAORnZbmg7WPiT85ZLu0bBMXM026xM+hDAZdCPG5ECJZCFHs9FuhMVMIcVQIsVcI4TTvXZuEVcPbw2Bzg342+yzvbn+XzrU7M6HNBJuurXAQG97RPvZ5QV8d5cAgDEyKnsQnAz4hOTeZMSvHsDFpo96ynJvgCKgTA4ccE0cvKmZs6cQe+nxgUAmPDwaaFd4mAZ9UXpZt8PIw0KxWNZtmulillZc2v4QQgjd7vIlBuPWbHPck+RDs/ga63A/V6+utptz0qNeDxUMWExEYwcO/Pcwnuz/BKvUf6OC0tLgFknZA1kW7b1XUbkQvgy7KcsAihIgEVkopr4ktCCFmA+ullAsLP08E+kgpz5e0ZqdOnWR8fHy5BV94+21MCYfKfP2xlGzSc810bFij3HsVx8Xci5zOOk2joEbU9KtpkzUVDiY5AYwZENFJy4JwUazSyqnMU6QaU6nuU51GwY3xFI4/iHN6zDlw9k+tx0u12nbd6mhyNplGMx0alGxvfFq2oPa//12hPYQQO6WUnYp7zBbuZT3g73XKSYVfK07IJCFEvBAiPiUlxQZbl46/tyfmAiv5Nsh0ySswkpSVRHWf6oQqY+6aGDMgNw2C67m0MQctBBMZ3IiGgQ3JMGWQkHaQXEue3rKcD68A8PTVfu92Jje/AH9v/f6ubLFzcTXuxbr9Uso5wBzQPPSKbFbeV7Xzx9N4fs425k/oTJ/m4RXZEtCqQe9Zcw+nMkP4fvj3yjt3RaxW+LQv5PjBo6vAy09vRTYhEvBM/pOn1j9FjvkUr/d4nUGRJUVJqyBr/g07PoPn1oOPfSo4jeYCnnjlZybf0ISbbtKnwNAWHnoS8PdAZARwzgbr2oSiWNaeM5VrpbsocRF7UvbwfJfnlTF3VfYtgfO7of/LbmPMi2gf3p7FQxYTVSOKZzc8y/Sd01XLgL/T4mYoMMHR3+y2xcHzmRRYpW4ZLmAbg74CGFeY7dINyCgtfu5Igv28aBcRzO+JyRVe42z2WWbsmkGver24pdEtNlSncBj5ufDba1C3PbS9XW81diHcP5x5N81jdPPRzNs/j8lrJ5NuTNdblnNQvxv4hcAh+6Uv/p6QjEFAbONQu+1RGmVJW1wIxAHNhRBJQoiJQogHhRAPFl6yCjgOHAU+BR6ym9oKMqBlLfacSSc5q/yjvqSUvLb1NQSCqd2mqi6Krsq2jyDzLNz4FhjcNzPJy8OLl7q9xGvdX2PnxZ2M+WkMhy6VPYnAbfHw1HLSD/8Mlny7bLE24SKdIkOoEeBtl/XLQql/2VLKO6WUdaSUXlLKCCnlXCnlLCnlrMLHpZTyYSllEyllWyll+VNX7Ez/lrUA7RW0vKw4toK483E80fEJ6lSrY2tpCkeQeQ42va9VDUb20FuNQ7it2W0sGLQAs9XM3avu5qfjjimscWpaDgNTBpzYYPOlky7ncuhCFgNaVvyczha4r6vyN1rWCaRedT/WltOgpxvTeS/+PWLCYhjdfLSd1Cnszq+vgNUCN76ptxKH0jasLYuHLKZVaCumbJrCf3f8t2rH1Zv0Be9AOPiDzZf+rdC2DCh0HvWiShh0IQT9W4az+WgKRnPZeyO/v+t9svOzmRo7VRUQuSqn4mDft9DjMQhppLcah1PTryaf3fQZY1uM5cuDXzLp10mk5dk/fc8p8fTReqQf+gkKbPvCtjbhIo1rBtA4TN9ZCFXGSg1oWQuj2cqWo2XrK/1n8p8sO7KMu1vdTVSNKDurU9gFawGsfhaCIqDnU3qr0Q0vgxcvdH2Bt3q+xd6UvYz5aQwHUg/oLUsfWg3T2iWf3mqzJbOMZrYdT2NAK329c6hCBr1r4xACvD3KFHYxW828se0NagfU5sF2D5Z6vcJJ2TkPLuyDm94Eb3+91ejOsCbDWDB4AQLBuNXjWH50ud6SHE/TAeDpBwdX2GzJTUdSMRdI+rfQN34OVcig+3h60DsqjN8SLmK1llzT9E3CNxy5fIQpXabg76UMgUuSkwq/v6nNCW11q95qnIbWoa1ZNGQR7cPbM3XLVN7c9ibmAscNUtcd7wBoNgASftQKzWzA2oSLBPt52ay9SGWoMgYdtLBLcpaJ/eeuX2SUkpvCx7s/pndEb/rV7+dAdQqb8stLYMqGm98DlWr6D0J8Q5g1cBbjW49nceJiJv4ykZRcx7TicApaDofsC1rDrkpSYJWsO5RMvxbheHrob071V+BA+rYIxyBg1b4L173mf7v+h9lq5vnOz6ucc1fl+AbYsxB6PA7hLfRW45R4Gjx5utPT/F/v/+PQpUOMXlmFRtxF3QQe3pBQ+bDLHyfSuJxrpr/O6YpFVCmDHhLgTf+WtVi84zR5+ddmu+xN2cuKYysY12ocDYIa6KBQUWnMRlj5JNRoBL2f0VuN0zOo0SC+uvkrfD19mfDzBBYdWuT+I+58g6BxXy19sZJhl/lbTlLD34v+LfQ/EIUqZtAB7uvZiMu5Zpb9mfSPr1ullf/88R/C/cKZFD1JJ3WKSrN5Olw6BkOmu12/FnsRVSOKRUMW0b1ud9764y1e2vISRkv5q6pdijYjIeNMpcIuJ1Nz+DXhInd1a4ift3O0La5yBr1LoxDa1gtm7uYT/zgcXXFsBfvT9vNExyfUQairkpwAm9/XerU0Uecf5SHIO4gP+n3A5HaTtXepq8eRlJVU+hNdleaDtZa6+7+r8BLztpzAy2Dg7tiGNhRWOaqcQRdCcF+vRhxPyWH9YS2FMdecy4xdM2gX1o4hjYforFBRIQossPwhrTXqTf/RW41LYhAGHop5iI/6f0RSdhJjfhrDlrNb9JZlH3yDoNmNcOB7rV6hnGTkmvk2PolhMXUJD/S1g8CKUeUMOsDNbetQO8iXzzadAGDegXmk5qXyXOfn1EGoqxL3AZzbpWW1VAvTW41L0zuiN4tvWUwt/1pMXjuZWXtmueeIuzYjIScZTm4u91O/2X6aPHMBE3s6V/VxlTToXh4GxveIZOuxNDYdP8r8/fMZFDmI6LBovaUpKkJKIqx7W2u+1HqE3mrcgvpB9fnq5q+4pfEtfLT7Ix79/VEyTJWbKeB0RN0E3tXKHXbJt1iZv/UEPZvW1G126PWokgYd4M7ODfD39uDVTdMokAU83uFxvSUpKoK1AH54WPvHvGWayjm3IX6efrzd821e7PoiW89tZfTK0RxMO6i3LNvh5acNkE5YUa6Wukt2nuFipomJvZzLO4cqbNCD/b24o7uBi3ILvWqNICIwQm9JioqwabqWqXDz/0E158gFdieEEIxpMYb5g+ZjsVq4e9XdLDuyTG9ZtqPNSMi7DMfXleny1GwT/12TSNdGIfSJcr7QXpU16ABnxLcYpB8793QoVxdGhZOQFA/r/6NltbQdpbcat6ZdWDu+HfotHWt15JWtrzB1y1Ty3GEgdeO+4Fu9zGGX/6w6RG6+hbdGtHHK87Yqa9DjzsWx/cIfjGoygaQ0+GjdUb0lKcqDKQu+uw+C6mkHoQq7E+IbwicDPuGB6Af44egP/GvVvziZcVJvWZXD0xtaDYeElVqriBLYdjyN73YlcX+vxjQNt8+g6cpSJQ26lJIZu2ZQJ6AOU3rcy4j29Zi14RjHUkr+hSqciNVTIP0U3DYb/KrrrabK4GHw4JH2j/DxgI9JyU1hzE9j+Pnkz3rLqhwxY8GcU2IrgHyLlZeW7yeihh+P9mvmQHHlo0oa9LWn13Ig7QAPxTyEt4c3/765JX5eHjz97R5MFhV6cXr2LoHdX0Gvp6Fhd73VVEl61uvJkqFLaFK9Cc9seIa3tr1FfoF9ZnXanfpdIaQx7P7mupdM+zWRo8nZvD68tdNUhRZHlTPoFquFmbtm0iS4CUMbDwUgLNCHd0dGs/tMOi8vP+D+vSxcmZRE+PFxaBALNzyvt5oqTe2A2sy/aT7jWo1jUeIi7lp1F2cyz+gtq/wIoXnpJzfB5VPXPPzD7rPM3nCcsV0b0M9JerZcjypn0FccW8HJzJM82uFRPAx/vdIObluHR/s1ZXH8Gb6Iu/aXqnACTNmw+G5tWMWoeeDhpbeiKo+XhxfPdn6WGX1nkJSdxB0r72DNyTV6yyo/0WMAAXsW/ePL+89m8Px3e+kcWYNXh7bWR1s5KJNBF0IMEkIkCiGOCiGmFPP4eCFEihBid+HtPttLrTymAhMf7/6Y6JrRxfY6f3JAFANahvP6yoNsPVa2UXUKByGl1kUx9TCM/AyC6uitSPE3+jXox5KhS2hcvTHPbniW1+Ned60GX9XrQ6PesPvrKx0YU7NNPPDlTmr4e/Pxvzri7en8/m+pCoUQHsBHwGCgFXCnEKJVMZcullLGFN4+s7FOm/Dd4e+4mHuRR9o/UmzKkcEgeH90DI1qBjDpi53En7ykg0pFsfwxWxv23PdFaNxHbzWKYqhXrR7zB81nQpsJLDm8hLGrxnIs/ZjesspOzL+0g/bTcaRmm7jrsz9IzTYx5+5OhAX66K2uTJTlJacLcFRKeVxKmQ8sAobbV5btMVqMfLbvMzrW6ki3Ot2ue12grxdfTuxCeKAP4z7frjx1Z+DY7/DzC9D8Fu0gVOG0eBm8eKrjU3wy4BPS8tIYs3IMSw4vcY1zqZZDwDuQvB1fMGbONk6m5TD3ns60jQjWW1mZKYtBrwf8/aQjqfBrVzNSCLFXCLFUCFG/uIWEEJOEEPFCiPiUFMeOvFp6eCkpeSk8HPNwqQUBdYL9WPRANyJq+DFh3g7WJ5Y+WFphJ1KPwpLxENZCS1E0OP/bXoWWBfPdsO9oH96e1+Ne5+kNTzt/LxjvAHKaDYUDy8lKT2PBhC70bFZTb1Xloiz/HcVZv6tfbn8EIqWU0cBaYEFxC0kp50gpO0kpO4WFOa5sNs+Sx2f7PqNr7a50rt25TM8JD/Rl0aRYmoRVY+KCeD7deNw1vAx3Ii8dFo4BgyfcuVBrjatwGWr61WTWwFk81fEp1p1Zx20rbmP7+e16y7ouW4+l8uChaPwwsiT2FF0bh+otqdyUxaAnAX/3uCOAc3+/QEqZJqU0FX76KdDRNvJsw7eJ35JmTOOhmIfK9byQAG8WP9CNgS1r8daqBB7+ZhfZJoudVCr+gdkIi/4Fl0/CHV9CjUi9FSkqgEEYmNBmAl/d/BX+nv7c98t9TI+f7lQ561JKZm84xl2f/cE5/5YYw2NocOwb7SDexSiLQd8BNBNCNBJCeANjgH+UVAkh/p5yMAxIsJ3EypFrzuXz/Z8TWyeWDrU6lPv5gb5efHJXB14Y3II1+y8w7IPNbD+hDkvtitUK3z8ApzbDiFkQ2UNvRYpK0jq0NYuHLGZk1EjmHZjH2J/GcuTyEb1lcT4jj4kL4vnP6kMMalObHx7piW/3B7VsqhMb9JZXbko16FJKC/AI8DOaof5WSnlACPG6EGJY4WWPCSEOCCH2AI8B4+0luLwsPbyUS8ZLTI6ZXOE1hBA8cEMTvr6vG/kFVu6YHcdLy/eRZTTbUKkC0LyiNVPg4HK48S3VdMuN8Pfy55XYV/iw34ek5KUwZuUYFhxYoMvwDKtV8mXcSQZO38jWY6m8PKQVH43tQDUfT62nvn8obP/U4boqi9ArLtypUycZHx9v1z1MBSYGfzeYRsGNmHvTXJusmZtvYdovh5m35QQ1q/nw5MAobu8YgaeHOqyzCRv+C+vegthH4Ka39FajsBOXjJd4deurrDuzjo61OvJGjzeoH1hsLoXNiTuWxjtrDrHnTDq9mtXkrVvb0iD0qjnCv74CW2fCE/sg2LlaawshdkopOxX3mFtboeVHlpOSl8Kk6Ek2W9Pf25OpQ1qx7KEe1A/x54Vl+7jxfxtZve/8P4ZOKyrApumaMW83Fga+obcahR0J8Q1hRt8ZvNnjTRIvJTJyxUi+TfzWrokH+89mMH7edu78dBvJmUam39GOL+7tcq0xB+h0r/ZuMX6e3fTYA7f10M1WM0OWDSHMP4wvB39pl97FUkp+PXiRd9cc4lhKDo3DApjUqzG3tq+Hr5fzNvBxSrZ+AL+8BG3v0OLmBvXzqypcyLnAy1teJu58HF1rd+XV7q/abOCMlJItR9OYvfEYm46kEuznxcN9mzAuNrL0/9FvxsDZeHhiP3g5zyDokjx0tzXoy48uZ+qWqXzU/yN6R/S22z4AlgIrq/ZfYM7GY+w/m0nNaj7c3imC0Z3qE1kzwK57uwVbZsKvU6H1bXDbp+DhqbcihYORUvLdke94L/49rNLK4x0e584Wd2IQFQsiZOSaWb77LAu3n+bQhSzCAn2Y0COSf3VtSLBfGXsAndgIC4Zqow07O083kypn0AusBQz/YTj+nv4sHrLYYZNFpJRsPZbGvC0nWZeYTIFV0q1xCLfG1GNQm9pU9/d2iA6XQUr47XXYPF07iLrtU9Vwq4pzIecCr8a9ypazW4gJi+HV7q/SpHqTMj3XZClg0+FUftx7jjX7L2CyWGlbL5i7ujXg1vb18PEs57s+KWHuQMi6CI/tcpq/zSpn0NecXMOzG55lep/pDGw40C57lMbFTCNLdyaxJP4MJ9Ny8TQIejaryY2tatOvRTi1g53nLZwuWAvgp6dh5zzoOB5uma7CLApAc4x+PP4j/93xX3LMOdzX9j7ua3sfPh7X9lPJMVnYdCSVtQkX+eXABTKNFoL9vBjarg5jOjegTb1Klu0nroGFo+HWT7QWu05AlTLoUkpGrxxNriWXH4b/8I8WuXogpWT/2UxW7j3HT/vOk3RZm8PYum4QvaPC6NGkJp0ia1StmHt+DiybBIdWQs+noP/LWk9qheJvpOWl8d8d/2XViVVEBkXyYrcX6VKrKwfPZ7LlaCqbj6byx/FL5BdYCfT1ZGCrWgyNrkuPpjVt1xlRSpjVCyxGePgPp3A6qpRBjzsXx6RfJ/FK7CuMinKuHGYpJUeSs/ktIZnfD13kz9PpWKwSb08DMfWr0zmyBp0ahtC+QXX3Dc9knNXK+S/u1/LMY8tXvauoWhjNBSza/yufHZxGhuUCIrsDWecHIS1BNA2vRt/mYfRrUYtOkTXwslfq8P5l+NH11wAAEohJREFUsHQC3D5fCw3qTJUy6Pf/cj9H04+yZuSaYt+iORPZJgs7Tlxiy9FUdpy6zIGzGVgKUx8bhvoTHVGdtvWCaFUnmJZ1Agmt5tzfT6kk7YRFYzUPfdRciLpJb0UKJyIvv4BDFzI5eD6T/Wcz2ZuUTuKFLO1/QpgJr78FY8BavA1eTGj9IJNi7sbL4IC4trUAPuoKnr7w4Cbd301WGYN+IO0AY1aO4YkOTzCx7USbru0IcvMt7D6Tzp4zGew5k87epHTOZfw1JCAs0IeoWtVoFh5Is1rVaFyzGo3DAggP9HHYwW+FkBJ2fAY//xuq1Yaxi6FWcS31FVWBTKOZ4yk5HE/J5lhKNocvZnPkYhanLuVeaZ8S6OtJu4jqREcEEx1RnY4NaxAW6MPJjJO8u+NdNp/dTJPgJjzX5Tm613XAXNndC2H5g3DbZxB9u/33K4EqY9Cf3fAsm89u5pdRvxDo7R6d+S7l5JNwPpOD5zJJvJjFkYtZHEnOJjf/r2HW/t4eNAwNoGGIPw1D/Ymo4UdEDX/q1fCjbnU/rZxZL4yZsOJRrZS/2U1ajrl/iH56FHbHZCngQoaRs+l5JF3Wbmcu5XIqLYdTabmk5fzVmMvDIGhUM+CKo9KqbhCt6gQRUcPvuk6KlJINSRt4d/u7JGUncUPEDTzd6WkaBTey3zdltcKcGyA3DR6J18Yg6kSVMOhnss4w5Psh3NP6Hp7q+JTN1nVGrFbJ+Uwjx1OyOZ6Sw8nCf5RTaTmcuZxHvuWfvTECfT2pG+xHrWBfagX6UDvYl/BAH8IKbzWr+RBazYcAbw/bevonN8PyhyAjSTv47P6Y6mfuwpgsBaRl55OabSI120RypomULBMXs4xcyDBxMdPI+QwjqdmmfzzPILQZAw1D/QtvATSuGUDjsGo0CPGv8AFmfkE+Xyd8zey9szFZTNze/HYeiH6AUD87tb09tRXmDYY+L0CfayZxOowqYdDf2f4Oiw8tZs3INdQKcO7J3PbEapWkZps4czmXpMt5nM8wci5d+3gxU7ulZJkorkuBr5eBEH9vQqp5U8O/6OZFdX9vqvt7Eez31y3Iz4sgXy8CfT3xv/qFID8XfnsN/pgFIY3h1lnQoKvjfgiK6yKlxGSxkmW0kGk0k5lnJtNoISPPrN1y87mcayY910x6bj6XcvO5lJPPpex8sq7TOrqGvxe1gnypFeRL7SBf6lb3o051X+pV96N+DX9qB/vadR5nWl4aH+/+mO+OfIePhw/j24znnlb34O9lBy/623vg8M/waLxuPV7c3qBn5mcyYMkABjQYwNu93rbJmu5MgVWSllPoYWWbrnhdadkmLuWYuZRj4lJOPul5Zi7n5JNpLLkHvEFANR9PAn296GfYxcPGT6ltvcDvwSP4pc4DePpWw9/bEz8vD/y9PfDz9sDXywM/L+2jr5dB++jpgY+XAR9PA96eBnw8PfD2NODtYcDLQzj3OYENkFJisUrMBVZMZiv5BVbyLVZMlgKMZismixWTuYA8s/a58cr9AnLztVtevuXK/Zx8CzkmC9mmgsKPFrKMZswFJf/P+3t7UN1PeyEPCfjrFhrgTc1AH0ICvP/xDq/cBTt24kTGCWbumsna02sJ8Q1hYpuJ3NH8Dnw9bVjzcfkUfNgZWg3ThpXrgNsb9Hn75zF953S+HfItLUNb2mRNxV9YCjSPLr3Qi9O8OjOZeZqByDZZ+P/27j46qvJO4Pj3NzNJAEkCEl4UCARBBFJUpC607qpBWcBuqdVugKXWU8+66nH3jz0ed/dQ93TbP/blbNdTi+1CsWurp9W29kBOi6uHtFW3COUtCNEqQV4aIbwTSCBhZu5v/3huyBAImSTzkrnz+5wz3HvnXmZ+T2bmd5/73HufJ9J8gHsOPMuMlo18Einne0OfYIvOoKU9xjk/ybTH+tdNakdiL4iEKAiHKAgJkXCISEiIhIVIKEQkLIRDQiQkhMTNhxPmQ+K6Qw4JhMQ93zEml+DWdTdEl6q6oboUPFXUn3oXp0rcS5h6EPM84ur+hrG4EvM8Yp4Si7vEHfOUaMxP3nGvX2MqFITF32lGuKYozNCiiD8foXiQe654UAFDiyKUDIpccpQ1bIg76iodXDBgEnRf7Ty2k5U7VrLp8CZGDR7FVz/1VR6Y8kDqEnvtN+Gd/4Tlv4DJ81Lzmr0Q6IQe9aIs+sUiyovLU9ZFrumFlmPuy73lBQgXwl3/AH/yOEQuv44+7unFGuX5C27aFvVoi7llVxt1Nc/E2umFmMeFuHIh5hGNJz6UWNwj6inxuKvdxjyPuOcSasxTPE+Jq5t2JN6455Kx4uaBi4k68degqpccFYj/j9C5M5COHUOo87lwSAgn7kxC4u98Onc6kZC/cwr7O6ewdB6NRDqPUgrD7uilKBKiyJ9ePKqJuCOeQYXuaCdt12HnqC1NW1i5YyXbj26nbHAZD894mC/d+KX+N8VcOAffvxvOnYTHN8LQzA2nCQFP6K/ve52n336alVUruXP8nSmIzCSl9bhrI3/3u+4uulv/yp0sKrk+25EZc5GqsvXIVla9t4rNhzdTUlhC9dRqlk1bRtngfgwA3bQbvl8Fk+5yl+FmsDkwsAldVVn6q6W0RltZ94V1fe6ZzfTCqf2wcSXseBli52H6F6Dqa1A2JduRGXNVO4/t5MXdL1J7sJZIKMJ9k+5j2U3L+t5Mu3kVvP40LPh3mPNYaoO9iqsl9Jzup7TuWB31J+p5Zs4zlszTyYvDnjddZ/973oRQBGZWw2f/DkZOzXZ0xiTl5pE38+zdz3LgzAFeev8lavbWsLZhLbeMvIXqm6q5d8K9vbu7/PZHoaHWdf18/S1QPid9wScpp2voT731FBsPbWTDgxvSc4lSPlOFpvdg189g12tw9pC7y3PWl91oLta0YnLcmQtnWNewjlf+8AoHzx6kpLCEz036HPdPuZ+pw6cmd1VV63F4YT60HIWH1sG429IedyCbXJpam1jw2gKWT1vOU59+KoWR5bF4DBp/Dx/9L/xhPZzY42rjN8yDW5fD1IUDpk9oY1LFU48tTVt4bc9rbDiwgagX5YbSG7hv0n0srFjY8+hJzZ+4G47aTsNXfgnXzUxrvIFM6M9tf441u9aw/ovrUzZcVd7xPDj+Iex7B/a/7e7sPH8KQgUw8bOufXz6YrtV3+SN022neWP/G6zft57tR7cDMO3aacwrn0dVeRWTh02+cs391AH4n0UQPQfVL7vfT5oELqG3x9uZ//P5zBw5k+9UfSfFkQWUF4eT+1y3tU273FiJn+yA9ma3vrQcKv4UpsyHG6pgUEl24zUmyw61HOLN/W9Se7CWumN1AIweMpo7xt7B3OvnMnv07Eu7GTixF15+wF048JknoeoZiKS+h9R+J3QRWQB8GwgDa1T137qsLwJ+BNwGnACqVXX/1V6zPwm9Zm8NK/5vBavvXc3c6+f26TUCKXoezhyC0wdd/ymn9sOJBvdFO9HgrkoBkLDr7XDsbBg3GybeAcMnZjNyYwa0o+eO8lbjW2z8ZCObDm+iJdoCQEVpBbNGzaKyrJIZI2YwefAYCmr/xY3ENWo63L3CdROdwqbKfiV0EQkDHwH3Ao3AFmCpqr6fsM0TwExVfUxElgD3q2r11V63rwldVVnyqyW0xdpYu3ht7t8OrgrxKMTb3TTW7hJvtM0l6Gir6z/8Qiu0NUP7GTc9f8r1/HbupDsh03LErUskYRg+AUZMdo/RM9xj5E1QMDg75TUmx0W9KPXH69l2ZBvbjmyj7mgdZ6NnASgIFTChZAIV4Wuo+GMd41pPMiZSzJipi7h2wl0MvW4m4dIJ/eqkrr8JfS7wdVX9c3/5nwBU9V8TtnnD3+ZdEYkATcBIvcqL9zWh121fw5d3fZuvtRVQHevnXi/p5qaE7S75P5qw7M+rlzDvL2vcTb24m/fi4MVcAtc4vRaKwODhMGRE56N4DAwd7a4+KR0Pw8ZDyVg7iWlMmnnq0Xi2kfoT9Xxw8gP2nd7Hx80f09jSiKeXdnchqhSr8tC1s/ibxS/16f36ex36WOCPCcuNQNeu8y5uo6oxEWkGRgDHuwTyKPAoQHl5eVLBX6ZwCJ8JFfMXwypcDbTfkqzhX3IkIF2e7+gMJNS5LCH/gYtTQm48wlDELYdC7lb5UAGEIxAucsuRQogMhoJBblp4TeejqAQGlbrada4fmRgTECEJUV5STnlJOQsrFl58PhqPcuTcEQ63Hqbp1F6aT+6hufkgzS1NTC6rTEssyST07voq6u02qOpqYDW4GnoS732ZWyqXsapyYIy+bYwx3SkIFzCueJy7Cm/MpzPynsk05DQC4xOWxwGHutvGb3IpBU6mIkBjjDHJSSahbwGmiEiFiBQCS4CaLtvUAF/x5x8Efn219nNjjDGp12OTi98m/iTwBu6yxR+oar2IfAPYqqo1wAvASyLSgKuZL0ln0MYYYy6XVOdcqroeWN/luX9OmG8DsjsUtjHG5DnrotAYYwLCEroxxgSEJXRjjAkIS+jGGBMQWettUUSOAQf6+N/L6HIXap7Ix3LnY5khP8udj2WG3pd7gqpecWTqrCX0/hCRrd31ZRBk+VjufCwz5Ge587HMkNpyW5OLMcYEhCV0Y4wJiFxN6KuzHUCW5GO587HMkJ/lzscyQwrLnZNt6MYYYy6XqzV0Y4wxXVhCN8aYgBjQCV1EFojIhyLSICL/eIX1RSLyqr9+s4hMzHyUqZdEuf9eRN4XkfdEpFZEJmQjzlTqqcwJ2z0oIioigbi8LZlyi8hf+p93vYj8ONMxploS3+9yEfmNiOzwv+OLshFnKonID0TkqIjs7ma9iMhz/t/kPRGZ1ac3UtUB+cB11bsXmAQUAjuB6V22eQL4b39+CfBqtuPOULnvBob484/nermTKbO/XTHwNrAJmJ3tuDP0WU8BdgDD/eVR2Y47A2VeDTzuz08H9mc77hSU+8+AWcDubtYvAl7Hjf42B9jcl/cZyDX024EGVf1YVS8ArwCLu2yzGPihP/9zYJ5Izg+22WO5VfU3qnrOX9yEG0UqlyXzWQN8E/gPoC2TwaVRMuX+a+B5VT0FoKpHMxxjqiVTZgVK/PlSLh8hLeeo6ttcfRS3xcCP1NkEDBOR63r7PgM5oV9pcOqx3W2jqjGgY3DqXJZMuRM9gtuz57IeyywitwLjVfWXmQwszZL5rG8EbhSR34nIJhFZkLHo0iOZMn8dWC4ijbhxGP42M6FlVW9/91eU1AAXWZKywalzTNJlEpHlwGzgzrRGlH5XLbOIhIBngYczFVCGJPNZR3DNLnfhjsTeEZFKVT2d5tjSJZkyLwVeVNVvichc3GholarqpT+8rElJLhvINfR8HZw6mXIjIvcAK4DPq2p7hmJLl57KXAxUAr8Vkf24NsaaAJwYTfY7vk5Vo6q6D/gQl+BzVTJlfgT4KYCqvgsMwnVgFWRJ/e57MpATer4OTt1juf3mh1W4ZJ7rbarQQ5lVtVlVy1R1oqpOxJ03+Lyqbs1OuCmTzHd8Le4kOCJShmuC+TijUaZWMmU+CMwDEJFpuIR+LKNRZl4N8JB/tcscoFlVD/f6VbJ99reHM8OLgI9wZ8VX+M99A/djBvdB/wxoAH4PTMp2zBkq9wbgCFDnP2qyHXO6y9xl298SgKtckvysBfgv4H1gF7Ak2zFnoMzTgd/hroCpA+ZnO+YUlPknwGEgiquNPwI8BjyW8Dk/7/9NdvX1+223/htjTEAM5CYXY4wxvWAJ3RhjAsISujHGBIQldGOMCQhL6MYYExCW0I0xJiAsoRtjTEBYQjd5T0Sm+32s13ezvlREPBEJwl25JsAsoRvjOjgD2NbN+lm4O/m2ZyYcY/rGEroxnQm9u75hbvOnltDNgGYJ3Ziea+iW0E1OsL5cTF4TkTBwBigCSrRzJKjEbT7CdVk7SV0XtsYMSFZDN/luOjAE+KCbZF4CTAZOWTI3A50ldJPvkj0huiMz4RjTd5bQTb7rSOh13axf5E+t/dwMeJbQTb7rOOF5vusKESkFlvmL3dXgjRkwLKGbvOWPQ3uzv7hURIoS1l0HvErnyOsfZjg8Y3rNErrJZ5W4YQwbgduBAyKyXkTexQ0FNgI3ZBjAGhF5KDthGpMcS+gmn3W0n9cC9+PGfKwCyoEXgHuA7+KaY0K4AY6NGbDsOnSTt0Tke7iBep9U1eezHY8x/WU1dJPPerrl35icYjV0k5dEpAA4i6vUlKhqW5ZDMqbfrIZu8tWncLf777ZkboLCaujGGBMQVkM3xpiAsIRujDEBYQndGGMCwhK6McYEhCV0Y4wJCEvoxhgTEJbQjTEmIP4fIZqwG6zII1MAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mus = arange(0,1,.01)\n", "\n", "# Plot for various values of a and b\n", "for ab in [(.1,.1),(8,4),(2,3), (1,1)]:\n", " plot(mus, vectorize(Beta)(mus,*ab), label=\"a=%s b=%s\" % ab)\n", "\n", "legend(loc=0)\n", "xlabel(r\"$\\mu$\", size=22)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now define the Binomial, which gives our likelihood (as a function of $\\mu$)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\operatorname{Bin}(m|N,\\mu) = $" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKgAAAAcCAYAAAD1EPTXAAAABHNCSVQICAgIfAhkiAAABYVJREFUeJztm1tsVEUYgL+2VLdZZI02IV641BsSIe22BkW81Fh48sFqbSRRo8XEKGriJSbqg0h8kMREE/DyojGKV0RUXtTEopAqokWsRkVCRA0UkSgVUQnq+vDPyY6zc7a7Z+ac3cbzJZtz+v/n/Oeff+bM/DNnCikpKSllGAQKwPWG/FrgT6A5cY/c8FqeRk9OpUQnD+wB+gx5JzACHEnco8poDZFHKU+YLWdORt6WfuTN+R34DJgFnAMMKdlmdW0YzwL7gGxcjtaILiQ+S0L0pyr9TcBhIKfpNgFPePbHV311Awst8qjl6VE2vXOpcug94CJgLrAd+AB4FzgPeaO+Bx4OsXE28A9whyHvA1YiBftVPWe1V++TYR0wCky26PqRijwa+A64RskbgYPADZ598VFfM4DHQ3Qu5VkFzKyoFFVwHzAGTNVkK5He8DhN9hzwfIiNd4ADQIsh34YE8yDwFRO3gc5DfL/XolsBDKvzR4A31PlsdU/esy8+6ms10uPacCnP6bZnuuagHcB64EdNNgN4DfhZk00HvrXcfwbSvb8C/GHoblf6KciQMVHZAnwN3Ag0GbouYKs6XwssAo5R8sPAFyE2lyEVXu7XbbnPtb7mAm1Ir2sjankAdgDTMBqx2UBvQwo3YDGQQ4biQU3WjuQrOnmLrB341GJzAGgAXrboNiinCxZdLak2RgAvIZXeY8jzFHucIeAXZBgeb4K0CumVyv22WO5zra+7kMYcRtTyBKwB7tYFk4wLOtVxmFI6kcYU6LJIUqwXpBVJrnVZG1JxtgL3AH9TGqB6ppoYBQyp40LgbXXehgyrQY9TQPLVPiSONvsB+9WvGlzrqxnoBZ4Ose9SnoBh4CEggyxJlTCCDLVmwwV5ewrAYvX3fKRx6cn/IkrXunqRHLPBsJcF/gI+r8DxbuonB60mRgE5Jdd7tT6kR8losouRWXQcEyTX+pqPlCGHHR/lySofuwOBPsRnkKFhBGk4Jl3qGLwJ7cgQ/Jt2TR7JM/SuvJ3ihEfnJCQnG63A8Xqh2hgFjCENYbpx7Zf8t6fYiMRzssWGK6711YHkqWMh9n2U5xCSH3fYlMFsM2wJYQey3GO+WVEJ3khb/mnSTbQedBfjTyb033j2XWK0G3ujniiswD5x8s124NHgD32YKpdb5ZD8ZRP+Ji3BrD1T9io3dhKSy4SwZxy9S4xaKF2pmEjkkJcvbsbQ0ohKG+gCpFfYatFFZZ86Hu/Rpsklnu1FjVEjcCzJ9EBxMYlkPrseQWuXeg4aBH+35aar1NFnXjQK/ET4om89EjVGs5DGuy0OpxLiELImHTdT1LOAYgNtBuao8/ONG/qBq9V5JTPuSikgSXQrcJpHu3HhEqNz1XFDDH4lxSjJNdCSDiCPNJi9SM72AvIJbCOyRLBX6ddSDLYPFiu7Sy26y4Bn1O8tdd1OTRb2rTguXGL0IjJBmpaUszHQi3wNsi2v+dqE0oD0nleaiiXqATcjFb9fXTiIbOa4U/39CXBitSUrw1FIxX5k0S2j/Ix7l0c/KiFqjHLI5Oj1JJ2NgalI+edYdD42oYB82i4AJ5iKx5RiXlTvHbiHeDZG+CZqjG5V913g3aPk2Uxxh5KOj00oIHm87RMtHyJDkLmjKAkyyNas9TV4djVEiVELsnT1aiweJc8AsnfXZA2la8hvAk8asveBB8vYfwrL3tkmZGgqt9Mkbi4E7qd+NyxHjdFsJFWZ6dmfWtGErFKYe1u/AW4xZD8A1xmyA8AVIbazwMdYctyzqJ/v3PVKGqMivUjaEhB8P1+gyVqReLVrsjYlOyXE7lLgcn9upvyfWU5xIui6CQVkUrTc9qD0n+ZSovAAxQ8srptQAM5UNlNSUlJSUlJSUlJqzr8rGcQSHZqrgwAAAABJRU5ErkJggg==\n", "text/latex": [ "$\\displaystyle \\mu^{m} \\left(1 - \\mu\\right)^{N - m} {\\binom{N}{m}}$" ], "text/plain": [ " m N - m ⎛N⎞\n", "μ ⋅(1 - μ) ⋅⎜ ⎟\n", " ⎝m⎠" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "N, m = sym.symbols('N m')\n", "\n", "Bin_norm = sym.binomial(N,m)\n", "Bin_f = mu**m * (1-mu) ** (N-m)\n", "\n", "Bin = sym.Lambda((m, N, mu), Bin_norm * Bin_f)\n", "\n", "maths(r\"\\operatorname{Bin}(m|N,\\mu) = \")\n", "display(Bin.expr)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Ignoring normalisation for now, we can find the posterior distribution\n", "by multiplying our Beta prior by the Binomial likelihood" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\operatorname{Beta}(\\mu|a,b) \\times \\operatorname{Bin}(m|N,\\mu) \\propto \\mu^{a + m - 1} \\left(1 - \\mu\\right)^{N + b - m - 1}$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p = Beta_f * Bin_f\n", "p = p.powsimp()\n", "\n", "maths(r\"\\operatorname{Beta}(\\mu|a,b) \\times \\operatorname{Bin}(m|N,\\mu) \\propto %s\" % sym.latex(p))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let $l = N-m$ (the number of tails, if we are talking about coin tosses) then we can write it as" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMsAAAAcCAYAAADPwEukAAAABHNCSVQICAgIfAhkiAAABXBJREFUeJztm3uMXUMcxz/bVdmHZEUrXq2sWH1kpfuoFCmbFaQRiVRUqwmJ6B+iKAkipNFVVIW/iOcfrEc1qyqNhnrEStsUDVWqbKilHq3Vx0ZZRdH1x3fO7nH2nHvuY+69596dT3Iy9/xmzsxv5je/ee0sOBwOh6Mo3Ae8XYByHgDWFaCcIcYUsrASpQ1YC+wGBoE5xVUn8TQDnxaonE8s5hdr52I7SyfQUWQd4qhFxr+h2IqUCLl04k7S7w9NOZQTRqydbTjLXcA2YADYCzwNVFnI188E5O1zgW7gIKrYZOBMYJORfWDS2mQdsBh4Jc30zwF7UOOXG9ORHRZExB8LnAgcBt5h2E5nWdbjeOA47DpLrJ1zdZYKoBK4DmgE5gOzgJtzzDdIswkXAncjB6kCngGWAbcBM5GhbJedCWcAVwLLgd8DcXOAR4CNwK+o071QUO1yZwuwBrgXOCokvsWEtwJLgVa0rOkCjrCoRwtyxB0W84wl1woMAkt8798BrwFTcsw3SBPqYPOAn43sLfM+Beg3svXACZbLzoRlSM/HQ+IWo3oMAD9iv40Kxf3AZmARqq+fZuAQcCnqCwC3o9mlHvjakg7NaDVz2FJ+BWEi8DDwGeqwA8DfwD0R6e80aQZ8aQ8FZOeGfLeKkaPwq8ATAdl6NOqF0YGcO9XTHvGtR6oN/iRkvKci4s8DTkOzcTulObN49CBnqAzIVwIrArJTUV0nheSTbX/oInxAgvzbeYhFJuE1IXF1qDN0m/dxQB/wElp6NQJTUQWjCjoGaPA9q5Gz+WXVId99xciN1w/A1QHZL8BlEWWPR6N5qqcm4luPVI243MSfH5MHJMtZMrG5xxLzzayAvActwfxcjuwStorJtj98CVwbWps82jlYgVYTbgnJoBWNil7cxWjfMM9kDuq8tcDWCCX6GV4yAfxm3lNNz7VodPLnOR5t5P2yU5Bxo8reZ558cQHwLzpkKCUysbnHJhNeCLxpftegmXOsL10FcBPwPPBPSP7Z9Ica5ERRm/u82TnMWf4EPg9JO92EH5twP9rkzUbLsIvQtHoA+MaijtNM6D+7bwX+Ar7wyZpN2d9aLBtUxwbfe70pqx/43shqjayHkRv7pJOJzT0+NGGbTzYNDZpXoZloP1oSTQAusaQraN8H6nM2ibWz/zSsCi2jthE+CngN540yrwNPAs8C76FRZQXy+MERX2dPEzr1GPDJWoDtaI3rT2e7bNAJ11aGZ6wHze+lvjQnofX7T5bLzjeZ2tzjAHKwk32yZqAXHWSsMnlWo2PjfuzRhJblBy3mCenZeYgZqKM9FpHZDnTSU2FXx7LgbNR2XWmmbye7PctO4jev/icu/1xsvotwBytb/MuwVGvXOrRv2Ij9kbsc+MOEtv8YG6QXjejpsjsmPhebVzNc71FBus4yE40uwbWrQ+wx4bg8l5POSVsmZGvzMcDR2N8fJhr/nsVruF0h6a4wYVijOrRX2Yuu35QS2dp8MnIkm9dNEo/nLGOB083vcwJp5qIrHGD/BKJcGAQ2oCPthpi0SSEXm3t3vd7Ng16JpwUZvA+tiV9E95g2oLPvPhO/GvuX4sqF+aiNro+In41u1XYCb5i0vT7ZQ3nWL0guNl+JNvcTC6VskliAGmYhMto+9PeCbnSkdot5/whdVnSM5EjUwTZHxHeQ+uRqZ941/D/Z2rwObezXFFLZJPEoargZxVakxLkDtWNLXMIEkK3NbzTfhd3ZGhW8j6bVsHs4jvSpQhcM1xZbkTTIxubV6Dj65bxoVAJUoul2e7EVKRPa0EXDJP/zV7Y2n4qWk/WW9SkZGknODVhHYXA2dzgcDofD4XA4HA6HwzE6+Q/uP4r5myRwtAAAAABJRU5ErkJggg==\n", "text/latex": [ "$\\displaystyle \\mu^{a + m - 1} \\left(1 - \\mu\\right)^{b + l - 1}$" ], "text/plain": [ " a + m - 1 b + l - 1\n", "μ ⋅(1 - μ) " ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "l = sym.Symbol('l')\n", "p = p.subs(N-m, l)\n", "p" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "So the posterior distribution has the same functional dependence on $\\mu$ as the prior, it is just another Beta distribution" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle p(\\mu|m,l,a,b) \\propto \\mu^{a + m - 1} \\left(1 - \\mu\\right)^{b + l - 1}$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "maths(\"p(\\mu|m,l,a,b) \\propto %s\" % sym.latex(p))" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle p(\\mu|m,l,a,b) = \\operatorname{Beta}(\\mu|a+m,b+l) =b$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPkAAAAnCAYAAADNe+7PAAAABHNCSVQICAgIfAhkiAAAB+hJREFUeJzt3XusHGUZx/FPQbAUjBdqQMKlSAWU9pxtxVijksYo/IEmmKhVEcXgXbzEVE2QaMVbCN6Dghr1EEVEUWMVseIlGBUqIloqJiog2mIJpaBWQDDUP57ZnDmzM2d3Z/bsnu55v8lkzrw78z7zPPPeZ/d3SCQSib2IjbgHV4z6RgbI83HRgPK5sOa14xjXBcM+c5TvWmyYo7xn4xN45QjsziUTuLGH89aaPeYT+H3NexjHuC4Y+q3kL8dm3IRNWFLT7vfwSVyHP+Pp+A7+ivU184Sf4d8Vnz0Wd+KYBvnPJVfgHSXpEzgWv8EfsbJm/pPYUvPa2eLajVHG/UO4egh2zsdVQ7DTjdIy1G8l3yQq5ErcgufVvJkVuBlrRAH6GF6F5+Csmnl24xz8UNx3m5NEg3MH9uBFc2S7F87DuXh0IX1C3N+JojC9q2b+K7C19t2VMyXiVtzyw/qyuA+Llvqjl37t/G4IdrpRWoYe0WcmZ+HF2B9H4LuFz3+MpTgo207L0k8TvTQ8SjQun8+OH8Bn8K/smn+W2K0qnCeLCtCNJXgtXlBIP1AUgi/jWz3kM5dswW14hYgHPBIHiFFP+5x1het6ifliLMJ/ctc1jWne/hmFtAeyfVXcmzAl/NrQw7ktXDYEO5P40oDzrENZGZrRkz8DD4uC0mapaJknRU/7FNH7TWIn/lAw8lwR2Nfg4uzvlunCRvQov644ruptVlRsvRbGUzPfflFIv0q0fN/uMZ+6dIttm41iStTmBDFEfzg7XqWzx+g15sVn1TSmbf6LHYXt3uyzqri/VxTI3bhLVJDFfdrtxuNxWGb/J7hPNOhrBmznUBxisD354aJsvAQ/NX3vx4mR9C+ztOuyc/MUy9CMSt4S8+PdubRVeFAMrVfgWtyP1+NgbKvhwAoz54ZH49bs7wkx3x80z8INInCjoFts22zG00TvTcTjidhPNApvwKdq2G8yH29CWdwXYV+8UTRiL8MpePuAba/K9uvFMHa1aMAu1/8Itpud+8TzHRStbP8mvF9U7MVixPlhvBPPFI1YMW7FMjTD2RZ+W2LsZjyEr4jh+eliHl230KwULSs8QbT8e3KfXVIzX2LoOCmG4dvE1OJaLNN/DzVIusW2zR2iQh8m5rArxTTievGs1ot49cuE6BHqUhXXbizTGfc9eF/u+HZcieMb3F8ZLdGIvjCzAe8WPeIy/GWAdraYHm0NgkkxfV0nFi3hR9nx8diVpV0j6lCeYhmawWadizpfE63H3s4mfK7LOb0svG1QvtCU39aWXNdrbJ+U5bGiy33MF6bw/Vk+L4v7Efi0GLHtEqObh/CBijzOyc7ZnTv3wULas0uuuwyXFtKOEfE9doB2Llf9PYa6eX4TXy2kbRTTsTzX4IOFtI4y1O7J98kS31O44Kn4bIUDexM7xaucplyIr3c552+F435i+7hsf1etu5t/FON+sBiV/FyMSraJHvB61avgF+MbuePzsV00FG22l1zXwhcLaavFwu6tnac3svPxshtvkOdk4RxiWlBcOyo7r6MMtSv5cWIlND+0WilavPnwaqApN+LMAeSzM9v6oZ/YrszOu9N4UIz7qWJuuc70FO1MMQ2o+sLPLtPDU+J9/S6zD7eXiB5tv1zaIrxNTDv/N0A7y1XXkTp5HihGHPl4LBULbPm0o8WrsmLcOspQe+GtPdE/WwTnZNM91v6z3NDewiY8WQQrz0GmV6OJuVoLRw7Qdj+xPUm8Ux4XinG/2/RrvuV4Cz6iunety4RoRM4Qi1bLReU+3Mz1gKa034wMcrF4ItvnRzarxVuM/CJtS8TttsL1HWUoX8mvFkHYio+KV0v3iAext3OTmBe/tJB+omgJ263hBdnf5w3Qdq+xPUAU/i8M0PaoKcb9B2KOfgl+JRq9S0VPOMg3Hy2x6HSumN9uEfFdY2bP2pRJ/Emsrg8yz7I3MVvNXKSd1Bm3WcvQJtGijjOniAey75Dt9hrbs8UK6rgxqrgvRGYtQzvE+8px5604asg2e43t68T8fRwZRdwXIpVl6BDR5Z8w1NtZGKTYJhKJRCKRSCQSiUQikUgsUBZl+1H9OiuRSCTmPfNdVqpIlcxUL5T5OgqJpSY+JBIdTJld4ugCzX4eO2wmxLe+ijJTU7pLOZX5eqX4Ft9ck/9SUZUPiUQtpkRPdWhhe4z4gcK9yn8u2MTehgHmV8YNeHOJ3So/qfZ1u5AbqsOU3n3dYab0VZkPiRLmSpJ53KiSOKqSN2I4Ekd1ZYI6JIIy+pVyGqXEUpUPiQKpkjejSlZqWBJHdWWCOiSCeqDM11FKLNXxYUEyyAexEFmmXFZqWBJHdWWCZpUIqmCZTl9HKbFUx4cFSerJm7HYtPxwnjKJo1fj7yXnFiWCTi9Jq5rzt4RufF5k4iihIJL/SeWRZv7u+P5s308vWObrKrEwd3surS37XFa26vpapmtex4cFSerJm1EmK9WvxFFdiSDqywTVkZkq83WUEkvjJpU1Z6RK3owyWal+JY7qSATRTCaojsxU0ddRSyyNm1TWnJGG680ok5UapsQR9WSC6shMFX0dtcTSuEllzRmpkjejTFZqWBJHdWWC6spMFX0dpcTSOEplJeYxe5u8UROZqfni67hKZc0Jo35Y48AtYi76D+X/rHG+0RL/DODuGtfOF1+b+JBIJBKJRCKRSCQSiUQiMRT+D4Nxnk6pkdtSAAAAAElFTkSuQmCC\n", "text/latex": [ "$\\displaystyle \\frac{\\mu^{a + m - 1} \\left(1 - \\mu\\right)^{b + l - 1} \\Gamma\\left(a + b + l + m\\right)}{\\Gamma\\left(a + m\\right) \\Gamma\\left(b + l\\right)}$" ], "text/plain": [ " a + m - 1 b + l - 1 \n", "μ ⋅(1 - μ) ⋅Γ(a + b + l + m)\n", "────────────────────────────────────────────\n", " Γ(a + m)⋅Γ(b + l) " ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "maths(\"p(\\mu|m,l,a,b) = \\operatorname{Beta}(\\mu|a+m,b+l) =b\")\n", "Beta(mu,a+m,b+l)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's see how this works graphically, with prior a=2, b=2 and likelihood N=1, m=1" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\mathbf{Prior}: \\operatorname{Beta}(\\mu|a=2,b=2) = 6 \\mu \\left(1 - \\mu\\right)$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEQCAYAAACgBo8fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXiU1dnH8e8dIFhFESVaKyigYEXEoiOtS20FX0Gq4NYWrBtiURR9cde6VKXKK7gjyOLWulNRi0vrikstqEEUQYoiKEZRggpqFVly3j/uUCMkZAgzc2ae+X2uKxdM5mGe+yHJL2fOcxYLISAiIoWvJHYBIiKSGQp0EZGEUKCLiCSEAl1EJCEU6CIiCdE41olbtmwZ2rRpE+v0IiIFadq0aYtDCGW1PRct0Nu0aUN5eXms04uIFCQze7+u59TlIiKSEAp0EZGEUKCLiCSEAl1EJCEU6CIiCaFAFxFJCAW6iEhCRBuHLhLdihWwYAEsXAiLF8Onn8IXX8C33/pHCFBaCk2bwiabwJZb+sfWW0ObNv45kTyiQJfkW7UKZs6E6dNhxgz/eOcdqKiAqqqGv25ZGeywA+y6K3TuDLvtBnvsARtvnLnaRdaDAl2Sp6oKXnsN/vEPePFFmDIFvvzSn/vBD6BTJ/jFL6BtW29pb7uth/OWW8Jmm3mLvLQUSkpg+XL/+Oorb8EvXgwffwzvvQfz5/svhgcfhPHj/fUbN4YuXWDffeHAA/08P/hBrP8JKTIWa8eiVCoVNPVfMmbZMnjySZg4Ef7+d6isBDNvPe+zj3/suae3qBs1yuy5Q/CQnzYN/vUveOkleOUVr2mjjeCXv4TDD4fDDoOWLTN7bik6ZjYthJCq9TkFuhSsqip47jm44w7429+8/7tFC+jVCw46yFvIZbWuYZR933wDL7zgv1wefRTefdd/kXTrBsccA0ccoa4ZaRAFuiTLJ594F8dtt3m3R/PmHpC/+Y0HZpMmsSv8vhDgjTdgwgS4/36YN8+7do46Ck4+2fveRdKkQJdkeP11uP56uPde79fu1g0GDPCujELpp66q8n79W2+FBx7wlvz++8OQIfCrX2W+O0gSZ12BrnHokv9eegl69vSbjQ88AAMHwpw58Mwz3sotlDAHv9H6i1/AX/4CH34Iw4fD3LnQp4/frL3rLli5MnaVUqAU6JK/pkyB7t19xMhrr8GwYT7UcORI6NAhdnUbrkULOOcc74K57z4fWXPMMfDjH3vgr1oVu0IpMAp0yT+zZ3s3yt57w6xZcO213ld+/vmw+eaxq8u8xo3ht7/1cfIPP+z3BI47zt+RPP6498GLpEGBLvnjs8/g1FO96+GZZ2DoUB8dcsYZxTErs6TEu15efdVvnn79tferH3CA/2ITqUe9gW5mt5nZIjObWc9xe5rZKjM7MnPlSVFYtcpHrXToAGPGwKBBHuQXXVQcQb6mkhIfsfPWW969NH26j4Q580xYujR2dZLH0mmh3wH0XNcBZtYIuAp4IgM1STGZNcv7yAcOhI4dva/8ppvijR/PJ6WlMHgwvP22j+a5/nrYeWd46KHYlUmeqjfQQwgvAJ/Vc9hpwERgUSaKkiLw7bdw6aXeT/zOO3DnnfD88xqTXZuWLWHsWJ99uvXWPuv0iCN8UTGRGja4D93MtgUOA8akcexAMys3s/LKysoNPbUUqhkzoGtXuOwy71qYPRuOPtqn6kvdUikP9WHD4LHHYJddfLKSSLVM3BS9HjgvhFDvGKsQwrgQQiqEkCrTW+ris2oVjBjha6p88gk88oiPu9b3QvqaNPHRPjNm+D2H3/7WfxkuWRK7MskDmQj0FHCfmb0HHAmMNrNDM/C6kiQLF8L//A+ce66P3HjzTTj44NhVFa4OHeCf/4TLL/cx7J07+8JgUtQ2ONBDCG1DCG1CCG2AB4BTQggPb3BlkhxPPQU/+QlMnepT3idOVKs8Exo3hosv9glYTZrAfvvBVVdt2BrvUtDSGbZ4LzAF2MnMKsxsgJmdbGYnZ788KWirVsEll0CPHh7g5eVwwgnqK8+0Pff00UGHH+7dMQcf7GP6pejUu8FFCKFfui8WQjh+g6qR5FiyBH73O5/pePzxMGqUlovNpubNfTLSL3/pC33tuacPb+zcOXZlkkOaKSqZN3OmB8pTT8HNN/sytwrz7DODU07x4Z/ffAN77eUhL0VDgS6Z9fjjHiRffQWTJ/t63+piya299vLdk7p0gb594Y9/1HowRUKBLpkzciQccgi0b+/rkeyzT+yKitc228Czz0L//j4S5ne/8y3xJNEU6LLhVq2C00/3j4MP9q3XWrWKXZWUlvqoomHDfFOQbt18k2tJLAW6bJhly3xyy8iRvnjUgw9Cs2axq5LVzHzkywMP+CJf++wD770XuyrJEgW6NNySJT4kceJEX7P8mmu0hVq+OuIIv0m9aJH3sb/xRuyKJAsU6NIwH3/sE1mmTIF77vE1yyW/7buvzy5t3Ni/di+8ELsiyTAFuqy/BQs8EN591xeJ6pf2VAWJbZddfImAH/3I92l9QiteJ4kCXdbP3Lnw85/74lpPPunrs0hhad3ax6p36AC9e/u2d5IICnRJ35w53jL/z398jLmGJRaurbbyr+Huu8ORR2oCUkIo0CU9c+b4tPJVq7x1t/vusSuSDdWihb/L2ntvH6eutdULngJd6jdnDuy/v4f5s896P6wkw6ab+uzevfeGo45SqBc4Bbqs27vv+oSUlSv9LbrCPHmaNftuyYajjvJhqFKQFOhStw8+gO7dff9PtcyTbXWod+3qo5b+/vfYFUkDKNCldp98AgccAJ9/7v2snTrFrkiybXX3S6dOvrb688/HrkjWkwJd1vb553DggVBR4T/gugFaPDbf3Memt23r6/K8+mrsimQ9KNDl+775xscmz57t45M1NLH4lJX5MgEtW0KvXvD227ErkjQp0OU7K1f6TbGXXoI779SkoWK27bbfzSLt0cM3+Za8l86eoreZ2SIzm1nH878zsxnVH/8ys90yX6ZkXQhw6qneKr/hBl9BUYpbhw7e5VZZCQcdBEuXxq5I6pFOC/0OoOc6np8P/CKE0BkYCozLQF2Sa1dcAePGwQUXwGmnxa5G8sWee/qSyLNm+YqNy5fHrkjWod5ADyG8ANS5hXgI4V8hhM+rH04FtLNBobnnHrj4Yjj6aA92kZoOPBBuuQWeeca3FNR2dnmrcYZfbwBQ5wBWMxsIDATYbrvtMnxqaZAXX/Rtyvbbz39otf+n1Oa443yS2dCh0K4dXHRR7IqkFhkLdDPbHw/0fes6JoQwjuoumVQqpV/zsb3zDhx6qA9Re+ghaNo0dkWSzy67DObP93dz7dr5DXTJKxkZ5WJmnYFbgD4hhE8z8ZqSZUuW+IbOJSV+42uLLWJXJPnOzN/F7bcfnHACvPJK7IpkDRsc6Ga2HfAgcEwIQQNWC8HKldC3L8yb5+t2tGsXuyIpFE2b+vfMNtv4u7sPP4xdkdSQzrDFe4EpwE5mVmFmA8zsZDM7ufqQS4AtgdFm9rqZlWexXsmEc87xMcY33+ytLZH10bIlPPIIfPmlh/rXX8euSKpZiHTHOpVKhfJyZX/O3X67v10eMgSuuy52NVLIHnkE+vTxd3t3360b6jliZtNCCKnantNM0WJSXg6DBvmiWyNGxK5GCt0hh/gw13vvheuvj12NoEAvHosW+Qp6P/wh3Hef7/wusqHOP98nHJ1zjq+XL1Ep0IvBypU+lb+y0ocnbrll7IokKcy8G2+nneA3v4EFC2JXVNQU6MXgD3+A557zqf1dusSuRpJm0029obB8ubfWv/02dkVFS4GedA8/7P3lp5wCxxwTuxpJqg4d4M9/9vs0Z50Vu5qipUBPsnnz4PjjIZWCa6+NXY0k3aGHepiPGgX33x+7mqKkQE+qZcvgyCO9j3PCBE3rl9wYNgz23htOPBHmzIldTdFRoCfVmWfC9Onwl7/4Wi0iudCkibfON9rIGxTffBO7oqKiQE+iiRN9FuhZZ/lYYZFcatXKd7yaOdMbFpIzCvSkee89GDDANya48srY1Uix6tnTx6aPGQMPPBC7mqKhQE+SFSt8SdMQfPJQaWnsiqSY/elP0LWr96e/917saoqCAj1JLr0Upkzx8eZaQVFiKy31ZQFCgH79fIKbZJUCPSleeMFHGPTvrw2eJX+0awdjx8LUqd5il6xSoCfBkiU+aahdO7jxxtjViHxf376+X+3Qof4OUrJGgZ4Ep57qGw3cfTc0axa7GpG13XQTbLedB/uXX8auJrEU6IXunnv845JL4Kc/jV2NSO2aN4e77vKbo6efHruaxFKgF7KKCl+jZa+9fAEukXy2zz7+fXrHHb7GkGRcOlvQ3WZmi8xsZh3Pm5ndaGZzzWyGme2e+TJlLSH4cLAVK3w2qNY3l0Jw8cW+4udJJ/lyzpJR6bTQ7wB6ruP5g4D21R8DgZs3vCyp17hxvi/oiBGw446xqxFJT2mpN0CWLIGTT/aGiWRMvYEeQngB+Gwdh/QB/hLcVGBzM9smUwVKLd5916f1H3CA/1CIFJJOnXzEy4MP+v0fyZhM9KFvC3xQ43FF9efWYmYDzazczMor9XarYaqqfJPnxo3httugRLdBpACddZavyjh4MHz0UexqEiMTaVDbVt+1vo8KIYwLIaRCCKmysrIMnLoIjR7tk4iuuw5at45djUjDNGrkN0eXLfONy9X1khGZCPQKoGaytAL0Kzcb5s/3TXl79PCNK0QKWfv2Pnt00iRfIkA2WCYCfRJwbPVol58BS0MICzPwulLT6lEtJSUwfrxvXCFS6IYMgZ/9DE47DT75JHY1BS+dYYv3AlOAncyswswGmNnJZrb6btzjwDxgLjAeOCVr1Raz8ePh2Wfh6qvV1SLJ0aiR3wv66ivvT5cNYiFS31UqlQrl5eVRzl1wPvwQdt7Z1zh/+mm1ziV5hg3zSUcPPgiHHRa7mrxmZtNCCKnantMQiUIweLAvPTpunMJckunss2G33XxdoqVLY1dTsBTo+e7BB32a9KWXwg47xK5GJDuaNIFbbvF+9PPPj11NwVKg57MlS7zF0qWL9maU5Eul/CbpmDHw4ouxqylICvR8dt55sGiR3xDVWi1SDC6/HNq0gYED4dtvY1dTcBTo+epf//I+8yFDYI89YlcjkhubbOIt9H//G4YPj11NwVGg56MVK3yNltat4bLLYlcjkls9evg2ildcAXPnxq6moCjQ89ENN8Cbb8LIkdqBSIrTtddC06a+3r+WBUibAj3fLFgAf/wj9O4NffrErkYkjh/9yFvoTz0F998fu5qCoUDPN6u359Jmz1LsBg3y+0dDhviIL6mXAj2fPPYY/O1v3kLffvvY1YjE1aiR3yBdtMh/JqReCvR8sWyZt85//GNvkYiIj00/6SS46SaYMSN2NXlPgZ4vhg+HefP8G7e0NHY1IvnjiiugRQufZKcbpOukQM8H8+f74kS//jV07x67GpH8ssUW/vPxz3/C3XfHriavKdDzwRlneH/htdfGrkQkPw0YAF27wjnnaPGudVCgx/bEE34j9KKLoFWr2NWI5KeSEhg1yhfvGjo0djV5S4Ee04oVfgN0xx29lS4idUulfIP0G26AOXNiV5OXFOgxjRrla1asnhUnIut2xRWw8cZafbQOaQW6mfU0szlmNtfM1lqs2My2M7PJZjbdzGaYWa/Ml5owlZW+xnmPHnDwwbGrESkMW28Nl1wCjz/uH/I96ewp2ggYBRwEdAT6mVnHNQ67CJgQQugC9AVGZ7rQxLnoIvjPf+C667QLkcj6OO006NDBuymXL49dTV5Jp4XeFZgbQpgXQlgO3AesuchIADar/ntz4KPMlZhAb7zha5wPHux7hYpI+kpLvSH09ts+b0P+K51A3xb4oMbjiurP1XQpcLSZVQCPA6fV9kJmNtDMys2svLKysgHlJkAI3v+3xRb+1lFE1l+vXtCzp494+fTT2NXkjXQCvbb+gDWna/UD7gghtAJ6AXea2VqvHUIYF0JIhRBSZWVl619tEjz6KDz7rPeft2gRuxqRwnX11fDFF9ozoIZ0Ar0CaF3jcSvW7lIZAEwACCFMATYCWmaiwERZscJ3N99pJ1+fQkQabpddfKu60aN9tJikFeivAu3NrK2ZleI3PSetccwCoDuAme2MB3qR9qmsw5gx3u939dW+y7mIbJjLLvNt6849N3YleaHeQA8hrAQGA08As/HRLLPM7HIz61192FnA783sDeBe4PgQtIrO93z+uXezdO8Ov/pV7GpEkmGrreDCC+GRR+CZZ2JXE53Fyt1UKhXKy8ujnDuKc8/1lvlrr8FPfhK7GpHkWLbMR4u1aAHl5b5MQIKZ2bQQQqq255J95fni/fd9B6Jjj1WYi2TaRhv5DNLp0+Gee2JXE5UCPRcuusj/1KJCItnRty/svrt3vyxbFruaaBTo2TZ9Otx1ly/C1bp1/ceLyPorKYERI3yT9ZEjY1cTjQI9m0Lw9Zu33BIuuCB2NSLJ1q2bTzi64oqinWykQM+mJ5/0O+8XXwzNm8euRiT5rroKvvwSrrwydiVRKNCzparKW+Vt2sDJJ8euRqQ4dOrkgw9GjfLulyKjQM+Wv/7V+8+HDtVa5yK5dNll3t1ZhEsCKNCzYcUKH9my667Qr1/sakSKy3bbwamnwh13wFtvxa4mpxTo2XDrrTB3rvfjNWoUuxqR4vOHP/iSAKuHDBcJBXqmff21v9Xbd19N8ReJpWVLH2H20EPw8suxq8kZBXqm3XgjfPwxDBumnYhEYjrjDF/rpYiGDCvQM2npUhg+3MfC7rtv7GpEiluzZt71Mnmy70FQBBTomXTttb6q4p/+FLsSEQHfd6BVK18SoAgWgFWgZ8rixR7oRx4JXbrErkZEwBfuuvhimDoVHnssdjVZp0DPlOHD4T//KcqxryJ5rX9/2GEHD/aqqtjVZJUCPRMWLvTdx48+Gjp2jF2NiNTUpIlvLvP66zBxYuxqskqBngnDhvlkoksvjV2JiNSmXz9vbF1yCaxaFbuarFGgb6iKChg71t/WtWsXuxoRqU2jRt7g+ve/4f77Y1eTNWkFupn1NLM5ZjbXzM6v45jfmNlbZjbLzIpn25Bhw/zu+YUXxq5ERNbliCN8OY7LL09sK73eQDezRsAo4CCgI9DPzDqucUx74AJgnxDCLsCQLNSafxYsgFtugRNOgO23j12NiKxLSQn88Y8wZw7ce2/sarIinRZ6V2BuCGFeCGE5cB/QZ41jfg+MCiF8DhBCWJTZMvPUlVd66/wPf4hdiYik47DDoHNnH422cmXsajIunUDfFvigxuOK6s/V1AHoYGYvmdlUM+tZ2wuZ2UAzKzez8srKyoZVnC/efx9uuw1OPNFXdxOR/FdS4n3pc+fC3XfHribj0gn02hYkWXPKVWOgPfBLoB9wi5ltvtY/CmFcCCEVQkiVlZWtb6355corfa0Wtc5FCsuhh/rkv6FDE9dKTyfQK4Cauxu3Aj6q5Zi/hRBWhBDmA3PwgE+mBQvg9tu9dd6qVexqRGR9mPnwxXffTVxfejqB/irQ3szamlkp0BeYtMYxDwP7A5hZS7wLZl4mC80rV13lf553Xtw6RKRhevf2vvQ//SlRI17qDfQQwkpgMPAEMBuYEEKYZWaXm1nv6sOeAD41s7eAycA5IYRkbrv94Yc+sqV/f/WdixSqkhJfCuDtt327yISwEGkFslQqFcrLy6Oce4P87//C6NHwzju+AbSIFKaqKh+XDvDmmx7yBcDMpoUQUrU9VxhXkC8WLoRx43xXcYW5SGFb3Up/6y148MHY1WSEAn19XH21r9mikS0iyfDrX8NOO/mIlwSsl65AT9fixTBmjC/ys8MOsasRkUxo1MgbaDNmwKOPxq5mgynQ03XDDb4BdBHtTyhSFPr18y7UK64o+Fa6Aj0dS5fCyJFw+OFa71wkaZo08SHIL79c8HuPKtDTMXq0h7r6zkWS6fjjYZttvJVewBTo9fn6a7juOujZE/bYI3Y1IpING20EZ58NkyfDlCmxq2kwBXp9xo+Hykqtdy6SdAMHwpZbFnQrXYG+LsuXwzXXwM9/DvvuG7saEcmmZs184uBjj/molwKkQF+Xe++FDz7QyBaRYjF4sAf78OGxK2kQBXpdqqp8Ea7Onb3/XESSr0ULOOkkuO8+mD8/djXrTYFel0cegdmz4fzzfblNESkOZ5zhywJcc03sStabAr02Ifjmz23b+tRgESke227r6zXdeissKqzdNBXotXnhBZ9kcPbZ0Lhx7GpEJNfOOQe+/RZuvDF2JetFgV6b//s/KCvzNc9FpPjstJPPDB81Cr78MnY1aVOgr+nNN+Ef/4DTT4cf/CB2NSISy3nnwZIlvqFNgVCgr+nqq2GTTeCUU2JXIiIx7bkn7LefzxRfsSJ2NWlJK9DNrKeZzTGzuWZ2/jqOO9LMgpnVuptG3quogHvugQEDYIstYlcjIrGde67PRZkwIXYlaak30M2sETAKOAjoCPQzs7WWHDSzTYHTgZczXWTO3HCDj3A544zYlYhIPjjoIF9hdcSIglhaN50WeldgbghhXghhOXAf0KeW44YCw4FlGawvd5YuhbFjfZiitpcTEfDx6GefDW+8AU89FbuaeqUT6NsCH9R4XFH9uf8ysy5A6xDCOrf8MLOBZlZuZuWVlZXrXWxWjR3rd7PPOSd2JSKST446ypfWHTEidiX1SifQa5sm+d/3HmZWAlwHnFXfC4UQxoUQUiGEVFlZWfpVZtvy5d7d0r077L577GpEJJ80beqLdj39NLz+euxq1imdQK8AWtd43Ar4qMbjTYFOwHNm9h7wM2BSQd0YnTABPvoIzqr3d5KIFKOBA33023XXxa5kndIJ9FeB9mbW1sxKgb7ApNVPhhCWhhBahhDahBDaAFOB3iGE8qxUnGkhwLXXws47Q48esasRkXzUogWccIKvwPrRR/UfH0m9gR5CWAkMBp4AZgMTQgizzOxyM+ud7QKz7vnnYfp0OPNMvwEiIlKbIUNg5Uq46abYldTJQqShOKlUKpSX50Ej/pBDfN2W99/XzFARWbcjjvBt6j74wLtgIjCzaSGEWru0i7tJOmcOPPqozwpVmItIfc48Ez7/HO64I3YltSruQL/+er+DrWn+IpKOvfeGn/7Us2PVqtjVrKV4A/3TT+HPf4ajj4attopdjYgUAjNvpc+d6+/u80zxBvr48fDNN36jQ0QkXYcfDq1b+9yVPFOcgb5ihd+pPuAA6NQpdjUiUkgaN/bNpCdP9iUB8khxBvrEifDhh2qdi0jDnHgibLxx3rXSizPQr78e2rf3ldRERNbXFlvAccf5ctt5tO9o8QX61Kk+7vz00zWRSEQa7vTTfd/RsWNjV/JfxZdo118PzZvD8cfHrkRECtmPf+zv8keP9mDPA8UV6BUV8MADviNRs2axqxGRQjdkCHz8Mfz1r7ErAYot0MeMgaoqv0MtIrKhDjgAdtoJRo6MXQlQTIG+bJn3dR1yCLRtG7saEUmCkhI47TR45RW/Nxe7nNgF5Mz998Pixf6fLyKSKcceC5tumhet9OII9BD8P3vnnX1XIhGRTNl0U+jf3zfK+fjjqKUUR6BPmQLTpnnr3GrbUU9EZAMMHuwz0CMPYSyOQB850ocqHnNM7EpEJIlWT1QcM8b3KI4k+YG+cKEPVTzhBA1VFJHsOe0073KZODFaCWkFupn1NLM5ZjbXzM6v5fkzzewtM5thZs+Y2faZL7WBxo3zbaO05rmIZFOPHrDjjjBqVLQS6g10M2sEjAIOAjoC/cys4xqHTQdSIYTOwAPA8EwX2iCr+7R69vT/aBGRbCkpgUGD4KWXoq3CmE4LvSswN4QwL4SwHLgP6FPzgBDC5BDC19UPpwKtMltmAz38sHe5nHpq7EpEpBj07+/bWUZqpacT6NsCH9R4XFH9uboMAP6+IUVlzKhR0KaNVlUUkdxo0QKOOgruvhuWLMn56dMJ9NrG+YVaDzQ7GkgBI+p4fqCZlZtZeWVlZfpVNsTMmfD88/4WqFGj7J5LRGS1U0+Fr7+OspF0OoFeAbSu8bgV8NGaB5nZAcCFQO8QQq1Lj4UQxoUQUiGEVFlZWUPqTd/o0b4B9AknZPc8IiI1dekCe+3lGVRVldNTpxPorwLtzaytmZUCfYFJNQ8wsy7AWDzM46/2/sUXcOed0LcvtGwZuxoRKTaDB8M778DTT+f0tPUGeghhJTAYeAKYDUwIIcwys8vNrHf1YSOAZsBfzex1M5tUx8vlxl13wVdfaaiiiMRxxBFQVgY335zT01oItXaHZ10qlQrl5eWZf+EQoHNnKC2F8nJN9ReROC64AIYPh/ffh1aZG/hnZtNCCKnankveTNGXXvIbooMGKcxFJJ6BA72BOX58zk6ZvEC/+WZft6Vfv9iViEgxa9vWJzWOH++THHMgWYFeWenrthx7LGyySexqRKTYDRrkkxsfeSQnp0tWoN9+u690dtJJsSsREYFevWC77XJ2czQ5gV5V5eu27Lcf7LJL7GpERHxS48CBPnzxnXeyfrrkBPpTT8G8ef4WR0QkXwwYAI0b52Tzi+QE+tixPonosMNiVyIi8p0f/hD69PGlAJYty+qpkhHoCxfCpEm+0lnTprGrERH5vpNOgk8/hYceyuppkhHot90Gq1bB738fuxIRkbV17w7t2mW926XwA72qysd5duvm+/qJiOSbkhJvcD7/PPz739k7TdZeOVeefNKn1mqooojks/79/eZoFmeOFn6gjxvni+AcemjsSkRE6rb11p5TWbw5WtiBXvNmaGlp7GpERNbtpJPgs89g4sSsvHxhB/rtt/vN0BNPjF2JiEj9unWD3XaDTz7Jyss3zsqr5kJVFdxyC+y/v26GikhhKCmB117zP7Px8ll51VyYPBnmz1frXEQKS5bCHAo50MeP9x22Dz88diUiInmhMAN98WKfcXXMMbDRRrGrERHJC2kFupn1NLM5ZjbXzM6v5fmmZnZ/9fMvm1mbTBf6PXfe6cvkqrtFROS/6g10M2sEjAIOAjoC/cys4xqHDQA+DyHsCFwHXJXpQv9r9aLDzugAAAZASURBVJZOP/0p7Lpr1k4jIlJo0mmhdwXmhhDmhRCWA/cBfdY4pg/w5+q/PwB0N8vShp5TpsDs2Vq3RURkDekE+rbABzUeV1R/rtZjQggrgaXAlmu+kJkNNLNyMyuvrKxsWMUAPXrAb3/b8H8vIpJA6QR6bS3t0IBjCCGMCyGkQgipsrKydOpb2957wz/+Ac2aNezfi4gkVDqBXgG0rvG4FfBRXceYWWOgOfBZJgoUEZH0pBPorwLtzaytmZUCfYFJaxwzCTiu+u9HAs+GENZqoYuISPbUO/U/hLDSzAYDTwCNgNtCCLPM7HKgPIQwCbgVuNPM5uIt877ZLFpERNaW1louIYTHgcfX+NwlNf6+DPh1ZksTEZH1UZgzRUVEZC0KdBGRhFCgi4gkhAJdRCQhLNboQjOrBN5v4D9vCSzOYDmFohivuxivGYrzuovxmmH9r3v7EEKtMzOjBfqGMLPyEEIqdh25VozXXYzXDMV53cV4zZDZ61aXi4hIQijQRUQSolADfVzsAiIpxusuxmuG4rzuYrxmyOB1F2QfuoiIrK1QW+giIrIGBbqISELkdaDn3ebUOZLGdZ9pZm+Z2Qwze8bMto9RZybVd801jjvSzIKZJWJ4WzrXbWa/qf56zzKze3JdY6al8f29nZlNNrPp1d/jvWLUmUlmdpuZLTKzmXU8b2Z2Y/X/yQwz271BJwoh5OUHvlTvu0A7oBR4A+i4xjGnAGOq/94XuD923Tm67v2Bjav/PqjQrzuda64+blPgBWAqkIpdd46+1u2B6UCL6sdbxa47B9c8DhhU/feOwHux687Ade8H7A7MrOP5XsDf8d3ffga83JDz5HMLPb82p86deq87hDA5hPB19cOp+C5ShSydrzXAUGA4sCyXxWVROtf9e2BUCOFzgBDCohzXmGnpXHMANqv+e3PW3iGt4IQQXmDdu7j1Af4S3FRgczPbZn3Pk8+BnrHNqQtMOtdd0wD8N3shq/eazawL0DqE8GguC8uydL7WHYAOZvaSmU01s545qy470rnmS4GjzawC34fhtNyUFtX6/tzXKq0NLiLJ2ObUBSbtazKzo4EU8IusVpR967xmMysBrgOOz1VBOZLO17ox3u3yS/yd2Itm1imEsCTLtWVLOtfcD7gjhHCNme2F74bWKYRQlf3yoslIluVzC71YN6dO57oxswOAC4HeIYRvc1RbttR3zZsCnYDnzOw9vI9xUgJujKb7Pf63EMKKEMJ8YA4e8IUqnWseAEwACCFMATbCF7BKsrR+7uuTz4FerJtT13vd1d0PY/EwL/Q+VajnmkMIS0MILUMIbUIIbfD7Br1DCOVxys2YdL7HH8ZvgmNmLfEumHk5rTKz0rnmBUB3ADPbGQ/0ypxWmXuTgGOrR7v8DFgaQli43q8S++5vPXeGewFv43fFL6z+3OX4DzP4F/qvwFzgFaBd7JpzdN1PA58Ar1d/TIpdc7aveY1jnyMBo1zS/FobcC3wFvAm0Dd2zTm45o7AS/gImNeBA2PXnIFrvhdYCKzAW+MDgJOBk2t8nUdV/5+82dDvb039FxFJiHzuchERkfWgQBcRSQgFuohIQijQRUQSQoEuIpIQCnQRkYRQoIuIJIQCXYqemXWsXmN9Vh3PNzezKjNLwqxcSTAFuogvcAYwrY7nd8dn8r2Wm3JEGkaBLvJdoNe1Nswe1X8q0CWvKdBF6m+hK9ClIGgtFylqZtYI+AJoCmwWvtsJquYxb+NL1rYLvoStSF5SC12KXUdgY2B2HWG+GbAj8LnCXPKdAl2KXbo3RKfnphyRhlOgS7FbHeiv1/F8r+o/1X8ueU+BLsVu9Q3Pb9Z8wsyaA0dVP6yrBS+SNxToUrSq96HdrfphPzNrWuO5bYD7+W7n9Tk5Lk9kvSnQpZh1wrcxrAC6Au+b2eNmNgXfCmxLfMswgFvM7Ng4ZYqkR4EuxWx1//kzwGH4no/dgO2AW4EDgNF4d0wJvsGxSN7SOHQpWmZ2M75R7+AQwqjY9YhsKLXQpZjVN+VfpKCohS5FycyaAF/ijZrNQgjLIpckssHUQpditSs+3X+mwlySQi10EZGEUAtdRCQhFOgiIgmhQBcRSQgFuohIQijQRUQSQoEuIpIQCnQRkYT4f7eWkXnbIyzZAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "prior = Beta(mu,2,2)\n", "\n", "plot(mus, (sym.lambdify(mu,prior))(mus), 'r')\n", "xlabel(\"$\\mu$\", size=22)\n", "\n", "maths(\"\\mathbf{Prior}: \\operatorname{Beta}(\\mu|a=2,b=2) = %s\" % sym.latex(prior))" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\mathbf{Likelihood}: \\operatorname{Bin}(m=1|N=1,\\mu) = \\mu$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEQCAYAAACgBo8fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAe9UlEQVR4nO3deXRV5bnH8e/j1C4qIAW8VVHBK94SJsWoOFuHFtAL7a11gZc6oTjhUK1XbdVarFpxwAnUMBTFCRTFqCAgg4ACEuapKCAIRSUqIooCgef+8UZNQ0JOwsnZ5+z9+6zFIuec3eTZgj+fvvsdzN0REZHct1vUBYiISHoo0EVEYkKBLiISEwp0EZGYUKCLiMTEHlH94EaNGnnTpk2j+vEiIjlp1qxZn7p744o+iyzQmzZtSlFRUVQ/XkQkJ5nZqso+05CLiEhMKNBFRGJCgS4iEhNVBrqZDTazdWa2sJLPzcweNrNlZjbfzNqlv0wREalKKh36EKDDTj7vCDQv/dUTeGzXyxIRkeqqMtDdfTLw+U4u6QI85cF0YB8z2y9dBYqISGrSMYZ+ALC6zOs1pe/twMx6mlmRmRUVFxen4UeLiMh30hHoVsF7Fe7J6+4F7p7v7vmNG1c4L15EJLY2bYIbb4RVlc4k3zXpCPQ1wIFlXjcB1qbh+4qIxMbEidC6NfTpA6NG1c7PSEegFwLnlc52aQ9scPeP0vB9RURy3oYN0LMnnHoq7LYbTJoEl19eOz+ryqX/ZvYccArQyMzWAH8B9gRw98eBUUAnYBmwCbiwdkoVEckthYUhvD/+GG64AW6/HerUqb2fV2Wgu3u3Kj534Mq0VSQikuPWrYOrr4Zhw8IwyyuvQH5+7f9crRQVEUkTd3j6aWjRAl56CXr3hqKizIQ5RLjboohInKxeDZddFh54tm8PgwZBXl5ma1CHLiKyC7Zvh8ceg5YtwwPPvn1h6tTMhzmoQxcRqbH334eLL4bJk+H006GgAJo1i64edegiItVUUhLmk7dpA/PmheGVsWOjDXNQhy4iUi3z5kGPHjBrFvz619CvH+y/f9RVBerQRURSsHkz3HprmLGyejUMHx5msmRLmIM6dBGRKk2bFrryJUvgvPPggQegYcOoq9qROnQRkUp89RVccw0cfzx8/TWMHg1PPpmdYQ7q0EVEKjRuXNiDZeVK6NUL7roL6taNuqqdU4cuIlLG+vVw0UXwy1/CXnvBlCnwyCPZH+agQBcR+d7IkWFB0FNPwc03hxktJ5wQdVWp05CLiCTeJ5/AVVfBCy/A4YfD669Duxw87l4duogklnvoxlu0CFvd3nUXvPtuboY5qEMXkYRatQouvRTGjIHjjgurPX/+86ir2jXq0EUkUbZvD6s7W7UKm2g98kh48JnrYQ7q0EUkQZYuDZtpTZ0aZrEUFMDBB0ddVfqoQxeR2Nu6Fe6+G9q2hUWLYMgQeOONeIU5qEMXkZibMycs258zB377W3j0UfjZz6KuqnaoQxeRWPr2W/jTn+Coo2DtWhgxAl58Mb5hDurQRSSGpk4NY+VLl8IFF4TNtBo0iLqq2qcOXURiY+PGsEDopJPCdrdjx8I//pGMMAcFuojExJgxYSpiv34h1BcsgDPOiLqqzFKgi0hO+/zzMKzSoQPUqROGWx56CPbeO+rKMk+BLiI5a8SIsJnW00/Dn/8cZrIcd1zUVUVHD0VFJOd8/DFceWU4Aq5duzCn/PDDo64qeurQRSRnuIdFQS1ahB0R//53mDFDYf4ddegikhNWrgwnCI0bByeeCAMHwmGHRV1VdlGHLiJZbds2ePjhMINl2rQwi2XSJIV5RdShi0jWWrIkLBB65x3o2BEefxwOOijqqrKXOnQRyTpbt8Kdd4ax8X/+E4YODWPmCvOdSynQzayDmS01s2VmdlMFnx9kZhPNbI6ZzTezTukvVUSSYNassP/KLbfAr38duvTu3cEs6sqyX5WBbma7A/2AjkAe0M3M8spddgsw3N2PALoC/dNdqIjE2zffwE03wTHHwLp18PLLMGwY7Ltv1JXljlQ69KOBZe6+wt23AM8DXcpd40C90q/rA2vTV6KIxN3kyWGv8nvugQsvhMWLQ3cu1ZNKoB8ArC7zek3pe2XdDnQ3szXAKOCqir6RmfU0syIzKyouLq5BuSISJ19+CVdcASefDCUl8OabMGAA7LNP1JXlplQCvaKRKy/3uhswxN2bAJ2AoWa2w/d29wJ3z3f3/MaNG1e/WhGJjdGjw1TExx+HP/whbKZ12mlRV5XbUgn0NcCBZV43YcchlR7AcAB3nwb8GGiUjgJFJF4++wx+/3vo1Anq1g1TEh94AH7yk6gry32pBPpMoLmZNTOzvQgPPQvLXfMhcBqAmbUgBLrGVETke+4wfHhYtv/883DbbTB7NrRvH3Vl8VHlwiJ3LzGzXsAYYHdgsLsvMrPeQJG7FwLXAwPM7A+E4ZgL3L38sIyIJNTatWEzrZEjIT8/jJW3aRN1VfGT0kpRdx9FeNhZ9r3byny9GDg+vaWJSK5zh8GD4frrwwlC994L114Le2iNeq3QP1YRqRUrVsAll8CECWEWy8CBcOihUVcVb1r6LyJptW0bPPggtG4NM2eGWSwTJijMM0EduoikzaJF0KNH2KP8zDNDmDdpEnVVyaEOXUR22ZYt0Ls3HHEELF8Ozz4Lr76qMM80degisktmzgxd+YIF0LVr2Ltc6wajoQ5dRGpk0ya44YYwj/zzz6GwEJ57TmEeJXXoIlJtkyaFgyeWLw8zWe69F+rXj7oqUYcuIinbsAEuvRR+8YvwesIEKChQmGcLBbqIpOS116BlyzCf/PrrYf78H4JdsoMCXUR2qrgYzj0X/vu/oUGDcFDzffdBnTpRVyblKdBFpELu4SFnXh68+CLcfns4Hu7oo6OuTCqjh6IisoM1a8LBE6++GgJ80KCwd7lkN3XoIvK97dvDQ86WLcOOiPffH/YrV5jnBnXoIgLAsmVhCuKkSeFh54AB8J//GXVVUh3q0EUSrqQkPORs3TocOFFQAOPHK8xzkTp0kQRbsCAs2585Ezp3hv794YDyR8BLzlCHLpJAmzfDX/4C7drBypXhSLiRIxXmuU4dukjCzJgBF10EixdD9+5h7/KGDaOuStJBHbpIQnz9NVx3HRx7LHz5Jbz+OgwdqjCPE3XoIgkwfnyYwfLBB2F++d13Q716UVcl6aYOXSTGvvgi7Ip4+unhYOa33oJ+/RTmcaVAF4mpV14Jy/aHDIEbb4R58+Ckk6KuSmqThlxEYmbdOrjqKhg+HNq0Ccv3jzwy6qokE9Shi8SEOzzzTOjKR46Ev/0NiooU5kmiDl0kBlavhssug1GjwpFwgwaFYJdkUYcuksO2b4fHHgubaU2aBA89BFOnKsyTSh26SI56770wg2XKlDCLpaAAmjWLuiqJkjp0kRxTUgJ9+kDbtmEvlsGDYexYhbmoQxfJKfPmhWX7s2fDb34T5pTvt1/UVUm2UIcukgM2b4Zbb4X8fPjXv8KRcC+9pDCXf6cOXSTLvfNO2OL2n/+E886Dvn3hpz+NuirJRil16GbWwcyWmtkyM7upkmvOMbPFZrbIzJ5Nb5kiyfPVV3DNNXDCCbBpE7zxBjz5pMJcKldlh25muwP9gDOANcBMMyt098VlrmkO3Awc7+7rzWzf2ipYJAnGjoWePeHDD3/YTKtu3airkmyXSod+NLDM3Ve4+xbgeaBLuWsuAfq5+3oAd1+X3jJFkmH9+vDQ81e/gh//GCZPhkcfVZhLalIJ9AOA1WVeryl9r6zDgMPM7G0zm25mHSr6RmbW08yKzKyouLi4ZhWLxNRLL4UFQU89BTffDHPnhuEWkVSlEuhWwXte7vUeQHPgFKAbMNDM9tnhf+Re4O757p7fuHHj6tYqEksffwxnnw2//S387GfhfM+77godukh1pBLoa4ADy7xuAqyt4JpX3H2ru38ALCUEvIhUwj085MzLg9degzvvhHffhSOOiLoyyVWpBPpMoLmZNTOzvYCuQGG5a0YCvwAws0aEIZgV6SxUJE5WrYJOneCCC6BFizC88qc/wZ57Rl2Z5LIqA93dS4BewBhgCTDc3ReZWW8z61x62RjgMzNbDEwEbnD3z2qraJFctX17WN3ZqlXYg+Xhh8PvP/951JVJHJh7+eHwzMjPz/eioqJIfrZIFJYuDZtpTZ0aZrE88QQcfHDUVUmuMbNZ7p5f0Wda+i9Sy7ZuDfPI27aFRYvCkXCjRyvMJf209F+kFs2ZE5btz5kTZrI88kiYySJSG9Shi9SCb78NDzmPOgo++ghGjIAXXlCYS+1Shy6SZlOnhrHypUvhwgvh/vuhQYOoq5IkUIcukiYbN0KvXnDiiaFDHzs2HD6hMJdMUaCLpMGYMWEqYv/+YYfEhQvhjDOirkqSRoEusgs++wzOPx86dIA6dcJwy4MPwt57R12ZJJECXaQG3MOpQXl58OyzcMstYbXnccdFXZkkmR6KilTTRx/BlVfCyy9Du3ZhrLxt26irElGHLpIyd/jHP0JXPno03HMPzJihMJfsoQ5dJAUrV4YThMaNC7NYBg6Eww6LuiqRf6cOXWQntm0LG2i1agXTpoVZLJMmKcwlO6lDF6nEkiVh2f60adCxIzz+OBx0UNRViVROHbpIOVu3hsMmDj88rPYcOhRef11hLtlPHbpIGbNmhUOa58+Hc84Jm2ntu2/UVYmkRh26CPDNN3DjjXDMMVBcHKYkDhumMJfcog5dEu+tt+CSS+D998OY+X33wT47HHEukv3UoUtiffklXH45nHIKlJTAm2+G6YgKc8lVCnRJpFGjoGXLcAzctdfCggVw2mlRVyWyaxTokiiffgrdu8OZZ0K9emFKYt++8JOfRF2ZyK5ToEsiuMPzz0OLFuFh5223wezZ4SGoSFzooajE3r/+BVdcAYWFkJ8fDp1o3TrqqkTSTx26xJY7DBgQNtMaOxb69AlDLApziSt16BJLK1aEqYgTJsDJJ4fZK4ceGnVVIrVLHbrEyrZt4SFnq1ZQVBT2X5kwQWEuyaAOXWJj4UK4+OKwR/lZZ8Fjj0GTJlFXJZI56tAl523ZAn/9azg9aPnycCRcYaHCXJJHHbrktJkzw2ZaCxfCueeGA5obN466KpFoqEOXnLRpE9xwA7RvD+vXh478mWcU5pJs6tAl50ycGGawLF8Ol14azvasXz/qqkSipw5dcsaGDSHATz01vJ44McxiUZiLBCkFupl1MLOlZrbMzG7ayXVnm5mbWX76ShSBV18NC4QGDoQ//jEcQHHKKVFXJZJdqgx0M9sd6Ad0BPKAbmaWV8F1dYGrgRnpLlKSq7g4POzs3Bl++tOw0vPee6FOnagrE8k+qXToRwPL3H2Fu28Bnge6VHDdHUAf4Ns01icJ5R6mH7ZoAS++CL17h+Phjj466spEslcqgX4AsLrM6zWl733PzI4ADnT313b2jcysp5kVmVlRcXFxtYuVZFizJnTk//u/YYXnnDlw662w115RVyaS3VIJdKvgPf/+Q7PdgL7A9VV9I3cvcPd8d89vrPllUs727VBQEA6emDAhLOF/++3wWkSqlsq0xTXAgWVeNwHWlnldF2gFTDIzgJ8BhWbW2d2L0lWoxNuyZWEq4qRJYRbLgAFwyCFRVyWSW1Lp0GcCzc2smZntBXQFCr/70N03uHsjd2/q7k2B6YDCXFJSUhIOZW7dOgytDBwYzvZUmItUX5UduruXmFkvYAywOzDY3ReZWW+gyN0Ld/4dRCq2YAH06BGW73fpAv37w/77R12VSO5KaaWou48CRpV777ZKrj1l18uSONu8Ge66K/xq0CAcCfe734FV9LRGRFKmpf+SUdOnh6588eJwWPODD0LDhlFXJRIPWvovGfH113DddXDccbBxI7z+OgwdqjAXSSd16FLrxo8PM1g++AAuvxz+/neoVy/qqkTiRx261JovvggnCJ1+OuyxB7z1VnjwqTAXqR0KdKkVr7wSNtMaMgRuvBHmzYOTToq6KpF405CLpNUnn8DVV8Pw4dC2bdgl8cgjo65KJBnUoUtauIeHnHl5MHIk3HFHmF+uMBfJHHXosss+/BAuuwxGj4Zjj4VBg8IuiSKSWerQpca2bw8POVu2DA88H3wQpkxRmItERR261Mh774UZLFOmhFksBQXQrFnUVYkkmzp0qZaSknAoc5s2YS+WwYNh7FiFuUg2UIcuKZs3Dy66CGbPht/8Bvr1g/32i7oqEfmOOnSp0rffwi23QH5+OE3ohRfgpZcU5iLZRh267NQ774Sx8iVL4Pzz4YEHwmHNIpJ91KFLhb76Cq65Bk44IWys9cYbYdWnwlwke6lDlx2MGwc9e8KqVXDllWHf8rp1o65KRKqiDl2+t349XHgh/PKX8KMfweTJ8MgjCnORXKFAFyA85MzLC8v3b74Z5s4Nwy0ikjs05JJwH38MvXrBiBFw+OHh4Il27aKuSkRqQh16QrnDk0+Grvy11+Duu+HddxXmIrlMHXoCrVoFl14KY8bA8ceHzbT+67+irkpEdpU69ATZvh0efTRspjV1avh68mSFuUhcqENPiKVLoUcPePvtMIuloAAOPjjqqkQkndShx9zWrWF8vG1bWLw4jJu/8YbCXCSO1KHH2Jw5YTOtuXPh7LPDEMt//EfUVYlIbVGHHkPffBPmkh91VJiWOGJE2FBLYS4Sb+rQY2bq1DBW/t57YdXn/fdDgwZRVyUimaAOPSY2bgwLhE48EbZsCYdODB6sMBdJEgV6DIweHaYi9u8PV18dThI644yoqxKRTFOg57DPPoPzzoNOnWDvvcOUxIceCl+LSPIo0HOQe3jImZcHzz0XThOaMweOPTbqykQkSikFupl1MLOlZrbMzG6q4PPrzGyxmc03s/FmplnOtWTtWvif/4FzzoEDD4SiIrjjjrDdrYgkW5WBbma7A/2AjkAe0M3M8spdNgfId/c2wItAn3QXmnTuYc+VvLywMOiee2D69LBgSEQEUuvQjwaWufsKd98CPA90KXuBu090902lL6cDTdJbZrKtWBEecl58cQjw+fPh//4P9tCkUxEpI5VAPwBYXeb1mtL3KtMDGF3RB2bW08yKzKyouLg49SoTats2ePBBaN0aZswIs1gmToTmzaOuTESyUSo9nlXwnld4oVl3IB84uaLP3b0AKADIz8+v8HtIsHhx6MinTYOOHeGJJ8KYuYhIZVLp0NcAZaOkCbC2/EVmdjrwZ6Czu29OT3nJs3Ur/O1vcMQRYbXn00+HU4QU5iJSlVQ69JlAczNrBvwL6AqcW/YCMzsCeALo4O7r0l5lQhQVhWX78+dD165hTvm++0ZdlYjkiio7dHcvAXoBY4AlwHB3X2Rmvc2sc+ll9wJ7Ay+Y2VwzK6y1imPom2/CQ85jjoFPP4WRI8P8coW5iFRHSvMk3H0UMKrce7eV+fr0NNeVGG+9FcbKly2DSy6BPn1gn32irkpEcpFWikbkyy/h8svhlFPC0XDjx4dThBTmIlJTCvQIjBoVNtMqKIDrrgubaZ16atRViUiuU6Bn0KefQvfucOaZUL8+vPNO2K+8Tp2oKxOROFCgZ4A7DBsWlu0PHw5/+QvMnh0egoqIpIsWj9eytWvDWHlhYTgSbtCgsPJTRCTd1KHXEncYODB05ePGwX33hVWfCnMRqS3q0GvB8uVhCuLEiWEWy4ABcOihUVclInGnDj2Ntm2Dvn1DFz5rVpjFMn68wlxEMkMdeposXBiW7b/7Lpx1Fjz2GDTRJsIikkHq0HfRli3w179Cu3Zh3/LnngsPQBXmIpJp6tB3wcyZcNFFoTs/99ywmVajRlFXJSJJpQ69BjZtgj/+Edq3h/Xr4dVX4ZlnFOYiEi116NU0cWLYTGvFCrj00nC2Z/36UVclIqIOPWUbNoQAP/VU2G23EOyPP64wF5HsoUBPwauvhgVCAweGoZZ588L8chGRbKJA34niYujWDTp3hoYNw0HN996rzbREJDsp0CvgDs8+Cy1awIgR0Lt3OB4uPz/qykREKqeHouWsXh0203r99bAb4qBBYe9yEZFspw691Pbt4SFny5bhgecDD8DbbyvMRSR3qEMH3n8/bKb11ltw2mlhD5ZDDom6KhGR6kl0h15SEra1bdMG5s4NwyvjxinMRSQ3JbZDnz8/bKZVVARdukD//rD//lFXJSJSc4nr0DdvDkfAHXkkrFoVjoZ7+WWFuYjkvkR16NOnh6588WL4/e/D3uUNG0ZdlYhIeiSiQ//6a7j2WjjuONi4MUxJfOophbmIxEvsO/Q33wwzWFauhCuugLvvhnr1oq5KRCT9Ytuhf/FFGF454wzYc88wJbFfP4W5iMRXLAN95MiwmdaTT8JNN4XNtE46KeqqRERqV6yGXD75BK66Cl54Adq2DbskHnlk1FWJiGRGLDp0dxg6NHTlr7wCd94ZjodTmItIkuR8h/7hh+HgiTfeCLNYBg4MuySKiCRNSh26mXUws6VmtszMbqrg8x+Z2bDSz2eYWdN0F1re9u3hIWfLljBlCjz8cPhdYS4iSVVloJvZ7kA/oCOQB3Qzs7xyl/UA1rv7oUBf4J50F1rW0qVw8snQqxcceywsXBjGzneLxQCSiEjNpBKBRwPL3H2Fu28Bnge6lLumC/Bk6dcvAqeZmaWvzB8MHhweeC5aBEOGwJgx0LRpbfwkEZHckkqgHwCsLvN6Tel7FV7j7iXABmCHdZhm1tPMisysqLi4uEYFH3YYnHVWWL5//vlQO//ZEBHJPak8FK0oMr0G1+DuBUABQH5+/g6fp+KEE8IvERH5d6l06GuAA8u8bgKsrewaM9sDqA98no4CRUQkNakE+kyguZk1M7O9gK5AYblrCoHzS78+G5jg7jXqwEVEpGaqHHJx9xIz6wWMAXYHBrv7IjPrDRS5eyEwCBhqZssInXnX2ixaRER2lNLCIncfBYwq995tZb7+FvhdeksTEZHq0MxtEZGYUKCLiMSEAl1EJCYU6CIiMWFRzS40s2JgVQ3/542AT9NYTq5I4n0n8Z4hmfedxHuG6t/3we7euKIPIgv0XWFmRe6eH3UdmZbE+07iPUMy7zuJ9wzpvW8NuYiIxIQCXUQkJnI10AuiLiAiSbzvJN4zJPO+k3jPkMb7zskxdBER2VGudugiIlKOAl1EJCayOtCz8XDq2pbCPV9nZovNbL6ZjTezg6OoM92quu8y151tZm5mOT+9LZV7NrNzSv+8F5nZs5musTak8Hf8IDObaGZzSv+ed4qiznQys8Fmts7MFlbyuZnZw6X/TOabWbsa/SB3z8pfhK16lwOHAHsB84C8ctdcATxe+nVXYFjUdWfgnn8B1Cn9+vJcv+dU77v0urrAZGA6kB913Rn4s24OzAEalL7eN+q6M3TfBcDlpV/nASujrjsN930S0A5YWMnnnYDRhNPf2gMzavJzsrlDz6rDqTOkynt294nuvqn05XTCCVK5LpU/a4A7gD7At5ksrpakcs+XAP3cfT2Au6/LcI21IZX7dqBe6df12fGEtJzj7pPZ+SluXYCnPJgO7GNm+1X352RzoKftcOockso9l9WD8F/1XFflfZvZEcCB7v5aJgurRan8WR8GHGZmb5vZdDPrkLHqak8q93070N3M1hDOYbgqM6VFqrr/7lcopQMuIpK2w6lzSMr3Y2bdgXzg5FqtKDN2et9mthvQF7ggUwVlQCp/1nsQhl1OIfw/sSlm1srdv6jl2mpTKvfdDRji7veb2bGE09Baufv22i8vMmnJsmzu0JN4OHUq94yZnQ78Gejs7pszVFttquq+6wKtgElmtpIwxliY4w9GU/37/Yq7b3X3D4ClhIDPZancdw9gOIC7TwN+TNjAKs5S+ne/Ktkc6Ek8nLrKey4deniCEOZxGFOFKu7b3Te4eyN3b+ruTQnPDjq7e1E05aZFKn+/RxIegmNmjQhDMCsyWmX6pXLfHwKnAZhZC0KgF2e0yswrBM4rne3SHtjg7h9V+7tE/fS3iifDnYD3CE/F/1z6Xm/Cv8wQ/qBfAJYB7wKHRF1zBu75TeATYG7pr8Koa87EfZe7dhI5PsslxT9rAx4AFgMLgK5R15yh+84D3ibMgJkL/DLqmtNwz88BHwFbCd14D+Ay4LIyf9b9Sv+ZLKjp328t/RcRiYlsHnIREZFqUKCLiMSEAl1EJCYU6CIiMaFAFxGJCQW6iEhMKNBFRGJCgS6JZ2Z5pXusL6rk8/pmtt3M4rIyV2JKgS4SNjkDmFXJ5+0IK/lmZ6YckZpRoIv8EOiV7Q1zZOnvCnTJagp0kao7dAW65ATt5SKJZma7A18CPwLq+Q+nQZW95j3CtrWHeNjGViQrqUOXpMsD6gBLKgnzesChwHqFuWQ7BbokXaoPROdkphyRmlOgS9J9F+hzK/m8U+nvGj+XrKdAl6T77oHnN+U/MLP6wLmlLyvr4EWyhgJdEqv0HNq2pS+7mdmPyny2HzCMH05eX5rh8kSqTYEuSdaKcIzhGuBoYJWZjTKzaYSjwBoSjgwDGGhm50VTpkhqFOiSZN+Nn48HfkM48/FU4CBgEHA60J8wHLMb4YBjkayleeiSWGb2GOGg3l7u3i/qekR2lTp0SbKqlvyL5BR16JJIZrYnsJHQ1NRz928jLklkl6lDl6RqTVjuv1BhLnGhDl1EJCbUoYuIxIQCXUQkJhToIiIxoUAXEYkJBbqISEwo0EVEYkKBLiISE/8PYzZsP1itUpoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "likelihood = Bin(1,1,mu)\n", "\n", "plot(mus, (sym.lambdify(mu,likelihood))(mus), 'b')\n", "xlabel(\"$\\mu$\", size=22)\n", "\n", "maths(\"\\mathbf{Likelihood}: \\operatorname{Bin}(m=1|N=1,\\mu) = %s\" % sym.latex(likelihood))" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\mathbf{Posterior (unnormalised)}: \\operatorname{Beta}(\\mu|2,2) \\times \\operatorname{Bin}(1|1,\\mu) = 6 \\mu^{2} \\left(1 - \\mu\\right)$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEQCAYAAACgBo8fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3dfZzVY/7H8denUkKF6sdSVKtl21qL+blZq6SoxERLSmm1VhvFbmgRa8nuIuuulci9isrN0lJyU6x7DRFFu+lG/YSEQrqbuX5/fCZmx2TO1DnnOud73s/HYx4zZ87XnM/XTO/5zvW9rs9lIQRERCT/1YpdgIiIpIcCXUQkIRToIiIJoUAXEUkIBbqISELUifXCTZo0CS1atIj18iIieem11177JITQtKrnogV6ixYtKCkpifXyIiJ5ycyWbO45DbmIiCSEAl1EJCEU6CIiCaFAFxFJCAW6iEhCKNBFRBJCgS4ikhDR5qGLSAasWgULFsDy5fDBB/Dpp7Bxo7/VqgU77wyNG8Muu8Dee8Nuu4FZ7KolTRToIvlq40aYNQuefhpefRXmzIElm11zUrWGDaFtW2jfHjp2hEMPhe23z0y9knEWa4OLoqKioJWiIjX01Vfw2GMweTI88QR88YV//sc/hn339bd99vEr7x/8wK/G69aF2rWhtNSv2Feu9Cv4d9+FefNg9mz/hbBxox979NFw8slwzDFQv37c85XvMLPXQghFVT6nQBfJcSHACy/AmDHw8MOwZg3suisUF0OnTn5l3bTK1h6p+/JLf41p0/yXxfLl0KABDBgAv/89tGyZnnORraZAF8lH69bB3XfDjTfCW29Bo0bQu7e/HXaYX3VnQmkpPPss3HknTJwIZWXQsydccgm0a5eZ15SUfV+ga5aLSK5ZswZuuAFatYLf/taD+9Zb4f/+D26+GQ4/PHNhDv61jzgCxo2DxYth2DB48kn42c/gtNO8DslJCnSRXLFxI9xyiwf5738PrVvDU0/B66/Db34T52bl7rvDlVfCwoVe0/jxXteVV3q9klMU6CK5YNo0v6E5aJAH5r/+Bc8842PkuTCtcOed4Zpr/EZqt25w4YVw0EHw5puxK5MKFOgiMS1dCscf7zNL1q+Hhx7yMD/ssNiVVa1lS3jwQXjgAVi2DIqK4C9/8XF2iU6BLhJDaSlcfz20aQPTp/sQxty5Hu65cEVenV/+0qc8/vKXcPHF0L07fPJJ7KoKngJdJNvee89vbA4d6lfic+fC+ef7HPB80rgx3Hcf3HQTzJgB++8Pr7wSu6qCpkAXyZYQYOxYHyufM8enJD72WH7P8TaDM86AF1/02TGHH+7DRhKFAl0kGz7/HE44wachHnIIvP029O+fH8MrqTjgAF9t+rOf+Xled53/ApOsUqCLZFpJiQ9HTJkCV1/tY+bNm8euKv2aNvWhl5494Zxz/E2hnlUKdJFMuuUW+PnP/Sboc8/Beed518Okql/fWwf87nd+03fwYM2AySJ1WxTJhPXr4ayzfMy8WzdfkLPzzrGryo5atXzIZdtt4aqrYMMG/8WW5F9kOUKBLpJuH33kww4vvugLcC6/PLNL9XORGVxxBWyzDfz5zz70cuutyblnkKMU6CLpNG+ez8n+6COYNAl69YpdUTxm/sts0/vGjf2KXTImpb+BzKyrmc03swVmdkEVz+9hZjPNbLaZzTGzo9NfqkiOmzHDx8u//tpXexZymFd02WVw5pkwciT87W+xq0m0agPdzGoDo4FuQBugj5m1qXTYxcDkEMJ+QG/gpnQXKpLTJkyArl2hWTNfXFNUZXfTwmQGo0b5L7hhw3z+vWREKlfoBwILQggLQwjrgYlAj0rHBKBh+ceNgA/SV6JIjrvhBujXz7dve/552HPP2BXlntq14Z57oHNn7xz57LOxK0qkVAJ9d2BphcfLyj9X0aVAPzNbBkwFzqrqC5nZQDMrMbOSFStWbEG5IjkkBLjoIm8r27Ond0zcccfYVeWuevW8qddee3kPmEWLYleUOKkEelW3pSuvFugD3BVCaAYcDYwzs+987RDC2BBCUQihqOnWbpklElNZmc+x/utf4fTTfe71ttvGrir3NWrkC6zKyuDYY2H16tgVJUoqgb4MqLisrRnfHVI5DZgMEEJ4CdgWaJKOAkVyTmmpDxuMGQN/+IPPsS60aYlbo3VruP9+761+yilaeJRGqQT6LKC1mbU0s7r4Tc8plY55H+gEYGY/xgNdYyqSPBs3eg+WO++EP/3J295qbnXNderki4+mTNHMlzSqNtBDCBuBIcB04B18NstcMxthZsXlh50LnG5mbwL3AaeGWLtPi2TKxo3Qty/ce68vmrn0UoX51hgyBE48EYYP95vJstUsVu4WFRWFkpKSKK8tUmOlpT48cN99fkV57rmxK0qGVau8U+PatfDGG9BEI7XVMbPXQghVzotVcwWR6pSWwoABHuZXXqkwT6dGjXw8/ZNPfChL4+lbRYEu8n3KymDgQBg3zpevn39+7IqSZ7/94NprfdrnmDGxq8lrCnSRzQnBr8bvuMP3zbz44tgVJdcZZ/hK22HDYP782NXkLQW6yOZcdpn39D77bBgxInY1yWYGt9/u/dT79/cb0FJjCnSRqlx/vQf6gAE+vU6zWTJvt918yOXVV30WkdSYAl2ksgkTYOhQX55+663amCGbevWCPn38L6I5c2JXk3f0kypS0fTpcOqpvnv9+PFaARrD3/8OO+3kq3FLS2NXk1cU6CKbzJrlV+U/+Qk8/LB6s8TSuLEPec2aBaNHx64mryjQRQAWLvSdhpo29elzjRrFrqiw9enjs16GD4f3349dTd5QoIusXOkbOZeWwuOPww9+ELsiMfMbpCH4bkfqJJISBboUtrVr4bjjYPFiH2bZe+/YFckmLVr4BtOPPebfG6mWAl0KV1mZT0t8/nnfTeeww2JXJJWddZbf0zjnHN+rVb6XAl0K12WXwcSJPuf5pJNiVyNVqVPHZ70sXqw2uylQoEthmjDB5zoPGKD+LLmuY0dvs3vFFbBkSexqcpoCXQrPiy/Cr38N7dvDzTdrFWg+uPpqfz9sWNw6cpwCXQrL++/7TdDmzeGhh6Bu3dgVSSr23BMuuMBb7T77bOxqcpYCXQrHV19Bjx6wbh08+qgvYJH8MWwYNGvm7zWNsUoKdCkMIfh4+Ztv+o3QffaJXZHUVP363pN+1iy/UpfvUKBLYfjLXzwERo70RUSSn045Bdq1gwsvhPXrY1eTcxToknyPPgp//CP066ft4/Jd7dr+S3nhQrjlltjV5BwFuiTb/PnQt69vRDx2rGa0JEGXLtCpk087XbUqdjU5RYEuybV6tc9oqVvXZ7TUrx+7IkkHM79K/+QT33xEvqFAl2QqK/O+5v/5D0yeDHvsEbsiSaf994eePT3QP/00djU5Q4EuyTRyJPzjH74gpWPH2NVIJlx2GXzxBVxzTexKcoYCXZLn6afhoou8P8vvfx+7GsmUtm19y7obbvDhF1GgS8IsXQq9e/s889tu003QpPvTn7wL48iRsSvJCQp0SY5167yJ07p1fhN0hx1iVySZ9uMfw8knw403wocfxq4mOgW6JMd558Err8Cdd2qjikJyySX+S/zaa2NXEp0CXZJh0iS/Shs61Dd6lsLRurXfLxkzpuBnvCjQJf/Nnw+/+Q0ccghcdVXsaiSGCy+EL7/0X+oFTIEu+W3NGh83r1fPr9K32SZ2RRJDu3ZQXOwzXr78MnY10SjQJb+dfTa8/bbvQNS8eexqJKYLL/QhlwLu8aJAl/w1fjzcfjsMH+79PaSwHXwwHHGELzRauzZ2NVEo0CU/vfsuDBoEhx0Gl14auxrJFcOHw/LlcM89sSuJQoEu+efrr32FYP36cN99vjO8CPgV+gEH+BTGsrLY1WRdSoFuZl3NbL6ZLTCzCzZzTC8zm2dmc83s3vSWKVLB0KHw1lswbhzsvnvsaiSXmME55/jMp6lTY1eTddUGupnVBkYD3YA2QB8za1PpmNbAhcChIYSfAGqgIZlx//1+0+v886Fr19jVSC468UTfe7QAFxqlcoV+ILAghLAwhLAemAj0qHTM6cDoEMJnACGEj9NbpgiwaJHPNz/oIN9bUqQq22zjs59mzoTZs2NXk1WpBPruwNIKj5eVf66iHwE/MrMXzOxlM6vy0snMBppZiZmVrFixYssqlsK0YYM33TLzcXPNN5fvc/rp3sunwFrrphLoVbWrC5Ue1wFaA4cDfYDbzGzH7/xHIYwNIRSFEIqaNm1a01qlkF18Mbz6qndQbNkydjWS63bc0f+amzQJli2LXU3WpBLoy4CKKzaaAR9UccwjIYQNIYRFwHw84EW23pNPenvUgQPhhBNiVyP54ne/85kuBdQOIJVAnwW0NrOWZlYX6A1MqXTMw0BHADNrgg/BLExnoVKgVqyA/v29Tar2j5SaaNHC95S97baCWWhUbaCHEDYCQ4DpwDvA5BDCXDMbYWbF5YdNB1aa2TxgJjAshLAyU0VLgQjB9wX97DOYOBG22y52RZJvBg+GlSt96KUAWAiVh8Ozo6ioKJSUlER5bckTo0b5n82jRsFZZ8WuRvJRCPCTn8D228OsWbGrSQszey2EUFTVc1opKrlpzhwYNgy6d4chQ2JXI/nKzH9+Skr8pnrCKdAl93z9NfTpAzvv7LsPaV9Q2RqnnAINGhTEzVEFuuSeYcNg3jy4+27Q9FbZWg0awK9+5ePoHyd7zaMCXXLLo4/C6NHer+Woo2JXI0kxeDCsX+8zXhJMgS6546OP4Ne/hn33hSuuiF2NJMk++3gnxltvTXQXRgW65IYQYMAA+OILuPde31JOJJ0GDoTFi+Gpp2JXkjEKdMkNo0fDtGlw9dXQpk31x4vU1HHHQZMmMHZs7EoyRoEu8c2b5zdCu3XzsU6RTKhXz2+OPvIIfPhh7GoyQoEuca1bB337+kyEO+7QFEXJrNNPh40b4a67YleSEQp0ieuSS+CNN3yz5113jV2NJN3ee0OHDj7bJYE3RxXoEs8zz/iY+cCBcOyxsauRQjFwILz3nm+AkTAKdInj88+9i+JeexXkVmESUc+evgr51ltjV5J2CnSJY/Bg+OADGD/eGyeJZMu22/p9m4cf9k6eCaJAl+ybONHnml9yCRx4YOxqpBCdeqrfkJ84MXYlaaVAl+xatgzOOAMOPhiGD49djRSq/faDdu0SN9tFgS7ZU1bmV0br18O4cVCnTuyKpFCZ+c/iq6/6OoiEUKBL9vz97/D0076V3F57xa5GCl3fvlC7tnf1TAgFumTHvHlw/vlwzDG+uEMktl128Q1Uxo3zxUYJoECXzFu//ttNBm67TatBJXeceiosXw5PPhm7krRQoEvmjRgBr7/u83532SV2NSLf6t4dGjf2nbESQIEumfXii97bfMAA73Ynkkvq1oXeveGf/4TVq2NXs9UU6JI5X37pq0H32AOuvz52NSJV69cP1q6Fhx6KXclWU6BL5px3Hixc6LMIGjaMXY1I1Q46CH74Q1+1nOcU6JIZjz0Gt9ziod6+fexqRDbPDE4+GWbM8HYUeUyBLun3ySdw2mm+Eu/yy2NXI1K9vn19G8Q8bwWgQJf0CgEGDfKmR+PHa29QyQ977w0HHAATJsSuZKso0CW9xo+HBx/0K/Of/jR2NSKp69fPp9e++27sSraYAl3S5/33YcgQOOwwOPfc2NWI1Ezv3lCrVl5fpSvQJT02Nd4qK/NZLbVrx65IpGZ23RU6dfLWziHErmaLKNAlPW64wbf0uv56aNkydjUiW6Z3b59q+9prsSvZIgp02Xrz5sGFF/q+oL/+dexqRLbcccd5W+fJk2NXskUU6LJ11q/3m0kNG3qvFjXekny2885w1FEe6Hk47KJAl61z+eUwezaMHavGW5IMvXrBkiW++UWeUaDLlnv5ZfjrX/1mqBpvSVIcd5w37Zo0KXYlNaZAly3z1VfeeKt5c78hKpIUjRpB165w//0+ayuPKNBlywwbBgsWqPGWJFOvXr6h+Usvxa6kRlIKdDPrambzzWyBmV3wPcedYGbBzIrSV6LknGnTYMwYXzzUoUPsakTSr7jY21bk2WyXagPdzGoDo4FuQBugj5m1qeK4BsDZwCvpLlJyyMqVPjWxbVs13pLkatAAjj4674ZdUrlCPxBYEEJYGEJYD0wEelRx3OXASGBtGuuTXLKp8dbKld6zZdttY1ckkjknnOD7jb6SP9eoqQT67sDSCo+XlX/uG2a2H9A8hPDo930hMxtoZiVmVrJixYoaFyuRTZgADzzge4Tuu2/sakQyq3t32GabvNrJKJVAr2qlyDcz7s2sFnAdUG03phDC2BBCUQihqGnTpqlXKfEtWQKDB8MvfuE3REWSrlEjOPJI7x6aJ4uMUgn0ZUDzCo+bARW39WgAtAWeMbPFwMHAFN0YTZCyMvjVr/z9Pfeo8ZYUjp49YdEiePPN2JWkJJVAnwW0NrOWZlYX6A1M2fRkCGFVCKFJCKFFCKEF8DJQHEIoyUjFkn3XXQfPPgujRqnxlhSW4mJvqZsnwy7VBnoIYSMwBJgOvANMDiHMNbMRZlac6QIlsrfeguHD4fjjfUWoSCFp2tSn5j74YOxKUmIh0thQUVFRKCnRRXxOW7sWDjwQPv7Yg133PaQQ3XgjnHUWvPMO7LNP7Gows9dCCFUOaWulqGzeRRd5kN9xh8JcCtemPkX/+EfcOlKgQJeqzZgB114LZ57pCyxEClWzZnDQQXkx7KJAl+/67DOf1bL33nD11bGrEYmvZ0/fxWjp0uqPjUiBLv8tBDjjDPjwQ18Nut12sSsSia+4fP7HP/8Zt45qKNDlv02Y4H2gL7sMirSUQATwm6E/+hFMmVL9sREp0OVbixd/uxr0/PNjVyOSW4qL/d7S6tWxK9ksBbq40lI45RT/eNw4rQYVqay4GDZsgOnTY1eyWQp0cVdeCc8/D6NHQ4sWsasRyT0//zk0bpzTwy4KdPH2oH/6E/TuDX37xq5GJDfVrg3HHAOPPeZX6jlIgV7ovvjCQ7xZM9+FyKpqrikigA+7fPYZvPBC7EqqpEAvdGef7d3kxo+HHXeMXY1IbjvqKN+aLkeHXRTohWzyZLjrLl/i/4tfxK5GJPftsAN06gSPPJKTPdIV6IVqyRIYONCXNP/xj7GrEckfxx4LCxfCu+/GruQ7FOiFaONGHzcvK4N77/VttkQkNZt6Gz32WNw6qqBAL0R//rPf1Ln5ZmjVKnY1Ivlljz2gXTsFuuSA55+Hyy+H/v3h5JNjVyOSn7p3939Lq1bFruS/KNALyaefeoi3bOlN+0Vkyxx9tA9dPvFE7Er+iwK9UIQAp53mXRQnToQGDWJXJJK/DjkEdtoJpk6NXcl/qRO7AMmSMWPg4Yfhb39TF0WRrVWnDnTp4oFeVuYbSeeA3KhCMmvOHDjnHOjWDYYOjV2NSDJ07+777b72WuxKvqFAT7ovv4STTvI/D++6K2euJETyXpcu3iojh2a76F930g0eDPPn+8YV//M/sasRSY6mTX1hXg6NoyvQk+zuu+Gee3wl6BFHxK5GJHm6d4dZs+Cjj2JXAijQk2vePDjzTDj8cLjkktjViCRTt27+/skn49ZRToGeRGvWQK9evsHzhAnafUgkU/bbz4deHn88diWApi0m0+DBfoX++OOw226xqxFJrlq1vKXuE0/kxPRFXaEnzZ13+myWiy/2HzQRyawuXWDFCnjjjdiVKNAT5a23/Oq8Y0ffUk5EMm/ThVMODLso0JNi9Wo44QRo1Mhb4mrcXCQ7dtnFx9KnT49diQI9ETb1aXnvPe/TsuuusSsSKSxdusCLL/qFVUQK9CS4/np44AG48kro0CF2NSKFp2tX7744Y0bUMhTo+e6552DYMOjZE849N3Y1IoXpkEN8v9HIwy4K9Hy2fLnPN2/VCu64w/tKiEj21a3rm0c//njUzaMV6Plq/Xq/Cbp6NTz0kN8MFZF4unSBxYvhP/+JVoICPV8NHeo3Ye68E9q2jV2NiGyavvjUU9FKSCnQzayrmc03swVmdkEVz59jZvPMbI6ZPW1me6a/VPnGXXfBTTfBeef5kIuIxNeqFbRokduBbma1gdFAN6AN0MfM2lQ6bDZQFEL4KfAAMDLdhUq5WbNg0CBfPHTFFbGrEZFNzKBzZ5/psnFjlBJSuUI/EFgQQlgYQlgPTAR6VDwghDAzhLCm/OHLQLP0limA3wQ97jifZz55sm+DJSK548gjYdWqaLsYpRLouwNLKzxeVv65zTkNmFbVE2Y20MxKzKxkxYoVqVcpsG4d/PKX8Pnn8Mgj0KRJ7IpEpLJN+w5EGnZJJdCrmgtX5bwcM+sHFAFXV/V8CGFsCKEohFDUtGnT1KssdCHAkCHw0ks+fr7vvrErEpGqNGnibQAi9UdPJdCXAc0rPG4GfFD5IDPrDFwEFIcQ1qWnPAFg1Ci47Ta46CI48cTY1YjI9+nc2WegffVV1l86lUCfBbQ2s5ZmVhfoDUypeICZ7Qfcgof5x+kvs4BNmwbnnAPHHw8jRsSuRkSq07kzbNjgq7izrNpADyFsBIYA04F3gMkhhLlmNsLMissPuxrYAbjfzN4wsymb+XJSE3PnwkknwU9/CuPGRW+eLyIpOOwwqFcvyrBLStMkQghTgamVPndJhY87p7kuWbECjj0Wtt8epkzx9yKS++rXh0MPjXJjVJd8uejrr6G42KcpPvIING9e/X8jIrmjc2eYMwc++iirL6tAzzVlZdC/P7zyCowfDwceGLsiEampzuWDFllup6tAzzXDh3tv85Ejfd65iOSf/faDhg1h5sysvqwCPZfcdBNcdZUv7Vdvc5H8VaeObzajQC9Q//iHLx469lj4+9/V21wk33XsCAsWwNKl1R+bJgr0XPDCC3DyyT5ePnGierSIJMGmNgBZvEpXoMc2b55flTdvDo8+CtttF7siEUmHdu2gcWMFesFYvNi7s9Wr51tXqeGWSHLUqgWHH+4zXbK0LZ0CPZYPP/QwX7MGnnjCm+OLSLJ07Ajvvw+LFmXl5RToMXz+OXTtCh98AFOn+p9mIpI8WR5HV6Bn2+rVHubz5vnMlkMOiV2RiGTKPvv4hjRZWmCkQM+mr76C7t19N5P77/92U1kRSSYzH0efOTMr4+gK9GzZ1J/lxRdhwgTo0aP6/0ZE8t8RR3hfpn//O+MvpUDPhjVrfGrizJm+41CvXrErEpFs6djR32dhHF2BnmmbhllmzPAwP+WU2BWJSDb98Iew227w7LMZfyktScykL7/0MH/+ed+gom/f2BWJSLaZQfv28K9/+Th6Btt66Ao9Uz791FtovvCCj5krzEUKV4cOPk35vfcy+jIK9ExYvty/gbNneyvc3r1jVyQiMXXo4O8zPOyiQE+3RYt8T8FFi3zR0HHHxa5IRGLbZx9o2tSHXTJIY+jpVFICxxwD69fD00/DQQfFrkhEcsGmcXRdoeeJqVN9AcG22/pcc4W5iFTUoQMsWeJvGaJAT4dbbvFFQ3vvDS+/7H9eiYhU1L69v8/gsIsCfWts2OC7DA0aBF26+J9Tu+4auyoRyUXt2sFOO2V02EWBvqU+/RS6dYPRo+G882DKFNhhh9hViUiuqlXLJ0wo0HPM669DURE895yv/rz6aqhdO3ZVIpLr2rf3fUY/+CAjX16BXhMhwNix8POfw8aN/pv2V7+KXZWI5ItN89EzNI6uQE/VqlXQvz/89rc+m+X11+Hgg2NXJSL55Gc/84vA3XfPyJfXPPRUPP889OsHy5bBiBEwfLiGWESk5urU8WHaDNEV+vdZu9bDu0MHD/DnnoM//lFhLiI5SVfom/Pcc3D66TB/PgwYADfcAA0axK5KRGSzdIVe2Sef+Dh5+/awbh1Mnw533KEwF5Gcp0DfZMMGGDUKWreG22+HoUPh7be176eI5A0FelkZTJrkq7h+9zv43/+FOXPg2mth++1jVycikrLCDfTSUu9Vvu++3q+8dm145BEfYmnTJnZ1IiI1VniBvmoVXHcd7LUXnHiiD7Xcd59flRcXZ3R7KBGRTCqMWS6lpb7j9rhx8OCDvnHzL34B11wDPXpoGqKIJEJKgW5mXYEbgNrAbSGEKys9Xw+4BzgAWAmcFEJYnN5Sa+iLL+Cpp2DaNHj0Ud8WrmFDH14ZNMh7sYiIJEi1gW5mtYHRwJHAMmCWmU0JIcyrcNhpwGchhL3MrDdwFXBSJgr+jtJSX8G5eDH8+9++a9CsWfDWW95vpUEDOPJIOOkkOPZYqF8/K2WJiGRbKlfoBwILQggLAcxsItADqBjoPYBLyz9+ALjRzCyEENJYq7v9dhg50odNvvrKr8RLS799fscd/er7D3/wID/0UNhmm7SXISKSa1IJ9N2BpRUeLwMq76/2zTEhhI1mtgpoDHxS8SAzGwgMBNhjjz22rOKmTb3BzQ47+LTChg1hzz2hZUto1QpatPC+wyIiBSaVQK9q2kflK+9UjiGEMBYYC1BUVLRlV+/Fxf4mIiL/JZVL2WVA8wqPmwGVu7N/c4yZ1QEaAZ+mo0AREUlNKoE+C2htZi3NrC7QG5hS6ZgpwKadHk4AZmRk/FxERDar2iGX8jHxIcB0fNriHSGEuWY2AigJIUwBbgfGmdkC/Mq8dyaLFhGR70ppHnoIYSowtdLnLqnw8VrgxPSWJiIiNaHpICIiCaFAFxFJCAW6iEhCKNBFRBLCYs0uNLMVwJIt/M+bUGkVaoEoxPMuxHOGwjzvQjxnqPl57xlCaFrVE9ECfWuYWUkIoeDaJRbieRfiOUNhnnchnjOk97w15CIikhAKdBGRhMjXQB8bu4BICvG8C/GcoTDPuxDPGdJ43nk5hi4iIt+Vr1foIiJSiQJdRCQhcjrQzayrmc03swVmdkEVz9czs0nlz79iZi2yX2V6pXDO55jZPDObY2ZPm9meMepMt+rOu8JxJ5hZMLO8n96WyjmbWa/y7/dcM7s32zVmQgo/43uY2Uwzm13+c350jDrTyczuMLOPzeztzTxvZjaq/P/JHDPbf4teKISQk294q973gFZAXeBNoE2lY84Ebi7/uDcwKXbdWTjnjsB25Va93WQAAAQ1SURBVB+fke/nnOp5lx/XAPgX8DJQFLvuLHyvWwOzgZ3KH/9P7LqzdN5jgTPKP24DLI5ddxrOuz2wP/D2Zp4/GpiG7/52MPDKlrxOLl+hf7M5dQhhPbBpc+qKegB3l3/8ANDJzKraDi9fVHvOIYSZIYQ15Q9fxneQynepfK8BLgdGAmuzWVyGpHLOpwOjQwifAYQQPs5yjZmQynkHoGH5x4347g5peSeE8C++fxe3HsA9wb0M7GhmP6jp6+RyoFe1OfXumzsmhLAR2LQ5db5K5ZwrOg3/rZ7vqj1vM9sPaB5CeDSbhWVQKt/rHwE/MrMXzOxlM+uateoyJ5XzvhToZ2bL8H0YzspOaVHV9N9+lVLa4CKStG1OnUdSPh8z6wcUAR0yWlF2fO95m1kt4Drg1GwVlAWpfK/r4MMuh+N/iT1nZm1DCJ9nuLZMSuW8+wB3hRCuMbND8N3Q2oYQyjJfXjRpybJcvkIvxM2pUzlnzKwzcBFQHEJYl6XaMqm6824AtAWeMbPF+BjjlDy/MZrqz/cjIYQNIYRFwHw84PNZKud9GjAZIITwErAt3sAqyVL6t1+dXA70QtycutpzLh96uAUP8ySMqUI15x1CWBVCaBJCaBFCaIHfOygOIZTEKTctUvn5fhi/CY6ZNcGHYBZmtcr0S+W83wc6AZjZj/FAX5HVKrNvCtC/fLbLwcCqEMLyGn+V2Hd/q7kzfDTwb/yu+EXlnxuB/2MG/0bfDywAXgVaxa45C+f8FPAR8Eb525TYNWfjvCsd+wx5Psslxe+1AdcC84C3gN6xa87SebcBXsBnwLwBHBW75jSc833AcmADfjV+GjAIGFThez26/P/JW1v6862l/yIiCZHLQy4iIlIDCnQRkYRQoIuIJIQCXUQkIRToIiIJoUAXEUkIBbqISEIo0KXgmVmb8h7rczfzfCMzKzOzpKzMlYRSoIt4kzOA1zbz/P74Sr7Xs1OOyJZRoIt8G+ib6w1zQPl7BbrkNAW6SPVX6Ap0yQvq5SIFzcxqA6uBekDD8O1uUBWP+TfetrZV8Da2IjlJV+hS6NoA2wHvbCbMGwJ7AZ8pzCXXKdCl0KV6Q3R2dsoR2XIKdCl0mwL9jc08f3T5e42fS85ToEuh23TD8+vKT5hZI+Dk8oebu4IXyRkKdClY5fvQ7lv+sI+Z1avw3A+ASXy78/r8LJcnUmMKdClkbfFtDJcBBwJLzGyqmb2EbwXWGN8yDOA2M+sfp0yR1CjQpZBtGj9/Gjge3/PxCGAP4HagM3ATPhxTC9/gWCRnaR66FCwzG4Nv1DskhDA6dj0iW0tX6FLIqlvyL5JXdIUuBcnMtgG+wC9qGoYQ1kYuSWSr6QpdClU7fLn/2wpzSQpdoYuIJISu0EVEEkKBLiKSEAp0EZGEUKCLiCSEAl1EJCEU6CIiCaFAFxFJiP8HIx/CWyMZNuMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "posterior = prior * likelihood\n", "\n", "plot(mus, (sym.lambdify(mu,posterior))(mus), 'r')\n", "xlabel(\"$\\mu$\", size=22)\n", "\n", "maths(r\"\\mathbf{Posterior (unnormalised)}: \\operatorname{Beta}(\\mu|2,2) \\times \\operatorname{Bin}(1|1,\\mu) = %s\" % sym.latex(posterior))" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\mathbf{Posterior}: p(\\mu|m,l,a,b) = p(\\mu|1,0,2,2) = \\operatorname{Beta}(\\mu,2+1,2+0) = 12 \\mu^{2} \\left(1 - \\mu\\right)$" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEQCAYAAAC+z7+sAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3debyWc/7H8denUogMyt6GUAi5J/xIslREIUtI+BkJYZjNLsz8xjJkyxLSyCi7slSWrFNxTqNSGUmWjpoplQxKqs/vj+/duNU5nbtz7vv+3sv7+Xicxzn3dV33uT9Xp/M+1/29vou5OyIiUrzqxC5ARESyS0EvIlLkFPQiIkVOQS8iUuQU9CIiRa5e7AIq07hxY2/RokXsMkRECsakSZO+cvcmle3Ly6Bv0aIF5eXlscsQESkYZvZ5VfvUdCMiUuQU9CIiRU5BLyJS5BT0IiJFTkEvIlLkFPQiIkVOQS8iUuTysh+9iGTYypXw2WcwZw7Mmxc+li2DFSvCvoYNYcstw0eLFrDrrrDRRrGrlgypNujNbAhwNDDf3feoZP/vgNNSvl9roIm7LzKzz4D/ACuBFe6eyFThIrIO//oXjBsHb74JkyfDtGnw/ffpP98MdtwREgno1AkOPRR23jlsl4Jj1S08YmYHA98Cj1QW9Gscewxwibsfmnz8GZBw96/Wp6hEIuEaGSuynqZPh8cfh2efDcEOsNlm0K4d7LUX7LknNG8O220H22wTruLr1Qvh/d13sHAhfPUVzJ4NM2aE7zF+PMydG77XTjvBqaeGj912i3eeUikzm1TVxXS1V/Tu/paZtUjztU4BhqdfmojUyuLFMHQoPPRQCPo6deDgg+HGG+Gww2CffaBu3eq/zyabhI/mzWHffX/a7g4ffwyvvQbPPAN/+hPccAO0bw+XXgo9e4Y/FpLXqr2iB0gG/QvruqI3s42BCmBnd1+U3PYpsBhw4H53H7yO5/cF+gI0a9Zs388/r3LaBhGZNg0GDoThw2HpUth/fzjtNDjhhHC1ni3z5oV3DffcE/4ANG8Ov/kNnHsu1K+fvdeVaq3rij6TvW6OAf6+OuSTDnT3dsCRwAXJZqBKuftgd0+4e6JJk0onYBORf/wDjj8+NMOMGAGnnw7vvw8TJkD//tkNeYBtt4Vf/xr++U947jlo2hQuughat4YnngjvACTvZDLoe7FGs427z01+ng88C7TP4OuJlI6PPw7NJPvuG26yXn01fPEF3H8/7L137uupUwd69IC334YxY0J7/8knw4EHhiYkySsZCXoz2wzoCIxM2dbQzDZd/TXQGZiWidcTKRkLF8LFF0ObNjB2LFx3HXz+OVx/fegKmQ+6dAnvKoYMCX+Q9tkn1Ld8eezKJKnaoDez4cAEYFczqzCzs82sn5n1SznsOOBld/8uZdvWwDtmNgV4D3jR3cdksniRouUOf/1r6N1y991w9tkwaxZcc03oSZNv6taFs84KvXV69oRrrw03bD/6KHZlQpo3Y3NN3SulpM2cGW5uvvEGHHAA3HcftG0bu6r1M3Jk+OP0ww/w4IOhWUeyKlc3Y0WkNlatgjvvDH3eJ08O7e/vvFN4IQ+h/f7990PtvXqFG8U//hi7qpKloBfJB3PmQOfOoT3+0ENDE0jfvuGmZ6Fq2jS8K7n0Uhg0CI46CpYsiV1VSSrg/0UiRWLUqHAV/+67MHgwvPBC6MZYDDbYAG69NdyofeMN6NAh/FGTnFLQi8SyfHkYbNSjB7RsGZprzjmnOOeTOessGD069Bjabz91wcwxBb1IDHPnQseOcNttcOGFYU6ZnXaKXVV2HX44/P3v4etOnWDq1Lj1lBAFvUiuTZgQBj598EEYTXrnndCgQeyqcmOPPcKMmvXrh3sRkyfHrqgkKOhFcmnIEDjkENh4Y5g4EU48MXZFudeqVQj7hg1D2E+ZEruioqegF8mFVavgsstC3/KOHaGsLFzdlqqddvop7Lt0CVMjS9Yo6EWybdkyOOUUuOkm6NcPXnoJttgidlXxtWgBL78c+td37hwWS5GsUNCLZNPixeEm5BNPwM03h+l9NX/7T1q3Dn/45s2Drl3Vzz5LFPQi2fLll2ERkLKyEPS/+11xdp2srf32C4uaTJ8eRtGuXBm7oqKjoBfJhpkzw5S9n30W+o+X4k3X9dGlSxg9O2YM/OEPsaspOnoPKZJpU6eG5hoIo0FTl+aTqvXtG1bOuvVW2H33MMhKMkJX9CKZNGlSGAzUoEFYlEMhv35uuy38kTz33DDeQDJCQS+SKRMmhH7hjRrBW2/BrrvGrqjw1KsX7mc0awYnnQRffRW7oqKgoBfJhPHjQxfBrbYKId+yZeyKCtfmm4ewnz8f+vQJYxCkVhT0IrX17ruha+B224VBQE2bxq6o8LVrBwMHhhvZt9wSu5qCp6AXqY2ysp+u5MeNC2EvmXHeeWFlqiuvDAuwSI0p6EVqasqUEPJbbBFCfvvtY1dUXMzC/PwtWsDpp8N//hO7ooKVzuLgQ8xsvplNq2L/IWa2xMwmJz+uSdnX1cw+MrNZZnZZJgsXiWrmzBDym2wCr78ebh5K5jVqFBZJ/+ILuOSS2NUUrHSu6IcCXas55m133zv5cT2AmdUFBgFHAm2AU8ysTW2KFckLX3wRugC6w6uvhitOyZ4DD4Tf/x4eeiisxiXrrdqgd/e3gEU1+N7tgVnuPtvdlwMjgB41+D4i+WP+/BDy33wTJuRSF8rcuO66sNziOefAggWxqyk4mWqjP8DMppjZaDPbPblteyB1cciK5LZKmVlfMys3s/IF+kFKPvrPf+DII6GiIkzEtffesSsqHfXrw7Bh8PXX0L9/7GoKTiaC/h9Ac3ffC7gLeC65vbLZm7yqb+Lug9094e6JJk2aZKAskQxavhyOPz7cgH3ySfif/4ldUenZc0+46qrQx/6FF2JXU1BqHfTu/o27f5v8+iVgAzNrTLiCT+1QvAMwt7avJ5Jzq1bBGWeE9viHHoJu3WJXVLr+8IcwD87556sXznqoddCb2TZmYe5VM2uf/J4LgTKglZm1NLP6QC9Ad1Kk8Pz+9zBiBNx4Ywh8iad+fXjggdB8dtVVsaspGNXOXmlmw4FDgMZmVgFcC2wA4O73AScA55nZCmAp0MvdHVhhZv2BsUBdYIi7T8/KWYhky6BBYTbFCy4IgS/xHXBAuKK/6y447TRo3z52RXnPQibnl0Qi4eXl5bHLkFL3/PNw7LGhqebZZ6Fu3dgVyWrffANt2oQRyWVl+tkAZjbJ3ROV7dPIWJHKTJoUVjtq1w6GD1eQ5JtGjcI7rfffhwcfjF1N3lPQi6ypogKOOQaaNAlX9Q0bxq5IKnPSSdCxY5gLZ1FNhvqUDgW9SKpvvw0h/+23oQvfNtvErkiqYgZ33hkWYL/mmuqPL2EKepHVVq2C3r3DUoAjRsAee8SuSKrTtm2Y5fLee8PPTSqloBdZ7corYeTIMA/6UUfFrkbSdf31YbGSiy8O8w/JWhT0IgCPPRb6yfftCxdeGLsaWR9bbAEDBoSF2EePjl1NXlL3SpHycujQIfTHfuWVMChHCsvy5WHEbIMGYZqKEuwlpe6VIlX5179CX/mtt4annlLIF6r69eHPf4bp08P89fIzCnopXcuXQ8+eodfGyJGhO6UUrp49Yb/9Qg+c77+PXU1eUdBL6br4Yhg/HoYODXOdS2EzCwuJf/kl3H577GryioJeStODD8J994XZEE88MXY1kikdOkD37nDzzeGdmgAKeilF774bJinr3Bn+9KfY1Uim3XADLFkSuskKoKCXUjN/PpxwAmy3XehSWYK9M4pe27bhXdrtt8PChbGryQsKeikdK1bAKafAV1/BM8/AllvGrkiy5dprwzQWf/lL7ErygoJeSsfVV8O4cWG4/D77xK5Gsmn33cPso3feGd7FlTgFvZSGkSN/Gvl65pmxq5FcuPZaWLYMbropdiXRKeil+H3yCfTpA4kE3HFH7GokV3bdNaxAde+9sGBB7GqiUtBLcVu2LNyYq1sXnnwSNtwwdkWSS5dfHv4PlPgfeAW9FLeLLw6rED3yCLRoEbsaybXWreH44+Huu0OXyxJVbdCb2RAzm29m06rYf5qZTU1+jDezvVL2fWZmH5jZZDPTLGWSW48+CoMHh0FRRx8duxqJ5YorQsjfe2/sSqJJ54p+KNB1Hfs/BTq6e1vgBmDwGvs7ufveVc2qJpIVH34I554LBx8Mf/xj7GokpnbtoGtXuO22kp0Dp9qgd/e3gCoXZHT38e6+eqzxRGCHDNUmUjPffx/WE23YMCzsXa9e7IoktiuuCDdkH3oodiVRZLqN/mwgdeZ/B142s0lm1nddTzSzvmZWbmblC0r8DrnU0sUXw7RpMGxYGAEr0qEDHHRQmPTsxx9jV5NzGQt6M+tECPo/pGw+0N3bAUcCF5jZwVU9390Hu3vC3RNNNF2s1NTf/hYmLLviCujSJXY1kk9+/3uYMyesO1BiMhL0ZtYWeBDo4e7/nVzC3ecmP88HngXaZ+L1RCo1c2Zol+/QAa67LnY1km+6dQt962+9teTWlq110JtZM+AZ4HR3n5myvaGZbbr6a6AzUGnPHZFa++GHMOS9QYMwWZna5WVNderAJZfApEnw9tuxq8mpdLpXDgcmALuaWYWZnW1m/cysX/KQa4AtgXvW6Ea5NfCOmU0B3gNedPcxWTgHkfC2/P33wyIiO6g/gFTh9NPDZHa33hq7kpyq9rLH3U+pZv+vgF9Vsn02oGV7JPtGjQqTV118MRxzTOxqJJ9tvDGcf37ocjtzJuyyS+yKckIjY6WwVVTAWWeFvtKavErSccEFYTHxElpuUEEvhWvlSujdO7TPDx8e2udFqrP11uH/zdChsKjKIUJFRUEvhevPf4Y334RBg0rmLbhkyEUXwdKl8PDDsSvJCQW9FKbx42HAgLBiVJ8+sauRQtO2beiGe8894Z1hkVPQS+FZsgROPRWaNg0TVZnFrkgKUf/+MHs2jCn+zoAKeiks7nDeeeEm7GOPwWabxa5ICtVxx8G224amvyKnoJfC8uij4cbrgAFwwAGxq5FCtsEG0K8fjB4Ns2bFriarFPRSOD75JHSN69AhrBwkUlvnnBNGUd9zT+xKskpBL4Xhxx/D+p916oSr+rp1Y1ckxWDbbeGEE2DIEPjuu9jVZI2CXgrDH/8I774L998PzZrFrkaKyfnnhxv8Tz4Zu5KsUdBL/hs/PgR9nz5w8smxq5Fic9BBsNtuYdnJIqWgl/z2zTehyaZ5c7jrrtjVSDEyg759YcIE+OCD2NVkhYJe8tuFF8IXX4TVoho1il2NFKvTTw/z3zzwQOxKskJBL/nriSfgkUfgyivhwANjVyPFrHFj6NkzXFAsXRq7moxT0Et+qqgIfZzbt4err45djZSCvn3h66+LcqlBBb3kn1Wr4Mwzw6yUjz4aBraIZFvHjtCqVVHelFXQS/654w547bUwX3irVrGrkVJhFgZQvfMOfPhh7GoySkEv+eWDD+Cyy6B7d/jVWguXiWRXnz5hMN5f/xq7koxS0Ev++OGHsCDEL34Rej9oVkrJta23hqOOCp0AVqyIXU3GpBX0ZjbEzOab2bQq9puZ3Wlms8xsqpm1S9l3hpl9nPw4I1OFSxG6+mqYOjUMR99qq9jVSKk680yYNw9eeSV2JRmT7hX9UKDrOvYfCbRKfvQF7gUwsy2Aa4H9gPbAtWa2eU2LlSL25pvwl7/AuedCt26xq5FSdvTRsOWWYanBIpFW0Lv7W8C6FlfsATziwUTgF2a2LdAFeMXdF7n7YuAV1v0HQ0rRkiWhbXSnnULYi8RUv35Y2Oa552Dx4tjVZESm2ui3B+akPK5Ibqtq+1rMrK+ZlZtZ+YIFCzJUlhSEiy6CL78Mg1U22SR2NSJw1lmwfDmMGBG7kozIVNBXdtfM17F97Y3ug9094e6JJk2aZKgsyXtPP/3T6Nf9949djUiw995hXdkiab7JVNBXAE1THu8AzF3HdpFww+vccyGRgKuuil2NyE/Mwk3Z996DGTNiV1NrmQr6UUCfZO+b/YEl7j4PGAt0NrPNkzdhOye3Salzh7PPDos9DBum0a+Sf049NSx087e/xa6k1tLtXjkcmADsamYVZna2mfUzs37JQ14CZgOzgAeA8wHcfRFwA1CW/Lg+uU1K3eDBYa3OW24Jc4GL5Jutt4YjjgiL0K9aFbuaWjH3SpvMo0okEl5eXh67DMmWWbNgr73CjJRjxoSrJpF8NGxY6BH29tthgZI8ZmaT3D1R2T79hklurVjx09zfDz+skJf8duyxsNFGBd98o98yya2bboKJE+Gee2D7SnvaiuSPTTeFHj3C2gjLl8eupsYU9JI7//gHDBgAvXrBKafErkYkPaedBosWwdjC7UeioJfcWLYsNNlstRUMGhS7GpH0dekSpkQo4OYbBb3kxhVXhP7IDz8MW2wRuxqR9G2wAZx8MowcGRarL0AKesm+ceNg4EC44ALo3Dl2NSLr77TTwrvSkSNjV1IjCnrJrq+/DiMMd9kFbr45djUiNXPAAdCsGTz+eOxKakRBL9l10UUwd27oj7zxxrGrEakZMzjxRHj55YKc0VJBL9nz9NMh4K+8Etq3j12NSO2cfDL8+GOYvrjAKOglO/71L01YJsUlkYCWLUOf+gKjoJfMcw8Le2vCMikmZnDSSfDqq7BwYexq1ouCXjLvwQfhxRfDKFhNWCbF5OSTwzQezz4bu5L1oqCXzPrkE7jkEjjsMOjfP3Y1Ipm1996w884F1/tGQS+Zs3JlmOlvgw3CyjyasEyKjVm4qh83DgpoyVP9Jkrm3HwzjB8fpjjYYYfY1Yhkx0knhfnpn3kmdiVpU9BLZrz/PlxzTfgl0IRlUsz23BNatVLQS4lZtgx694YmTeDee8PbW5FiZQbHHx+abwpk8JSCXmrv8ss1YZmUlp49Q++b55+PXUla0l0ztquZfWRms8zsskr2DzSzycmPmWb2dcq+lSn7RmWyeMkDr70Gt98eeth06RK7GpHcSCTCfagCab6pV90BZlYXGAQcAVQAZWY2yt1nrD7G3S9JOf5CYJ+Ub7HU3ffOXMmSNxYvDhOW7bZb6DMvUipWN98MHgzffgubbBK7onVK54q+PTDL3We7+3JgBNBjHcefAgzPRHGS5y64IEx1oAnLpBT17BnuT40eHbuSaqUT9NsDc1IeVyS3rcXMmgMtgXEpmzc0s3Izm2hmx9a4Uskvjz0Gw4fDtdeGt7EipebAA0MHhAJovqm26QaorAuFV3FsL+Apd1+Zsq2Zu881sx2BcWb2gbt/staLmPUF+gI0a9YsjbIkmi++gPPPD3N0X7bWLRuR0lC3Lhx7bLjgWbYMNtwwdkVVSueKvgJomvJ4B2BuFcf2Yo1mG3efm/w8G3iDn7ffpx432N0T7p5o0qRJGmVJFKtWhXb5lStDk029dK4VRIrU8ceHNvpXX41dyTqlE/RlQCsza2lm9QlhvlbvGTPbFdgcmJCybXMza5D8ujFwIDBjzedKARk4EF5/He64A3baKXY1InEdeig0apT3c9RXeznm7ivMrD8wFqgLDHH36WZ2PVDu7qtD/xRghLunNuu0Bu43s1WEPyo3pvbWkQIzZUpY5PvYY+Gss2JXIxJf/fpw5JGhP/2qVXk7v5P9PJfzQyKR8PLy8thlSKqlS+GXvwzzcH/wATRuHLsikfwwfDiceipMmAD77x+tDDOb5O6V9ozIzz8/kn8uuwymTw+zUirkRX7StWu4VzVyZOxKqqSgl+qNHQt33hkW+tboV5Gf23xzOPhgGJW/A/8V9LJuCxaEXja77w433hi7GpH81KNHmO9p1qzYlVRKQS9Vc4ezz4ZFi8IAqY02il2RSH465pjwOU+v6hX0UrX77w+9CW66Cdq2jV2NSP5q2TLMU6+gl4Ly4Ydw6aWhTf6ii2JXI5L/evSAd94JPdPyjIJe1vbDD6G7WMOGWvtVJF3du4cR4y+9FLuSteg3WNZ2+eUweXJYSGSbbWJXI1IY9t0Xtt0WXnghdiVrUdDLz40ZE6Y56N8fjj46djUihaNOnTBKduxY+PHH2NX8jIJefvLvf8MZZ4SbSrfcErsakcLTrRssWQLjx8eu5GcU9BKsnpXym29gxIi8nnJVJG8dfjhssAG8+GLsSn5GQS/BwIGh2ea226BNm9jViBSmRo2gQwcFveShsrIwl81xx0G/frGrESls3bqFUbKffRa7kv9S0Je6b76BXr1Cb4EHHwyLHotIzXXrFj7nUTdLBX0pc4fzzoPPPw9TrW6xReyKRArfLruERXnyqPlGQV/KHn44zGEzYEBY6FhEas8MjjoKxo2D77+PXQ2goC9d06eHvvKHHhoGSIlI5nTrFhYMf+ON2JUACvrS9N13cNJJsOmm8Le/hdXsRSRzOnaEjTfOm+abateMlSJ04YVh0rKXX9YUByLZsOGG4d3y2LGxKwHSvKI3s65m9pGZzTKzyyrZf6aZLTCzycmPX6XsO8PMPk5+nJHJ4qUGHnkktM1ffnkY3CEi2dGlC3zySV4sRlLtFb2Z1QUGAUcAFUCZmY1y9xlrHPq4u/df47lbANcCCcCBScnnLs5I9bJ+pk8PvWw6doTrrotdjUhx69o1fB47FnbeOWop6VzRtwdmuftsd18OjAB6pPn9uwCvuPuiZLi/AnStWalSK99+CyeeCJtsEnra1FOrnUhW7bwz7LhjXjTfpBP02wNzUh5XJLetqaeZTTWzp8ys6Xo+FzPra2blZla+YMGCNMqStLnD+efDP/8ZQn677WJXJFIaunQJ3SyXL49aRjpBX9lQSV/j8fNAC3dvC7wK/HU9nhs2ug9294S7J5o0aZJGWZK2Bx6AYcNCf/nDDotdjUjp6No19HL7+9+jlpFO0FcATVMe7wDMTT3A3Re6+w/Jhw8A+6b7XMmy8vLQy6ZzZ7jyytjViJSWTp1CM+mYMVHLSCfoy4BWZtbSzOoDvYCfrYBrZtumPOwOfJj8eizQ2cw2N7PNgc7JbZILixbBCSeELpTqLy+Se5tuCgcdFL2dvtqgd/cVQH9CQH8IPOHu083sejPrnjzsIjObbmZTgIuAM5PPXQTcQPhjUQZcn9wm2bZqFfTuDfPmwVNPQePGsSsSKU1dusCUKeF3MRJzr7TJPKpEIuHl5eWxyyhs110X2uTvvVdTD4vENHky7LMPDB0aVnDLEjOb5O6JyvZpCoRi9MILIeT79IFzz41djUhp22sv2HrrqM03Cvpi8/HHocmmXTu47z7NLy8SmxkccQS8+mpoUo1AQV9Mvv02rBJVrx488wxstFHsikQEwnQjCxbAtGlRXl5BXyzcw+LeH34YFhFp3jx2RSKy2up5pV55JcrLK+iLxZ/+BE8/DTffHN4mikj+2H57aN06NN9EoKAvBqNGwdVXh7b5Sy+NXY2IVObww+Gtt+CHH6o/NsMU9IVuxowQ8IkEDB6sm68i+erww8PSghMm5PylFfSFbOFCOOaYsJLNs8/q5qtIPuvYMYxOj9B8o6AvVMuXQ8+e8OWX8NxzsMMOsSsSkXXZbDNo315BL2laPe3wm2/CkCGw//6xKxKRdBxxBJSVweLcrr2koC9EAwfCQw/BVVfBqafGrkZE0nX44WHQ1Btv5PRlFfSF5tln4be/DbNSajlAkcKy337QsGHOm28U9IXkvffgtNNCO98jj0Ad/fhECkr9+uGmrIJeKvXZZ6GHzTbbhH7z6mEjUpg6dYKZM2Fu7tZgUtAXgoULw5Jky5fDiy/CVlvFrkhEaurQQ8Pn11/P2Usq6PPd0qXQvXu4oh85MgyjFpHCtdde8ItfhEXDc6Rezl5J1t/KlaFXzYQJ8MQTcPDBsSsSkdqqWxcOOURX9ELoK3/BBWEw1B13hF42IlIcOnWCTz8N79RzQEGfr66+Gu6/H/7wB7jwwtjViEgm5bidPq2gN7OuZvaRmc0ys8sq2X+pmc0ws6lm9pqZNU/Zt9LMJic/RmWy+KJ1661h2uFzzoE//zl2NSKSabvvDk2a5Czoq22jN7O6wCDgCKACKDOzUe4+I+Ww94GEu39vZucBNwMnJ/ctdfe9M1x38RoyJAyIOvHEsLC3ZqMUKT5mofnm9ddDM22Wf8/TuaJvD8xy99nuvhwYAfRIPcDdX3f375MPJwKaYasmhg8PV/GdO8OwYeGmjYgUp06doKICZs3K+kulE/TbA3NSHlckt1XlbGB0yuMNzazczCaa2bFVPcnM+iaPK1+wYEEaZRWZZ56B00+Hgw4K0xw0aBC7IhHJpk6dwuccNN+kE/SVvafwSg806w0kgFtSNjdz9wRwKnC7me1U2XPdfbC7J9w90aRJkzTKKiIvvAC9eoWpDV54IcwvLyLFbZddYLvt8iboK4CmKY93ANYau2tmhwNXAt3d/b9rZbn73OTn2cAbwD61qLf4PP98mFe+bVsYPRo23TR2RSKSC2u202dROkFfBrQys5ZmVh/oBfys94yZ7QPcTwj5+SnbNzezBsmvGwMHAqk3cUvbyJE/hfwrr4SFCUSkdHTsCP/+d5j7JouqDXp3XwH0B8YCHwJPuPt0M7vezLonD7sF2AR4co1ulK2BcjObArwO3LhGb53S9cwzYRBUu3Yh5DffPHZFIpJrq0e7v/VWVl/GPMtvGWoikUh4eXl57DKyZ9gwOOus0CY/Zgw0ahS7IhGJwR223TYsSPLoo7X6VmY2KXk/dC0aGZtrd90FffqEuS5eflkhL1LKzELzzZtvZrWdXkGfK+5www1w0UVw7LGhd80mm8SuSkRiO/jg0J8+i/PeKOhzYcUKOO88uOaacDX/5JOw4YaxqxKRfNCxY/j85ptZewkFfbZ99x0cd1yYoOzyy2HoUKin2aFFJKlNG9hyy6wGvRInm+bNgx49YNKkMG9Nv36xKxKRfFOnDnTokNWeN7qiz5bycvjlL2HGjDCnvEJeRKrSsSPMnh3a6rNAQZ8Njz8e/kLXqwfjx4dFvUVEqpLldnoFfSb9+CP85jdh3ppEAt57L4x6FRFZl7Ztw8j4LAW92ugzZe5cOPlkeOcd6N8/LB5Sv37sqkSkEBhL3Y4AAAhCSURBVNStG2auzVI7vYI+E0aPhjPPDD1shg8PV/QiIuujd2+YOhVWrsz4WhRquqmNpUvDeq5HHQVbbQVlZQp5EamZXr3g//4vKwsOKehrqqwM9t0X7r4bfv3r8Lh169hViYisRUG/vr77Di69FPbfH5YsgbFjYeBAjXQVkbyloE+Xe5ifZs89Q7D37Rv6yHfuHLsyEZF1UtCn48MP4cgjQ3/4Bg1CF6h779VCISJSEBT06zJnThjRuueeMHEi3H57uCu+erEAEZECoO6VlZkzB266CR54IDTZnHsuDBgApbZouYgUBQV9qgkTwlX700+HBQH+93/hyiuhWbPYlYmI1JiCft48eOyxsLzflCmh3f2SS8Lo1ubNY1cnIlJrabXRm1lXM/vIzGaZ2WWV7G9gZo8n979rZi1S9l2e3P6RmXXJXOk15A7TpsHNN0OnTrDDDvDb34bpCgYNCrPH3XKLQl5Eika1V/RmVhcYBBwBVABlZjbK3WekHHY2sNjddzazXsBNwMlm1gboBewObAe8ama7uPvKTJ9Ipb75Bj79NHxMnRomGSsrg/nzw/62bcNiIL17w2675aQkEZFcS6fppj0wy91nA5jZCKAHkBr0PYABya+fAu42M0tuH+HuPwCfmtms5PebkJny17DvvrB4MXz7bRjY9P33P+0zCyNXjzwyTB7UtWu4mhcRKXLpBP32wJyUxxXAflUd4+4rzGwJsGVy+8Q1nrt9ZS9iZn2BvgDNanrzs3XrEOgNG4aPbbaBFi2gZUvYZRdo1Khm31dEpIClE/RWyTZP85h0nhs2ug8GBgMkEolKj6nWo4/W6GkiIsUsnZuxFUDTlMc7AHOrOsbM6gGbAYvSfK6IiGRROkFfBrQys5ZmVp9wc3XUGseMAs5Ifn0CMM7dPbm9V7JXTkugFfBeZkoXEZF0VNt0k2xz7w+MBeoCQ9x9upldD5S7+yjgIWBY8mbrIsIfA5LHPUG4cbsCuCBnPW5ERAQACxfe+SWRSHh5eXnsMkRECoaZTXL3RGX7NKmZiEiRU9CLiBQ5Bb2ISJFT0IuIFLm8vBlrZguAz2v49MbAVxkspxCU4jlDaZ53KZ4zlOZ5r+85N3f3ShfNyMugrw0zK6/qznOxKsVzhtI871I8ZyjN887kOavpRkSkyCnoRUSKXDEG/eDYBURQiucMpXnepXjOUJrnnbFzLro2ehER+blivKIXEZEUCnoRkSJXsEFfmwXLC1Ua53ypmc0ws6lm9pqZFcUK59Wdd8pxJ5iZm1nBd8NL55zN7KTkz3u6mT2W6xozLY3/383M7HUzez/5f/yoGHVmkpkNMbP5Zjativ1mZncm/02mmlm7Gr2QuxfcB2G65E+AHYH6wBSgzRrHnA/cl/y6F/B47LpzcM6dgI2TX59X6Oec7nknj9sUeIuwdGUidt05+Fm3At4HNk8+3ip23Tk458HAecmv2wCfxa47A+d9MNAOmFbF/qOA0YTV+vYH3q3J6xTqFf1/Fyx39+XA6gXLU/UA/pr8+ingsOSC5YWq2nN299fdffWK6BMJK3oVunR+1gA3ADcDy3JZXJakc87nAIPcfTGAu8/PcY2Zls45O7B64efNKILV6tz9LcIaHlXpATziwUTgF2a27fq+TqEGfWULlq+56PjPFiwHVi9YXqjSOedUZxOuBApdtedtZvsATd39hVwWlkXp/Kx3AXYxs7+b2UQz65qz6rIjnXMeAPQ2swrgJeDC3JQW1fr+3lcqncXB81FtFiwvVGmfj5n1BhJAx6xWlBvrPG8zqwMMBM7MVUE5kM7Puh6h+eYQwju3t81sD3f/Osu1ZUs653wKMNTdbzWzAwir2u3h7quyX140GcmxQr2ir82C5YUqrYXWzexw4Eqgu7v/kKPasqm6894U2AN4w8w+I7RjjirwG7Lp/v8e6e4/uvunwEeE4C9U6Zzz2cATAO4+AdiQMPFXMUvr9746hRr0tVmwvFBVe87JJoz7CSFf6G22q63zvN19ibs3dvcW7t6CcG+iu7sX8lqU6fz/fo5w8x0za0xoypmd0yozK51z/gI4DMDMWhOCfkFOq8y9UUCfZO+b/YEl7j5vfb9JQTbdeC0WLC9UaZ7zLcAmwJPJ+85fuHv3aEVnQJrnXVTSPOexQGczmwGsBH7n7gvjVV07aZ7zb4AHzOwSQvPFmQV+8YaZDSc0vzVO3nu4FtgAwN3vI9yLOAqYBXwPnFWj1ynwfycREalGoTbdiIhImhT0IiJFTkEvIlLkFPQiIkVOQS8iUuQU9CIiRU5BLyJS5BT0IlUwszbJ+e2nV7F/MzNbZWbFMgpZipSCXqRqq+fLmVTF/naESaf+kZtyRGpGQS9StdVBX9W8OfsmPyvoJa8p6EWqVt0VvYJeCoLmuhGphJnVBb4BGgCNUlbuSj1mJmFq4B2TUwWL5CVd0YtUrg2wMfBhFSHfCNgZWKyQl3ynoBepXLo3Yt/PTTkiNaegF6nc6qCfXMX+o5Kf1T4veU9BL1K51Tdal665w8w2A05NPqzqil8kbyjoRdaQXGN4r+TDU8ysQcq+bYHHge2Tmz7KcXki601BL7K2PQjrkVYA7YHPzewlM5sAfAJsCfyYPPZBM+sTp0yR9CjoRda2un3+NeA4YB5wKNCMsBbx4cA9hGadOoSFrUXylvrRi6zBzO4F+gH93X1Q7HpEaktX9CJrq27qA5GCoit6kRRmtgHwH8JFUCN3Xxa5JJFa0xW9yM/tSZj2YJpCXoqFruhFRIqcruhFRIqcgl5EpMgp6EVEipyCXkSkyCnoRUSKnIJeRKTIKehFRIrc/wOSBhjpd/mBUgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "posterior = Beta(mu, 2+1, 2+0)\n", "\n", "plot(mus, (sym.lambdify(mu,posterior))(mus), 'r')\n", "xlabel(\"$\\mu$\", size=22)\n", "\n", "maths(r\"\\mathbf{Posterior}: p(\\mu|m,l,a,b) = p(\\mu|1,0,2,2) = \\operatorname{Beta}(\\mu,2+1,2+0) = %s\" % sym.latex(posterior))" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.2" } }, "nbformat": 4, "nbformat_minor": 1 }