Saddlepoint approximations for the computation of survival and hazard functions are introduced in the context of parametric survival analysis. Although these approximations are computationally fast, accurate, and relatively straightforward to implement, their use in survival analysis has been lacking. We approximate survival functions using the Lugannani and Rice saddlepoint approximation to the distribution function or by numerically integrating the saddlepoint density approximation. The hazard function is approximated using the saddlepoint density and distribution functions. The approximations are especially useful for consideration of survival and hazard functions for waiting times in complicated models. Examples include total or partial waiting times for a disease that progresses through various stages (convolutions of distributions).
We present methodology giving highly accurate approximations for Bayesian predictive densities and distribution functions of first passage times between states of a semi-Markov process with a finite number of states. When the states describe a degenerative disorder with an absorbing end state, such predictive distributions are the survival distributions of a patient. We illustrate these methods with a variety of examples, including data from the San Francisco AIDS study. We achieve our approximations using a three-step sequence. First, we introduce advanced concepts of flowgraph theory, which allow us to compute the moment generating function of the first passage time given the model parameters. Next, we use saddlepoint approximations to convert this into a density or distribution function conditional on the model parameter. Finally, we use Monte Carlo methods to remove dependence on the model parameter. These methods apply quite generally to all finite-state semi-Markov models in discrete or continuous time. Currently, there are no competing alternative methods that can achieve the saddlepoint accuracy of these computations.
Modeling recurrent event data is of current interest in statistics and engineering. This article proposes a framework for incorporating covariates in flowgraph models, with application to recurrent event data in systems reliability settings. A flowgraph is a generalized transition graph (GTG) originally developed to model total system waiting times for semi-Markov processes. The focus of flowgraph models is expanded by linking covariates into branch transition models, enriching the toolkit of available data analysis methods for complex stochastic systems. This article takes a Bayesian approach to the analysis of flowgraph models. Potential applications are not limited to engineering systems, but also extend to survival analysis.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.