Multivariate count data generalized linear models : Three approaches based on the Sarmanov distribution