mstatecox: A Package for Simulating Transition Probabilities from Semiparametric Multistate Survival Models

Multistate duration models are a valuable tool used in multiple fields to examine how subjects move through a series of discrete phases/stages. The models themselves may be estimated using common statistical software, but their broader adoption has been limited due to a lack of software to substantively interpret their results. Transition probabilities are the common post-estimation quantity for interpreting multistate duration model results. De Wreede, Fiocco, and Putter’s mstate package provides R with the functionality to estimate these quantities from semiparametric multistate models, yet no Stata equivalent exists for semiparametric models. We introduce a new set of Stata commands to meet this need. Our mstatecox suite calculates transition probabilities from semiparametric multistate duration models via simulation. It can accommodate any configuration of stages, and also has the ability to accommodate time-interacted covariates. We demonstrate our package’s functionality using de Wreede, Fiocco, and Putter’s (2011) European Registry of Blood and Marrow Transplantation example dataset from R-mstate.