Cluster approximations for infection dynamics on random networks

The effects of heterogeneity on stochastic cycles in epidemics

Models of biological processes are often subject to different sources of noise. Developing an understanding of the combined effects of different types of uncertainty is an open challenge. In this paper, we study a variant of the susceptible-infective-recovered model of epidemic spread, which combines both agent-to-agent heterogeneity and intrinsic noise. We focus on epidemic cycles, driven by the stochasticity of infection and recovery events, and study in detail how heterogeneity in susceptibilities and propensities to pass on the disease affects these quasi-cycles. While the system can only be described by a large hierarchical set of equations in the transient regime, we derive a reduced closed set of equations for population-level quantities in the stationary regime. We analytically obtain the spectra of quasi-cycles in the linear-noise approximation. We find that the characteristic frequency of these cycles is typically determined by population averages of susceptibilities and infectivities, but that their amplitude depends on higher-order moments of the heterogeneity. We also investigate the synchronisation properties and phase lag between different groups of susceptible and infected individuals.

SIR model with local and global infective contacts: A deterministic approach and applications

An epidemic model with births and deaths is considered on a two-dimensional L×L lattice. Each individual can have global infective contacts according to the standard susceptible-infected-recovered (SIR) model rules or local infective contacts with their nearest neighbors. We propose a deterministic approach to this model and, for the parameters corresponding to pertussis and rubella in the prevaccine era, verify that there is a close agreement with the stochastic simulations when epidemic spread or endemic stationarity is considered. We also find that our approach captures the characteristic features of the dynamic behavior of the system after a sudden decrease in global contacts that may arise as a consequence of health care measures. By using the deterministic approach, we are able to characterize the exponential growth of the epidemic behavior and analyze the stability of the system at the stationary values. Since the deterministic approximation captures the essential features of the disease transmission dynamics of the stochastic model, it provides a useful tool for performing systematic studies as a function of the model parameters. We give an example of this potentiality by analyzing the likelihood of the endemic state to become extinct when the weight of the global contacts is drastically reduced.

Phase lag in epidemics on a network of cities

We study the synchronization and phase lag of fluctuations in the number of infected individuals in a network of cities between which individuals commute. The frequency and amplitude of these oscillations is known to be very well captured by the van Kampen system-size expansion, and we use this approximation to compute the complex coherence function that describes their correlation. We find that, if the infection rate differs from city to city and the coupling between them is not too strong, these oscillations are synchronized with a well-defined phase lag between cities. The analytic description of the effect is shown to be in good agreement with the results of stochastic simulations for realistic population sizes.

Stochastic formulation of ecological models and their applications

The increasing use of computer simulation by theoretical ecologists started a move away from models formulated at the population level towards individual-based models. However, many of the models studied at the individual level are not analysed mathematically and remain defined in terms of a computer algorithm. This is not surprising, given that they are intrinsically stochastic and require tools and techniques for their study that may be unfamiliar to ecologists. Here, we argue that the construction of ecological models at the individual level and their subsequent analysis is, in many cases, straightforward and leads to important insights. We discuss recent work that highlights the importance of stochastic effects for parameter ranges and systems where it was previously thought that such effects would be negligible.

Stochastic oscillations in models of epidemics on a network of cities

We carry out an analytic investigation of stochastic oscillations in a susceptible-infected-recovered model of disease spread on a network of n cities. In the model a fraction f(jk) of individuals from city k commute to city j, where they may infect, or be infected by, others. Starting from a continuous-time Markov description of the model the deterministic equations, which are valid in the limit when the population of each city is infinite, are recovered. The stochastic fluctuations about the fixed point of these equations are derived by use of the van Kampen system-size expansion. The fixed point structure of the deterministic equations is remarkably simple: A unique nontrivial fixed point always exists and has the feature that the fraction of susceptible, infected, and recovered individuals is the same for each city irrespective of its size. We find that the stochastic fluctuations have an analogously simple dynamics: All oscillations have a single frequency, equal to that found in the one-city case. We interpret this phenomenon in terms of the properties of the spectrum of the matrix of the linear approximation of the deterministic equations at the fixed point.

Outbreak properties of epidemic models: the roles of temporal forcing and stochasticity on pathogen invasion dynamics

Despite temporally forced transmission driving many infectious diseases, analytical insight into its role when combined with stochastic disease processes and non-linear transmission has received little attention. During disease outbreaks, however, the absence of saturation effects early on in well-mixed populations mean that epidemic models may be linearised and we can calculate outbreak properties, including the effects of temporal forcing on fade-out, disease emergence and system dynamics, via analysis of the associated master equations. The approach is illustrated for the unforced and forced SIR and SEIR epidemic models. We demonstrate that in unforced models, initial conditions (and any uncertainty therein) play a stronger role in driving outbreak properties than the basic reproduction number R0, while the same properties are highly sensitive to small amplitude temporal forcing, particularly when R0 is small. Although illustrated for the SIR and SEIR models, the master equation framework may be applied to more realistic models, although analytical intractability scales rapidly with increasing system dimensionality. One application of these methods is obtaining a better understanding of the rate at which vector-borne and waterborne infectious diseases invade new regions given variability in environmental drivers, a particularly important question when addressing potential shifts in the global distribution and intensity of infectious diseases under climate change.

Stochastic effects in a seasonally forced epidemic model

The interplay of seasonality, the system's nonlinearities and intrinsic stochasticity, is studied for a seasonally forced susceptible-exposed-infective-recovered stochastic model. The model is explored in the parameter region that corresponds to childhood infectious diseases such as measles. The power spectrum of the stochastic fluctuations around the attractors of the deterministic system that describes the model in the thermodynamic limit is computed analytically and validated by stochastic simulations for large system sizes. Size effects are studied through additional simulations. Other effects such as switching between coexisting attractors induced by stochasticity often mentioned in the literature as playing an important role in the dynamics of childhood infectious diseases are also investigated. The main conclusion is that stochastic amplification, rather than these effects, is the key ingredient to understand the observed incidence patterns.