We develop a new methodology for the efficient computation of epidemic final size distributions for a broad class of Markovian models. We exploit a particular representation of the stochastic epidemic process to derive a method which is both computat
ionally efficient and numerically stable. The algorithms we present are also physically transparent and so allow us to extend this method from the basic SIR model to a model with a phase-type infectious period and another with waning immunity. The underlying theory is applicable to many Markovian models where we wish to efficiently calculate hitting probabilities.
We calculate both the exponential and pre-factor contributions in a WKB approximation of the master equation for a stochastic SIR model with highly oscillatory dynamics. Fixing the basic parameters of the model we investigate how the outbreak distrib
ution changes with the population size. We show that this distribution rapidly becomes highly non-Gaussian, acquiring large tails indicating the presence of rare, but large outbreaks, as the population is made smaller. The analytic results are found to be in excellent agreement with simulations until the systems become so small that the dynamics are dominated by fade-out of the disease.
We study the stochastic susceptible-infected-recovered (SIR) model with time-dependent forcing using analytic techniques which allow us to disentangle the interaction of stochasticity and external forcing. The model is formulated as a continuous time
Markov process, which is decomposed into a deterministic dynamics together with stochastic corrections, by using an expansion in inverse system size. The forcing induces a limit cycle in the deterministic dynamics, and a complete analysis of the fluctuations about this time-dependent solution is given. This analysis is applied when the limit cycle is annual, and after a period-doubling when it is biennial. The comprehensive nature of our approach allows us to give a coherent picture of the dynamics which unifies past work, but which also provides a systematic method for predicting the periods of oscillations seen in whooping cough and measles epidemics.
