I am currently teaching a class on Advanced Numerical Methods in Biomedical Engineering (BME 502) at Stony Brook University and I wanted to include some basic Bayesian Analysis. I chose to go through some of the examples in “Bayesian Data Analysis” by A. Gelman, J.B. Carlin, H.S. Stern, D.B. Dunson, A.Vehtari and D.B. Rubin. Since my class is taught in Python/Numpy/Scipy/Matplotlib, I translated his examples into an iPython notebook.
1000 draws from the posterior LD50 distribution for the bioassay data set from Gelman’s “Bayesian Data Analysis” book.