Abstract:A model fit by general estimating equations (GEE) has been used extensively for the analysis of longitudinal data in medical studies. To some extent, GEE tries to minimize a quadratic form of the residuals, and therefore is not robust in the sense that it, like least squares estimates, is sensitive to heavy-tailed distributions, contaminated distributions and extreme values. This paper describes the family of truncated robust estimating equations and its properties for the analysis of quantitative longitudinal data. Like GEE, the robust estimating equations aim to assess the covariate effects in the generalized linear model in the complete population of observations, but in a manner that is more robust to the influence of aberrant observations. A simulation study has been conducted to compare the finite-sample performance of GEE and the robust estimating equations under a variety of error distributions and data structures. It shows that the parameter estimates based on GEE and the robust estimating equations are approximately unbiased and the type I errors of Wald tests do not tend to be inflated. GEE is slightly more efficient with pure normal data, but the efficiency of GEE declines much more quickly than the robust estimating equations when the data become contaminated or have heavy tails, which makes the robust estimating equations advantageous with non-normal data. Both GEE and the robust estimating equations are applied to a longitudinal analysis of renal function in the Diabetes Control and Complications Trial (DCCT). For this application, GEE seems to be sensitive to the working correlation specification in that different working correlation structures may lead to different conclusions about the effect of intensive diabetes treatment. On the other hand, the robust estimating equations consistently conclude that the treatment effect is highly significant no matter which working correlation structure is used. The DCCT Research Group also demonstrated a significant effect using a mixed-effects longitudinal model.

geessbin: an R package for analyzing small-sample binary data using modified generalized estimating equations with bias-adjusted covariance estimators

A bias-reduced generalized estimating equation approach for proportional odds models with small-sample longitudinal ordinal data

Application of robust estimating equations to the analysis of quantitative longitudinal data

Sample Size and Power Calculations with Correlated Binary Data

Doubly Robust-Based Generalized Estimating Equations for the Analysis of Longitudinal Ordinal Missing Data

On GEE for Mean-Variance-Correlation Models: Variance Estimation and Model Selection

Robust Estimating Equations and Bias Correction of Correlation Parameters for Longitudinal Data

An efficient and robust method for analyzing population pharmacokinetic data in genome-wide pharmacogenomic studies: a generalized estimating equation approach

Generalized Estimating Equations in Longitudinal Studies: A Non-Parametric Alternative for Two-Way Repeated Measures Mixed ANOVA

The hidden factor: accounting for covariate effects in power and sample size computation for a binary trait

A Comparison of Variance Estimators for Logistic Regression Models Estimated Using Generalized Estimating Equations (GEE) in the Context of Observational Health Services Research

The Testing in GEE Method without PA Condition

Robust Estimating Functions and Bias Correction for Longitudinal Data Analysis

Outlier detection method in GEEs

Efficient Estimation for Marginal Generalized Partially Linear Single-Index Models with Longitudinal Data

Bias Correction Estimator for A Dynamic Panel Data Model with Fixed Effects Using an Iterated Bootstrap

Assessing the use of GEE methods for analyzing binary outcomes in family studies: the Strong Heart Family Study

Empirical Sandwich Variance Estimator for Iterated Conditional Expectation g‐Computation

Small sample adjustment for inference without assuming orthogonality in a mixed model for repeated measures analysis

On Binscatter

Robust estimation in joint mean–covariance regression model for longitudinal data