Estimating summary functionals in multistate models with an application to hospital infection data

Inhomogeneous Markov chains are a relevant framework for analysing event histories. The fundamental estimator in the presence of incompletely observed data is the Nelson-Aalen estimator of the multivariate cumulative hazards. It may be summarised in terms of probability estimates via the empirical transition matrix. The empirical transition matrix has only slowly entered applications, one reason being previous lack of software. In a number of applications, further summary measures are desired. We illustrate how they may be computed from the empirical transition matrix and why bootstrapping their variance works. In contrast, computing such summaries outside the present framework has typically led to biased results. As an example, we consider in more detail hospital stay following infectious complication. This summary quantity is often considered by clinical decision makers, but reliable estimates require modelling the timing of infection as in the present set-up. In this context, we also derive new summary measures that further distinguish between patients discharged and patients deceased.