“…Several studies recognizing the importance of unobserved heterogeneity have developed multivariate approaches that account for the potential dependency across count variables. The various model structures developed from multivariate models include multivariate Poisson regression model (Ye et al, 2009), multivariate Poisson lognormal model (Serhiyenko et al, 2016), multinomialgeneralized Poisson model (Chiou and Fu, 2013), multivariate Poisson gamma mixture count model (Mothafer et al, 2016), multivariate Poisson lognormal spatial and temporal model (Aguero-Valverde et al, 2016;Cheng et al, 2017), Integrated Nested Laplace Approximation Multivariate Poisson Lognormal model (Wang et al, 2017), Bayesian latent class flexible mixture multivariate model (Heydari et al, 2017) and multivariate random-parameters zeroinflated negative binomial model (Anastasopoulos, 2016).…”