Consistency Inference Property of QIC in Selecting the True Working Correlation Structure for Generalized Estimating Equations

Abstract

The generalized estimating equations (GEE) is one of the statistical approaches for the analysis of longitudinal data with correlated response. A working correlation structure for the repeated measures of the outcome variable of a subject needs to be specified by this method and the GEE estimator for the regression parameter will be the most efficient if the working correlation matrix is correctly specified. Hence it is desirable to choose a working correlation matrix that is the closest to the underlying structure among a set of working correlation structures. The quasi-likelihood Information criteria (QIC) was proposed for the selection of the working correlation structure and the best subset of explanatory variables in GEE. However, its success rate in selecting the true correlation structure has been established to be about 29.4%. Likewise, past studies have shown that its bias increases with the number of parameters. By considering longitudinal data with binary response, we establish numerically through simulations the consistency property of QIC in selecting the true working correlation structure and the conditions for its consistency. Further, we propose a modified QIC that penalizes for the number of parameter estimates in the original QIC and numerically establish that the penalization enhances the consistency of QIC in selecting the true working correlation structure. The results indicate that QIC selects the true correlation structure with probability approaching one if only parsimonious structures are considered otherwise the selection rates are less than 50% regardless of the increase in the sample size, measurements per subject and level of correlation. Further, we established that the probability of selecting the true correlation structure R 0 almost surely converges to one when we penalize for the number of correlation parameters estimated.

Profile of health facilities and determinants of the family planning services supply to unmarried adolescents in Burkina Faso, Ghana, and Niger

“…The choice of this model is justified by the fact that it considers the intra-subject correlation, which, if not considered, leads to invalid statistical inferences. The intra-subject correlations are considered by the specification of a correlation structure, "A working correlation structure" 21 , allowing sample population-level results.…”

“…The empirical validation of these models was carried out using the global validity test (Wald Khi-deux test) at the 5% threshold and using the quasi-likelihood under the Independence model Criterion (QIC) with different correlation matrices. [21][22][23] After numerous simulations, we selected the GEE models with a "compound symmetry" correlation matrix with the lowest QIC values. [21][22][23] The Odds Ratio (OR) is used to interpret the results of analysis produced using Stata software.…”

“…[21][22][23] After numerous simulations, we selected the GEE models with a "compound symmetry" correlation matrix with the lowest QIC values. [21][22][23] The Odds Ratio (OR) is used to interpret the results of analysis produced using Stata software.…”