AM_mcmc_fit function

Performs a Gibbs sampling