Dynamic predictions from longitudinal CD4 count measures and time to death of HIV/AIDS patients using a Bayesian joint model

A Bayesian joint modeling approach to dynamic prediction of HIV progression and mortality allows individualized predictions to be made for HIV patients, based on monitoring of their CD4 counts. This study aims to provide predictions of patient-specific trajectories of HIV disease progression and survival. Longitudinal data on 254 HIV/AIDS patients who received ART between 2009 and 2014, and who had at least one CD4 count observed, were employed in a Bayesian joint model of disease progression. Different forms of association structure that relate the longitudinal CD4 biomarker and time to death were assessed; and predictions were averaged over the different models using Bayesian model averaging. The individual follow-up times ranged from 1 to 120 months, with a median of 22 months and IQR 7–39 months. The estimates of the association structure parameters from two of the three models considered indicated that the HIV mortality hazard at any time point is associated with the rate of change in the underlying value of the CD4 count. Model averaging the dynamic predictions resulted in only one of the hypothesized association structures having non-zero weight in almost all time points for each individual, with the exception of twelve patients, for whom other association structures were preferred at a few time points. The predictions were found to be different when we averaged them over models than when we derived them from the highest posterior weight model alone. The model with highest posterior weight for almost all time points for each individual gave an estimate of the association parameter of −0.02 implying that for a unit increase in the CD4 count, the hazard of HIV mortality decreases by a factor (hazard ratio) of 0.98. Functional status and alcohol intake are important contributing factors that affect the mean square root of CD4 measurements.