Python (PyMC) adaptation of the R code from "Doing Bayesian Data Analysis"