Second look at the spread of epidemics on networks

In an important paper, M.E.J. Newman claimed that a general network-based stochastic Susceptible-Infectious-Removed (SIR) epidemic model is isomorphic to a bond percolation model, where the bonds are the edges of the contact network and the bond occupation probability is equal to the marginal probability of transmission from an infected node to a susceptible neighbor. In this paper, we show that this isomorphism is incorrect and define a semi-directed random network we call the epidemic percolation network that is exactly isomorphic to the SIR epidemic model in any finite population. In the limit of a large population, (i) the distribution of (self-limited) outbreak sizes is identical to the size distribution of (small) out-components, (ii) the epidemic threshold corresponds to the phase transition where a giant strongly-connected component appears, (iii) the probability of a large epidemic is equal to the probability that an initial infection occurs in the giant in-component, and (iv) the relative final size of an epidemic is equal to the proportion of the network contained in the giant out-component. For the SIR model considered by Newman, we show that the epidemic percolation network predicts the same mean outbreak size below the epidemic threshold, the same epidemic threshold, and the same final size of an epidemic as the bond percolation model. However, the bond percolation model fails to predict the correct outbreak size distribution and probability of an epidemic when there is a nondegenerate infectious period distribution. We confirm our findings by comparing predictions from percolation networks and bond percolation models to the results of simulations. In an appendix, we show that an isomorphism to an epidemic percolation network can be defined for any time-homogeneous stochastic SIR model.Comment: 29 pages, 5 figure

The epidemiology and transmissibility of Zika virus in Girardot and San Andres Island, Colombia, September 2015 to January 2016

Transmission of Zika virus (ZIKV) was first detected in Colombia in September 2015. As of April 2016, Colombia had reported over 65,000 cases of Zika virus disease (ZVD). We analysed daily surveillance data of ZVD cases reported to the health authorities of San Andres and Girardot, Colombia, between September 2015 and January 2016. ZVD was laboratory-confirmed by reverse transcription-polymerase chain reaction (RT-PCR) in the serum of acute cases within five days of symptom onset. We use daily incidence data to estimate the basic reproductive number (R0) in each population. We identified 928 and 1,936 reported ZVD cases from San Andres and Girardot, respectively. The overall attack rate for reported ZVD was 12.13 cases per 1,000 residents of San Andres and 18.43 cases per 1,000 residents of Girardot. Attack rates were significantly higher in females in both municipalities (p < 0.001). Cases occurred in all age groups with highest rates in 20 to 49 year-olds. The estimated R0 for the Zika outbreak was 1.41 (95% confidence interval (CI): 1.15-1.74) in San Andres and 4.61 (95% CI: 4.11-5.16) in Girardot. Transmission of ZIKV is ongoing in the Americas. The estimated R0 from Colombia supports the observed rapid spread

Effects of epidemic threshold definition on disease spread statistics

We study the statistical properties of the SIR epidemics in heterogeneous networks, when an epidemic is defined as only those SIR propagations that reach or exceed a minimum size s_c. Using percolation theory to calculate the average fractional size of an epidemic, we find that the strength of the spanning link percolation cluster $P_{\infty}$ is an upper bound to . For small values of s_c, $P_{\infty}$ is no longer a good approximation, and the average fractional size has to be computed directly. The value of s_c for which $P_{\infty}$ is a good approximation is found to depend on the transmissibility T of the SIR. We also study Q, the probability that an SIR propagation reaches the epidemic mass s_c, and find that it is well characterized by percolation theory. We apply our results to real networks (DIMES and Tracerouter) to measure the consequences of the choice s_c on predictions of average outcome sizes of computer failure epidemics.Comment: 12 pages, 8 figure

Commentary on the use of the reproduction number R during the COVID-19 pandemic

Since the beginning of the COVID-19 pandemic, the reproduction number R has become a popular epidemiological metric used to communicate the state of the epidemic. At its most basic, R is defined as the average number of secondary infections caused by one primary infected individual. R seems convenient, because the epidemic is expanding if R>1 and contracting if R<1. The magnitude of R indicates by how much transmission needs to be reduced to control the epidemic. Using R in a naïve way can cause new problems. The reasons for this are threefold: (1) There is not just one definition of R but many, and the precise definition of R affects both its estimated value and how it should be interpreted. (2) Even with a particular clearly defined R, there may be different statistical methods used to estimate its value, and the choice of method will affect the estimate. (3) The availability and type of data used to estimate R vary, and it is not always clear what data should be included in the estimation. In this review, we discuss when R is useful, when it may be of use but needs to be interpreted with care, and when it may be an inappropriate indicator of the progress of the epidemic. We also argue that careful definition of R, and the data and methods used to estimate it, can make R a more useful metric for future management of the epidemic

Dynamics of multi-stage infections on networks

