Abstract:Studies of health outcomes commonly involve binary measures, assessed multiple times. Although generalized linear mixed (GLIMMIX) models are well suited for analyzing these data, there does not exist a formal statistic to assess the goodness of fit (GOF) for GLIMMIX models. We developed an extension of the Hosmer-Lemeshow GOF test used for logistic regression that can be applied to GLIMMIX models. The correlation among repeated measurements of the binary outcome variable was accounted for by a random effect in the GLIMMIX model. The principles of Hosmer-Lemeshow method were followed. The linear unbiased estimate of dependent variables were transformed to the original probability, sorted from least to largest, and divided into deciles. A Chi-square statistic and corresponding p-value with eight degrees of freedom, was calculated based upon the expected and observed numbers among deciles. The proposed GOF test was validated by a simulation study with 1000 runs generated from logistic regression models with and without random effects. The results were compared with the Hosmer-Lemeshow GOF test in situations where the latter is appropriate. The proposed method was used in the analysis of a comparative effectiveness (CE) study of ophthalmologic treatments for open-angle glaucoma patients. When there was no random effect, the proposed GOF test results from the GLIMMIX procedure were almost identical to those of Hosmer-Lemeshow GOF test from the logistic procedure . With a random effect built in a correctly specified model, the goodness of fit rejection rate was 5.1%, which is close to the nominal level 5%. The proposed test did not indicate lack of fit for the regression models in the CE study. The proposed GOF test provides an assessment of model fit for models with binary outcomes and repeated measurements for predictor variables.

On Fitting Generalized Linear Mixed-Effects Models for Binary Responses Using Different Statistical Packages.

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

General Design Bayesian Generalized Linear Mixed Models

Multivariate probit linear mixed models for multivariate longitudinal binary data

A Semiparametric Bayesian Approach to Binomial Distribution Logistic Mixed-Effects Models for Longitudinal Data

Flexible Bayesian Modeling for Longitudinal Binary and Ordinal Responses

Fitting Longitudinal Mixed Effects Logistic Models in S-Plus

A Semiparametric Bayesian Approach to Generalized Partial Linear Mixed Models for Longitudinal Data.

Extension of the Hosmer-Lemeshow Goodness of Fit Statistic to Linear Models with Repeated Measurements

On the impact of parametric assumptions and robust alternatives for longitudinal data analysis.

Latent Gaussian Copula Models for Longitudinal Binary Data.

Statistical Inference in Generalized Linear Mixed Models by Joint Modelling Mean and Covariance of Non-Normal Random Effects

GENERALIZED LINEAR MIXED MODELS AND ITS APPLICATION

Modeling Heterogeneity and Missing Data of Multiple Longitudinal Outcomes in Electronic Health Records

Efficient Semiparametric Mean-Association Estimation for Longitudinal Binary Responses

Does the Package Matter? A Comparison of Five Common Multilevel Modeling Software Packages

Statistical methods for longitudinal and clustered designs with binary responses

Robust data integration from multiple external sources for generalized linear models with binary outcomes

glmmPen: High Dimensional Penalized Generalized Linear Mixed Models

A comprehensive analysis comparing linear and generalized linear models in detecting adaptive SNPs