More than 720,000 patients with end-stage renal disease in the US require life-sustaining dialysis treatment that is predominantly received at local dialysis facilities. In this population of typically older patients with a high morbidity burden, hospitalization is frequent at a rate of about twice per patient-year. Aside from frequent hospitalizations, which is a major source of mortality, overall mortality in dialysis patients is higher than other comparable populations, including Medicare patients with cancer. Thus, understanding patient- and facility-level risk factors that jointly contribute to longitudinal hospitalizations and mortality is of interest. Towards this objective, we propose a novel methodology to jointly model hospitalization, a binary longitudinal outcome, and survival, based on multilevel data from the United States Renal Data System (USRDS), with repeated observations over time nested in patients and patients nested in dialysis facilities. In our approach, the outcomes are modeled through a common set of multilevel random effects. In order to accommodate the USRDS data structure, we depart from the literature on joint modeling of longitudinal and survival data by including multilevel random effects and multilevel covariates, at both the patient and facility levels. An approximate EM algorithm is developed for estimation where fully exponential Laplace approximations are utilized to address computational challenges.

Prof. Esra Kurum