bivariate negative binomial regression model