The Negative Binomial Distribution (NBD) Regression is used when modeling an overdispersed count variable.
This article describes how to create a Negative Binomial Distribution (NBD) regression output as shown below. The example below models a survey respondent’s fastfood consumption based on characteristics like age, gender, and work status.
Requirements
 A count Outcome variable to be predicted. Ideally, a Numeric variable. The count variable must only include positive integers. When using stacked data the Outcome variable should be a single question in a Multi type structure (eg. NumericMulti).
 Continuous, categorical, or binary Predictor variables will be considered as predictors of the outcome variable. When using stacked data the Predictor(s) need to be a single question in a Grid type structure (Binary  Grid).
Method
 Go to Anything > Advanced Analysis > Regression > NBD Regression.
 In the object inspector, go to the Inputs tab.
 In the Outcome dropdown, select the numeric variable to be predicted by the predictor variables.
 Select the predictor variable(s) from the Predictor(s) dropdown.
 OPTIONAL: Select the desired Output type:
 Summary: The default; as shown in the example above.
 Detail: Typical R output, some additional information compared to Summary, but without the pretty formatting.
 ANOVA: Analysis of variance table containing the results of Chisquared likelihood ratio tests for each predictor.
 Relative Importance Analysis: The results of a relative importance analysis.
 Effects Plot Plots the relationship between each of the Predictors and the Outcome.
 OPTIONAL: Select the desired Missing Data treatment. (See Missing Data Options).
 OPTIONAL: Select Variable names to display variable names in the output instead of labels.
 OPTIONAL: Select Correction. Choose between None (the default), False Discovery Rate, Bonferroni. This is the multiple comparisons correction applied when computing the pvalues of the posthoc comparisons.
 OPTIONAL: Select Crosstab Interaction. Optional variable to test for interaction with other variables in the model. The interaction variable is treated as a categorical variable.
 OPTIONAL: Specify the Automated outlier removal percentage (between 0 and 50, including 0 but not 50) to remove possible outliers. See below for details about fitted models and residuals.
 OPTIONAL: Select Stack data to stack the input data prior to analysis. Stacking can be desirable when each individual in the data set has multiple cases and an aggregate model is desired. See requirements below.
 OPTIONAL: Update Random seed. This is used to initialize the (pseudo)random number generator for the model fitting algorithm. Different seeds may lead to slightly different answers, but should normally not make a large difference.
Additional Properties
When using this feature you can obtain additional information that is stored by the R code which produces the output.
 To do so, select Calculation > Custom.
 In the R CODE, paste: item = YourReferenceName
 Replace YourReferenceName with the reference name of your output. Find this by selecting the output and then going to Properties > GENERAL > Name from the object inspector.
 Below the first line of code, you can paste in snippets from below or type in str(item) to see a list of available information.
For a more in depth discussion on extracting information from objects in R, checkout our blog post here.
Properties which may be of interest are:
 Summary outputs from the regression model:

 item$summary$coefficients # summary regression outputs
Technical Details
If Crosstab Interaction is selected, coefficients in the table are computed by creating separate regressions for each level of the interaction variable. To evaluate whether a coefficient is significantly higher (blue) or lower (red), we perform a ttest of the coefficient compared to the coefficient using the remaining data as described in Driver Analysis. Pvalues are corrected for multiple comparisons across the whole table (excluding the NET column). The Pvalue in the subtitle is calculated using a the likelihood ratio test between the pooled model with no interaction variable, and a model where all predictors interact with the interaction variable.
If a zerovalue Automated outlier removal percentage is indicated, then no outlier removal is performed and a standard regression output for the entire (possibly filtered) dataset is applied. If a nonzero value is selected for this option then the regression model is fitted twice. The first regression model uses the entire dataset (after filters have been applied) and identifies the observations that generate the largest residuals. The user specified percent of cases in the data that have the largest residuals are then removed. The regression model is refitted on this reduced dataset and output returned. The specific residual used varies depending on the regression type. For NBD regression, a studentized deviance residual in an unweighted regression and the Pearson residual in a weighted regression. The studentized deviance residual computes the contribution the fitted point has to the likelihood and standardizes (adjusts) based on the influence of the point and an externally adjusted variance calculation (see rstudent function in R and Davison and Snell (1991)^{} for more details). The Pearson residual in the weighted case computes the distance between the observed and fitted value and adjusts appropriately for the provided survey weights. See rstudent function in R and Davison and Snell (1991) for more details of the specifics of the calculations.
If Stack data is selected, then the Outcome needs to be a single v . Similarly, the Predictor(s) need to be a single v . In the process of stacking, the data reduction is inspected. Any constructed NETs are removed unless comprised of source values that are mutually exclusive to other codes, such as the result of merging two categories.
References
Davison, A. C. and Snell, E. J. (1991) Residuals and diagnostics. In: Statistical Theory and Modelling. In Honour of Sir David Cox, FRS, eds. Hinkley, D. V., Reid, N. and Snell, E. J., Chapman & Hall.
See Also
How to Run Binary Logit Regression
How to Run a Generalized Linear Model
How to Run a Multinomial Logit Regression
How to Run Ordered Logit Regression
How to Run QuasiPoisson Regression
How to Run a Stepwise Regression
How to Create Regression Multicollinearity Table (VIF)
How to Create a PredictionAccuracy Table
How to Create a GoodnessofFit Plot
How to Save Predicted Values of Regression Models
How to Save Fitted Values of Regression Models
How to Test Residual Normality (ShapiroWilk) of Regression Models
How to Test Residual Serial Correlation (DurbinWatson) of Regression Models
How to Save Residuals of Regression Models
Comments
0 comments
Article is closed for comments.