test the mcmc code separately using mcmc_time.R