Joint analysis of left-censored longitudinal biomarker and binary outcome via latent class modeling

Joint latent class modeling is an appealing approach for evaluating the association between a longitudinal biomarker and clinical outcome when the study population is heterogeneous. The link between the biomarker trajectory and the risk of event is reflected by the latent classes, which accommodate the underlying population heterogeneity. The estimation of joint latent class models may be complicated by the censored data in the biomarker measurements due to detection limits. We propose a modified likelihood function under the parametric assumption of biomarker distribution and develop a Monte Carlo expectation-maximization algorithm for joint analysis of a biomarker and a binary outcome. We conduct simulation studies to demonstrate the satisfactory performance of our Monte Carlo expectation-maximization algorithm and the superiority of our method to the naive imputation method for handling censored biomarker data. In addition, we apply our method to the Genetic and Inflammatory Markers of Sepsis study to investigate the role of inflammatory biomarker profile in predicting 90-day mortality for patients hospitalized with community-acquired pneumonia.