Joint modelling of longitudinal ordinal and multi-state data

Abstract

Joint modeling of longitudinal and survival data is increasingly used in biomedical studies. However, existing joint models are not applicable to model the longitudinal ordinal responses with non-ignorable missing values caused by the occurrence of events in a multi-state process. In this article, we introduce a joint model for longitudinal ordinal measurements and multi-state data. Our proposed joint model consists of two sub-models: a proportional odds sub-model for longitudinal ordinal measurements and a multi-state sub-model with transition-specific proportional hazards for times of transitions between different health states, both linked by shared random effects. The model parameters were estimated employing the maximum likelihood method for a piecewise constant baseline hazard function. The proposed joint model is evaluated in a simulation study and, as an illustration, it is fitted to real data from people with human immunodeficiency virus.

Keywords

Joint model longitudinal ordinal data multi-state shared random effects human immunodeficiency virus/acquired immunodeficiency syndrome

Get full access to this article

View all access options for this article.

References

Little

Rubin

. Statistical analysis with missing data. New York: John Wiley, 2002.

Varin

Czado

. A mixed autoregressive probit model for ordinal longitudinal data. Biostatistics 2010; 11: 127–138.

Zhang

. Bayesian analysis of multivariate longitudinal ordinal data using multiple multivariate probit models. Am J Theoretic Appl Stat 2023; 12: 1–12.

Armero

Forné

Rué

, et al. Bayesian joint ordinal and survival modeling for breast cancer risk assessment. Stat Med 2016; 35: 5267–5282.

World Health

Organization

. Interim WHO clinical staging of HVI/AIDS and HIV/AIDS case definitions for surveillance: African Region. Switzerland: World Health Organization, 2005.

Hamidi

Tapak

Poorolajal

, et al. Identifying risk factors for progression to AIDS and mortality post-HIV infection using illness-death multistate model. Clin Epidemiol Glob Health 2017; 5: 163–168.

Tapak

Kosorok

Sadeghifar

, et al. Multistate recursively imputed survival trees for time-to-event data analysis: an application to AIDS and mortality post-HIV infection data. BMC Med Res Methodol 2018; 18: 1–12.

Ferrer

Rondeau

Dignam

, et al. Joint modelling of longitudinal and multi-state processes: application to clinical progressions in prostate cancer. Stat Med 2016; 35: 3933–3948.

Dantan

Joly

Dartigues

J-F

, et al. Joint model with latent state for longitudinal and multistate data. Biostatistics 2011; 12: 723–736.

10.

Mwanyekange

Mwalili

Ngesa

. Bayesian joint models for longitudinal and multi-state survival data. Int J Stat Probab 2019; 8: 34

11.

Furgal

. Bayesian Models for Joint Longitudinal and Multi-State Survival Data 2021.

12.

Alafchi

Mahjub

Tapak

, et al. Two-stage joint model for multivariate longitudinal and multistate processes, with application to renal transplantation data. J Probab Stat 2021; 2021: 1–10.

13.

Alam

Maity

Sinha

, et al. Joint modeling of longitudinal continuous, longitudinal ordinal, and time-to-event outcomes. Lifetime Data Anal 2021; 27: 64–90.

14.

Baghfalaki

Kalantari

Ganjali

, et al. Bayesian joint modeling of ordinal longitudinal measurements and competing risks survival data for analysing Tehran lipid and glucose study. J Biopharm Stat 2020; 30: 689–703.

15.

Elashoff

, et al. Joint modeling of longitudinal ordinal data and competing risks survival times and analysis of the NINDS rt-PA stroke trial. Stat Med 2010; 29: 546–557.

16.

Elashoff

, et al. Joint analysis of bivariate longitudinal ordinal outcomes and competing risks survival times with nonparametric distributions for random effects. Stat Med 2012; 31: 1707–1721.

17.

Hogan

Laird

. Model-based approaches to analysing incomplete longitudinal and failure time data. Stat Med 1997; 16: 259–272.

18.

Schluchter

. Methods for the analysis of informatively censored longitudinal data. Stat Med 1992; 11: 1861–1870.

19.

Weiser

. Methods for Multivariate Quadrature 2016 [Available from: https://cran.r-project.org/web/packages/mvQuad/index.html.

20.

Beger

. Create Panels of Independent States 2021 [Available from: https://cran.r-project.org/web/packages/states/index.html.

21.

Archer

Hedeker

Nordgren

, et al. mixor: an R package for longitudinal and clustered ordinal response modeling. 2015.

22.

de Wreede

Fiocco

Putter

. Mstate: an R package for the analysis of competing risks and multi-state models. J Stat Softw 2011; 38: 1–30.

23.

Liu

, et al. Analysis of longitudinal and survival data: joint modeling, inference methods, and issues. J Probab Stat 2012; 1: 640153.

24.

Beyersmann

Allignol

Schumacher

. Competing risks and multistate models with R. New York, Dordrecht, Heidelberg, London: Springer Science & Business Media, 2011.

25.

Crowther

Lambert

. Simulating biologically plausible complex survival data. Stat Med 2013; 32: 4118–4134.

26.

Brent

Richard

. Algorithms for minimization without derivatives. Mineola, New York: Dover Publications, 1973.

27.

Edathodu

Ali

Alrajhi

. CD4 validation for the World Health Organization classification and clinical staging of HIV/AIDS in a developing country. Int J Infect Dis 2009; 13: 243–246.

28.

Alafchi

Tapak

Doosti

, et al. A two-step algorithm to estimate variable importance for multi-state data: an application to COVID-19. Comput Model Eng Sci 2023; 135: 2046–2063.

29.

Abbring

Van den Berg

. The identifiability of the mixed proportional hazards competing risks model. J R Stat Soc Ser B Stat Methodol 2003; 65: 701–710.

30.

Cullagh

. Regression models for ordinal data (with discussion). J Roy Statist Soc, B 1980; 42: 109–142.

31.

Peterson

Harrell Jr

. Partial proportional odds models for ordinal response variables. J R Stat Soc: Ser C (Appl Stat) 1990; 39: 205–217.

32.

Asanjarani

Liquet

Nazarathy

. Estimation of semi-Markov multi-state models: a comparison of the sojourn times and transition intensities approaches. Int J Biostat 2021; 18: 243–262.

33.

Gillaizeau

Giral

Dantan

, et al. Multi-state analysis of kidney transplant recipients outcome: a semi-Markov model for studying the role of pre-transplant sensitization against Angiotensin II Type 1 receptor. Journal de la Société Française de Statistique 2014; 155: 117–133.

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.02 MB

0.04 MB