“…First-order estimating equations using conditional residuals in a model for the association between clustered outcomes can be moderately efficient in relation to second-order equations with less burden computationally (Lipsitz and Fitzmaurice, 1996). In the case of multivariate binary data, alternating logistics regressions (ALR), which specify within-cluster associations in term of pairwise odds ratios (Carey, Zeger, and Diggle, 1993), have been applied in diverse settings, including studies of multilevel geographical clustering of drug (Petronis and Anthony, 2003) and alcohol use (Reboussin et al, 2012), surveillance of occupational illnesses with workplace clustering (Preisser, Arcury, and Quandt, 2003) and clusterrandomized trials with longitudinal data (Yi and Cook, 2002;Yi, Cook, and Chen, 2010). Heagerty and Zeger (1996) defined estimating equations for the associations of correlated ordinal data in an extension of Carey et al's method.…”