Functions for hierarchical Bayesian estimation: a flexible approach
