fit_mcmc function

Fit the base imputation model using a Bayesian approach