A Python package for defining hierarchical Bayesian models, with routines to optimise and sample from the posterior distribution.