How to fit Poisson data and do MCMC and stuff