In recent work with Martin Bootsma, we revisit stochastic epidemic simulation using the Sellke construction. We show how heaps enable efficient simulation even in complex epidemiological scenarios, while stochastic coupling enables counterfactual simulation on a per-person basis
