Re-factor MCMC.* in order to avoid using system()