Methods for Assessing Longitudinal Biomarkers of Time-to-Event Outcomes in CKD: A Simulation Study

Identifying novel biomarkers is critical to advancing diagnosis and treatment of CKD, but relies heavily on the statistical methods used. Inappropriate methods can lead to both false positive and false negative associations between biomarkers and outcomes. This study assessed accuracy of methods using computer simulations and compared biomarker association estimates in the NEPhrotic syndrome sTUdy NEtwork (NEPTUNE), a prospective cohort study of patients with glomerular disease.

We compared three methods for analyzing repeatedly measured biomarkers in proportional hazards models: (1) time-invariant average, that averages values over all follow-up and uses the average as a baseline covariate, (2) time-varying last observation carried forward (LOCF), that assumes the covariate is unchanged until the next observed value, and (3) time-varying cumulative average, that updates the average using values at or before each measurement.

Under both true mechanisms of LOCF and cumulative average, simulation results showed the time-invariant average method often gave extremely inaccurate results. When LOCF was the true association mechanism, the cumulative average method often gave overestimated association estimates that were further away from the null. When cumulative average was the true mechanism, LOCF always underestimated the associations, i.e., closer to the null. In NEPTUNE, compared with the LOCF or cumulative average methods, hazard ratios estimated from the time-invariant average method were always higher.

Different analytic methods resulted in markedly different results. Using the time-invariant average produces inaccurate association estimates, whereas other methods can estimate additive (cumulative average) or instantaneous (LOCF) associations depending on the hypothesized underlying association mechanism and research question.