Fitting Bayesian Poisson regression
