Functions for mcmc, in progress
Functions for mcmc, in progress