Diagnostic Tools for a Multivariate Negative Binomial Model for Fitting Correlated Data with Overdispersion