Multivariate correlated data with various types are encountered in many research areas. Such data are generally more difficult to analyze due to the scarcity of appropriate statistic models. We demonstrate that the multivariate negative binomial model can be easily modified to provide asymptotically legitimate inferences about the regression parameters for correlated data of mixing types. Contrast is also made with models proposed in the literatures.