Robust likelihood inference for multivariate correlated count data