Tooling-up for infectious disease transmission modelling

JAGS model specification for spatiotemporal epidemiological modelling

Bayesian inference in modelling infectious diseases using Bayesian inference using Gibbs Sampling (BUGS) is notable in the last two decades in parallel with the advancements in computing and model development. The ability of BUGS to easily implement the Markov chain Monte Carlo (MCMC) method brought Bayesian analysis to the mainstream of infectious disease modelling. However, with the existing software that runs MCMC to make Bayesian inferences, it is challenging, especially in terms of computational complexity, when infectious disease models become more complex with spatial and temporal components, in addition to the increasing number of parameters and large datasets. This study investigates two alternative subscripting strategies for creating models in Just Another Gibbs Sampler (JAGS) environment and their performance in terms of run times. Our results are useful for practitioners to ensure the efficiency and timely implementation of Bayesian spatiotemporal infectious disease modelling.

Modeling COVID-19 spread using multi-agent simulation with small-world network approach

BACKGROUND

The rapid global spread of COVID-19 has seriously impacted people's daily lives and the social economy while also posing a threat to their lives. The analysis of infectious disease transmission is of significant importance for the rational allocation of epidemic prevention and control resources, the management of public health emergencies, and the improvement of future public health systems.

METHODS

We propose a spatiotemporal COVID-19 transmission model with a neighborhood as an agent unit and an urban spatial network with long and short edge connections. The spreading model includes a network of defined agent attributes, transformation rules, and social relations and a small world network representing agents' social relations. Parameters for each stage are fitted by the Runge-Kutta method combined with the SEIR model. Using the NetLogo development platform, accurate dynamic simulations of the spatial and temporal evolution of the early epidemic were achieved.

RESULTS

Experimental results demonstrate that the fitted curves from the four stages agree with actual data, with only a 12.27% difference between the average number of infected agents and the actual number of infected agents after simulating 1 hundred times. Additionally, the model simulates and compares different "city closure" scenarios. The results showed that implementing a 'lockdown' 10 days earlier would lead to the peak number of infections occurring 7 days earlier than in the normal scenario, with a reduction of 40.35% in the total number of infections.

DISCUSSION

Our methodology emphasizes the crucial role of timely epidemic interventions in curbing the spread of infectious diseases, notably in the predictive assessment and evaluation of lockdown strategies. Furthermore, this approach adeptly forecasts the influence of varying intervention timings on peak infection rates and total case numbers, accurately reflecting real-world virus transmission patterns. This highlights the importance of proactive measures in diminishing epidemic impacts. It furnishes a robust framework, empowering policymakers to refine epidemic response strategies based on a synthesis of predictive modeling and empirical data.

Scenario prediction of public health emergencies using infectious disease dynamics model and dynamic Bayes

This study was aimed to discuss the predictive value of infectious disease dynamics model (IDD model) and dynamic Bayesian network (DBN) for scenario deduction of public health emergencies (PHEs). Based on the evolution law of PHEs and the meta-scenario representation of basic knowledge, this study established a DBN scenario deduction model for scenario deduction and evolution path analysis of PHEs. At the same time, based on the average field dynamics model of the SIR network, the dimensionality reduction process was performed to calculate the epidemic scale and epidemic time based on the IDD model, so as to determine the calculation methods of threshold value and epidemic time under emergency measures (quarantine). The Corona Virus Disease (COVID) epidemic was undertaken as an example to analyze the results of DBN scenario deduction, and the infectious disease dynamics model was used to analyze the number of reproductive numbers, peak arrival time, epidemic time, and latency time of the COVID epidemic. It was found that after the M1 measure was used to process the S1 state, the state probability and the probability of being true (T) were the highest, which were 91.05 and 90.21, respectively. In the sixth stage of the development of the epidemic, the epidemic had developed to level 5, the number of infected people was about 26, and the estimated loss was about 220 million yuan. The comprehensive cumulative foreground (CF) values of O1  ∼  O3 schemes were -1.34, -1.21, and -0.77, respectively, and the final CF values were -1.35, 0.01, and -0.08, respectively. The final CF value of O2 was significantly higher than the other two options. The household infection probability was the highest, which was 0.37 and 0.35 in Wuhan and China, respectively. Under the measures of home quarantine, the numbers of confirmed cases of COVID in China and Wuhan were 1.503 (95% confidential interval (CI) = 1.328  ∼  1.518) and 1.729 (95% CI = 1.107  ∼  1.264), respectively, showing good fits with the real data. On the 21st day after the quarantine measures were taken, the number of COVID across the country had an obvious peak, with the confirmed cases of 24495, and the model prediction value was 24085 (95% CI = 23988  ∼  25056). The incubation period 1/q was shortened from 8 days to 3 days, and the number of confirmed cases showed an upward trend. The peak period of confirmed cases was advanced, shortening the overall epidemic time. It showed that the prediction results of scenario deduction based on DBN were basically consistent with the actual development scenario and development status of the epidemic. It could provide corresponding decisions for the prevention and control of COVID based on the relevant parameters of the infectious disease dynamic model, which verified the rationality and feasibility of the scenario deduction method proposed in this study.

