Network-based prediction of COVID-19 epidemic spreading in Italy

Inference on spatiotemporal dynamics for coupled biological populations

Mathematical models in ecology and epidemiology must be consistent with observed data in order to generate reliable knowledge and evidence-based policy. Metapopulation systems, which consist of a network of connected sub-populations, pose technical challenges in statistical inference owing to nonlinear, stochastic interactions. Numerical difficulties encountered in conducting inference can obstruct the core scientific questions concerning the link between the mathematical models and the data. Recently, an algorithm has been proposed that enables computationally tractable likelihood-based inference for high-dimensional partially observed stochastic dynamic models of metapopulation systems. We use this algorithm to build a statistically principled data analysis workflow for metapopulation systems. Via a case study of COVID-19, we show how this workflow addresses the limitations of previous approaches. The COVID-19 pandemic provides a situation where mathematical models and their policy implications are widely visible, and we revisit an influential metapopulation model used to inform basic epidemiological understanding early in the pandemic. Our methods support self-critical data analysis, enabling us to identify and address model weaknesses, leading to a new model with substantially improved statistical fit and parameter identifiability. Our results suggest that the lockdown initiated on 23 January 2020 in China was more effective than previously thought.

Quality assessment and community detection methods for anonymized mobility data in the Italian Covid context

We discuss how to assess the reliability of partial, anonymized mobility data and compare two different methods to identify spatial communities based on movements: Greedy Modularity Clustering (GMC) and the novel Critical Variable Selection (CVS). These capture different aspects of mobility: direct population fluxes (GMC) and the probability for individuals to move between two nodes (CVS). As a test case, we consider movements of Italians before and during the SARS-Cov2 pandemic, using Facebook users' data and publicly available information from the Italian National Institute of Statistics (Istat) to construct daily mobility networks at the interprovincial level. Using the Perron-Frobenius (PF) theorem, we show how the mean stochastic network has a stationary population density state comparable with data from Istat, and how this ceases to be the case if even a moderate amount of pruning is applied to the network. We then identify the first two national lockdowns through temporal clustering of the mobility networks, define two representative graphs for the lockdown and non-lockdown conditions and perform optimal spatial community identification on both graphs using the GMC and CVS approaches. Despite the fundamental differences in the methods, the variation of information (VI) between them assesses that they return similar partitions of the Italian provincial networks in both situations. The information provided can be used to inform policy, for example, to define an optimal scale for lockdown measures. Our approach is general and can be applied to other countries or geographical scales.

Emergent encoding of dispersal network topologies in spatial metapopulation models

We address a generalization of the concept of metapopulation capacity for trees and networks acting as the template for ecological interactions. The original measure had been derived from an insightful phenomenological model and is based on the leading eigenvalue of a suitable landscape matrix. It yields a versatile predictor of metapopulation persistence through a threshold value of the eigenvalue determined by ecological features of the focal species. Here, we present an analytical solution to a fundamental microscopic model that incorporates key ingredients of metapopulation dynamics and explicitly distinguishes between individuals comprising the "settled population" and "explorers" seeking colonization. Our approach accounts for general network characteristics (in particular graph-driven directional dispersal which is known to significantly constrain many ecological estimates) and yields a generalized version of the original model, to which it reduces for particular cases. Through examples, including real landscapes used as the template, we compare the predictions from our approach with those of the standard model. Results suggest that in several cases of practical interest, differences are significant. We also examine, with both models, how changes in habitat fragmentation, including removal, addition, or alteration in size, affect metapopulation persistence. The current approach demonstrates a high level of flexibility, enabling the incorporation of diverse "microscopic" elements and their impact on the resulting biodiversity landscape pattern.

Modeling Consequences of COVID-19 and Assessing Its Epidemiological Parameters: A System Dynamics Approach

In 2020, coronavirus (COVID-19) was declared a global pandemic and it remains prevalent today. A necessity to model the transmission of the virus has emerged as a result of COVID-19's exceedingly contagious characteristics and its rapid propagation throughout the world. Assessing the incidence of infection could enable policymakers to identify measures to halt the pandemic and gauge the required capacity of healthcare centers. Therefore, modeling the susceptibility, exposure, infection, and recovery in relation to the COVID-19 pandemic is crucial for the adoption of interventions by regulatory authorities. Fundamental factors, such as the infection rate, mortality rate, and recovery rate, must be considered in order to accurately represent the behavior of the pandemic using mathematical models. The difficulty in creating a mathematical model is in identifying the real model variables. Parameters might vary significantly across models, which can result in variations in the simulation results because projections primarily rely on a particular dataset. The purpose of this work was to establish a susceptible-exposed-infected-recovered (SEIR) model describing the propagation of the COVID-19 outbreak throughout the Kingdom of Saudi Arabia (KSA). The goal of this study was to derive the essential COVID-19 epidemiological factors from actual data. System dynamics modeling and design of experiment approaches were used to determine the most appropriate combination of epidemiological parameters and the influence of COVID-19. This study investigates how epidemiological variables such as seasonal amplitude, social awareness impact, and waning time can be adapted to correctly estimate COVID-19 scenarios such as the number of infected persons on a daily basis in KSA. This model can also be utilized to ascertain how stress (or hospital capacity) affects the percentage of hospitalizations and the number of deaths. Additionally, the results of this study can be used to establish policies or strategies for monitoring or restricting COVID-19 in Saudi Arabia.

Preventive control strategy on second wave of Covid-19 pandemic model incorporating lock-down effect

