stochastic MCMC epidemiological compartment model using Guillespie tau leap