The analysis of binary longitudinal data with time-dependent covariates

• Published in: Statistics in medicine Vol. 31; no. 10; pp. 931 - 948
• Main Authors: Guerra, Matthew W., Shults, Justine, Amsterdam, Jay, Ten-Have, Thomas
• Format: Journal Article
• Language: English
• Published: Chichester, UK John Wiley & Sons, Ltd 10.05.2012; Wiley Subscription Services, Inc
• Subjects: alternating logistic regression; Antidepressive Agents, Second-Generation - therapeutic use; binary data; Computer Simulation; Cyclohexanols - therapeutic use; Data Interpretation, Statistical; Depression - drug therapy; Estimating techniques; generalized estimating equations; Humans; Likelihood Functions; Lithium - therapeutic use; Logistic Models; longitudinal binary data; Longitudinal Studies; maximum likelihood; Maximum likelihood method; Measurement; Regression Analysis; Simulation; Venlafaxine Hydrochloride
• ISSN: 0277-6715; 1097-0258; 1097-0258
• DOI: 10.1002/sim.4465

We consider longitudinal studies with binary outcomes that are measured repeatedly on subjects over time. The goal of our analysis was to fit a logistic model that relates the expected value of the outcomes with explanatory variables that are measured on each subject. However, additional care must be taken to adjust for the association between the repeated measurements on each subject. We propose a new maximum likelihood method for covariates that may be fixed or time varying. We also implement and make comparisons with two other approaches: generalized estimating equations, which may be more robust to misspecification of the true correlation structure, and alternating logistic regression, which models association via odds ratios that are subject to less restrictive constraints than are correlations. The proposed estimation procedure will yield consistent and asymptotically normal estimates of the regression and correlation parameters if the correlation on consecutive measurements on a subject is correctly specified. Simulations demonstrate that our approach can yield improved efficiency in estimation of the regression parameter; for equally spaced and complete data, the gains in efficiency were greatest for the parameter associated with a time‐by‐group interaction term and for stronger values of the correlation. For unequally spaced data and with dropout according to a missing‐at‐random mechanism, MARK1ML with correctly specified consecutive correlations yielded substantial improvements in terms of both bias and efficiency. We present an analysis to demonstrate application of the methods we consider. We also offer an R function for easy implementation of our approach. Copyright © 2012 John Wiley & Sons, Ltd.