There are few computational and methodological tools available for the analysis of general multi-state models with interval censoring. Here, we propose a general framework for parametric inference with interval censored multi-state data. Our framework can accommodate any parametric model for the transition times, and covariates may be included in various ways. We present a general method for constructing the likelihood, which we have implemented in a ready-to-use R package, smms, available on GitHub. The R package also computes the required high-dimensional integrals in an efficient manner. Further, we explore connections between our modelling framework and existing approaches: our models fall under the class of semi-Markovian multi-state models, but with a different, and sparser parameterisation than what is often seen. We illustrate our framework through a dataset monitoring heart transplant patients. Finally, we investigate the effect of some forms of misspecification of the model assumptions through simulations.
KalbfleischJLawlessJF. The analysis of panel data under a Markov assumption. J Am Stat Assoc1985; 80: 863–871.
2.
JacksonCH. Multi-state models for panel data: the MSM package for R. J Stat Softw2011; 38: 1–28.
3.
de WreedeLCFioccoMPutterH. Mstate: an R package for the analysis of competing risks and multi-state models. J Stat Softw2011; 38: 1–30.
4.
KròlASaint-PierreP. SemiMarkov: an R package for parametric estimation in multi-state semi-Markov models. J Stat Softw2015; 66: 1–16.
5.
JacksonCH. flexsurv: a platform for parametric survival modeling in R. J Stat Softw2016; 70: 1–38.
6.
KangMLagakosSW. Statistical methods for panel data from a semi-Markov process, with application to HPV. Biostatistics2006; 8: 252–264.
7.
FoucherYGiralMSoulillouJet al. A flexible semi-Markov model for interval-censored data and goodness-of-fit testing. Stat Methods Med Res2010; 19: 127–145.
8.
KapetanakisVMatthewsFEvan den HoutA. A semi-Markov model for stroke with piecewise-constant hazards in the presence of left, right and interval censoring. Stat Med2013; 32: 697–713.
9.
WeiSKryscioRJ. Semi-Markov models for interval censored transient cognitive states with back transitions and a competing risk. Stat Methods Med Res2016; 25: 2909–2924.
10.
AralisHBrookmeyerR. A stochastic estimation procedure for intermittently-observed semi-Markov multistate models with back transitions. Stat Methods Med Res2019; 28: 770–787.
11.
WilliamsJPStorlieCBTherneauTMet al. A Bayesian approach to multistate hidden Markov models: application to dementia progression. J Am Stat Assoc2020; 115: 16–31.
12.
GentlemanRLawlessJLindseyJet al. Multi-state Markov models for analysing incomplete disease history data with illustrations for HIV disease. Stat Med1994; 13: 805–821.
13.
JepsenPVilstrupHAndersenPK. The clinical course of cirrhosis: the importance of multistate models and competing risks analysis. Hepatology2015; 62: 292–302.
14.
CommengesD. Inference for multi-state models from interval-censored data. Stat Methods Med Res2002; 11: 167–182.
15.
AsanjaraniALiquetBNazarathyY. Estimation of semi-Markov multi-state models: a comparison of the sojourn times and transition intensities approaches. Int J Biostat2021; 18: 1–20.
16.
AastveitME. Multi-State Models for Interval-Censored Data with Transition Times from Gamma Processes. Master’s Thesis, 2021.
17.
HougaardP. Multi-state models: a review. Lifetime Data Anal1999; 5: 239–264.
18.
HuangJWellnerJA. Interval censored survival data: a review of recent progress. In Proceedings of the first Seattle symposium in Biostatistics. Springer, 1997, pp. 123–169.
19.
KimJS. Maximum likelihood estimation for the proportional hazards model with partly interval-censored data. J R Stat Soc: Ser B (Stat Methodol)2003; 65: 489–502.
20.
ClaeskensGHjortNL. Model Selection and Model Averaging. Cambridge: Cambridge University Press, 2008
21.
AalenOBorganØGjessingH. Survival and event history analysis: a process point of view. New York: Springer Science & Business Media, 2008.
22.
DattaSSattenGA. Validity of the Aalen–Johansen estimators of stage occupation probabilities and Nelson–Aalen estimators of integrated transition hazards for non-Markov models. Stat Probab Lett2001; 55: 403–411.
CsardiGNepuszT. The igraph software package for complex network research. Int J Complex Syst2006; 1695: 1–9. https://igraph.org.
25.
NarasimhanBJohnsonSGHahnTet al. Cubature: Adaptive Multivariate Integration over Hypercubes, 2022. https://bnaras.github.io/cubature/. R package version 2.0.4.4.
26.
PengRD. R Programming for Data Science. BC, Canada: Leanpub Victoria, 2016.
