When modelling multivariate binomial data, it often occurs that it is necessary to take into consideration both clustering and overdispersion, the former arising from the dependence between data, and the latter due to the additional variability in the data not prescribed by the distribution. If interest lies in accommodating both phenomena at the same time, we can use separate sets of random effects that capture the within-cluster association and the extra variability. In particular, the random effects for overdispersion can be included in the model either additively or multiplicatively. For this purpose, we propose a series of Bayesian hierarchical models that deal simultaneously with both phenomena. The proposed models are applied to bivariate repeated prevalence data for hepatitis C virus (HCV) and human immunodeficiency virus (HIV) infection in injecting drug users in Italy from 1998 to 2007.
BreslowNEClaytonDG (1993) Approximate inference in generalized linear mixed models. Journal of the American Statistical Association, 88, 9–25.
7.
CareyVZegerSLDiggleP (1993) Modelling multivariate binary data with alternating logistic regression. Biometrika, 80, 517–26.
8.
CeleuxGForbesFRobertCPTitteringtonDM (2006) Deviance information criteria for missing data models. Bayesian Analysis, 1, 651–74.
9.
ChenJJAhnH (1997) Marginal models with multiplicative variance components for overdispersed binomial data. Journal of Agricultural, Biological, and Environmental Statistics, 2, 440–50.
10.
ClaytonD (1996) Generalized linear mixed models. In GilksWRRichardsonSSpiegelhalterDJ (eds), Markov Chain Monte Carlo in practice. London: Chapman & Hall.
11.
CoutinhoRA (1998). HIV and hepatitis C among injecting drug users. BMJ (Clinical research ed.), 317(7156), 424–25.
12.
Del FavaEKasimAUsmanMShkedyZHensNAertsMBollaertsKScalia TombaGVickermanPSuttonAJWiessingLKretzschmarM (2011) Joint modeling of HCV and HIV infections among injecting drug users in italy using repeated cross-sectional prevalence data. Statistical Communications in Infectious Diseases, 3, 1–24.
13.
GelmanARubinDB (1992) Inference from iterative simulation using multiple sequences. Statistical Science, 7, 457–511.
14.
GelmanAMengXLSternH (1996) Posterior predictive assessment of model fitness via realized discrepancies. Statistica Sinica, 6, 733–807.
15.
GelmanA (2006) Prior distributions for variance parameters in hierarchical models. Bayesian Analysis, 1, 515–33.
16.
GewekeJ (1992) Evaluating the accuracy of sampling-based approaches to calculating posterior moments. In BernadoJMBergerJODawidAPSmithAFM (eds), Bayesian Statistics 4. Oxford, UK: Clarendon Press.
17.
HindeJDemétrioCGB (1998) Overdis-persion: models and estimation. Computational Statistics and Data Analysis, 27, 151–70.
18.
LesaffreELawsonAB (2012) Bayesian biostatistics. Hoboken, NJ: John Wiley & Sons.
19.
McLachlanGJ (1997) On the EM algorithm for overdispersed count data. Statistical Methods in Medical Research, 6, 76–98.
20.
MolenberghsGVerbekeG (2005) Models for discrete longitudinal data. Berlin: Springer-Verlag.
21.
MolenberghsGVerbekeGDemétrioCGB (2007) An extended random-effects approach to modeling repeated, overdispersed count data. Lifetime Data Analysis, 13, 513–31.
22.
MolenberghsGVerbekeGDemétrioCGBVieiraAMC (2010) A family of generalized linear models for repeated measures with normal and conjugate random effects. Statistical Science, 25, 325–47.
23.
NeelonBO’MalleyANormandS (2010) A Bayesian model for repeated measures zero-inflated count data with application to outpatient psychiatric service use. Statistical Modelling, 10, 421–39.
24.
O’MalleyAJZaslavskyAM (2008) Domain-level covariance analysis for multilevel survey data with structured nonresponse. Journal of the American Statistical Association, 103, 1405–18.
25.
PlummerM (2003) JAGS: a program for analysis of Bayesian graphical models using Gibbs sampling. Proceedings of the 3rd International Workshop on Distributed Statistical Computing (DSC 2003), March20–22, Vienna, Austria. ISSN 1609–395X.
26.
PlummerM (2011) JAGS Version 3.3. 0 user manual. International Agency for Research on Cancer.
27.
PlummerMBestNCowlesKVinesK (2006) CODA: convergence diagnosis and output analysis for MCMC, R News, 6, 7–11.
28.
PlummerM (2008) Penalized loss functions for Bayesian model comparison. Biostatistics, 9, 523–39.
29.
SmithB.J. (2007) Boa: An R package for MCMC output convergence assessment and posterior inference. Journal of Statistical Software, 21, 1–37.
30.
QaqishBFLiangK-Y (1992) Marginal models for correlated binary responses with multiple classes and multiple levels of nesting. Biometrics, 48, 939–50.
31.
SkellamJG (1948) A probability distribution derived from the binomial distribution by regarding the probability of success as variable between the sets of trials. Journal of the Royal Statistical Society, Series B, 10, 257–61.
32.
SpiegelhalterDJBestNJCarlinBPVan der LindeA (2002) Bayesian measures of model complexity and fit. Journal of the Royal Statistical Society, series B, 64, 583–640.
