Joint modeling of multivariate longitudinal mixed measurements and time to event data using a Bayesian approach

In many longitudinal studies multiple characteristics of each individual, along with time to occurrence of an event of interest, are often collected. In such data set, some of the correlated characteristics may be discrete and some of them may be continuous. In this paper, a joint model for analysing multivariate longitudinal data comprising mixed continuous and ordinal responses and a time to event variable is proposed. We model the association structure between longitudinal mixed data and time to event data using a multivariate zero-mean Gaussian process. For modeling discrete ordinal data we assume a continuous latent variable follows the logistic distribution and for continuous data a Gaussian mixed effects model is used. For the event time variable, an accelerated failure time model is considered under different distributional assumptions. For parameter estimation, a Bayesian approach using Markov Chain Monte Carlo is adopted. The performance of the proposed methods is illustrated using some simulation studies. A real data set is also analyzed, where different model structures are used. Model comparison is performed using a variety of statistical criteria.

[1]  R. Little,et al.  Maximum likelihood estimation for mixed continuous and categorical data with missing values , 1985 .

[2]  Ingram Olkin,et al.  Multivariate Correlation Models with Mixed Discrete and Continuous Variables , 1961 .

[3]  Xin Huang,et al.  A joint model of longitudinal and competing risks survival data with heterogeneous random effects and outlying longitudinal measurements. , 2010, Statistics and its interface.

[4]  J. Ibrahim,et al.  A Bayesian semiparametric joint hierarchical model for longitudinal and survival data. , 2003, Biometrics.

[5]  Yan Wang,et al.  Jointly Modeling Longitudinal and Event Time Data With Application to Acquired Immunodeficiency Syndrome , 2001 .

[6]  Anastasios A. Tsiatis,et al.  A semiparametric estimator for the proportional hazards model with longitudinal covariates measured with error , 2001 .

[7]  M Ganjali,et al.  A MODEL FOR MIXED CONTINUOUS AND DISCRETE RESPONSES WITH POSSIBILITY OF MISSING RESPONSES , 2003 .

[8]  B P Carlin,et al.  Response of CD4 lymphocytes and clinical consequences of treatment using ddI or ddC in patients with advanced HIV infection. , 1996, Journal of acquired immune deficiency syndromes and human retrovirology : official publication of the International Retrovirology Association.

[9]  David E. Booth,et al.  Analysis of Incomplete Multivariate Data , 2000, Technometrics.

[10]  Rupert G. Miller,et al.  Survival Analysis , 2022, The SAGE Encyclopedia of Research Design.

[11]  Chi-Hong Tseng,et al.  Joint analysis of bivariate longitudinal ordinal outcomes and competing risks survival times with nonparametric distributions for random effects , 2012, Statistics in medicine.

[12]  D. Dunson,et al.  Bayesian latent variable models for mixed discrete outcomes. , 2005, Biostatistics.

[13]  Joseph G Ibrahim,et al.  Joint Models for Multivariate Longitudinal and Multivariate Survival Data , 2006, Biometrics.

[14]  Jeremy M. G. Taylor,et al.  A Stochastic Model for Analysis of Longitudinal AIDS Data , 1994 .

[15]  J. Ware,et al.  Random-effects models for longitudinal data. , 1982, Biometrics.

[16]  Geert Verbeke,et al.  A Two‐Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with Excess Zeros , 2008, Biometrics.

[17]  Juvêncio Santos Nobre,et al.  Residual Analysis for Linear Mixed Models , 2007, Biometrical journal. Biometrische Zeitschrift.

[18]  J. Ibrahim,et al.  A Flexible B‐Spline Model for Multiple Longitudinal Biomarkers and Survival , 2005, Biometrics.

[19]  Jayant V. Deshpande,et al.  Life-Time Data - Statistical Models and Methods , 2006, Series on Quality, Reliability and Engineering Statistics.

[20]  Dimitris Rizopoulos,et al.  A Bayesian semiparametric multivariate joint model for multiple longitudinal outcomes and a time‐to‐event , 2011, Statistics in medicine.

[21]  Joseph G. Ibrahim,et al.  A new joint model for longitudinal and survival data with a cure fraction , 2004 .

[22]  N M Laird,et al.  Regression models for mixed discrete and continuous responses with potentially missing values. , 1997, Biometrics.

[23]  Jie Zhang,et al.  Regression models for mixed Poisson and continuous longitudinal data , 2007, Statistics in medicine.

[24]  R. Hashemi,et al.  A Latent Process Model for Joint Modeling of Events and Marker , 2003, Lifetime data analysis.

[25]  Menggang Yu,et al.  JOINT LONGITUDINAL-SURVIVAL-CURE MODELS AND THEIR APPLICATION TO PROSTATE CANCER , 2004 .

[26]  M. Ganjali,et al.  Bayesian latent variable model for mixed continuous and ordinal responses with possibility of missing responses , 2011 .

[27]  Joseph G. Ibrahim,et al.  BAYESIAN METHODS FOR JOINT MODELING OF LONGITUDINAL AND SURVIVAL DATA WITH APPLICATIONS TO CANCER VACCINE TRIALS , 2004 .

[28]  M. Ganjali,et al.  A latent variable model for mixed continuous and ordinal responses with nonignorable missing responses: Assessing the local influence via covariance structure , 2010 .

[29]  R. Weiss,et al.  Residual plots for repeated measures. , 1992, Statistics in medicine.

[30]  David A. Schoenfeld,et al.  A Random-Effects Model for Multiple Characteristics with Possibly Missing Data , 1997 .

[31]  Geert Verbeke,et al.  Pairwise Fitting of Mixed Models for the Joint Modeling of Multivariate Longitudinal Profiles , 2006, Biometrics.

[32]  A. Gelfand,et al.  Bayesian Model Choice: Asymptotics and Exact Calculations , 1994 .

[33]  Dimitris Rizopoulos,et al.  JM: An R package for the joint modelling of longitudinal and time-to-event data , 2010 .

[34]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[35]  Amanda G Chetwynd,et al.  Joint modelling of repeated measurements and time‐to‐event outcomes: The fourth Armitage lecture , 2008, Statistics in medicine.

[36]  S. Oman.  Checking the assumptions in mixed-model analysis of variance: a residual analysis approach , 1995 .

[37]  E. Acuña,et al.  A Meta analysis study of outlier detection methods in classification , 2004 .

[38]  Gang Li,et al.  Robust Joint Modeling of Longitudinal Measurements and Competing Risks Failure Time Data , 2009, Biometrical journal. Biometrische Zeitschrift.

[39]  N M Laird,et al.  Model-based approaches to analysing incomplete longitudinal and failure time data. , 1997, Statistics in medicine.

[40]  Yueh-Yun Chi,et al.  Bayesian approaches to joint longitudinal and survival models accommodating both zero and nonzero cure fractions , 2007 .

[41]  Jason Roy,et al.  Analysis of Multivariate Longitudinal Outcomes With Nonignorable Dropouts and Missing Covariates , 2002 .

[42]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[43]  C. McCulloch.  Joint modelling of mixed outcome types using latent variables , 2008, Statistical methods in medical research.

[44]  M. May.  Bayesian Survival Analysis. , 2002 .

[45]  R Henderson,et al.  Joint modelling of longitudinal measurements and event time data. , 2000, Biostatistics.

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

[47]  F. Hsieh,et al.  Joint modelling of accelerated failure time and longitudinal data , 2005 .