A Bayesian approach to joint analysis of multivariate longitudinal data and parametric accelerated failure time

Impairment caused by Parkinson's disease (PD) is multidimensional (e.g., sensoria, functions, and cognition) and progressive. Its multidimensional nature precludes a single outcome to measure disease progression. Clinical trials of PD use multiple categorical and continuous longitudinal outcomes to assess the treatment effects on overall improvement. A terminal event such as death or dropout can stop the follow‐up process. Moreover, the time to the terminal event may be dependent on the multivariate longitudinal measurements. In this article, we consider a joint random‐effects model for the correlated outcomes. A multilevel item response theory model is used for the multivariate longitudinal outcomes and a parametric accelerated failure time model is used for the failure time because of the violation of proportional hazard assumption. These two models are linked via random effects. The Bayesian inference via MCMC is implemented in ‘BUGS’ language. Our proposed method is evaluated by a simulation study and is applied to DATATOP study, a motivating clinical trial to determine if deprenyl slows the progression of PD. © 2013 The authors. Statistics in Medicine published by John Wiley & Sons, Ltd.

[1]  S. Luo,et al.  Joint modeling of multivariate longitudinal measurements and survival data with applications to Parkinson’s disease , 2016, Statistical methods in medical research.

[2]  A. Dreher.  Modeling Survival Data Extending The Cox Model , 2016 .

[3]  Cécile Proust-Lima,et al.  Joint latent class models for longitudinal and time-to-event data: A review , 2014, Statistical methods in medical research.

[4]  Pan Tong,et al.  integIRTy: a method to identify genes altered in cancer by accounting for multiple mechanisms of regulation using item response theory , 2012, Bioinform..

[5]  Wen Chung 王文中 Wang,et al.  The Generalized Multilevel Facets Model for Longitudinal Data , 2012 .

[6]  S Bandyopadhyay,et al.  A review of multivariate longitudinal data analysis , 2011, Statistical methods in medical research.

[7]  Inês Sousa,et al.  A review on joint modelling of longitudinal measurements and time-to-event , 2011 .

[8]  S. McKay Curtis,et al.  BUGS Code for Item Response Theory , 2010 .

[9]  J. Fox.  Bayesian Item Response Modeling: Theory and Applications , 2010 .

[10]  Peter Müller,et al.  A Bayesian Semiparametric Survival Model with Longitudinal Markers , 2010, Biometrics.

[11]  Robert E Weiss,et al.  Bayesian methods for data analysis. , 2010, American journal of ophthalmology.

[12]  J. Varni,et al.  Bayesian Estimation of Graded Response Multilevel Models Using Gibbs Sampling: Formulation and Illustration , 2010 .

[13]  Jean-Paul Fox,et al.  Bayesian Item Response Modeling , 2010 .

[14]  P. Huang,et al.  Using global statistical tests in long‐term Parkinson's disease clinical trials , 2009, Movement disorders : official journal of the Movement Disorder Society.

[15]  Hanneke Geerlings,et al.  Analysis of longitudinal randomized clinical trials using item response models. , 2009, Contemporary clinical trials.

[16]  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..

[17]  J. Jankovic,et al.  A randomized, double‐blind, placebo‐controlled, delayed start study to assess rasagiline as a disease modifying therapy in Parkinson's disease (the ADAGIO study): Rationale, design, and baseline characteristics , 2008, Movement disorders : official journal of the Movement Disorder Society.

[18]  Cécile Proust-Lima,et al.  The International Journal of Biostatistics Pattern Mixture Models and Latent Class Models for the Analysis of Multivariate Longitudinal Data with Informative Dropouts , 2011 .

[19]  Peng Huang,et al.  A rank‐based sample size method for multiple outcomes in clinical trials , 2008, Statistics in medicine.

[20]  Peng Huang,et al.  The Use of a Global Statistical Approach for the Design and Data Analysis of Clinical Trials with Multiple Primary Outcomes , 2008 .

[21]  S. Rabe-Hesketh,et al.  Latent Variable Modelling: A Survey * , 2007 .

[22]  Jean-Paul Fox,et al.  Multilevel IRT Modeling in Practice with the Package mlirt , 2007 .

[23]  Jean-Paul Fox,et al.  Modeling of Responses and Response Times with the Package cirt , 2007 .

