Numerical Integration of the Master Equation in Some Models of Stochastic Epidemiology

The processes by which disease spreads in a population of individuals are inherently stochastic. The master equation has proven to be a useful tool for modeling such processes. Unfortunately, solving the master equation analytically is possible only in limited cases (e.g., when the model is linear), and thus numerical procedures or approximation methods must be employed. Available approximation methods, such as the system size expansion method of van Kampen, may fail to provide reliable solutions, whereas current numerical approaches can induce appreciable computational cost. In this paper, we propose a new numerical technique for solving the master equation. Our method is based on a more informative stochastic process than the population process commonly used in the literature. By exploiting the structure of the master equation governing this process, we develop a novel technique for calculating the exact solution of the master equation – up to a desired precision – in certain models of stochastic epidemiology. We demonstrate the potential of our method by solving the master equation associated with the stochastic SIR epidemic model. MATLAB software that implements the methods discussed in this paper is freely available as Supporting Information S1

Modelling the Proportion of Influenza Infections within Households during Pandemic and Non-Pandemic Years

Background: The key epidemiological difference between pandemic and seasonal influenza is that the population is largely susceptible during a pandemic, whereas, during non-pandemic seasons a level of immunity exists. The population-level efficacy of household-based mitigation strategies depends on the proportion of infections that occur within households. In general, mitigation measures such as isolation and quarantine are more effective at the population level if the proportion of household transmission is low. Methods/Results: We calculated the proportion of infections within households during pandemic years compared with non-pandemic years using a deterministic model of household transmission in which all combinations of household size and individual infection states were enumerated explicitly. We found that the proportion of infections that occur within households was only partially influenced by the hazard h of infection within household relative to the hazard of infection outside the household, especially for small basic reproductive numbers. During pandemics, the number of within-household infections was lower than one might expect for a given h because many of the susceptible individuals were infected from the community and the number of susceptible individuals within household was thus depleted rapidly. In addition, we found that for the value of h at which 30% of infections occur within households during non-pandemic years, a similar 31% of infections occur within households during pandemic years. Interpretation: We suggest that a trade off between the community force of infection and the number of susceptible individuals in a household explains an apparent invariance in the proportion of infections that occur in households in our model. During a pandemic, although there are more susceptible individuals in a household, the community force of infection is very high. However, during non-pandemic years, the force of infection is much lower but there are fewer susceptible individuals within the household. © 2011 Kwok et al.published_or_final_versio

A critical review of mathematical models and data used in diabetology

The literature dealing with mathematical modelling for diabetes is abundant. During the last decades, a variety of models have been devoted to different aspects of diabetes, including glucose and insulin dynamics, management and complications prevention, cost and cost-effectiveness of strategies and epidemiology of diabetes in general. Several reviews are published regularly on mathematical models used for specific aspects of diabetes. In the present paper we propose a global overview of mathematical models dealing with many aspects of diabetes and using various tools. The review includes, side by side, models which are simple and/or comprehensive; deterministic and/or stochastic; continuous and/or discrete; using ordinary differential equations, partial differential equations, optimal control theory, integral equations, matrix analysis and computer algorithms

Spreading of infection on temporal networks: an edge-centered perspective

We discuss a continuous-time extension of the contact-based (CB) model, as proposed in [Koher et al. Phys. Rev. X 9, 031017 (2019)], for infections with permanent immunity on temporal networks. At the core of our methodology is a fundamental change to an edge-centered perspective, which allows for an accurate model on temporal networks, where the underlying time-aggregated graph has a tree structure. From the continuous-time CB model, we derive the infection propagator for the low prevalence limit and propose a novel spectral criterion to estimate the epidemic threshold. In addition, we explore the relation between the continuous-time CB model and the previously proposed edge-based compartmental model, as well as the message-passing framework