Changed parameterisation of the Poisson model
The shape of the gamma distribution for the dispersion parameter of the Negative Binomial is now treated as free parameter