Modeling Longitudinal Data with Nonparametric Multiplicative Random Effects Jointly with Survival Data
- 28 June 2008
- journal article
- Published by Oxford University Press (OUP) in Biometrics
- Vol. 64 (2), 546-556
- https://doi.org/10.1111/j.1541-0420.2007.00896.x
Abstract
Summary In clinical studies, longitudinal biomarkers are often used to monitor disease progression and failure time. Joint modeling of longitudinal and survival data has certain advantages and has emerged as an effective way to mutually enhance information. Typically, a parametric longitudinal model is assumed to facilitate the likelihood approach. However, the choice of a proper parametric model turns out to be more elusive than models for standard longitudinal studies in which no survival endpoint occurs. In this article, we propose a nonparametric multiplicative random effects model for the longitudinal process, which has many applications and leads to a flexible yet parsimonious nonparametric random effects model. A proportional hazards model is then used to link the biomarkers and event time. We use B‐splines to represent the nonparametric longitudinal process, and select the number of knots and degrees based on a version of the Akaike information criterion (AIC). Unknown model parameters are estimated through maximizing the observed joint likelihood, which is iteratively maximized by the Monte Carlo Expectation Maximization (MCEM) algorithm. Due to the simplicity of the model structure, the proposed approach has good numerical stability and compares well with the competing parametric longitudinal approaches. The new approach is illustrated with primary biliary cirrhosis (PBC) data, aiming to capture nonlinear patterns of serum bilirubin time courses and their relationship with survival time of PBC patients.Keywords
This publication has 27 references indexed in Scilit:
- Dynamic Analysis of Recurrent Event Data Using the Additive Hazard ModelBiometrical Journal, 2006
- Functional Data Analysis for Sparse Longitudinal DataJournal of the American Statistical Association, 2005
- Functional Data AnalysisPublished by Wiley ,2005
- Ascent-Based Monte Carlo Expectation– MaximizationJournal of the Royal Statistical Society Series B: Statistical Methodology, 2005
- Modeling the Relationship of Survival to Longitudinal Data Measured with Error. Applications to Survival and CD4 Counts in Patients with AIDSJournal of the American Statistical Association, 1995
- Modeling the Relationship of Survival to Longitudinal Data Measured with Error. Applications to Survival and CD4 Counts in Patients with AIDSJournal of the American Statistical Association, 1995
- Primary biliary cirrhosis: Prediction of short-term survival based on repeated patient visitsHepatology, 1994
- Missing Data, Imputation, and the BootstrapJournal of the American Statistical Association, 1994
- Missing Data, Imputation, and the BootstrapJournal of the American Statistical Association, 1994
- An optimal selection of regression variablesBiometrika, 1981