Code from Kruschke, Doing Bayesian Data Analysis