{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Conjugate distributions" ] }, { "cell_type": "code", "execution_count": 83, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from numpy.random import beta\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Beta-Bernoulli\n", "\n", "Let's look at the shape of the Beta distribution" ] }, { "cell_type": "code", "execution_count": 84, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAO3ElEQVR4nO3cf4xlZX3H8fdHVrQ/1FV2IGR36WhcE4mJSia4jUmrrjGAjcsfYDC1rmbTjRYbG5tUWv9oa/sHNqm0Jsa6KcbFVIXSWjaW/qALxNZ00UUQRWpZKYXJbti1wFZDtKLf/nEfdNid2Tk7c+8M8+z7ldzcc57z3DnfZ+7MZ84895yTqkKS1JdnrXYBkqTxM9wlqUOGuyR1yHCXpA4Z7pLUoXWrXQDAhg0banp6erXLkKQ15c477/xOVU3Nt+0ZEe7T09McOHBgtcuQpDUlyX8vtM1pGUnqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6tAz4gpVrR3TV/39quz3wavfvCr7ldYqj9wlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShzwVUtJpb7VO8YXJnebrkbskdWjQkXuSB4HvAj8CnqyqmSQvAq4HpoEHgbdW1WNJAvw5cAnwBPDOqvrq+EvX6eR0vHjqdByzxudUjtxfX1WvqqqZtn4VsK+qtgD72jrAxcCW9tgFfHxcxUqShlnOnPt24HVteQ9wO/CB1n5dVRWwP8n6JOdW1eHlFKqfWs35wdON32utVUPDvYB/TlLAJ6pqN3DOU4FdVYeTnN36bgQenvPa2db2tHBPsovRkT3nnXfekgfQ4wchkrRcQ8P9tVV1qAX4LUn+4yR9M09bndAw+gOxG2BmZuaE7ZKkpRs0515Vh9rzEeDzwIXAI0nOBWjPR1r3WWDznJdvAg6Nq2BJ0uIWPXJP8nPAs6rqu235TcCHgL3ADuDq9nxTe8le4L1JPge8BjjmfLu0djjV2Ych0zLnAJ8fneHIOuAzVfWPSb4C3JBkJ/AQcHnrfzOj0yAPMjoV8l1jr1qSdFKLhntVPQC8cp72/wG2zdNewJVjqU6StCReoSpJHTLcJalDhrskdchwl6QOGe6S1CHv574M3ndE0jOV4S7pGcMDpvFxWkaSOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDg0O9yRnJLkryRfa+ouT3JHk/iTXJzmztT+nrR9s26cnU7okaSGncuT+PuC+OesfBq6pqi3AY8DO1r4TeKyqXgpc0/pJklbQoHBPsgl4M/CXbT3AG4AbW5c9wKVteXtbp23f1vpLklbI0CP3PwN+B/hxWz8LeLyqnmzrs8DGtrwReBigbT/W+j9Nkl1JDiQ5cPTo0SWWL0maz6LhnuRXgCNVdefc5nm61oBtP22o2l1VM1U1MzU1NahYSdIw6wb0eS3wliSXAM8Fns/oSH59knXt6HwTcKj1nwU2A7NJ1gEvAB4de+WSpAUteuReVb9bVZuqahq4Ari1qn4VuA24rHXbAdzUlve2ddr2W6vqhCN3SdLkLOc89w8A709ykNGc+rWt/VrgrNb+fuCq5ZUoSTpVQ6ZlfqKqbgdub8sPABfO0+f7wOVjqE2StEReoSpJHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6tGi4J3luki8n+VqSe5P8YWt/cZI7ktyf5PokZ7b257T1g2379GSHIEk63pAj9x8Ab6iqVwKvAi5KshX4MHBNVW0BHgN2tv47gceq6qXANa2fJGkFLRruNfK9tvrs9ijgDcCNrX0PcGlb3t7Wadu3JcnYKpYkLWrQnHuSM5LcDRwBbgG+DTxeVU+2LrPAxra8EXgYoG0/Bpw1z9fcleRAkgNHjx5d3igkSU8zKNyr6kdV9SpgE3Ah8PL5urXn+Y7S64SGqt1VNVNVM1NTU0PrlSQNcEpny1TV48DtwFZgfZJ1bdMm4FBbngU2A7TtLwAeHUexkqRhhpwtM5VkfVv+GeCNwH3AbcBlrdsO4Ka2vLet07bfWlUnHLlLkiZn3eJdOBfYk+QMRn8MbqiqLyT5JvC5JH8M3AVc2/pfC3w6yUFGR+xXTKBuSdJJLBruVXUP8Op52h9gNP9+fPv3gcvHUp0kaUm8QlWSOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1aNFwT7I5yW1J7ktyb5L3tfYXJbklyf3t+YWtPUk+muRgknuSXDDpQUiSnm7IkfuTwG9X1cuBrcCVSc4HrgL2VdUWYF9bB7gY2NIeu4CPj71qSdJJLRruVXW4qr7alr8L3AdsBLYDe1q3PcClbXk7cF2N7AfWJzl37JVLkhZ0SnPuSaaBVwN3AOdU1WEY/QEAzm7dNgIPz3nZbGs7/mvtSnIgyYGjR4+eeuWSpAUNDvckPw/8DfBbVfW/J+s6T1ud0FC1u6pmqmpmampqaBmSpAEGhXuSZzMK9r+qqr9tzY88Nd3Sno+09llg85yXbwIOjadcSdIQQ86WCXAtcF9VfWTOpr3Ajra8A7hpTvs72lkzW4FjT03fSJJWxroBfV4L/Brw9SR3t7bfA64GbkiyE3gIuLxtuxm4BDgIPAG8a6wVS5IWtWi4V9W/Mf88OsC2efoXcOUy65IkLYNXqEpShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOLRruST6Z5EiSb8xpe1GSW5Lc355f2NqT5KNJDia5J8kFkyxekjS/IUfunwIuOq7tKmBfVW0B9rV1gIuBLe2xC/j4eMqUJJ2KRcO9qr4IPHpc83ZgT1veA1w6p/26GtkPrE9y7riKlSQNs9Q593Oq6jBAez67tW8EHp7Tb7a1nSDJriQHkhw4evToEsuQJM1n3B+oZp62mq9jVe2uqpmqmpmamhpzGZJ0eltquD/y1HRLez7S2meBzXP6bQIOLb08SdJSLDXc9wI72vIO4KY57e9oZ81sBY49NX0jSVo56xbrkOSzwOuADUlmgd8HrgZuSLITeAi4vHW/GbgEOAg8AbxrAjVLkhaxaLhX1dsW2LRtnr4FXLncoiRJy+MVqpLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDEwn3JBcl+VaSg0mumsQ+JEkLG3u4JzkD+BhwMXA+8LYk5497P5KkhU3iyP1C4GBVPVBV/wd8Dtg+gf1IkhawbgJfcyPw8Jz1WeA1x3dKsgvY1Va/l+RbS9zfBuA7S3ztWuWYTw+O+TSQDy9rzL+w0IZJhHvmaasTGqp2A7uXvbPkQFXNLPfrrCWO+fTgmE8PkxrzJKZlZoHNc9Y3AYcmsB9J0gImEe5fAbYkeXGSM4ErgL0T2I8kaQFjn5apqieTvBf4J+AM4JNVde+49zPHsqd21iDHfHpwzKeHiYw5VSdMh0uS1jivUJWkDhnuktShNRPui93SIMlzklzftt+RZHrlqxyvAWN+f5JvJrknyb4kC57zulYMvXVFksuSVJI1f9rckDEneWt7r+9N8pmVrnHcBvxsn5fktiR3tZ/vS1ajznFJ8skkR5J8Y4HtSfLR9v24J8kFy95pVT3jH4w+mP028BLgTOBrwPnH9fkN4C/a8hXA9atd9wqM+fXAz7bl95wOY279ngd8EdgPzKx23SvwPm8B7gJe2NbPXu26V2DMu4H3tOXzgQdXu+5ljvmXgAuAbyyw/RLgHxhdJ7QVuGO5+1wrR+5DbmmwHdjTlm8EtiWZ74KqtWLRMVfVbVX1RFvdz+iagrVs6K0r/gj4E+D7K1nchAwZ868DH6uqxwCq6sgK1zhuQ8ZcwPPb8gtY49fKVNUXgUdP0mU7cF2N7AfWJzl3OftcK+E+3y0NNi7Up6qeBI4BZ61IdZMxZMxz7WT0l38tW3TMSV4NbK6qL6xkYRM05H1+GfCyJF9Ksj/JRStW3WQMGfMfAG9PMgvcDPzmypS2ak71931Rk7j9wCQMuaXBoNserCGDx5Pk7cAM8MsTrWjyTjrmJM8CrgHeuVIFrYAh7/M6RlMzr2P039m/JnlFVT0+4domZciY3wZ8qqr+NMkvAp9uY/7x5MtbFWPPr7Vy5D7klgY/6ZNkHaN/5U72b9Az3aDbOCR5I/BB4C1V9YMVqm1SFhvz84BXALcneZDR3OTeNf6h6tCf7Zuq6odV9V/AtxiF/Vo1ZMw7gRsAqurfgecyuqlYr8Z+25a1Eu5DbmmwF9jRli8Dbq32ScUateiY2xTFJxgF+1qfh4VFxlxVx6pqQ1VNV9U0o88Z3lJVB1an3LEY8rP9d4w+PCfJBkbTNA+saJXjNWTMDwHbAJK8nFG4H13RKlfWXuAd7ayZrcCxqjq8rK+42p8in8KnzZcA/8noU/YPtrYPMfrlhtGb/9fAQeDLwEtWu+YVGPO/AI8Ad7fH3tWuedJjPq7v7azxs2UGvs8BPgJ8E/g6cMVq17wCYz4f+BKjM2nuBt602jUvc7yfBQ4DP2R0lL4TeDfw7jnv8cfa9+Pr4/i59vYDktShtTItI0k6BYa7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6tD/A1xhMaiKhNxjAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "a, b = 1, 1\n", "densities = plt.hist([beta(a,b) for _ in range(5000)], bins=10)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Suppose we have observed a sequence of $n$ coin flips with $x$. Then the posterior -- what we know believe about $p$, given our prior and the new evidence -- also follows a Beta distribution.\n", "\n", "Here is the derivation (from the [Wiki page](https://en.wikipedia.org/wiki/Conjugate_prior)), where $s$ denotes \"successes\" (heads, here) and $f$ failures (so $n-x$)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "![alt text](https://wikimedia.org/api/rest_v1/media/math/render/svg/2fcd7af2c9aacbb32c7cc5ded34c8204f554484a \"Derivation\")\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is a Beta distribution with parameters $\\alpha + s, \\beta + f$." ] }, { "cell_type": "code", "execution_count": 85, "metadata": {}, "outputs": [], "source": [ "def posterior(a, b, x, n):\n", " samples = [beta(a + x, b + n - x) for _ in range(10000)]\n", " densities = plt.hist(samples, bins=10)\n", " return np.array(samples)" ] }, { "cell_type": "code", "execution_count": 88, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAPPElEQVR4nO3df4ylVX3H8fdHVrT+BNnF0N1tB+PaSkwayQSxJta6xgg2LH9Ag6l1JZtuYtFaMa3Y/kGjSYP9RUtisFuXujQWodSUjaU1hB+xbbrEQSzyo4YtWpiylbHA9gexSv32j3vQYZndvTvPzB0u5/1KJvd5znOee85hZj/33PPc+5CqQpLUh+etdQckSZNj6EtSRwx9SeqIoS9JHTH0Jakj69a6A0eyfv36mpmZWetuSNJUueOOO75dVRuWOvasDv2ZmRnm5ubWuhuSNFWS/Ovhjrm8I0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHXlWfyNXktbSzCV/vWZtf/Oyd67K8zrTl6SOONPXilirGdFqzYak5ypDX1NtLd9+rxVf6DTEUZd3klyV5JEkdy8qe0WSm5Lc3x5PbOVJckWS/UnuSnL6onO2t/r3J9m+OsORJB3JOGv6nwHecUjZJcDNVbUFuLntA5wFbGk/O4ErYfQiAVwKvAE4A7j0qRcKSdLkHHV5p6q+lGTmkOJtwFva9h7gNuAjrfzqqipgX5ITkpzS6t5UVY8CJLmJ0QvJNYNHIOk5r8dlvNWy3E/vvLKqDgC0x5Nb+UbgoUX15lvZ4cqfIcnOJHNJ5hYWFpbZPUnSUlb6Qm6WKKsjlD+zsGoXsAtgdnZ2yTrj8hMlkvR0yw39byU5paoOtOWbR1r5PLB5Ub1NwMOt/C2HlN+2zLZ1GL4FlnQ0yw39vcB24LL2eMOi8vcn+Ryji7YH2wvDF4HfXnTx9u3AR5ffbUlrwYnF9Dtq6Ce5htEsfX2SeUafwrkMuC7JDuBB4PxW/UbgbGA/8ARwIUBVPZrk48CXW72PPXVR97nIfxiSnq3G+fTOuw5zaOsSdQu46DDPcxVw1TH1TpK0orz3jiR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakj/p+zpCnjN741hDN9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUkUGhn+RDSe5JcneSa5K8MMmpSW5Pcn+Sa5Mc3+q+oO3vb8dnVmIAkqTxLTv0k2wEfgWYrarXAccBFwCfAC6vqi3AY8COdsoO4LGqejVweasnSZqgocs764AfSbIOeBFwAHgrcH07vgc4t21va/u041uTZGD7kqRjsOzQr6p/A34PeJBR2B8E7gAer6onW7V5YGPb3gg81M59stU/6dDnTbIzyVySuYWFheV2T5K0hCHLOycymr2fCvwo8GLgrCWq1lOnHOHYDwuqdlXVbFXNbtiwYbndkyQtYcjyztuAb1TVQlV9D/g88NPACW25B2AT8HDbngc2A7TjLwceHdC+JOkYDQn9B4Ezk7yorc1vBe4FbgXOa3W2Aze07b1tn3b8lqp6xkxfkrR6hqzp387oguxXgK+159oFfAS4OMl+Rmv2u9spu4GTWvnFwCUD+i1JWoZ1R69yeFV1KXDpIcUPAGcsUfc7wPlD2pMkDeM3ciWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6sig0E9yQpLrk/xzkvuSvDHJK5LclOT+9nhiq5skVyTZn+SuJKevzBAkSeMaOtP/I+Bvq+ongZ8C7gMuAW6uqi3AzW0f4CxgS/vZCVw5sG1J0jFadugneRnwZmA3QFV9t6oeB7YBe1q1PcC5bXsbcHWN7ANOSHLKsnsuSTpmQ2b6rwIWgD9NcmeSTyd5MfDKqjoA0B5PbvU3Ag8tOn++lT1Nkp1J5pLMLSwsDOieJOlQQ0J/HXA6cGVVvR74H364lLOULFFWzyio2lVVs1U1u2HDhgHdkyQdakjozwPzVXV727+e0YvAt55atmmPjyyqv3nR+ZuAhwe0L0k6RssO/ar6d+ChJD/RirYC9wJ7ge2tbDtwQ9veC7ynfYrnTODgU8tAkqTJWDfw/A8An01yPPAAcCGjF5LrkuwAHgTOb3VvBM4G9gNPtLqSpAkaFPpV9VVgdolDW5eoW8BFQ9qTJA3jN3IlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SerI4NBPclySO5N8oe2fmuT2JPcnuTbJ8a38BW1/fzs+M7RtSdKxWYmZ/geB+xbtfwK4vKq2AI8BO1r5DuCxqno1cHmrJ0maoEGhn2QT8E7g020/wFuB61uVPcC5bXtb26cd39rqS5ImZOhM/w+BXwe+3/ZPAh6vqifb/jywsW1vBB4CaMcPtvpPk2RnkrkkcwsLCwO7J0labNmhn+TngEeq6o7FxUtUrTGO/bCgaldVzVbV7IYNG5bbPUnSEtYNOPdNwDlJzgZeCLyM0cz/hCTr2mx+E/Bwqz8PbAbmk6wDXg48OqB9SdIxWvZMv6o+WlWbqmoGuAC4pap+AbgVOK9V2w7c0Lb3tn3a8Vuq6hkzfUnS6lmNz+l/BLg4yX5Ga/a7W/lu4KRWfjFwySq0LUk6giHLOz9QVbcBt7XtB4AzlqjzHeD8lWhPkrQ8fiNXkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOLDv0k2xOcmuS+5Lck+SDrfwVSW5Kcn97PLGVJ8kVSfYnuSvJ6Ss1CEnSeIbM9J8EPlxVrwXOBC5KchpwCXBzVW0Bbm77AGcBW9rPTuDKAW1LkpZh2aFfVQeq6itt+7+A+4CNwDZgT6u2Bzi3bW8Drq6RfcAJSU5Zds8lScdsRdb0k8wArwduB15ZVQdg9MIAnNyqbQQeWnTafCs79Ll2JplLMrewsLAS3ZMkNYNDP8lLgL8EfrWq/vNIVZcoq2cUVO2qqtmqmt2wYcPQ7kmSFhkU+kmezyjwP1tVn2/F33pq2aY9PtLK54HNi07fBDw8pH1J0rEZ8umdALuB+6rqDxYd2gtsb9vbgRsWlb+nfYrnTODgU8tAkqTJWDfg3DcBvwh8LclXW9lvAJcB1yXZATwInN+O3QicDewHngAuHNC2JGkZlh36VfX3LL1OD7B1ifoFXLTc9iRJw/mNXEnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOjLx0E/yjiRfT7I/ySWTbl+SejbR0E9yHPBJ4CzgNOBdSU6bZB8kqWeTnumfAeyvqgeq6rvA54BtE+6DJHVr3YTb2wg8tGh/HnjD4gpJdgI72+5/J/n6gPbWA98ecP606W284Jh70d2Y84lBY/7xwx2YdOhnibJ62k7VLmDXijSWzFXV7Eo81zTobbzgmHvhmFfOpJd35oHNi/Y3AQ9PuA+S1K1Jh/6XgS1JTk1yPHABsHfCfZCkbk10eaeqnkzyfuCLwHHAVVV1zyo2uSLLRFOkt/GCY+6FY14hqaqj15IkPSf4jVxJ6oihL0kdmfrQP9ptHZK8IMm17fjtSWYm38uVNcaYL05yb5K7ktyc5LCf2Z0W496+I8l5SSrJ1H+8b5wxJ/n59ru+J8mfT7qPK22Mv+0fS3Jrkjvb3/fZa9HPlZLkqiSPJLn7MMeT5Ir23+OuJKcPbrSqpvaH0cXgfwFeBRwP/BNw2iF1fhn4VNu+ALh2rfs9gTH/LPCitv2+Hsbc6r0U+BKwD5hd635P4Pe8BbgTOLHtn7zW/Z7AmHcB72vbpwHfXOt+Dxzzm4HTgbsPc/xs4G8YfcfpTOD2oW1O+0x/nNs6bAP2tO3rga1JlvqS2LQ46pir6taqeqLt7mP0fYhpNu7tOz4O/A7wnUl2bpWMM+ZfAj5ZVY8BVNUjE+7jShtnzAW8rG2/nCn/nk9VfQl49AhVtgFX18g+4IQkpwxpc9pDf6nbOmw8XJ2qehI4CJw0kd6tjnHGvNgORjOFaXbUMSd5PbC5qr4wyY6tonF+z68BXpPkH5LsS/KOifVudYwz5t8C3p1kHrgR+MBkurZmjvXf+1FN+jYMK+2ot3UYs840GXs8Sd4NzAI/s6o9Wn1HHHOS5wGXA++dVIcmYJzf8zpGSzxvYfRu7u+SvK6qHl/lvq2Wccb8LuAzVfX7Sd4I/Fkb8/dXv3trYsXza9pn+uPc1uEHdZKsY/SW8Ehvp57txrqVRZK3Ab8JnFNV/zuhvq2Wo435pcDrgNuSfJPR2ufeKb+YO+7f9g1V9b2q+gbwdUYvAtNqnDHvAK4DqKp/BF7I6GZsz1UrfuuaaQ/9cW7rsBfY3rbPA26pdoVkSh11zG2p448ZBf60r/PCUcZcVQeran1VzVTVDKPrGOdU1dzadHdFjPO3/VeMLtqTZD2j5Z4HJtrLlTXOmB8EtgIkeS2j0F+YaC8nay/wnvYpnjOBg1V1YMgTTvXyTh3mtg5JPgbMVdVeYDejt4D7Gc3wL1i7Hg835ph/F3gJ8BftmvWDVXXOmnV6oDHH/Jwy5pi/CLw9yb3A/wG/VlX/sXa9HmbMMX8Y+JMkH2K0zPHeaZ7EJbmG0fLc+nad4lLg+QBV9SlG1y3OBvYDTwAXDm5ziv97SZKO0bQv70iSjoGhL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjry/1v46pyDSvMQAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "samples = posterior(.01, 1, 1, 1)" ] }, { "cell_type": "code", "execution_count": 77, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.6677448353537242" ] }, "execution_count": 77, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.mean(samples)" ] }, { "cell_type": "code", "execution_count": 78, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAASKUlEQVR4nO3df4xlZ13H8feHVqooSGUHUne7TiFbY2l0gUmtMWBJEdpqWjCIu1FasHEBWyNCjIv+UYJpgkAlErG4yKatkZZqxW5sEWtFqoYFplDKtlCZlkKH3XRXihVTrXb5+sc9C5fd+XHn3jt3On3er+Rmzn3Oc875Pt3pZ84859wzqSokSW140loXIEmaHENfkhpi6EtSQwx9SWqIoS9JDTl+rQtYzoYNG2p6enqty5CkdeP222//96qaWmjd4z70p6enmZ2dXesyJGndSPKVxdY5vSNJDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ153H8iV5LW0vTOm9bkuPe//edWZb+e6UtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGLBv6SXYnOZhkX1/bh5Lc0b3uT3JH1z6d5L/71r2vb5sXJPl8krkk70mS1RmSJGkxg9ynfxXwx8A1Rxqq6peOLCe5Ani4r/+9VbV1gf1cCewA9gI3A+cAH1l5yZKkYS17pl9VtwEPLbSuO1t/FXDtUvtIchLwtKr6RFUVvR8gL195uZKkUYw6p/9C4MGq+lJf2ylJPpvk40le2LVtBOb7+sx3bQtKsiPJbJLZQ4cOjViiJOmIUUN/O999ln8A2FxVzwPeBHwwydOAhebva7GdVtWuqpqpqpmpqQX/oLskaQhDP3snyfHALwAvONJWVY8Cj3bLtye5FziV3pn9pr7NNwH7hz22JGk4o5zpvwT4YlV9e9omyVSS47rlZwNbgPuq6gDwzSRndtcBLgRuHOHYkqQhDHLL5rXAJ4AfTTKf5OJu1TaOvYD7IuDOJJ8D/gp4fVUduQj8BuDPgDngXrxzR5Imbtnpnaravkj7axZouwG4YZH+s8DpK6xPktbs8cZPRH4iV5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWrIIH8YfXeSg0n29bW9NcnXktzRvc7rW/eWJHNJ7knysr72c7q2uSQ7xz8USdJyBjnTvwo4Z4H2d1fV1u51M0CS04BtwHO7bf4kyXFJjgPeC5wLnAZs7/pKkibo+OU6VNVtSaYH3N8FwHVV9Sjw5SRzwBndurmqug8gyXVd37tXXLEkaWijzOlfmuTObvrnxK5tI/BAX5/5rm2xdknSBA0b+lcCzwG2AgeAK7r2LNC3lmhfUJIdSWaTzB46dGjIEiVJRxsq9Kvqwao6XFXfAt7Pd6Zw5oGT+7puAvYv0b7Y/ndV1UxVzUxNTQ1ToiRpAUOFfpKT+t6+AjhyZ88eYFuSE5KcAmwBPgV8GtiS5JQkT6Z3sXfP8GVLkoax7IXcJNcCZwEbkswDlwFnJdlKb4rmfuB1AFV1V5Lr6V2gfQy4pKoOd/u5FPgocBywu6ruGvtoJK2q6Z03rXUJGtEgd+9sX6D5A0v0vxy4fIH2m4GbV1SdJGms/ESuJDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1JBlQz/J7iQHk+zra3tnki8muTPJh5M8vWufTvLfSe7oXu/r2+YFST6fZC7Je5JkdYYkSVrMIGf6VwHnHNV2C3B6Vf048G/AW/rW3VtVW7vX6/varwR2AFu619H7lCStsuOX61BVtyWZPqrt7/ve7gVeudQ+kpwEPK2qPtG9vwZ4OfCRFdYrNW96501rXYLWsXHM6f8q3x3epyT5bJKPJ3lh17YRmO/rM9+1SZImaNkz/aUk+T3gMeAvuqYDwOaq+nqSFwB/k+S5wELz97XEfnfQmwpi8+bNo5QoSeoz9Jl+kouAnwd+uaoKoKoeraqvd8u3A/cCp9I7s9/Ut/kmYP9i+66qXVU1U1UzU1NTw5YoSTrKUKGf5Bzgd4Dzq+qRvvapJMd1y8+md8H2vqo6AHwzyZndXTsXAjeOXL0kaUWWnd5Jci1wFrAhyTxwGb27dU4AbunuvNzb3anzIuBtSR4DDgOvr6qHul29gd6dQN9H7xqAF3ElacIGuXtn+wLNH1ik7w3ADYusmwVOX1F1kqSx8hO5ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQ0Z6nr7UMv+CldYjz/QlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWrIQKGfZHeSg0n29bX9UJJbknyp+3pi154k70kyl+TOJM/v2+airv+Xklw0/uFIkpYy6Jn+VcA5R7XtBG6tqi3Ard17gHOBLd1rB3Al9H5IAJcBPwmcAVx25AeFJGkyBgr9qroNeOio5guAq7vlq4GX97VfUz17gacnOQl4GXBLVT1UVd8AbuHYHySSpFU0ypz+s6rqAED39Zld+0bggb5+813bYu2SpAlZjQu5WaCtlmg/dgfJjiSzSWYPHTo01uIkqWWjhP6D3bQN3deDXfs8cHJfv03A/iXaj1FVu6pqpqpmpqamRihRktRvlNDfAxy5A+ci4Ma+9gu7u3jOBB7upn8+Crw0yYndBdyXdm2SpAkZ6NHKSa4FzgI2JJmndxfO24Hrk1wMfBX4xa77zcB5wBzwCPBagKp6KMnvA5/u+r2tqo6+OCxJWkUDhX5VbV9k1dkL9C3gkkX2sxvYPXB1kqSx8hO5ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1ZKD79KXHs+mdN611CdK64Zm+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhoydOgn+dEkd/S9/jPJG5O8NcnX+trP69vmLUnmktyT5GXjGYIkaVBDP2Wzqu4BtgIkOQ74GvBh4LXAu6vqXf39k5wGbAOeC/ww8A9JTq2qw8PWIElamXFN75wN3FtVX1mizwXAdVX1aFV9GZgDzhjT8SVJAxjX8/S3Adf2vb80yYXALPDmqvoGsBHY29dnvms7RpIdwA6AzZs3j6lErSafaS+tDyOf6Sd5MnA+8Jdd05XAc+hN/RwArjjSdYHNa6F9VtWuqpqpqpmpqalRS5QkdcYxvXMu8JmqehCgqh6sqsNV9S3g/XxnCmceOLlvu03A/jEcX5I0oHGE/nb6pnaSnNS37hXAvm55D7AtyQlJTgG2AJ8aw/ElSQMaaU4/yVOAnwVe19f8jiRb6U3d3H9kXVXdleR64G7gMeAS79yRpMkaKfSr6hHgGUe1vXqJ/pcDl49yTEnS8PxEriQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktSQkf5GLkCS+4FvAoeBx6pqJskPAR8Cpun9cfRXVdU3kgT4I+A84BHgNVX1mVFr0HdM77xprUuQ9Dg2rjP9F1fV1qqa6d7vBG6tqi3Ard17gHOBLd1rB3DlmI4vSRrAak3vXABc3S1fDby8r/2a6tkLPD3JSatUgyTpKOMI/QL+PsntSXZ0bc+qqgMA3ddndu0bgQf6tp3v2r5Lkh1JZpPMHjp0aAwlSpJgDHP6wE9X1f4kzwRuSfLFJfpmgbY6pqFqF7ALYGZm5pj1kqThjHymX1X7u68HgQ8DZwAPHpm26b4e7LrPAyf3bb4J2D9qDZKkwYwU+km+P8lTjywDLwX2AXuAi7puFwE3dst7gAvTcybw8JFpIEnS6ht1eudZwId7d2JyPPDBqvq7JJ8Grk9yMfBV4Be7/jfTu11zjt4tm68d8fiSpBUYKfSr6j7gJxZo/zpw9gLtBVwyyjElScPzE7mS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNGcdfztJRpnfetNYlSNKCPNOXpIYY+pLUEENfkhpi6EtSQwx9SWrI0KGf5OQkH0vyhSR3JfnNrv2tSb6W5I7udV7fNm9JMpfkniQvG8cAJEmDG+WWzceAN1fVZ5I8Fbg9yS3dundX1bv6Oyc5DdgGPBf4YeAfkpxaVYdHqEGStAJDn+lX1YGq+ky3/E3gC8DGJTa5ALiuqh6tqi8Dc8AZwx5fkrRyY5nTTzINPA/4ZNd0aZI7k+xOcmLXthF4oG+zeRb5IZFkR5LZJLOHDh0aR4mSJMYQ+kl+ALgBeGNV/SdwJfAcYCtwALjiSNcFNq+F9llVu6pqpqpmpqamRi1RktQZKfSTfA+9wP+LqvprgKp6sKoOV9W3gPfznSmceeDkvs03AftHOb4kaWVGuXsnwAeAL1TVH/a1n9TX7RXAvm55D7AtyQlJTgG2AJ8a9viSpJUb5e6dnwZeDXw+yR1d2+8C25NspTd1cz/wOoCquivJ9cDd9O78ucQ7dyRpsoYO/ar6Fxaep795iW0uBy4f9piSpNH4iVxJaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWrIKI9heNyb3nnTWpcgSY8rnulLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDJh76Sc5Jck+SuSQ7J318SWrZREM/yXHAe4FzgdOA7UlOm2QNktSySZ/pnwHMVdV9VfW/wHXABROuQZKaNekHrm0EHuh7Pw/85NGdkuwAdnRv/yvJPSs4xgbg34eucH1qcczQ5rhbHDM0OO78wUhj/pHFVkw69LNAWx3TULUL2DXUAZLZqpoZZtv1qsUxQ5vjbnHM0Oa4V2vMk57emQdO7nu/Cdg/4RokqVmTDv1PA1uSnJLkycA2YM+Ea5CkZk10eqeqHktyKfBR4Dhgd1XdNebDDDUttM61OGZoc9wtjhnaHPeqjDlVx0ypS5KeoPxEriQ1xNCXpIas29Bf7nEOSU5I8qFu/SeTTE++yvEaYMxvSnJ3kjuT3Jpk0Xt115NBH92R5JVJKsm6v7VvkDEneVX3731Xkg9OusZxG+D7e3OSjyX5bPc9ft5a1DlOSXYnOZhk3yLrk+Q93X+TO5M8f+SDVtW6e9G7CHwv8GzgycDngNOO6vPrwPu65W3Ah9a67gmM+cXAU7rlN6z3MQ867q7fU4HbgL3AzFrXPYF/6y3AZ4ETu/fPXOu6JzDmXcAbuuXTgPvXuu4xjPtFwPOBfYusPw/4CL3POJ0JfHLUY67XM/1BHudwAXB1t/xXwNlJFvpw2Hqx7Jir6mNV9Uj3di+9z0Gsd4M+uuP3gXcA/zPJ4lbJIGP+NeC9VfUNgKo6OOEax22QMRfwtG75B3kCfManqm4DHlqiywXANdWzF3h6kpNGOeZ6Df2FHuewcbE+VfUY8DDwjIlUtzoGGXO/i+mdIax3y447yfOAk6vqbydZ2Coa5N/6VODUJP+aZG+ScyZW3eoYZMxvBX4lyTxwM/AbkyltTa30//tlTfoxDOMyyOMcBnrkwzoy8HiS/AowA/zMqlY0GUuOO8mTgHcDr5lUQRMwyL/18fSmeM6i9xvdPyc5var+Y5VrWy2DjHk7cFVVXZHkp4A/78b8rdUvb82MPcfW65n+II9z+HafJMfT+3VwqV+jHu8GeoRFkpcAvwecX1WPTqi21bTcuJ8KnA78U5L76c177lnnF3MH/f6+sar+r6q+DNxD74fAejXImC8Grgeoqk8A30vvQWxPZGN/dM16Df1BHuewB7ioW34l8I/VXRlZp5YdczfN8af0An+9z/EeseS4q+rhqtpQVdNVNU3vWsb5VTW7NuWOxSDf339D78I9STbQm+65b6JVjtcgY/4qcDZAkh+jF/qHJlrl5O0BLuzu4jkTeLiqDoyyw3U5vVOLPM4hyduA2araA3yA3q9/c/TO8LetXcWjG3DM7wR+APjL7pr1V6vq/DUregwGHPcTyoBj/ijw0iR3A4eB366qr69d1aMZcMxvBt6f5LfoTXG8Zp2fyJHkWnpTdBu6axWXAd8DUFXvo3ft4jxgDngEeO3Ix1zn/80kSSuwXqd3JElDMPQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQ/4fIrgUSTCZZoMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "samples = posterior_under_prior(1, .01, 1, 2)" ] }, { "cell_type": "code", "execution_count": 79, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAASVUlEQVR4nO3df4xlZ13H8feHFuoPQKo7kLrbdQrZGkujC0xqDQGrRegP0xaDuptoCzYuYGtUjHGRP0ogTeqPStKIhUU2bQ30B1TohhaxVrRqWGAKddkWKtOy0mE33YVqxVSrLV//uGfhsjs/7sy9c6ezz/uV3My53/Occ56nO/3MmeeceyZVhSSpDc9Y7Q5IksbH0Jekhhj6ktQQQ1+SGmLoS1JDjl/tDixm3bp1NTk5udrdkKQ145577vl6VU3Mte5pH/qTk5NMT0+vdjckac1I8m/zrXN6R5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGvK0/0TuMCa3374qx9131fmrclxJWoxn+pLUEENfkhpi6EtSQwx9SWrIoqGfZGeSg0n29tVuTnJv99qX5N6uPpnkv/vWvadvm5cl+UKSmSTXJMnKDEmSNJ9B7t65Dvgz4IbDhar65cPLSa4GHutr/2BVbZ5jP9cC24DdwB3AOcDHl95lSdJyLRr6VXV3ksm51nVn678E/OxC+0hyEvDcqvpU9/4G4CKO0dBfrVtFwdtFJS1s2Dn9VwCPVNWX+2qnJPl8kn9I8oquth6Y7Wsz29UkSWM07IeztgI39r0/AGysqm8keRnw0SQvBuaav6/5dppkG72pIDZu3DhkFyVJhy37TD/J8cAvADcfrlXVE1X1jW75HuBB4FR6Z/Yb+jbfAOyfb99VtaOqpqpqamJizr/tK0lahmGmd14FfKmqvj1tk2QiyXHd8guBTcBDVXUA+GaSM7vrABcDtw1xbEnSMgxyy+aNwKeAH00ym+TSbtUWvntqB+CVwJ4k/wJ8GHhTVT3arXsz8BfADL3fAI7Ji7iS9HQ2yN07W+epv36O2q3ArfO0nwZOX2L/JEkj5CdyJakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUkEVDP8nOJAeT7O2rvT3J15Lc273O61v31iQzSR5I8pq++jldbSbJ9tEPRZK0mEHO9K8Dzpmj/q6q2ty97gBIchqwBXhxt82fJzkuyXHAu4FzgdOArV1bSdIYHb9Yg6q6O8nkgPu7ELipqp4AvpJkBjijWzdTVQ8BJLmpa3v/knssSVq2Yeb0L0+yp5v+ObGrrQce7msz29Xmq88pybYk00mmDx06NEQXJUn9lhv61wIvAjYDB4Cru3rmaFsL1OdUVTuqaqqqpiYmJpbZRUnSkRad3plLVT1yeDnJ+4CPdW9ngZP7mm4A9nfL89UlSWOyrDP9JCf1vX0tcPjOnl3AliQnJDkF2AR8BvgssCnJKUmeRe9i767ld1uStByLnuknuRE4C1iXZBa4AjgryWZ6UzT7gDcCVNV9SW6hd4H2SeCyqnqq28/lwCeA44CdVXXfyEcjSVrQIHfvbJ2j/P4F2l8JXDlH/Q7gjiX1TpI0Usua09fT1+T221fluPuuOn9VjitpaXwMgyQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhiwa+kl2JjmYZG9f7Y+TfCnJniQfSfK8rj6Z5L+T3Nu93tO3zcuSfCHJTJJrkmRlhiRJms8gZ/rXAeccUbsTOL2qfhz4V+CtfeserKrN3etNffVrgW3Apu515D4lSSts0T+MXlV3J5k8ovY3fW93A69baB9JTgKeW1Wf6t7fAFwEfHyJ/dXT1Gr9QXbwj7JLSzGKOf1f47vD+5Qkn0/yD0le0dXWA7N9bWa72pySbEsynWT60KFDI+iiJAmGDP0kbwOeBD7QlQ4AG6vqJcBbgA8meS4w1/x9zbffqtpRVVNVNTUxMTFMFyVJfRad3plPkkuAnwfOrqoCqKongCe65XuSPAicSu/MfkPf5huA/cs9tiRpeZZ1pp/kHOD3gQuq6vG++kSS47rlF9K7YPtQVR0AvpnkzO6unYuB24buvSRpSRY9009yI3AWsC7JLHAFvbt1TgDu7O683N3dqfNK4B1JngSeAt5UVY92u3ozvTuBvpfeNQAv4krSmA1y987WOcrvn6ftrcCt86ybBk5fUu8kSSPlJ3IlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDVn2s3ekp4vVeqyzj3TWWuSZviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGDBT6SXYmOZhkb1/tB5PcmeTL3dcTu3qSXJNkJsmeJC/t2+aSrv2Xk1wy+uFIkhYy6Jn+dcA5R9S2A3dV1Sbgru49wLnApu61DbgWej8kgCuAnwTOAK44/INCkjQeA4V+Vd0NPHpE+ULg+m75euCivvoN1bMbeF6Sk4DXAHdW1aNV9e/AnRz9g0SStIKGmdN/QVUdAOi+Pr+rrwce7ms329Xmq0uSxmQlLuRmjlotUD96B8m2JNNJpg8dOjTSzklSy4YJ/Ue6aRu6rwe7+ixwcl+7DcD+BepHqaodVTVVVVMTExNDdFGS1G+Y0N8FHL4D5xLgtr76xd1dPGcCj3XTP58AXp3kxO4C7qu7miRpTAb6y1lJbgTOAtYlmaV3F85VwC1JLgW+Cvxi1/wO4DxgBngceANAVT2a5J3AZ7t276iqIy8OS5JW0EChX1Vb51l19hxtC7hsnv3sBHYO3DtJ0kj5iVxJaoihL0kNGWh6R9LRJrffvirH3XfV+atyXB0bPNOXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhiw79JP8aJJ7+17/meS3k7w9ydf66uf1bfPWJDNJHkjymtEMQZI0qGX/ucSqegDYDJDkOOBrwEeANwDvqqo/6W+f5DRgC/Bi4IeBv01yalU9tdw+SJKWZlTTO2cDD1bVvy3Q5kLgpqp6oqq+AswAZ4zo+JKkAYwq9LcAN/a9vzzJniQ7k5zY1dYDD/e1me1qR0myLcl0kulDhw6NqIuSpKFDP8mzgAuAD3Wla4EX0Zv6OQBcfbjpHJvXXPusqh1VNVVVUxMTE8N2UZLUGcWZ/rnA56rqEYCqeqSqnqqqbwHv4ztTOLPAyX3bbQD2j+D4kqQBjSL0t9I3tZPkpL51rwX2dsu7gC1JTkhyCrAJ+MwIji9JGtCy794BSPJ9wM8Bb+wr/1GSzfSmbvYdXldV9yW5BbgfeBK4zDt3JGm8hgr9qnoc+KEjar+6QPsrgSuHOabUusntt6/asfdddf6qHVuj4SdyJakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYMHfpJ9iX5QpJ7k0x3tR9McmeSL3dfT+zqSXJNkpkke5K8dNjjS5IGN6oz/Z+pqs1VNdW93w7cVVWbgLu69wDnApu61zbg2hEdX5I0gJWa3rkQuL5bvh64qK9+Q/XsBp6X5KQV6oMk6QjHj2AfBfxNkgLeW1U7gBdU1QGAqjqQ5Pld2/XAw33bzna1A/07TLKN3m8CbNy4cQRdlDQKk9tvX5Xj7rvq/FU57rFoFKH/8qra3wX7nUm+tEDbzFGrowq9Hxw7AKampo5aL0lanqGnd6pqf/f1IPAR4AzgkcPTNt3Xg13zWeDkvs03APuH7YMkaTBDhX6S70/ynMPLwKuBvcAu4JKu2SXAbd3yLuDi7i6eM4HHDk8DSZJW3rDTOy8APpLk8L4+WFV/neSzwC1JLgW+Cvxi1/4O4DxgBngceMOQx5ckLcFQoV9VDwE/MUf9G8DZc9QLuGyYY0qSls9P5EpSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGjKKP4wuSStqcvvtq3bsfVedv2rHXgme6UtSQwx9SWrIskM/yclJPpnki0nuS/JbXf3tSb6W5N7udV7fNm9NMpPkgSSvGcUAJEmDG2ZO/0ngd6vqc0meA9yT5M5u3buq6k/6Gyc5DdgCvBj4YeBvk5xaVU8N0QdJ0hIs+0y/qg5U1ee65W8CXwTWL7DJhcBNVfVEVX0FmAHOWO7xJUlLN5I5/SSTwEuAT3ely5PsSbIzyYldbT3wcN9ms8zzQyLJtiTTSaYPHTo0ii5KkhhB6Cd5NnAr8NtV9Z/AtcCLgM3AAeDqw03n2Lzm2mdV7aiqqaqampiYGLaLkqTOUKGf5Jn0Av8DVfVXAFX1SFU9VVXfAt7Hd6ZwZoGT+zbfAOwf5viSpKUZ5u6dAO8HvlhVf9pXP6mv2WuBvd3yLmBLkhOSnAJsAj6z3ONLkpZumLt3Xg78KvCFJPd2tT8AtibZTG/qZh/wRoCqui/JLcD99O78ucw7dyRpvJYd+lX1T8w9T3/HAttcCVy53GNKkobjJ3IlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDfFv5ErSAlbr7/Ou1N/m9Uxfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqyNhDP8k5SR5IMpNk+7iPL0ktG2voJzkOeDdwLnAasDXJaePsgyS1bNxn+mcAM1X1UFX9L3ATcOGY+yBJzRr3o5XXAw/3vZ8FfvLIRkm2Adu6t/+V5IFlHm8d8PVlbrtWtThmcNwtaWLM+cOjSksZ94/Mt2LcoZ85anVUoWoHsGPogyXTVTU17H7WkhbHDI57tfsxTi2OGUY37nFP78wCJ/e93wDsH3MfJKlZ4w79zwKbkpyS5FnAFmDXmPsgSc0a6/ROVT2Z5HLgE8BxwM6qum8FDzn0FNEa1OKYwXG3pMUxw4jGnaqjptQlSccoP5ErSQ0x9CWpIcdE6C/2aIckJyS5uVv/6SST4+/laA0w5rckuT/JniR3JZn3vt21ZNDHeCR5XZJKsuZv7RtkzEl+qfv3vi/JB8fdx5UwwPf4xiSfTPL57vv8vNXo5ygl2ZnkYJK986xPkmu6/yZ7krx0yQepqjX9ondB+EHghcCzgH8BTjuizW8A7+mWtwA3r3a/xzDmnwG+r1t+81of86Dj7to9B7gb2A1MrXa/x/BvvQn4PHBi9/75q93vMY17B/Dmbvk0YN9q93sE434l8FJg7zzrzwM+Tu8zT2cCn17qMY6FM/1BHu1wIXB9t/xh4Owkc31QbK1YdMxV9cmqerx7u5veZyLWukEf4/FO4I+A/xln51bIIGP+deDdVfXvAFV1cMx9XAmDjLuA53bLP8Ax8JmfqrobeHSBJhcCN1TPbuB5SU5ayjGOhdCf69EO6+drU1VPAo8BPzSW3q2MQcbc71J6Zwdr3aLjTvIS4OSq+tg4O7aCBvm3PhU4Nck/J9md5Jyx9W7lDDLutwO/kmQWuAP4zfF0bVUt9f/9o4z7MQwrYZBHOwz0+Ic1ZODxJPkVYAr46RXt0XgsOO4kzwDeBbx+XB0ag0H+rY+nN8VzFr3f6P4xyelV9R8r3LeVNMi4twLXVdXVSX4K+Mtu3N9a+e6tmqGz7Fg40x/k0Q7fbpPkeHq/Ci70K9TT3UCPs0jyKuBtwAVV9cSY+raSFhv3c4DTgb9Pso/enOeuNX4xd9Dv79uq6v+q6ivAA/R+CKxlg4z7UuAWgKr6FPA99B5Kdiwb+lE2x0LoD/Joh13AJd3y64C/q+6qyBq16Ji7aY730gv8Y2GOFxYZd1U9VlXrqmqyqibpXcu4oKqmV6e7IzHI9/dH6V24J8k6etM9D421l6M3yLi/CpwNkOTH6IX+obH2cvx2ARd3d/GcCTxWVQeWsoM1P71T8zzaIck7gOmq2gW8n96vfjP0zvC3rF6PhzfgmP8YeDbwoe6a9Ver6oJV6/QIDDjuY8qAY/4E8Ook9wNPAb9XVd9YvV4Pb8Bx/y7wviS/Q2+K4/Vr/GSOJDfSm6Zb112ruAJ4JkBVvYfetYvzgBngceANSz7GGv9vJElagmNhekeSNCBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXk/wH0zA2I5mgRuQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "samples = posterior_under_prior(.01, 1, 1, 2)" ] }, { "cell_type": "code", "execution_count": 91, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAARkklEQVR4nO3df4xl5V3H8fenS0GjjSwyNLis7la3ScFEqBNKbDS1VdjSxKXRGvqHXRt0NYGoiTHZ+g9KJUGjkjRWdBs2bo0Wib9YWyKuWKya0DJYpCxImC5YpruBsUtbmyq6+PWP+2y87N6ZufPrzu4+71dyc8/5nufc8zxz4DNnzjn3bKoKSVIfXrPRHZAkTY6hL0kdMfQlqSOGviR1xNCXpI6ct9EdWMzFF19c27Zt2+huSNJZ5dFHH/33qpoateyMDv1t27YxMzOz0d2QpLNKkn9baJmndySpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSNn9DdypaVs2/uJDdv2c3e8a8O2La2UR/qS1JElQz/JNyT5TJJ/SXI4ya+2+vYkn07yTJI/SXJ+q1/Q5mfb8m1Dn/WBVn86yXXrNShJ0mjjHOm/DLy9qr4HuBLYmeQa4NeBO6tqB/AScFNrfxPwUlV9F3Bna0eSy4EbgSuAncDvJtm0loORJC1uydCvga+12de2VwFvB/601Q8AN7TpXW2etvwdSdLq91TVy1X1LDALXL0mo5AkjWWsc/pJNiV5DHgROAR8HvhyVZ1oTeaALW16C/A8QFv+FeBbh+sj1hne1p4kM0lm5ufnlz8iSdKCxgr9qnqlqq4ELmNwdP6mUc3aexZYtlD91G3tq6rpqpqemhr5bwBIklZoWXfvVNWXgYeAa4ALk5y85fMy4GibngO2ArTl3wIcH66PWEeSNAHj3L0zleTCNv2NwA8BTwGfBH6sNdsN3NemD7Z52vK/q6pq9Rvb3T3bgR3AZ9ZqIJKkpY3z5axLgQPtTpvXAPdW1ceTPAnck+TXgM8Cd7f2dwN/mGSWwRH+jQBVdTjJvcCTwAng5qp6ZW2HI0lazJKhX1WPA1eNqB9hxN03VfVfwHsW+KzbgduX301J0lrwG7mS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOLBn6SbYm+WSSp5IcTvLzrf4rSb6Y5LH2un5onQ8kmU3ydJLrhuo7W202yd71GZIkaSHnjdHmBPCLVfXPSV4HPJrkUFt2Z1X95nDjJJcDNwJXAN8G/G2SN7bFHwZ+GJgDHklysKqeXIuBSJKWtmToV9Ux4Fib/o8kTwFbFlllF3BPVb0MPJtkFri6LZutqiMASe5pbQ19SZqQZZ3TT7INuAr4dCvdkuTxJPuTbG61LcDzQ6vNtdpC9VO3sSfJTJKZ+fn55XRPkrSEsUM/yTcDfwb8QlV9FbgL+E7gSgZ/CfzWyaYjVq9F6q8uVO2rqumqmp6amhq3e5KkMYxzTp8kr2UQ+H9UVX8OUFUvDC3/CPDxNjsHbB1a/TLgaJteqC5JmoBx7t4JcDfwVFX99lD90qFm7waeaNMHgRuTXJBkO7AD+AzwCLAjyfYk5zO42HtwbYYhSRrHOEf6bwV+Avhcksda7ZeB9ya5ksEpmueAnwGoqsNJ7mVwgfYEcHNVvQKQ5BbgAWATsL+qDq/hWCRJSxjn7p1/ZPT5+PsXWed24PYR9fsXW0+StL78Rq4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHVkrAeuSUvZtvcTG92FiduoMT93x7s2ZLs6N3ikL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkeWDP0kW5N8MslTSQ4n+flWvyjJoSTPtPfNrZ4kH0oym+TxJG8e+qzdrf0zSXav37AkSaOMc6R/AvjFqnoTcA1wc5LLgb3Ag1W1A3iwzQO8E9jRXnuAu2DwSwK4FXgLcDVw68lfFJKkyVgy9KvqWFX9c5v+D+ApYAuwCzjQmh0AbmjTu4CP1sDDwIVJLgWuAw5V1fGqegk4BOxc09FIkha1rHP6SbYBVwGfBl5fVcdg8IsBuKQ12wI8P7TaXKstVD91G3uSzCSZmZ+fX073JElLGDv0k3wz8GfAL1TVVxdrOqJWi9RfXajaV1XTVTU9NTU1bvckSWMYK/STvJZB4P9RVf15K7/QTtvQ3l9s9Tlg69DqlwFHF6lLkiZknLt3AtwNPFVVvz206CBw8g6c3cB9Q/X3tbt4rgG+0k7/PABcm2Rzu4B7batJkibkvDHavBX4CeBzSR5rtV8G7gDuTXIT8AXgPW3Z/cD1wCzwdeD9AFV1PMkHgUdau9uq6viajEKSNJYlQ7+q/pHR5+MB3jGifQE3L/BZ+4H9y+mgJGnt+I1cSeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SerIkqGfZH+SF5M8MVT7lSRfTPJYe10/tOwDSWaTPJ3kuqH6zlabTbJ37YciSVrKOEf6fwDsHFG/s6qubK/7AZJcDtwIXNHW+d0km5JsAj4MvBO4HHhvaytJmqDzlmpQVZ9Ksm3Mz9sF3FNVLwPPJpkFrm7LZqvqCECSe1rbJ5fdY0nSiq3mnP4tSR5vp382t9oW4PmhNnOttlD9NEn2JJlJMjM/P7+K7kmSTrXS0L8L+E7gSuAY8FutnhFta5H66cWqfVU1XVXTU1NTK+yeJGmUJU/vjFJVL5ycTvIR4ONtdg7YOtT0MuBom16oLkmakBUd6Se5dGj23cDJO3sOAjcmuSDJdmAH8BngEWBHku1JzmdwsffgyrstSVqJJY/0k3wMeBtwcZI54FbgbUmuZHCK5jngZwCq6nCSexlcoD0B3FxVr7TPuQV4ANgE7K+qw2s+GknSosa5e+e9I8p3L9L+duD2EfX7gfuX1TtJ0pryG7mS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOLBn6SfYneTHJE0O1i5IcSvJMe9/c6knyoSSzSR5P8uahdXa39s8k2b0+w5EkLWacI/0/AHaeUtsLPFhVO4AH2zzAO4Ed7bUHuAsGvySAW4G3AFcDt578RSFJmpwlQ7+qPgUcP6W8CzjQpg8ANwzVP1oDDwMXJrkUuA44VFXHq+ol4BCn/yKRJK2zlZ7Tf31VHQNo75e0+hbg+aF2c622UP00SfYkmUkyMz8/v8LuSZJGWesLuRlRq0Xqpxer9lXVdFVNT01NrWnnJKl3561wvReSXFpVx9rpmxdbfQ7YOtTuMuBoq7/tlPpDK9y2FrFt7yc2uguSzmArPdI/CJy8A2c3cN9Q/X3tLp5rgK+00z8PANcm2dwu4F7bapKkCVrySD/JxxgcpV+cZI7BXTh3APcmuQn4AvCe1vx+4HpgFvg68H6Aqjqe5IPAI63dbVV16sVhSdI6WzL0q+q9Cyx6x4i2Bdy8wOfsB/Yvq3eSpDXlN3IlqSMrvZAraYNs5MX65+5414ZtW2vDI31J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOrKq0E/yXJLPJXksyUyrXZTkUJJn2vvmVk+SDyWZTfJ4kjevxQAkSeNbiyP9H6yqK6tqus3vBR6sqh3Ag20e4J3AjvbaA9y1BtuWJC3Depze2QUcaNMHgBuG6h+tgYeBC5Ncug7blyQtYLWhX8DfJHk0yZ5We31VHQNo75e0+hbg+aF151rtVZLsSTKTZGZ+fn6V3ZMkDTtvleu/taqOJrkEOJTkXxdpmxG1Oq1QtQ/YBzA9PX3acknSyq3qSL+qjrb3F4G/AK4GXjh52qa9v9iazwFbh1a/DDi6mu1LkpZnxaGf5JuSvO7kNHAt8ARwENjdmu0G7mvTB4H3tbt4rgG+cvI0kCRpMlZzeuf1wF8kOfk5f1xVf53kEeDeJDcBXwDe09rfD1wPzAJfB96/im1LklZgxaFfVUeA7xlR/xLwjhH1Am5e6fYkSavnN3IlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6stp/OUsjbNv7iY3ugiSNZOhLGttGHdA8d8e7NmS75yJP70hSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqyMRDP8nOJE8nmU2yd9Lbl6SeTfTZO0k2AR8GfhiYAx5JcrCqnlyP7fngM+ncsJH/L59rz/2Z9JH+1cBsVR2pqv8G7gF2TbgPktStST9lcwvw/ND8HPCW4QZJ9gB72uzXkjw9ob4t18XAv290J9aJYzs7ObZ1kF+fyGbWenzfsdCCSYd+RtTqVTNV+4B9k+nOyiWZqarpje7HenBsZyfHdvaa5PgmfXpnDtg6NH8ZcHTCfZCkbk069B8BdiTZnuR84Ebg4IT7IEndmujpnao6keQW4AFgE7C/qg5Psg9r6Iw/BbUKju3s5NjOXhMbX6pq6VaSpHOC38iVpI4Y+pLUEUN/hHEfFZHkx5JUkumh2gfaek8nuW4yPR7fSseWZFuS/0zyWHv93uR6PZ6lxpbkJ5PMD43hp4aW7U7yTHvtnmzPl7bKsb0yVD/jbpwY57/JJD+e5Mkkh5P88VD9rN5vrc1CY1uf/VZVvoZeDC4wfx54A3A+8C/A5SPavQ74FPAwMN1ql7f2FwDb2+ds2ugxrdHYtgFPbPQYVjM24CeB3xmx7kXAkfa+uU1v3ugxrcXY2rKvbfQYVjm2HcBnT+4T4JJzaL+NHNt67jeP9E837qMiPgj8BvBfQ7VdwD1V9XJVPQvMts87U6xmbGe61Tzi4zrgUFUdr6qXgEPAznXq50qcy48vGWdsPw18uO0bqurFVj8X9ttCY1s3hv7pRj0qYstwgyRXAVur6uPLXXeDrWZsANuTfDbJ3yf5/nXs50qM+7P/0SSPJ/nTJCe/KHjW77dm1NgAviHJTJKHk9ywrj1dvnHG9kbgjUn+qY1h5zLW3UirGRus036b9GMYzgaLPioiyWuAOxn8Ob2sdc8AqxnbMeDbq+pLSb4X+MskV1TVV9elp8s3zs/+r4CPVdXLSX4WOAC8fcx1N9JqxgaD/XY0yRuAv0vyuar6/Dr2dznGGdt5DE6DvI3Bt/j/Icl3j7nuRlrx2Krqy6zTfvNI/3RLPSridcB3Aw8leQ64BjjYLnie6Y+ZWPHY2imrLwFU1aMMzlW+cSK9Hs+SP/uq+lJVvdxmPwJ877jrbrDVjI2qOtrejwAPAVetZ2eXaZyf/RxwX1X9Tztt+jSDoDzr9xsLj2399ttGX+w4014MfvMeYXAh9uTFlysWaf8Q/3+x8wpefSH3CGfWhdzVjG3q5FgYXJj6InDRRo9pOWMDLh2afjfwcJu+CHiWwcXAzW36XBnbZuCCNn0x8AwjLt6f4WPbCRwYGsPzwLeeI/ttobGt237z9M4paoFHRSS5DZipqgVvnWrt7gWeBE4AN1fVKxPp+BhWMzbgB4DbkpwAXgF+tqqOr3+vxzPm2H4uyY8w2DfHaaexqup4kg8yeDYUwG3nytiANwG/n+R/Gfxlf0et0z9atBJjju0B4NokTzL4b++Xqv3VeQ7st5FjS/J9rNN+8zEMktQRz+lLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktSR/wPmAbzUepH7HgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "samples = posterior_under_prior(1, .001, 100, 200)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Normal-Normal\n", "\n", "There is a nice interactive demo here: https://stephens999.github.io/fiveMinuteStats/shiny_normal_example.html" ] }, { "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.3" } }, "nbformat": 4, "nbformat_minor": 2 }