This paper investigates the dynamics of infectious diseases with a nonexponentially distributed infectious period. This is achieved by considering a multistage infection model on networks. Using pairwise approximation with a standard closure, a number of important characteristics of disease dynamics are derived analytically, including the final size of an epidemic and a threshold for epidemic outbreaks, and it is shown how these quantities depend on disease characteristics, as well as the number of disease stages. Stochastic simulations of dynamics on networks are performed and compared to output of pairwise models for several realistic examples of infectious diseases to illustrate the role played by the number of stages in the disease dynamics. These results show that a higher number of disease stages results in faster epidemic outbreaks with a higher peak prevalence and a larger final size of the epidemic. The agreement between the pairwise and simulation models is excellent in the cases we consider

A monotonic relationship between the variability of the infectious period and final size in pairwise epidemic modelling

For a recently derived pairwise model of network epidemics with non-Markovian recovery, we prove that under some mild technical conditions on the distribution of the infectious periods, smaller variance in the recovery time leads to higher reproduction number, and consequently to a larger epidemic outbreak, when the mean infectious period is fixed. We discuss how this result is related to various stochastic orderings of the distributions of infectious periods. The results are illustrated by a number of explicit stochastic simulations, suggesting that their validity goes beyond regular networks

Effects of Heterogeneous and Clustered Contact Patterns on Infectious Disease Dynamics

The spread of infectious diseases fundamentally depends on the pattern of contacts between individuals. Although studies of contact networks have shown that heterogeneity in the number of contacts and the duration of contacts can have far-reaching epidemiological consequences, models often assume that contacts are chosen at random and thereby ignore the sociological, temporal and/or spatial clustering of contacts. Here we investigate the simultaneous effects of heterogeneous and clustered contact patterns on epidemic dynamics. To model population structure, we generalize the configuration model which has a tunable degree distribution (number of contacts per node) and level of clustering (number of three cliques). To model epidemic dynamics for this class of random graph, we derive a tractable, low-dimensional system of ordinary differential equations that accounts for the effects of network structure on the course of the epidemic. We find that the interaction between clustering and the degree distribution is complex. Clustering always slows an epidemic, but simultaneously increasing clustering and the variance of the degree distribution can increase final epidemic size. We also show that bond percolation-based approximations can be highly biased if one incorrectly assumes that infectious periods are homogeneous, and the magnitude of this bias increases with the amount of clustering in the network. We apply this approach to model the high clustering of contacts within households, using contact parameters estimated from survey data of social interactions, and we identify conditions under which network models that do not account for household structure will be biased

Decay Properties of p pi- Systems Produced in Neutron Dissociation at Fermilab Energies

The authors have examined the decay distributions of (p{pi}{sup -}) systems produced in the reaction n + p {yields} (p{pi}{sup -}) + p for neutron momenta between 120 GeV/c and 300 GeV/c. Preliminary analysis of decay moments indicates the presence of large helicity-flip amplitudes even for small (p{pi}{sup -}) mass values, and does not support the hypothesis that the helicity non-flip (p{pi}{sup -}) states are produced peripherally in impact parameter. These results are in approximate agreement with predictions of the Deck mechanism. The experiment was performed at the M-3 neutral beam of Fermilab

Mean-field models for non-Markovian epidemics on networks

This paper introduces a novel extension of the edge-based compartmental model to epidemics where the transmission and recovery processes are driven by general independent probability distributions. Edge-based compartmental modelling is just one of many different approaches used to model the spread of an infectious disease on a network; the major result of this paper is the rigorous proof that the edge-based compartmental model and the message passing models are equivalent for general independent transmission and recovery processes. This implies that the new model is exact on the ensemble of configuration model networks of infinite size. For the case of Markovian transmission themessage passing model is re-parametrised into a pairwise-like model which is then used to derive many well-known pairwise models for regular networks, or when the infectious period is exponentially distributed or is of a fixed length

Erratic Flu Vaccination Emerges from Short-Sighted Behavior in Contact Networks

The effectiveness of seasonal influenza vaccination programs depends on individual-level compliance. Perceptions about risks associated with infection and vaccination can strongly influence vaccination decisions and thus the ultimate course of an epidemic. Here we investigate the interplay between contact patterns, influenza-related behavior, and disease dynamics by incorporating game theory into network models. When individuals make decisions based on past epidemics, we find that individuals with many contacts vaccinate, whereas individuals with few contacts do not. However, the threshold number of contacts above which to vaccinate is highly dependent on the overall network structure of the population and has the potential to oscillate more wildly than has been observed empirically. When we increase the number of prior seasons that individuals recall when making vaccination decisions, behavior and thus disease dynamics become less variable. For some networks, we also find that higher flu transmission rates may, counterintuitively, lead to lower (vaccine-mediated) disease prevalence. Our work demonstrates that rich and complex dynamics can result from the interaction between infectious diseases, human contact patterns, and behavior