{
"metadata": {
"language": "Julia",
"name": "",
"signature": "sha256:4a4b90c4d7b08c3f7ad49111ea29f8fe4161d3bb8517528d70647fd4490ac22d"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Based on https://probmods.org/hierarchical-models.html"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"using Stochy, Gadfly, Stochy.GadflySupport"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"posterior = @pp pmcmc(100,100) do\n",
" local\n",
" colors = [:black, :blue, :green, :orange, :red],\n",
" prototype = 5 * ~Dir(5, 1.0),\n",
" bag1 = ~Dir(prototype),\n",
" bag2 = ~Dir(prototype),\n",
" bag3 = ~Dir(prototype),\n",
" bagn = ~Dir(prototype)\n",
" \n",
" observe(Discrete(colors, bag1), :blue, :blue, :black, :blue, :blue, :blue)\n",
" observe(Discrete(colors, bag2), :blue, :green, :blue, :blue, :blue, :red)\n",
" observe(Discrete(colors, bag3), :blue, :blue, :blue, :blue, :blue, :orange)\n",
" \n",
" # Posterior predictive.\n",
" ~Discrete(colors, bagn)\n",
"end;"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 6
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"bar(posterior)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 7,
"svg": [
"\n",
"\n"
],
"text": [
"Plot(...)"
]
}
],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}