1
|
Ahmadi A, Foster JM, Protas B. Data-driven optimal closures for mean-cluster models: Beyond the classical pair approximation. Phys Rev E 2022; 106:025313. [PMID: 36109923 DOI: 10.1103/physreve.106.025313] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 07/25/2022] [Indexed: 06/15/2023]
Abstract
This study concerns the mean-clustering approach to modeling the evolution of lattice dynamics. Instead of tracking the state of individual lattice sites, this approach describes the time evolution of the concentrations of different cluster types. It leads to an infinite hierarchy of ordinary differential equations which must be closed by truncation using a so-called closure condition. This condition approximates the concentrations of higher-order clusters in terms of the concentrations of lower-order ones. The pair approximation is the most common form of closure. Here, we consider its generalization, termed the "optimal approximation," which we calibrate using a robust data-driven strategy. To fix attention, we focus on a recently proposed structured lattice model for a nickel-based oxide, similar to that used as cathode material in modern commercial Li-ion batteries. The form of the obtained optimal approximation allows us to deduce a simple sparse closure model. In addition to being more accurate than the classical pair approximation, this "sparse approximation" is also physically interpretable which allows us to a posteriori refine the hypotheses underlying construction of this class of closure models. Moreover, the mean-cluster model closed with this sparse approximation is linear and hence analytically solvable such that its parametrization is straightforward, although it offers a good approximation of the actual time evolution of the cluster concentrations on short timescales only. On the other hand, parametrization of the mean-cluster model closed with the pair approximation is shown to lead to an ill-posed inverse problem.
Collapse
Affiliation(s)
- Avesta Ahmadi
- School of Computational Science & Engineering, McMaster University, Hamilton, Ontario, Canada L8S 4L8
| | - Jamie M Foster
- School of Mathematics & Physics, University of Portsmouth, Portsmouth, Hampshire PO1 2UP, United Kingdom
| | - Bartosz Protas
- Department of Mathematics & Statistics, McMaster University, Hamilton, Ontario, Canada L8S 4L8
| |
Collapse
|
2
|
Wang H, Moore JM, Small M, Wang J, Yang H, Gu C. Epidemic dynamics on higher-dimensional small world networks. APPLIED MATHEMATICS AND COMPUTATION 2022; 421:126911. [PMID: 35068617 PMCID: PMC8759951 DOI: 10.1016/j.amc.2021.126911] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Revised: 12/26/2021] [Accepted: 12/29/2021] [Indexed: 06/14/2023]
Abstract
Dimension governs dynamical processes on networks. The social and technological networks which we encounter in everyday life span a wide range of dimensions, but studies of spreading on finite-dimensional networks are usually restricted to one or two dimensions. To facilitate investigation of the impact of dimension on spreading processes, we define a flexible higher-dimensional small world network model and characterize the dependence of its structural properties on dimension. Subsequently, we derive mean field, pair approximation, intertwined continuous Markov chain and probabilistic discrete Markov chain models of a COVID-19-inspired susceptible-exposed-infected-removed (SEIR) epidemic process with quarantine and isolation strategies, and for each model identify the basic reproduction number R 0 , which determines whether an introduced infinitesimal level of infection in an initially susceptible population will shrink or grow. We apply these four continuous state models, together with discrete state Monte Carlo simulations, to analyse how spreading varies with model parameters. Both network properties and the outcome of Monte Carlo simulations vary substantially with dimension or rewiring rate, but predictions of continuous state models change only slightly. A different trend appears for epidemic model parameters: as these vary, the outcomes of Monte Carlo change less than those of continuous state methods. Furthermore, under a wide range of conditions, the four continuous state approximations present similar deviations from the outcome of Monte Carlo simulations. This bias is usually least when using the pair approximation model, varies only slightly with network size, and decreases with dimension or rewiring rate. Finally, we characterize the discrepancies between Monte Carlo and continuous state models by simultaneously considering network efficiency and network size.
Collapse
Affiliation(s)
- Haiying Wang
- Business School, University of Shanghai for Science and Technology, 334 Jungong Road, Shanghai, 200093, China
| | - Jack Murdoch Moore
- School of Physics Science and Engineering, Tongji University, 1239 Siping Road, Shanghai, 200092, Western Australia, China
| | - Michael Small
- Complex Systems Group, Department of Mathematics and Statistics, University of Western Australia, 35 Stirling Highway, Crawley, 6009, Australia
- Mineral Resources, CSIRO, 26 Dick Perry Ave, Kensington, 6151, Western Australia, Australia
| | - Jun Wang
- School of Economics and Management, Beihang University, 37 Xueyuan Road, Beijing, 100191, China
| | - Huijie Yang
- Business School, University of Shanghai for Science and Technology, 334 Jungong Road, Shanghai, 200093, China
| | - Changgui Gu
- Business School, University of Shanghai for Science and Technology, 334 Jungong Road, Shanghai, 200093, China
| |
Collapse
|
3
|
Kuga K, Tanimoto J. Effects of void nodes on epidemic spreads in networks. Sci Rep 2022; 12:3957. [PMID: 35273312 PMCID: PMC8913681 DOI: 10.1038/s41598-022-07985-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 02/22/2022] [Indexed: 11/17/2022] Open
Abstract
We present the pair approximation models for susceptible–infected–recovered (SIR) epidemic dynamics in a sparse network based on a regular network. Two processes are considered, namely, a Markovian process with a constant recovery rate and a non-Markovian process with a fixed recovery time. We derive the implicit analytical expression for the final epidemic size and explicitly show the epidemic threshold in both Markovian and non-Markovian processes. As the connection rate decreases from the original network connection, the epidemic threshold in which epidemic phase transits from disease-free to endemic increases, and the final epidemic size decreases. Additionally, for comparison with sparse and heterogeneous networks, the pair approximation models were applied to a heterogeneous network with a degree distribution. The obtained phase diagram reveals that, upon increasing the degree of the original random regular networks and decreasing the effective connections by introducing void nodes accordingly, the final epidemic size of the sparse network is close to that of the random network with average degree of 4. Thus, introducing the void nodes in the network leads to more heterogeneous network and reduces the final epidemic size.
Collapse
Affiliation(s)
- Kazuki Kuga
- Faculty of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka, 816-8580, Japan.
| | - Jun Tanimoto
- Faculty of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka, 816-8580, Japan.,Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka, 816-8580, Japan
| |
Collapse
|
4
|
Yang S, Senapati P, Wang D, Bauch CT, Fountoulakis K. Targeted pandemic containment through identifying local contact network bottlenecks. PLoS Comput Biol 2021; 17:e1009351. [PMID: 34460813 PMCID: PMC8432902 DOI: 10.1371/journal.pcbi.1009351] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 09/10/2021] [Accepted: 08/13/2021] [Indexed: 01/24/2023] Open
Abstract
Decision-making about pandemic mitigation often relies upon simulation modelling. Models of disease transmission through networks of contacts–between individuals or between population centres–are increasingly used for these purposes. Real-world contact networks are rich in structural features that influence infection transmission, such as tightly-knit local communities that are weakly connected to one another. In this paper, we propose a new flow-based edge-betweenness centrality method for detecting bottleneck edges that connect nodes in contact networks. In particular, we utilize convex optimization formulations based on the idea of diffusion with p-norm network flow. Using simulation models of COVID-19 transmission through real network data at both individual and county levels, we demonstrate that targeting bottleneck edges identified by the proposed method reduces the number of infected cases by up to 10% more than state-of-the-art edge-betweenness methods. Furthermore, the proposed method is orders of magnitude faster than existing methods. During the COVID-19 pandemic decision makers frequently face questions like where to impose a lockdown, which traffic to close, and whom to quarantine, all required to be carried out at minimal costs. Establishing cost-effective pandemic control policies requires identifying good targets. New computational models from network theory and epidemic simulations over real contact networks provide a valuable tool for finding the right bottlenecks to target upon. Here we study a computationally efficient network centrality measure that enables us to detect local transmission bottlenecks, i.e., contact edges that are especially important for the spread of disease among small communities or local network structures inside large networks. We find that pandemic intervention strategies that target at local network structures significantly outperform interventions that solely focus on the entire network structure as a whole, which are traditionally believed to be the most effective.
Collapse
Affiliation(s)
- Shenghao Yang
- School of Computer Science, University of Waterloo, Waterloo, Ontario, Canada
- * E-mail:
| | - Priyabrata Senapati
- School of Computer Science, University of Waterloo, Waterloo, Ontario, Canada
| | - Di Wang
- Google Research, Mountain View, California, United States
| | - Chris T. Bauch
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada
| | - Kimon Fountoulakis
- School of Computer Science, University of Waterloo, Waterloo, Ontario, Canada
| |
Collapse
|
5
|
Kuga K, Tanaka M, Tanimoto J. Pair approximation model for the vaccination game: predicting the dynamic process of epidemic spread and individual actions against contagion. Proc Math Phys Eng Sci 2021; 477:20200769. [PMID: 35153542 PMCID: PMC8317980 DOI: 10.1098/rspa.2020.0769] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2020] [Accepted: 01/12/2021] [Indexed: 12/31/2022] Open
Abstract
We successfully establish a theoretical framework of pairwise approximation for the vaccination game in which both the dynamic process of epidemic spread and individual actions in helping prevent social behaviours are quantitatively evaluated. In contrast with mean-field approximation, our model captures higher-order effects from neighbours by using an underlying network that shows how the disease spreads and how individual decisions evolve over time. This model considers not only imperfect vaccination but also intermediate protective measures other than vaccines. Our analytical predictions are validated by multi-agent simulation results that estimate random regular graphs at varying degrees.
Collapse
Affiliation(s)
- Kazuki Kuga
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka 816-8580, Japan
| | - Masaki Tanaka
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka 816-8580, Japan
| | - Jun Tanimoto
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka 816-8580, Japan.,Faculty of Engineering Sciences, Kyushu University, Kasuga-koen, Kasuga-shi, Fukuoka 816-8580, Japan
| |
Collapse
|
6
|
Rodriguez-Brenes IA, Wodarz D, Komarova NL. Beyond the pair approximation: Modeling colonization population dynamics. Phys Rev E 2021; 101:032404. [PMID: 32289892 DOI: 10.1103/physreve.101.032404] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Accepted: 01/02/2020] [Indexed: 11/07/2022]
Abstract
The process of range expansion (colonization) is one of the basic types of biological dynamics, whereby a species grows and spreads outwards, occupying new territories. Spatial modeling of this process is naturally implemented as a stochastic cellular automaton, with individuals occupying nodes on a rectangular grid, births and deaths occurring probabilistically, and individuals only reproducing onto unoccupied neighboring spots. In this paper we derive several approximations that allow prediction of the expected range expansion dynamics, based on the reproduction and death rates. We derive several approximations, where the cellular automaton is described by a system of ordinary differential equations that preserves correlations among neighboring spots (up to a distance). This methodology allows us to develop accurate approximations of the population size and the expected spatial shape, at a fraction of the computational time required to simulate the original stochastic system. In addition, we provide simple formulas for the steady-state population densities for von Neumann and Moore neighborhoods. Finally, we derive concise approximations for the speed of range expansion in terms of the reproduction and death rates, for both types of neighborhoods. The methodology is generalizable to more complex scenarios, such as different interaction ranges and multiple-species systems.
Collapse
Affiliation(s)
| | - Dominik Wodarz
- Department of Population Health and Disease Prevention, University of California, Irvine, California 92617, USA
| | - Natalia L Komarova
- Department of Mathematics, University of California Irvine, Irvine, California 92697, USA
| |
Collapse
|
7
|
Horrocks J, Bauch CT. Algorithmic discovery of dynamic models from infectious disease data. Sci Rep 2020; 10:7061. [PMID: 32341374 PMCID: PMC7184751 DOI: 10.1038/s41598-020-63877-w] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Accepted: 04/07/2020] [Indexed: 11/09/2022] Open
Abstract
Theoretical models are typically developed through a deductive process where a researcher formulates a system of dynamic equations from hypothesized mechanisms. Recent advances in algorithmic methods can discover dynamic models inductively-directly from data. Most previous research has tested these methods by rediscovering models from synthetic data generated by the already known model. Here we apply Sparse Identification of Nonlinear Dynamics (SINDy) to discover mechanistic equations for disease dynamics from case notification data for measles, chickenpox, and rubella. The discovered models provide a good qualitative fit to the observed dynamics for all three diseases, However, the SINDy chickenpox model appears to overfit the empirical data, and recovering qualitatively correct rubella dynamics requires using power spectral density in the goodness-of-fit criterion. When SINDy uses a library of second-order functions, the discovered models tend to include mass action incidence and a seasonally varying transmission rate-a common feature of existing epidemiological models for childhood infectious diseases. We also find that the SINDy measles model is capable of out-of-sample prediction of a dynamical regime shift in measles case notification data. These results demonstrate the potential for algorithmic model discovery to enrich scientific understanding by providing a complementary approach to developing theoretical models.
Collapse
Affiliation(s)
- Jonathan Horrocks
- Department of Applied Mathematics, University of Waterloo, Waterloo, N2L 3G1, Canada
| | - Chris T Bauch
- Department of Applied Mathematics, University of Waterloo, Waterloo, N2L 3G1, Canada.
| |
Collapse
|
8
|
Hamada M, Takasu F. Equilibrium properties of the spatial SIS model as a point pattern dynamics - How is infection distributed over space? J Theor Biol 2019; 468:12-26. [PMID: 30738052 DOI: 10.1016/j.jtbi.2019.02.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2018] [Revised: 02/01/2019] [Accepted: 02/06/2019] [Indexed: 10/27/2022]
Abstract
We revisit the classical epidemiological SIS model as a stochastic point pattern dynamics with special focus on its spatial distribution at equilibrium. In this model, each point on a continuous space is either susceptible S or infectious I, and infection occurs with an infection kernel as a function of distance from I to S. This stochastic process has been mathematically described by the hierarchical dynamics of the probabilities that a point, a pair made by two points, and a triplet made by three points, etc., is in a specific configuration of status. Using a simple closure thereby triplet probabilities that appear in the dynamics are approximated, we show that the average singlet probabilities and the pair probabilities that describe spatial distributions of Ss and Is at equilibrium can be explicitly derived using the infection kernel; Is are spatially clustered in the same order of the infection kernel. The results highlight the advantage of point pattern approach to model spatial population dynamics in general ecology where local interactions among individuals likely depend on distance between them.
Collapse
Affiliation(s)
- Miki Hamada
- Graduate School of Humanities and Sciences, Nara Women's University, Kita-Uoya Nishimachi, Nara 630-8506, Japan.
| | - Fugo Takasu
- Department of Environmental Science, Nara Women's University, Kita-Uoya Nishimachi, Nara 630-8506, Japan.
| |
Collapse
|
9
|
Lee MJ, Lee DS. Understanding the temporal pattern of spreading in heterogeneous networks: Theory of the mean infection time. Phys Rev E 2019; 99:032309. [PMID: 30999425 PMCID: PMC7217551 DOI: 10.1103/physreve.99.032309] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Indexed: 11/12/2022]
Abstract
For a reliable prediction of an epidemic or information spreading pattern in complex systems, well-defined measures are essential. In the susceptible-infected model on heterogeneous networks, the cluster of infected nodes in the intermediate-time regime exhibits too large fluctuation in size to use its mean size as a representative value. The cluster size follows quite a broad distribution, which is shown to be derived from the variation of the cluster size with the time when a hub node was first infected. On the contrary, the distribution of the time taken to infect a given number of nodes is well concentrated at its mean, suggesting the mean infection time is a better measure. We show that the mean infection time can be evaluated by using the scaling behaviors of the boundary area of the infected cluster and use it to find a nonexponential but algebraic spreading phase in the intermediate stage on strongly heterogeneous networks. Such slow spreading originates in only small-degree nodes left susceptible, while most hub nodes are already infected in the early exponential-spreading stage. Our results offer a way to detour around large statistical fluctuations and quantify reliably the temporal pattern of spread under structural heterogeneity.
Collapse
Affiliation(s)
- Mi Jin Lee
- Department of Physics, Inha University, Incheon 22212, Korea
| | - Deok-Sun Lee
- Department of Physics, Inha University, Incheon 22212, Korea
| |
Collapse
|
10
|
SIR dynamics in random networks with communities. J Math Biol 2018; 77:1117-1151. [PMID: 29752517 DOI: 10.1007/s00285-018-1247-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2017] [Revised: 03/26/2018] [Indexed: 10/16/2022]
Abstract
This paper investigates the effects of the community structure of a network on the spread of an epidemic. To this end, we first establish a susceptible-infected-recovered (SIR) model in a two-community network with an arbitrary joint degree distribution. The network is formulated as a probability generating function. We also obtain the sufficient conditions for disease outbreak and extinction, which involve the first-order and second-order moments of the degree distribution. As an example, we then study the effect of community structure on epidemic spread in a complex network with a Poisson joint degree distribution. The numerical solutions of the SIR model well agree with stochastic simulations based on the Monte Carlo method, confirming that the model is reliable and accurate. Finally, by strengthening the community structure in the simulation, i.e. fixing the total degree distribution and reducing the number ratio of the external edges, we can increase or decrease the final cumulative epidemic incidence depending on the transmissibility of the virus between humans and the community structure at that point. Why community structure can affect disease dynamics in a complicated way is also discussed. In any case, for large-scale epidemics, strengthening the community structure to reduce the size of disease is undoubtedly an effective way.
Collapse
|
11
|
Arias JH, Gómez-Gardeñes J, Meloni S, Estrada E. Epidemics on plants: Modeling long-range dispersal on spatially embedded networks. J Theor Biol 2018; 453:1-13. [PMID: 29738720 DOI: 10.1016/j.jtbi.2018.05.004] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2018] [Revised: 05/03/2018] [Accepted: 05/04/2018] [Indexed: 10/17/2022]
Abstract
Here we develop an epidemic model that accounts for long-range dispersal of pathogens between plants. This model generalizes the classical compartmental models-Susceptible-Infected-Susceptible (SIS) and Susceptible-Infected-Recovered (SIR)-to take into account those factors that are key to understand epidemics in real plant populations. These ingredients are the spatial characteristics of the plots and fields in which plants are embedded and the effect of long-range dispersal of pathogens. The spatial characteristics are included through the use of random rectangular graphs which allow to consider the effects of the elongation of plots and fields, while the long-range dispersal is implemented by considering transformations, such as the Mellin and Laplace transforms, of a generalization of the adjacency matrix of the geometric graph. Our results point out that long-range dispersal favors the propagation of pathogens while the elongation of plant plots increases the epidemic threshold and decreases dramatically the number of affected plants. Interestingly, our model is able of reproducing the existence of patchy regions of infected plants and the absence of a clear propagation front centered in the initial infected plants, as it is observed in real plant epidemics.
Collapse
Affiliation(s)
- Juddy H Arias
- Department of Mathematics, Universidad del Valle, Colombia
| | - Jesus Gómez-Gardeñes
- GOTHAM Lab, Institute for Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, Zaragoza, Spain; Department of Condensed Matter Physics, University of Zaragoza, Zaragoza, Spain
| | - Sandro Meloni
- Institute for Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, Zaragoza 50018, Spain; Department of Theoretical Physics, University of Zaragoza, Zaragoza 50009, Spain
| | - Ernesto Estrada
- Department of Mathematics & Statistics, University of Strathclyde, 26 Richmond Street, Glasgow, G11XH, UK.
| |
Collapse
|
12
|
Dynamics analysis of SIR epidemic model with correlation coefficients and clustering coefficient in networks. J Theor Biol 2018; 449:1-13. [PMID: 29649430 DOI: 10.1016/j.jtbi.2018.04.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2017] [Revised: 04/02/2018] [Accepted: 04/04/2018] [Indexed: 11/22/2022]
Abstract
In this paper, the correlation coefficients between nodes in states are used as dynamic variables, and we construct SIR epidemic dynamic models with correlation coefficients by using the pair approximation method in static networks and dynamic networks, respectively. Considering the clustering coefficient of the network, we analytically investigate the existence and the local asymptotic stability of each equilibrium of these models and derive threshold values for the prevalence of diseases. Additionally, we obtain two equivalent epidemic thresholds in dynamic networks, which are compared with the results of the mean field equations.
Collapse
|
13
|
Ringa N, Bauch CT. Spatially-implicit modelling of disease-behaviour interactions in the context of non-pharmaceutical interventions. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2018; 15:461-483. [PMID: 29161845 DOI: 10.3934/mbe.2018021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Pair approximation models have been used to study the spread of infectious diseases in spatially distributed host populations, and to explore disease control strategies such as vaccination and case isolation. Here we introduce a pair approximation model of individual uptake of non-pharmaceutical interventions (NPIs) for an acute self-limiting infection, where susceptible individuals can learn the NPIs either from other susceptible individuals who are already practicing NPIs ("social learning"), or their uptake of NPIs can be stimulated by being neighbours of an infectious person ("exposure learning"). NPIs include individual measures such as hand-washing and respiratory etiquette. Individuals can also drop the habit of using NPIs at a certain rate. We derive a spatially defined expression of the basic reproduction number R0 and we also numerically simulate the model equations. We find that exposure learning is generally more efficient than social learning, since exposure learning generates NPI uptake in the individuals at immediate risk of infection. However, if social learning is pre-emptive, beginning a sufficient amount of time before the epidemic, then it can be more effective than exposure learning. Interestingly, varying the initial number of individuals practicing NPIs does not significantly impact the epidemic final size. Also, if initial source infections are surrounded by protective individuals, there are parameter regimes where increasing the initial number of source infections actually decreases the infection peak (instead of increasing it) and makes it occur sooner. The peak prevalence increases with the rate at which individuals drop the habit of using NPIs, but the response of peak prevalence to changes in the forgetting rate are qualitatively different for the two forms of learning. The pair approximation methodology developed here illustrates how analytical approaches for studying interactions between social processes and disease dynamics in a spatially structured population should be further pursued.
Collapse
Affiliation(s)
- Notice Ringa
- Botswana International University of Science and Technology, Department of Mathematics and Statistical Sciences, Private Bag 16, Palapye, Botswana
| | - Chris T Bauch
- University of Waterloo, Department of Applied Mathematics, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
14
|
Li J, Li W, Jin Z. The epidemic model based on the approximation for third-order motifs on networks. Math Biosci 2018; 297:12-26. [PMID: 29330075 DOI: 10.1016/j.mbs.2018.01.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2017] [Revised: 01/06/2018] [Accepted: 01/08/2018] [Indexed: 10/18/2022]
Abstract
The spread of an infectious disease may depend on the structure of the network. To study the influence of the structure parameters of the network on the spread of the epidemic, we need to put these parameters into the epidemic model. The method of moment closure introduces structure parameters into the epidemic model. In this paper, we present a new moment closure epidemic model based on the approximation of third-order motifs in networks. The order of a motif defined in this paper is determined by the number of the edges in the motif, rather than by the number of nodes in the motif as defined in the literature. We provide a general approach to deriving a set of ordinary differential equations that describes, to a high degree of accuracy, the spread of an infectious disease. Using this method, we establish a susceptible-infected-recovered (SIR) model. We then calculate the basic reproduction number of the SIR model, and find that it decreases as the clustering coefficient increases. Finally, we perform some simulations using the proposed model to study the influence of the clustering coefficient on the final epidemic size, the maximum number of infected, and the peak time of the disease. The numerical simulations based on the SIR model in this paper fit the stochastic simulations based on the Monte Carlo method well at different levels of clustering. Our results show that the clustering coefficient poses impediments to the spread of disease under an SIR model.
Collapse
Affiliation(s)
- Jinxian Li
- School of Mathematical Sciences, Shanxi University, Taiyuan 030006, PR China; Shanxi Key Laboratory of Mathematical Techniques and Big Data Analysis on Disease Control and Prevention, Taiyuan 030006, PR China
| | - Weiqiang Li
- School of Mathematical Sciences, Shanxi University, Taiyuan 030006, PR China
| | - Zhen Jin
- Complex System Research Center, Shanxi University, Taiyuan 030006, PR China; Shanxi Key Laboratory of Mathematical Techniques and Big Data Analysis on Disease Control and Prevention, Taiyuan 030006, PR China.
| |
Collapse
|
15
|
Chen S, Wang K, Sun M, Fu X. Spread of competing viruses on heterogeneous networks. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:rsta.2016.0284. [PMID: 28507229 PMCID: PMC5434075 DOI: 10.1098/rsta.2016.0284] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 12/23/2016] [Indexed: 05/08/2023]
Abstract
In this paper, we propose a model where two strains compete with each other at the expense of common susceptible individuals on heterogeneous networks by using pair-wise approximation closed by the probability-generating function (PGF). All of the strains obey the susceptible-infected-recovered (SIR) mechanism. From a special perspective, we first study the dynamical behaviour of an SIR model closed by the PGF, and obtain the basic reproduction number via two methods. Then we build a model to study the spreading dynamics of competing viruses and discuss the conditions for the local stability of equilibria, which is different from the condition obtained by using the heterogeneous mean-field approach. Finally, we perform numerical simulations on Barabási-Albert networks to complement our theoretical research, and show some dynamical properties of the model with competing viruses.This article is part of the themed issue 'Mathematical methods in medicine: neuroscience, cardiology and pathology'.
Collapse
Affiliation(s)
- Shanshan Chen
- Department of Mathematics, Shanghai University, Shanghai 200444, People's Republic of China
| | - Kaihua Wang
- College of Mathematics and Statistics, Hainan Normal University, Haikou 571158, People's Republic of China
| | - Mengfeng Sun
- Department of Mathematics, Shanghai University, Shanghai 200444, People's Republic of China
| | - Xinchu Fu
- Department of Mathematics, Shanghai University, Shanghai 200444, People's Republic of China
| |
Collapse
|
16
|
LUO XIAOFENG, CHANG LILI, JIN ZHEN. DEMOGRAPHICS INDUCE EXTINCTION OF DISEASE IN AN SIS MODEL BASED ON CONDITIONAL MARKOV CHAIN. J BIOL SYST 2017. [DOI: 10.1142/s0218339017500085] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Demographics have significant effects on disease spread in populations and the topological evolution of the underlying networks that represent the populations. In the context of network-based epidemic modeling, Markov chain-based approach and pairwise approximation are two powerful tools — the former can capture stochastic effects of disease transmission dynamics and the latter can characterize the dynamical correlations in each pair of connected individuals. However, to our best knowledge, the study on epidemic spreading in networks relying on these two techniques is still lacking. To fill this gap, in this paper, a deterministic pairwise susceptible–infected–susceptible (SIS) epidemic model with demographics on complex networks with arbitrary degree distributions is studied based on a continuous time conditional Markov chain. This deterministic model is rigorously derived — using the moment generating function — from the Kolmogorov differential equations for the evolution of individuals and pairs. It is found that demographics will induce the extinction of the disease by reducing the basic reproduction number or lowering the epidemic prevalence after the disease prevails. Moreover, due to the demographical effects, the resulting network tends to a homogeneous network with a degree distribution similar to Poisson distribution, irrespective of the initial network structure. Additionally, we find excellent agreement between numerical solutions and individual-based stochastic simulations using both Erdös–Renyi (ER) random and Barabási–Albert (BA) scale-free initial networks. Our results may provide new insights on the understanding of the influence of demographics on epidemic dynamics and network evolution.
Collapse
Affiliation(s)
- XIAOFENG LUO
- School of Computer and Information Technology, Shanxi University, Taiyuan, Shanxi 030006, P. R. China
| | - LILI CHANG
- Complex System Research Center, Shanxi University, Taiyuan, Shanxi 030006, P. R. China
| | - ZHEN JIN
- Complex System Research Center, Shanxi University, Taiyuan, Shanxi 030006, P. R. China
| |
Collapse
|
17
|
Dynamics of epidemic diseases on a growing adaptive network. Sci Rep 2017; 7:42352. [PMID: 28186146 PMCID: PMC5301221 DOI: 10.1038/srep42352] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2016] [Accepted: 01/08/2017] [Indexed: 12/03/2022] Open
Abstract
The study of epidemics on static networks has revealed important effects on disease prevalence of network topological features such as the variance of the degree distribution, i.e. the distribution of the number of neighbors of nodes, and the maximum degree. Here, we analyze an adaptive network where the degree distribution is not independent of epidemics but is shaped through disease-induced dynamics and mortality in a complex interplay. We study the dynamics of a network that grows according to a preferential attachment rule, while nodes are simultaneously removed from the network due to disease-induced mortality. We investigate the prevalence of the disease using individual-based simulations and a heterogeneous node approximation. Our results suggest that in this system in the thermodynamic limit no epidemic thresholds exist, while the interplay between network growth and epidemic spreading leads to exponential networks for any finite rate of infectiousness when the disease persists.
Collapse
|
18
|
|
19
|
Sharkey KJ, Wilkinson RR. Complete hierarchies of SIR models on arbitrary networks with exact and approximate moment closure. Math Biosci 2015; 264:74-85. [PMID: 25829147 DOI: 10.1016/j.mbs.2015.03.008] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2014] [Revised: 03/20/2015] [Accepted: 03/23/2015] [Indexed: 11/19/2022]
Abstract
We first generalise ideas discussed by Kiss et al. (2015) to prove a theorem for generating exact closures (here expressing joint probabilities in terms of their constituent marginal probabilities) for susceptible-infectious-removed (SIR) dynamics on arbitrary graphs (networks). For Poisson transmission and removal processes, this enables us to obtain a systematic reduction in the number of differential equations needed for an exact 'moment closure' representation of the underlying stochastic model. We define 'transmission blocks' as a possible extension of the block concept in graph theory and show that the order at which the exact moment closure representation is curtailed is the size of the largest transmission block. More generally, approximate closures of the hierarchy of moment equations for these dynamics are typically defined for the first and second order yielding mean-field and pairwise models respectively. It is frequently implied that, in principle, closed models can be written down at arbitrary order if only we had the time and patience to do this. However, for epidemic dynamics on networks, these higher-order models have not been defined explicitly. Here we unambiguously define hierarchies of approximate closed models that can utilise subsystem states of any order, and show how well-known models are special cases of these hierarchies.
Collapse
Affiliation(s)
- Kieran J Sharkey
- Department of Mathematical Sciences, University of Liverpool, Peach Street, Liverpool, L69 7ZL, United Kingdom.
| | - Robert R Wilkinson
- Department of Mathematical Sciences, University of Liverpool, Peach Street, Liverpool, L69 7ZL, United Kingdom
| |
Collapse
|
20
|
Ringa N, Bauch CT. Impacts of constrained culling and vaccination on control of foot and mouth disease in near-endemic settings: a pair approximation model. Epidemics 2014; 9:18-30. [PMID: 25480131 DOI: 10.1016/j.epidem.2014.09.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2014] [Revised: 09/19/2014] [Accepted: 09/21/2014] [Indexed: 11/18/2022] Open
Abstract
Many countries have eliminated foot and mouth disease (FMD), but outbreaks remain common in other countries. Rapid development of international trade in animals and animal products has increased the risk of disease introduction to FMD-free countries. Most mathematical models of FMD are tailored to settings that are normally disease-free, and few models have explored the impact of constrained control measures in a 'near-endemic' spatially distributed host population subject to frequent FMD re-introductions from nearby endemic wild populations, as characterizes many low-income, resource-limited countries. Here we construct a pair approximation model of FMD and investigate the impact of constraints on total vaccine supply for prophylactic and ring vaccination, and constraints on culling rates and cumulative culls. We incorporate natural immunity waning and vaccine waning, which are important factors for near-endemic populations. We find that, when vaccine supply is sufficiently limited, the optimal approach for minimizing cumulative infections combines rapid deployment of ring vaccination during outbreaks with a contrasting approach of careful rationing of prophylactic vaccination over the year, such that supplies last as long as possible (and with the bulk of vaccines dedicated toward prophylactic vaccination). Thus, for optimal long-term control of the disease by vaccination in near-endemic settings when vaccine supply is limited, it is best to spread out prophylactic vaccination as much as possible. Regardless of culling constraints, the optimal culling strategy is rapid identification of infected premises and their immediate contacts at the initial stages of an outbreak, and rapid culling of infected premises and farms deemed to be at high risk of infection (as opposed to culling only the infected farms). Optimal culling strategies are similar when social impact is the outcome of interest. We conclude that more FMD transmission models should be developed that are specific to the challenges of FMD control in near-endemic, low-income countries.
Collapse
Affiliation(s)
- N Ringa
- Department of Mathematics and Statistics, University of Guelph, 50 Stone Rd E, Guelph, ON N1G 2W1, Canada.
| | - C T Bauch
- Department of Mathematics and Statistics, University of Guelph, 50 Stone Rd E, Guelph, ON N1G 2W1, Canada; Department of Applied Mathematics, University of Waterloo, 200 University Avenue, West Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
21
|
Llensa C, Juher D, Saldaña J. On the early epidemic dynamics for pairwise models. J Theor Biol 2014; 352:71-81. [DOI: 10.1016/j.jtbi.2014.02.037] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2013] [Revised: 02/21/2014] [Accepted: 02/26/2014] [Indexed: 12/01/2022]
|
22
|
Ringa N, Bauch CT. Dynamics and control of foot-and-mouth disease in endemic countries: a pair approximation model. J Theor Biol 2014; 357:150-9. [PMID: 24853274 DOI: 10.1016/j.jtbi.2014.05.010] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2013] [Revised: 04/07/2014] [Accepted: 05/05/2014] [Indexed: 10/25/2022]
Abstract
Previous mathematical models of spatial farm-to-farm transmission of foot and mouth disease (FMD) have explored the impacts of control measures such as culling and vaccination during a single outbreak in a country normally free of FMD. As a result, these models do not include factors that are relevant to countries where FMD is endemic in some regions, like long-term waning natural and vaccine immunity, use of prophylactic vaccination and disease re-importations. These factors may have implications for disease dynamics and control, yet few models have been developed for FMD-endemic settings. Here we develop and study an SEIRV (susceptible-exposed-infectious-recovered-vaccinated) pair approximation model of FMD. We focus on long term dynamics by exploring characteristics of repeated outbreaks of FMD and their dependence on disease re-importation, loss of natural immunity, and vaccine waning. We find that the effectiveness of ring and prophylactic vaccination strongly depends on duration of natural immunity, rate of vaccine waning, and disease re-introduction rate. However, the number and magnitude of FMD outbreaks are generally more sensitive to the duration of natural immunity than the duration of vaccine immunity. If loss of natural immunity and/or vaccine waning happen rapidly, then multiple epidemic outbreaks result, making it difficult to eliminate the disease. Prophylactic vaccination is more effective than ring vaccination, at the same per capita vaccination rate. Finally, more frequent disease re-importation causes a higher cumulative number of infections, although a lower average epidemic peak. Our analysis demonstrates significant differences between dynamics in FMD-free settings versus FMD-endemic settings, and that dynamics in FMD-endemic settings can vary widely depending on factors such as the duration of natural and vaccine immunity and the rate of disease re-importations. We conclude that more mathematical models tailored to FMD-endemic countries should be developed that include these factors.
Collapse
Affiliation(s)
- N Ringa
- Department of Mathematics and Statistics, University of Guelph, 50 Stone Rd E, Guelph, Canada ON N1G 2W1.
| | - C T Bauch
- Department of Mathematics and Statistics, University of Guelph, 50 Stone Rd E, Guelph, Canada ON N1G 2W1; Department of Applied Mathematics, University of Waterloo, 200 University Avenue West Waterloo, Canada ON N2L 3G1
| |
Collapse
|
23
|
Droste F, Do AL, Gross T. Analytical investigation of self-organized criticality in neural networks. J R Soc Interface 2012; 10:20120558. [PMID: 22977096 DOI: 10.1098/rsif.2012.0558] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Dynamical criticality has been shown to enhance information processing in dynamical systems, and there is evidence for self-organized criticality in neural networks. A plausible mechanism for such self-organization is activity-dependent synaptic plasticity. Here, we model neurons as discrete-state nodes on an adaptive network following stochastic dynamics. At a threshold connectivity, this system undergoes a dynamical phase transition at which persistent activity sets in. In a low-dimensional representation of the macroscopic dynamics, this corresponds to a transcritical bifurcation. We show analytically that adding activity-dependent rewiring rules, inspired by homeostatic plasticity, leads to the emergence of an attractive steady state at criticality and present numerical evidence for the system's evolution to such a state.
Collapse
Affiliation(s)
- Felix Droste
- Bernstein Center for Computational Neuroscience, Haus 2, Philippstrasse 13, Berlin, Germany.
| | | | | |
Collapse
|
24
|
Messinger SM, Ostling A. The influence of host demography, pathogen virulence, and relationships with pathogen virulence on the evolution of pathogen transmission in a spatial context. Evol Ecol 2012. [DOI: 10.1007/s10682-012-9594-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
25
|
Outbreak analysis of an SIS epidemic model with rewiring. J Math Biol 2012; 67:411-32. [DOI: 10.1007/s00285-012-0555-4] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2012] [Revised: 05/18/2012] [Indexed: 11/24/2022]
|
26
|
The effect of predation on the prevalence and aggregation of pathogens in prey. Biosystems 2011; 105:300-6. [DOI: 10.1016/j.biosystems.2011.05.012] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2010] [Revised: 05/09/2011] [Accepted: 05/26/2011] [Indexed: 11/21/2022]
|
27
|
Wells CR, Tchuenche JM, Meyers LA, Galvani AP, Bauch CT. Impact of imitation processes on the effectiveness of ring vaccination. Bull Math Biol 2011; 73:2748-72. [PMID: 21409511 DOI: 10.1007/s11538-011-9646-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2010] [Accepted: 02/18/2011] [Indexed: 11/25/2022]
Abstract
Ring vaccination can be a highly effective control strategy for an emerging disease or in the final phase of disease eradication, as witnessed in the eradication of smallpox. However, the impact of behavioural dynamics on the effectiveness of ring vaccination has not been explored in mathematical models. Here, we analyze a series of stochastic models of voluntary ring vaccination. Contacts of an index case base vaccinating decisions on their own individual payoffs to vaccinate or not vaccinate, and they can also imitate the behaviour of other contacts of the index case. We find that including imitation changes the probability of containment through ring vaccination considerably. Imitation can cause a strong majority of contacts to choose vaccination in some cases, or to choose non-vaccination in other cases-even when the equivalent solution under perfectly rational (non-imitative) behaviour yields mixed choices. Moreover, imitation processes can result in very different outcomes in different stochastic realizations sampled from the same parameter distributions, by magnifying moderate tendencies toward one behaviour or the other: in some realizations, imitation causes a strong majority of contacts not to vaccinate, while in others, imitation promotes vaccination and reduces the number of secondary infections. Hence, the effectiveness of ring vaccination can depend significantly and unpredictably on imitation processes. Therefore, our results suggest that risk communication efforts should be initiated early in an outbreak when ring vaccination is to be applied, especially among subpopulations that are heavily influenced by peer opinions.
Collapse
Affiliation(s)
- Chad R Wells
- Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario, Canada.
| | | | | | | | | |
Collapse
|
28
|
Danon L, Ford AP, House T, Jewell CP, Keeling MJ, Roberts GO, Ross JV, Vernon MC. Networks and the epidemiology of infectious disease. Interdiscip Perspect Infect Dis 2011; 2011:284909. [PMID: 21437001 PMCID: PMC3062985 DOI: 10.1155/2011/284909] [Citation(s) in RCA: 183] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2010] [Accepted: 12/24/2010] [Indexed: 11/23/2022] Open
Abstract
The science of networks has revolutionised research into the dynamics of interacting elements. It could be argued that epidemiology in particular has embraced the potential of network theory more than any other discipline. Here we review the growing body of research concerning the spread of infectious diseases on networks, focusing on the interplay between network theory and epidemiology. The review is split into four main sections, which examine: the types of network relevant to epidemiology; the multitude of ways these networks can be characterised; the statistical methods that can be applied to infer the epidemiological parameters on a realised network; and finally simulation and analytical methods to determine epidemic dynamics on a given network. Given the breadth of areas covered and the ever-expanding number of publications, a comprehensive review of all work is impossible. Instead, we provide a personalised overview into the areas of network epidemiology that have seen the greatest progress in recent years or have the greatest potential to provide novel insights. As such, considerable importance is placed on analytical approaches and statistical methods which are both rapidly expanding fields. Throughout this review we restrict our attention to epidemiological issues.
Collapse
Affiliation(s)
- Leon Danon
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK
| | - Ashley P. Ford
- Department of Statistics, University of Warwick, Coventry CV4 7AL, UK
| | - Thomas House
- Warwick Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK
| | - Chris P. Jewell
- Department of Statistics, University of Warwick, Coventry CV4 7AL, UK
| | - Matt J. Keeling
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK
- Warwick Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK
| | - Gareth O. Roberts
- Department of Statistics, University of Warwick, Coventry CV4 7AL, UK
| | - Joshua V. Ross
- School of Mathematical Sciences, University of Adelaide, SA 5005, Australia
| | - Matthew C. Vernon
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK
| |
Collapse
|
29
|
Sharkey KJ. Deterministic epidemic models on contact networks: correlations and unbiological terms. Theor Popul Biol 2011; 79:115-29. [PMID: 21354193 DOI: 10.1016/j.tpb.2011.01.004] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2010] [Revised: 01/26/2011] [Accepted: 01/27/2011] [Indexed: 11/28/2022]
Abstract
The relationship between system-level and subsystem-level master equations is investigated and then utilised for a systematic and potentially automated derivation of the hierarchy of moment equations in a susceptible-infectious-removed (SIR) epidemic model. In the context of epidemics on contact networks we use this to show that the approximate nature of some deterministic models such as mean-field and pair-approximation models can be partly understood by the identification of implicit anomalous terms. These terms describe unbiological processes which can be systematically removed up to and including the nth order by nth order moment closure approximations. These terms lead to a detailed understanding of the correlations in network-based epidemic models and contribute to understanding the connection between individual-level epidemic processes and population-level models. The connection with metapopulation models is also discussed. Our analysis is predominantly made at the individual level where the first and second order moment closure models correspond to what we term the individual-based and pair-based deterministic models, respectively. Matlab code is included as supplementary material for solving these models on transmission networks of arbitrary complexity.
Collapse
Affiliation(s)
- Kieran J Sharkey
- Department of Mathematical Sciences, The University of Liverpool, Peach Street, Liverpool, L69 7ZL, United Kingdom.
| |
Collapse
|
30
|
Rozhnova G, Nunes A. Cluster approximations for infection dynamics on random networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:051915. [PMID: 20365014 DOI: 10.1103/physreve.80.051915] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/22/2009] [Indexed: 05/29/2023]
Abstract
In this paper, we consider a simple stochastic epidemic model on large regular random graphs and the stochastic process that corresponds to this dynamics in the standard pair approximation. Using the fact that the nodes of a pair are unlikely to share neighbors, we derive the master equation for this process and obtain from the system size expansion the power spectrum of the fluctuations in the quasistationary state. We show that whenever the pair approximation deterministic equations give an accurate description of the behavior of the system in the thermodynamic limit, the power spectrum of the fluctuations measured in long simulations is well approximated by the analytical power spectrum. If this assumption breaks down, then the cluster approximation must be carried out beyond the level of pairs. We construct an uncorrelated triplet approximation that captures the behavior of the system in a region of parameter space where the pair approximation fails to give a good quantitative or even qualitative agreement. For these parameter values, the power spectrum of the fluctuations in finite systems can be computed analytically from the master equation of the corresponding stochastic process.
Collapse
Affiliation(s)
- G Rozhnova
- Centro de Física Teórica e Computacional and Departamento de Física, Faculdade de Ciências da Universidade de Lisboa, P-1649-003 Lisboa Codex, Portugal
| | | |
Collapse
|
31
|
A Motif-Based Approach to Network Epidemics. Bull Math Biol 2009; 71:1693-706. [DOI: 10.1007/s11538-009-9420-z] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2008] [Accepted: 04/02/2009] [Indexed: 10/20/2022]
|
32
|
Long-range correlations improve understanding of the influence of network structure on contact dynamics. Theor Popul Biol 2008; 73:383-94. [DOI: 10.1016/j.tpb.2007.12.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2006] [Revised: 10/31/2007] [Accepted: 12/11/2007] [Indexed: 11/22/2022]
|
33
|
Trapman P. Reproduction numbers for epidemics on networks using pair approximation. Math Biosci 2007; 210:464-89. [PMID: 17681553 DOI: 10.1016/j.mbs.2007.05.011] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2006] [Revised: 05/08/2007] [Accepted: 05/29/2007] [Indexed: 11/19/2022]
Abstract
One way to describe the spread of an infection on a network is by using the method of pair approximation. This method is a deterministic pair-based variant of the usual methods used to describe the progress of an epidemic in randomly mixing populations. However, although the ideas of pair approximation are intuitively clear, it is not straightforward to make all assumptions used explicit. Furthermore, in literature problems arise in defining basic quantities like the basic reproduction number R(0) and the real-time epidemic growth rate parameter r. We formulate the pair approximations and the needed assumptions explicitly. We discuss problems inherent to this method. Furthermore, we define a new reproduction number, similar to R(0) and a new real-time growth rate parameter similar to r. We illustrate the methods of the paper by an example for which we can compare the approximation of the reproduction number with exact results.
Collapse
Affiliation(s)
- Pieter Trapman
- Faculty of Veterinary Medicine, Utrecht University, The Netherlands.
| |
Collapse
|