Dynamic survival prediction of end-stage kidney disease using random survival forests for competing risk analysis

A static predictive model relying solely on baseline clinicopathological data cannot capture the heterogeneity in predictor trajectories observed in the progression of chronic kidney disease (CKD). To address this, we developed and validated a dynamic survival prediction model using longitudinal clinicopathological data to predict end-stage kidney disease (ESKD), with death as a competing risk.

We trained a sequence of random survival forests using a landmarking approach and optimized the model with a pre-specified prediction horizon of 5 years. The predicted cumulative incidence function (CIF) values were used to generate a personalized dynamic prediction plot.

The model was developed using baseline demographics and 13 longitudinal clinicopathological variables from 4,950 patients. Variable importance analysis for ESKD and death informed the creation of a sequence of reduced models that utilized six key variables: age, serum albumin, bicarbonate, chloride, eGFR, and hemoglobin. The models demonstrated robust predictive performance, with a median concordance index of 84.84% for ESKD and 84.1% for death. The median integrated Brier scores were 0.03 for ESKD and 0.038 for death across all landmark times. External validation with 8,729 patients confirmed these results.

We successfully developed and validated a dynamic survival prediction model using common longitudinal clinicopathological data. This model predicts ESKD with death as a competing risk and aims to assist clinicians in dialysis planning for patients with CKD.