Using observed incidence to calibrate the transmission level of a mathematical model for Plasmodium vivax dynamics including case management and importation

In this work, we present a simple and flexible model for Plasmodium vivax dynamics which can be easily combined with routinely collected data on local and imported case counts to quantify transmission intensity and simulate control strategies. This model extends the model from White et al. (2016) by including case management interventions targeting liver-stage or blood-stage parasites, as well as imported infections. The endemic steady state of the model is used to derive a relationship between the observed incidence and the transmission rate in order to calculate reproduction numbers and simulate intervention scenarios. To illustrate its potential applications, the model is used to calculate local reproduction numbers in Panama and identify areas of sustained malaria transmission that should be targeted by control interventions.

A unifying nonlinear probabilistic epidemic model in space and time

Covid-19 epidemic dramatically relaunched the importance of mathematical modelling in supporting governments decisions to slow down the disease propagation. On the other hand, it remains a challenging task for mathematical modelling. The interplay between different models could be a key element in the modelling strategies. Here we propose a continuous space-time non-linear probabilistic model from which we can derive many of the existing models both deterministic and stochastic as for example SI, SIR, SIR stochastic, continuous-time stochastic models, discrete stochastic models, Fisher-Kolmogorov model. A partial analogy with the statistical interpretation of quantum mechanics provides an interpretation of the model. Epidemic forecasting is one of its possible applications; in principle, the model can be used in order to locate those regions of space where the infection probability is going to increase. The connection between non-linear probabilistic and non-linear deterministic models is analyzed. In particular, it is shown that the Fisher-Kolmogorov equation is connected to linear probabilistic models. On the other hand, a generalized version of the Fisher-Kolmogorov equation is derived from the non-linear probabilistic model and is shown to be characterized by a non-homogeneous time-dependent diffusion coefficient (anomalous diffusion) which encodes information about the non-linearity of the probabilistic model.

Modeling and Simulation of Dissemination of Cultivated Land Protection Policies in China

Cultivated land protection is the top priority of the national economy in China and the livelihood of people. Cultivated land protection policies (CLPP) play an important role in the protection of cultivated land. However, the process of dissemination of CLPP on social networks of farmers has problems, such as distortion of policy content, single dissemination channels, low level of farmers’ knowledge, and low dissemination efficiency. For revealing the characteristics of the dissemination of CLPP in the farmers’ social networks (FSN), this study combines the Suspected–Exposed–Infected–Recovered–Suspected (SEIRS) epidemic model to construct a model of CLPP dissemination suitable for FSN. In addition, a numerical simulation of the dissemination process of CLPP is conducted on the FSN, and the influence of the structural characteristics of the FSN and different model parameters on the dissemination of CLPP is analyzed. Results show that (1) the dissemination rate between farmers in FSN has a significant impact on the scale and speed of CLPP. A greater initial dissemination rate corresponds to faster speed and larger scale of CLPP dissemination. (2) A greater node degree in FSN means stronger dissemination ability for CLPP. Therefore, identifying structural holes (opinion leaders) in FSN can effectively promote the dissemination of CLPP. (3) The SEIRS model can dynamically describe the evolution law of CLPP dissemination process over time through the four states of farmer nodes of suspected, exposed, infected, and recovered. Numerical simulation results show that the immune degradation rate is proportional to CLPP. However, the direct immunization rate is inversely proportional. The increase in immune degradation rate can reduce the number of recovered farmers and improve the efficiency of CLPP dissemination. On the basis of the abovementioned conclusions, this study draws policy recommendations to increase the scale and speed of CLPP dissemination in China.