Adding a negative binomial regression model