A Tutorial on Discrete Event Simulation Models Using a Cost-Effectiveness Analysis Example in R

Discrete event simulation (DES) is a flexible and computationally efficient approach for modeling diverse processes; however, DES remains underutilized in health care and medical decision making due to a lack of reliable and reproducible implementations. We developed an open-source DES framework to simulate individual-level state-transition models (iSTMs) in continuous time accounting for treatment effects, time dependence on state residence, and age-dependent mortality. Our DES implementation employs a modular and easily adaptable structure, with each module corresponding to a unique transition between health states. To simulate the evolution of the process (i.e., individual state transitions), we adapted the next-reaction algorithm from the stochastic chemical reactions literature. Simulation-time dependence (age-dependent mortality) and state residence time dependence (transition from sick to sicker) are seamlessly incorporated into the DES framework via validated nonparametric and parametric sampling routines (e.g., inversion method) of event times. Treatment effects are integrated as scaling factors of the hazard functions (proportional hazards). We illustrate the framework's benefits by implementing the Sick-Sicker Model in R and conduct a cost-effectiveness analysis and probabilistic analysis. We also obtain epidemiological outcomes of interest from the DES output, such as disease prevalence, survival probabilities, and distributions of state-specific dwell times. Our DES framework offers a reliable and accessible alternative that enables the simulation of more realistic dynamics of state-transition processes at potentially lower implementation and computational costs than discrete-time iSTMs.HighlightsDiscrete event simulation (DES) is a flexible and efficient approach to simulate diverse processes in model-based decision analysis.The tutorial presents an open-source DES framework to simulate individual-level state-transition models (iSTMs) in continuous time.The modular structure of our DES framework accommodates treatment effects, time-dependent transitions, and age-dependent mortality using validated sampling methods.The coded example in R uses the Sick-Sicker Model to compute a cost-effectiveness analysis, epidemiological outcomes, probabilistic analysis, and value-of-information analysis.