[24]  I. Shoulson.  CEP-1347 Treatment Fails to Favorably Modify the Progression of Parkinson's Disease (PRECEPT) Study , 2006, Neurology.

[25]  G. Molenberghs,et al.  Models for Discrete Longitudinal Data , 2005 .

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

[27]  Stuart Lipsitz,et al.  Adjusting O'Brien's Test to Control Type I Error for the Generalized Nonparametric Behrens–Fisher Problem , 2005, Biometrics.

[28]  Stan Lipovetsky,et al.  Generalized Latent Variable Modeling: Multilevel,Longitudinal, and Structural Equation Models , 2005, Technometrics.

[29]  J-P Fox,et al.  Multilevel IRT using dichotomous and polytomous response data. , 2005, The British journal of mathematical and statistical psychology.

[30]  Svend Kreiner,et al.  Latent Regression in Loglinear Rasch Models , 2004 .

[31]  C. Tanner,et al.  Levodopa and the progression of Parkinson's disease. , 2004, The New England journal of medicine.

[32]  Jean-Paul Fox,et al.  Applications of Multilevel IRT Modeling , 2004 .

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

[34]  C. Parvin.  An Introduction to Multivariate Statistical Analysis, 3rd ed. T.W. Anderson. Hoboken, NJ: John Wiley & Sons, 2003, 742 pp., $99.95, hardcover. ISBN 0-471-36091-0. , 2004 .

[35]  Garrett M. Fitzmaurice,et al.  Analysis of longitudinal multiple‐source binary data using generalized estimating equations , 2004 .

[36]  E. B. Andersen,et al.  Latent Regression Analysis based on the Rating Scale Model , 2004 .

[37]  Anastasios A. Tsiatis,et al.  Joint Modeling of Longitudinal and Time-to-Event Data : An Overview , 2004 .

[38]  J. Nutt,et al.  A controlled trial of rasagiline in early Parkinson disease: the TEMPO Study. , 2002, Archives of neurology.

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

[40]  C. McCulloch,et al.  Latent Class Models for Joint Analysis of Longitudinal Biomarker and Event Process Data , 2002 .

[41]  J. Douglas,et al.  Item response models for joint analysis of quality of life and survival , 2002, Statistics in medicine.

[42]  S J Pocock,et al.  Impact of missing data due to drop‐outs on estimators for rates of change in longitudinal studies: a simulation study , 2001, Statistics in medicine.

[43]  Kimberly S. Maier.  A Rasch Hierarchical Measurement Model , 2001 .

[44]  J. Fox,et al.  Bayesian estimation of a multilevel IRT model using gibbs sampling , 2001 .

[45]  Akihito Kamata,et al.  Item Analysis by the Hierarchical Generalized Linear Model. , 2001 .

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

[47]  S. Pillemer,et al.  Global statistical tests for comparing multiple outcomes in rheumatoid arthritis trials. MIRA Trial Group. , 1999, Arthritis and rheumatism.

[48]  Ira Shoulson,et al.  DATATOP: A decade of neuroprotective inquiry , 1998 .

[49]  I. Shoulson.  DATATOP: a decade of neuroprotective inquiry. Parkinson Study Group. Deprenyl And Tocopherol Antioxidative Therapy Of Parkinsonism. , 1998, Annals of neurology.

[50]  Raymond J. Adams,et al.  Multilevel Item Response Models: An Approach to Errors in Variables Regression , 1997 .

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

[52]  N M Laird,et al.  Mixture models for the joint distribution of repeated measures and event times. , 1997, Statistics in medicine.

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

[54]  N L Geller,et al.  Use of a global test for multiple outcomes in stroke trials with application to the National Institute of Neurological Disorders and Stroke t-PA Stroke Trial. , 1996, Stroke.

[55]  D. Thomas,et al.  Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. , 1996, Statistics in medicine.

[56]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[57]  Robert J. Mislevy,et al.  Estimation of Latent Group Effects , 1985 .

[58]  P. O'Brien.  Procedures for comparing samples with multiple endpoints. , 1984, Biometrics.

[59]  W. Cleveland.  Robust Locally Weighted Regression and Smoothing Scatterplots , 1979 .

[60]  Z. Schwab,et al.  Projection technique for evaluating surgery in Parkinson’s disease , 1969 .

[61]  T. W. Anderson,et al.  An Introduction to Multivariate Statistical Analysis , 1959 .