Helper functions for Bayesian analyses
