MCMC for all species; and some output