A fast EM algorithm for fitting joint models of a binary response and multiple longitudinal covariates subject to detection limits

Joint modeling techniques have become a popular strategy for studying the association between a response and one or more longitudinal covariates. Motivated by the GenIMS study, where it is of interest to model the event of survival using censored longitudinal biomarkers, a joint model is proposed for describing the relationship between a binary outcome and multiple longitudinal covariates subject to detection limits. A fast, approximate EM algorithm is developed that reduces the dimension of integration in the E-step of the algorithm to one, regardless of the number of random effects in the joint model. Numerical studies demonstrate that the proposed approximate EM algorithm leads to satisfactory parameter and variance estimates in situations with and without censoring on the longitudinal covariates. The approximate EM algorithm is applied to analyze the GenIMS data set.

[1]  Raymond J. Carroll,et al.  Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process , 1988 .

[2]  Marie Davidian,et al.  Conditional Estimation for Generalized Linear Models When Covariates Are Subject‐Specific Parameters in a Mixed Model for Longitudinal Measurements , 2004, Biometrics.

[3]  Anthony N. Pettitt,et al.  Censored observations, repeated measures and mixed effects models: An approach using the EM algorithm and normal errors , 1986 .

[4]  F. James,et al.  Monte Carlo theory and practice , 1980 .

[5]  New York Dover,et al.  ON THE CONVERGENCE PROPERTIES OF THE EM ALGORITHM , 1983 .

[6]  Daowen Zhang,et al.  Likelihood and pseudo-likelihood methods for semiparametric joint models for a primary endpoint and longitudinal data , 2007, Comput. Stat. Data Anal..

[7]  Marie Davidian,et al.  Mixed model analysis of censored longitudinal data with flexible random-effects density. , 2012, Biostatistics.

[8]  Cécile Proust-Lima,et al.  Joint modelling of multivariate longitudinal outcomes and a time-to-event: A nonlinear latent class approach , 2009, Comput. Stat. Data Anal..

[9]  B. Reboussin,et al.  Mixed Effects Logistic Regression Models for Longitudinal Ordinal Functional Response Data with Multiple‐Cause Drop‐Out from the Longitudinal Study of Aging , 2000, Biometrics.

[10]  Xihong Lin,et al.  A penalized likelihood approach to joint modeling of longitudinal measurements and time-to-event data , 2008 .

[11]  Gina D'Angelo,et al.  An index approach for the Cox model with left censored covariates , 2008, Statistics in medicine.

[12]  J R Landis,et al.  Mixed effects logistic regression models for longitudinal binary response data with informative drop-out. , 1998, Biometrics.

[13]  Jessica Brooks,et al.  Understanding the Inflammatory Cytokine Response in Pneumonia and Sepsis , 2008 .

[14]  Hulin Wu,et al.  Joint inference for nonlinear mixed-effects models and time to event at the presence of missing data. , 2007, Biostatistics.

[15]  Geert Molenberghs,et al.  Analysis of non‐ignorable missing and left‐censored longitudinal data using a weighted random effects tobit model , 2011, Statistics in medicine.

[16]  Lisa A. Weissfeld,et al.  Joint modeling of censored longitudinal and event time data , 2013 .

[17]  Ryan C. May,et al.  Estimation Methods for Data Subject to Detection Limits , 2012 .

[18]  Naisyin Wang,et al.  Joint Models for a Primary Endpoint and Multiple Longitudinal Covariate Processes , 2007, Biometrics.

[19]  Geert Verbeke,et al.  Fully exponential Laplace approximations for the joint modelling of survival and longitudinal data , 2009 .

[20]  C Y Wang,et al.  Regression Analysis When Covariates Are Regression Parameters of a Random Effects Model for Observed Longitudinal Measurements , 2000, Biometrics.

[21]  Dimitris Rizopoulos.  Event Time Event Time , 2012 .

[22]  Hsun-Chih Kuo,et al.  The joint model of the logistic model and linear random effect model - An application to predict orthostatic hypertension for subacute stroke patients , 2011, Comput. Stat. Data Anal..

[23]  Ying Yuan,et al.  Mixed‐Effect Hybrid Models for Longitudinal Data with Nonignorable Dropout , 2009, Biometrics.

[24]  Mohsen Mohammadzadeh,et al.  Asymptotic normality of posterior distributions for generalized linear mixed models , 2012, J. Multivar. Anal..

[25]  Brian D. M. Tom,et al.  Bias in 2-part mixed models for longitudinal semicontinuous data , 2009, Biostatistics.

[26]  Andrew D. Martin,et al.  MCMCpack: Markov chain Monte Carlo in R , 2011 .

[27]  Abdus Sattar,et al.  A Parametric Survival Model When a Covariate is Subject to Left-Censoring. , 2012, Journal of biometrics & biostatistics.

[28]  Dennis R. Helsel,et al.  Statistics for Censored Environmental Data Using Minitab and R , 2012 .

[29]  G. C. Wei,et al.  A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms , 1990 .

[30]  JolyPierre,et al.  Joint modelling of multivariate longitudinal outcomes and a time-to-event , 2009 .

[31]  M. Wulfsohn,et al.  A joint model for survival and longitudinal data measured with error. , 1997, Biometrics.

[32]  Dimitris Rizopoulos,et al.  Joint Models for Longitudinal and Time-to-Event Data: With Applications in R , 2012 .

[33]  Robert H. Lyles,et al.  Random regression models for human immunodeficiency virus ribonucleic acid data subject to left censoring and informative drop‐outs , 2000 .

[34]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[35]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[36]  Dimitris Rizopoulos,et al.  Fast fitting of joint models for longitudinal and event time data using a pseudo-adaptive Gaussian quadrature rule , 2012, Comput. Stat. Data Anal..

[37]  Charles E McCulloch,et al.  Maximum likelihood estimation in the joint analysis of time‐to‐event and multiple longitudinal variables , 2002, Statistics in medicine.

[38]  Paul W. Bernhardt,et al.  Flexible modeling of survival data with covariates subject to detection limits via multiple imputation , 2014, Comput. Stat. Data Anal..

[39]  Lang Wu,et al.  A Joint Model for Nonlinear Mixed-Effects Models With Censoring and Covariates Measured With Error, With Application to AIDS Studies , 2002 .

[40]  John A Kellum,et al.  Understanding the inflammatory cytokine response in pneumonia and sepsis: results of the Genetic and Inflammatory Markers of Sepsis (GenIMS) Study. , 2007, Archives of internal medicine.

[41]  Paul W. Bernhardt,et al.  Statistical Methods for Generalized Linear Models with Covariates Subject to Detection Limits , 2015, Statistics in biosciences.

[42]  R. Fernando,et al.  Mixed-model analysis of a censored normal distribution with reference to animal breeding. , 1987, Biometrics.