mgibbs function

Modified Gibbs sampler for Bayesian inference in linear mixed models