This study presents an optimal control strategy through a mathematical model of the Covid-19 outbreak without lock-down. The pandemic model analyses the lock-down effect without control strategy based on the current scenario of second wave data to control the rapid spread of the virus. The pandemic model has been discussed with respect to the basic reproduction number and stability analysis of disease-free and endemic equilibrium. A new optimal control problem with treatment is framed to minimize the vulnerable situation of the second wave. This system is applied to study the effects of vaccines and treatment controls. Numerical solutions and the graphical presentation of the results predict the fate of India’s second wave situation on account of the control strategy. Lastly, a comparative study with control and without control has been analysed for the exposed phase, infective phase, and recovery phase to understand the effectiveness of the controls. This model is used to estimate the total number of infected and active cases, deaths, and recoveries in order to control the disease using this system and studying the effects of vaccines and treatment controls.

Optimal control design incorporating vaccination and treatment on six compartment pandemic dynamical system

In this paper, a mathematical model of the COVID-19 pandemic with lockdown that provides a more accurate representation of the infection rate has been analyzed. In this model, the total population is divided into six compartments: the susceptible class, lockdown class, exposed class, asymptomatic infected class, symptomatic infected class, and recovered class. The basic reproduction number (R0) is calculated using the next-generation matrix method and presented graphically based on different progression rates and effective contact rates of infective individuals. The COVID-19 epidemic model exhibits the disease-free equilibrium and endemic equilibrium. The local and global stability analysis has been done at the disease-free and endemic equilibrium based on R0. The stability analysis of the model shows that the disease-free equilibrium is both locally and globally stable when R0<1, and the endemic equilibrium is locally and globally stable when R0>1 under some conditions. A control strategy including vaccination and treatment has been studied on this pandemic model with an objective functional to minimize. Finally, numerical simulation of the COVID-19 outbreak in India is carried out using MATLAB, highlighting the usefulness of the COVID-19 pandemic model and its mathematical analysis.

Prediction of COVID-19 Infected Population for Indian States through a State Interaction Network-based SEIR Epidemic Model

Objective of this present study is to predict the COVID-19 trajectories in terms of infected population of Indian states. In this work, a state interaction network of sixteen Indian states with highest number of infected caseload is considered, based on networked Susceptible-Exposed-Infected-Recovered (SEIR) epidemic model. An intervention term has been introduced in order to capture the effect of lockdown with different stringencies at different periods of time. The model has been fitted using least absolute shrinkage and selection operator (LASSO). Machine learning methods have been used to train the parameters of the model, cross-validate the data, and predict the parameters. The predictions of infected population for each of the sixteen states have been shown using data considered from January 1, 2021 till writing this manuscript on June 25, 2021. Finally, the effectiveness of the model is manifested by the calculated mean error and confidence interval.

Comparing the accuracy of several network-based COVID-19 prediction algorithms

Researchers from various scientific disciplines have attempted to forecast the spread of coronavirus disease 2019 (COVID-19). The proposed epidemic prediction methods range from basic curve fitting methods and traffic interaction models to machine-learning approaches. If we combine all these approaches, we obtain the Network Inference-based Prediction Algorithm (NIPA). In this paper, we analyse a diverse set of COVID-19 forecast algorithms, including several modifications of NIPA. Among the algorithms that we evaluated, the original NIPA performed best at forecasting the spread of COVID-19 in Hubei, China and in the Netherlands. In particular, we show that network-based forecasting is superior to any other forecasting algorithm.

Country-specific optimization strategy for testing through contact tracing can help maintain a low reproduction number ([Formula: see text]) during unlock

In response to the COVID19 pandemic, many countries have implemented lockdowns in multiple phases to ensure social distancing and quarantining of the infected subjects. Subsequent unlocks to reopen the economies started next waves of infection and imposed an extra burden on quarantine to keep the reproduction number ([Formula: see text]) < 1. However, most countries could not effectively contain the infection spread, suggesting identification of the potential sources weakening the effect of lockdowns could help design better informed lockdown-unlock cycles in the future. Here, through building quantitative epidemic models and analyzing the metadata of 50 countries from across the continents we first found that the estimated value of [Formula: see text], adjusted w.r.t the distribution of medical facilities and virus clades correlates strongly with the testing rates in a country. Since the testing capacity of a country is limited by its medical resources, we investigated if a cost-benefit trade-off can be designed connecting testing rate and extent of unlocking. We present a strategy to optimize this trade-off in a country specific manner by providing a quantitative estimate of testing and quarantine rates required to allow different extents of unlocks while aiming to maintain [Formula: see text]. We further show that a small fraction of superspreaders can dramatically increase the number of infected individuals even during strict lockdowns by strengthening the positive feedback loop driving infection spread. Harnessing the benefit of optimized country-specific testing rates would critically require minimizing the movement of these superspreaders via strict social distancing norms, such that the positive feedback driven switch-like exponential spread phase of infection can be avoided/delayed.

Dynamical analysis of novel COVID-19 epidemic model with non-monotonic incidence function

In this study, we developed and analyzed a mathematical model for explaining the transmission dynamics of COVID-19 in India. The proposed SI u I k R model is a modified version of the existing SIR model. Our model divides the infected class I of SIR model into two classes: I u (unknown infected class) and I k (known infected class). In addition, we consider R a recovered and reserved class, where susceptible people can hide them due to fear of the COVID-19 infection. Furthermore, a non-monotonic incidence function is deemed to incorporate the psychological effect of the novel coronavirus diseases on India's community. The epidemiological threshold parameter, namely the basic reproduction number, has been formulated and presented graphically. With this threshold parameter, the local and global stability analysis of the disease-free equilibrium and the endemic proportion equilibrium based on disease persistence have been analyzed. Lastly, numerical results of long-run prediction using MATLAB show that the fate of this situation is very harmful if people are not following the guidelines issued by the authority.