Joint modeling of multiple repeated measures and survival data using multidimensional latent trait linear mixed model

Abstract

Impairment caused by Amyotrophic lateral sclerosis (ALS) is multidimensional (e.g. bulbar, fine motor, gross motor) and progressive. Its multidimensional nature precludes a single outcome to measure disease progression. Clinical trials of ALS use multiple 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 develop a joint model consisting of a multidimensional latent trait linear mixed model (MLTLMM) for the multiple longitudinal outcomes, and a proportional hazards model with piecewise constant baseline hazard for the event time data. Shared random effects are used to link together two models. The model inference is conducted using a Bayesian framework via Markov chain Monte Carlo simulation implemented in Stan language. Our proposed model is evaluated by simulation studies and is applied to the Ceftriaxone study, a motivating clinical trial assessing the effect of ceftriaxone on ALS patients.

Keywords

Amyotrophic lateral sclerosis informative dropout longitudinal data Markov chain Monte Carlo mixed model

Get full access to this article

View all access options for this article.

References

Wijesekera

Leigh

. Amyotrophic lateral sclerosis. Orphanet J Rare Dis 2009; 4: 3.

Berry

Shefner

Conwit

, et al. Design and initial results of a multi-phase randomized trial of ceftriaxone in amyotrophic lateral sclerosis. PLoS One 2013; 8: e61177.

Prize4Life. Pooled Resource Open-Access ALS Clinical Trials Database; 2013. Data Sets, https://nctu.partners.org/ProACT (accessed 12 March 2018).

Abe

Itoyama

Sobue

, et al. Confirmatory double-blind, parallel-group, placebo-controlled study of efficacy and safety of edaravone (MCI-186) in amyotrophic lateral sclerosis patients. Amyotroph Lateral Scler Frontotemporal Degener 2014; 15: 610–617.

Bacci

Staniewska

Coyne

, et al. Item response theory analysis of the amyotrophic lateral sclerosis functional rating scale-revised in the pooled resource open-access ALS clinical trials database. Amyotroph Lateral Scler Frontotemporal Degener 2016; 17: 157–167.

The ACTS Study Group. The amyotrophic lateral sclerosis functional rating scale. Assessment of activities of daily living in patients with amyotrophic lateral sclerosis. Arch Neurol 1996; 53: 141–147.

Clarke

Hickey

O'Boyle

, et al. Assessing individual quality of life in amyotrophic lateral sclerosis. Qual Life Res 2001; 10: 149–158.

Coco

Marchese

La Bella

, et al. The amyotrophic lateral sclerosis functional rating scale predicts survival time in amyotrophic lateral sclerosis patients on invasive mechanical ventilation. Chest 2007; 132: 64–69.

Bedlack

Vaughan

Wicks

, et al. How common are ALS plateaus and reversals?. Neurology 2016; 86: 808–812.

10.

Balsis

Unger

Benge

, et al. Gaining precision on the Alzheimer's disease assessment scale-cognitive: a comparison of item response theory-based scores and total scores. Alzheimer's Dementia 2012; 8: 288–294.

11.

Grimby

Tennant

Tesio

. The use of raw scores from ordinal scales: time to end malpractice?. J Rehabil Med 2012; 44: 97–98.

12.

Andrich

. Rasch models for measurement, Newbury Park, USA: Sage, 1988. .

13.

Samejima

. Estimation of latent ability using a response pattern of graded scores. ETS Res Bull Series 1968; 1968: i–169.

14.

Fox

Glas

. Bayesian estimation of a multilevel IRT model using Gibbs sampling. Psychometrika 2001; 66: 271–288.

15.

Glas

Geerlings

van de Laar

, et al. Analysis of longitudinal randomized clinical trials using item response models. Contemp Clin Trials 2009; 30: 158–170.

16.

Luo

Wang

. Bayesian hierarchical model for multiple repeated measures and survival data: an application to Parkinson's disease. Stat Med 2014; 33: 4279–4291.

17.

Wang

Luo

. Multidimensional latent trait linear mixed model: An application in clinical studies with multivariate longitudinal outcomes. Stat Med 2017; 36: 3244–3256.

18.

Henderson

Diggle

Dobson

. Joint modelling of longitudinal measurements and event time data. Biostatistics 2000; 1: 465–480.

19.

Faucett

Thomas

. Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. Stat Med 1996; 15: 1663–1685.

20.

Wulfsohn

Tsiatis

. A joint model for survival and longitudinal data measured with error. Biometrics 1997; 53: 330–339.

21.

Rizopoulos

. Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data. Biometrics 2011; 67: 819–829.

22.

Wang

Luo

. Dynamic prediction for multiple repeated measures and event time data: an application to Parkinson's disease. Ann Appl Stat 2017; 11: 1787.

23.

Lawless

Zhan

. Analysis of interval-grouped recurrent-event data using piecewise constant rate functions. Can J Stat 1998; 26: 549–565.

24.

Feng

Wolfe

Port

. Frailty survival model analysis of the National Deceased Donor Kidney Transplant Dataset using Poisson variance structures. J Am Stat Assoc 2005; 100: 728–735.

25.

Liu

Huang

. Joint analysis of correlated repeated measures and recurrent events processes in the presence of death, with application to a study on acquired immune deficiency syndrome. J Royal Stat Soc: Ser C 2009; 58: 65–81.

26.

Stan Development Team. Stan modeling language: user's guide and reference manual, version 2.17.0, 2017, http://mc-stan.org/ (accessed 14 July 2017).

27.

Gelman

Carlin

Stern

, et al. Bayesian data analysis, Boca Raton, USA: CRC Press, 2013. .

28.

Celeux

Forbes

Robert

, et al. Deviance information criteria for missing data models. Bayesian Analys 2006; 1: 651–673.

29.

Carlin

Louis

. Bayesian methods for data analysis, Boca Raton, USA: CRC Press, 2009. .

30.

Geskus

. Which individuals make dropout informative?. Stat Meth Med Res 2014; 23(1): 91–106.

31.

Luo

. Joint modeling of multivariate longitudinal measurements and survival data with applications to Parkinson's disease. Stat Meth Med Res 2016; 25: 1346–1358.

32.

Brown

. Assessing the association between trends in a biomarker and risk of event with an application in pediatric HIV/AIDS. Ann Appl Stat 2009; 3: 1163–1182.

33.

Rizopoulos

Hatfield

Carlin

, et al. Combining dynamic predictions from joint models for longitudinal and time-to-event data using Bayesian model averaging. J Am Stat Assoc 2014; 109: 1385–1397.

34.

Harrell

. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis, New York, USA: Springer, 2015. .

35.

Jacqmin-Gadda

Sibillot

Proust

, et al. Robustness of the linear mixed model to misspecified error distribution. Computat Stat Data Analys 2007; 51: 5142–5154.

36.

Rizopoulos

Verbeke

Molenberghs

. Shared parameter models under random effects misspecification. Biometrika 2008; 95: 63–74.

37.

McCulloch

Neuhaus

, et al. Misspecifying the shape of a random effects distribution: Why getting it wrong may not matter. Stat Sci 2011; 26: 388–402.

38.

Escobar

. Estimating normal means with a Dirichlet process prior. J Am Stat Assoc 1994; 89: 268–277.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.08 MB