Estimate the number of lives saved by a SARS-CoV-2 vaccination campaign in six states in the United States with a simple model.

Objectives: Vaccination and the emergence of the highly transmissible Omicron variant changed the fate of the COVID-19 pandemic. It is very challenging to estimate the number of lives saved by vaccination given the multiple doses of vaccination, the time-varying nature of transmissibility, the waning of immunity, and the presence of immune evasion. Methods: We established a S-SV-E-I-T-D-R model to simulate the number of lives saved by vaccination in six states in the United States (U.S.) from March 5, 2020, to March 23, 2023. The cumulative number of deaths were estimated under three vaccination scenarios based on two assumptions. Additionally, immune evasion by the Omicron and loss of protection afforded by vaccination or infection were considered. Results: The number of deaths averted by COVID-19 vaccinations (including three doses) ranged from 0.154-0.295% of the total population across six states. The number of deaths averted by the third dose ranged from 0.008-0.017% of the total population. Conclusions: Our estimate of death averted by COVID-19 vaccination in the U.S. was largely in line with an official estimate (at a level of 0.15-0.20% of the total population). We found that the additional contribution of the third dose was small but significant.

Global stability and optimal control of an age-structured SVEIR epidemic model with waning immunity and relapses.

The efficacy of vaccination, incomplete treatment and disease relapse are critical challenges that must be faced to prevent and control the spread of infectious diseases. Age heterogeneity is also a crucial factor for this study. In this paper, we investigate a new age-structured SVEIR epidemic model with the nonlinear incidence rate, waning immunity, incomplete treatment and relapse. Next, the asymptotic smoothness, the uniform persistence and the existence of interior global attractor of the solution semi-flow generated by the system are given. We define the basic reproduction number R 0 and prove the existence of the equilibria of the model. And we study the global asymptotic stability of the equilibria. Then the parameters of the model are estimated using tuberculosis data in China. The sensitivity analysis of R 0 is derived by the Partial Rank Correlation Coefficient method. These main theoretical results are applied to analyze and predict the trend of tuberculosis prevalence in China. Finally, the optimal control problem of the model is discussed. We choose to take strengthening treatment and controlling relapse as the control parameters. The necessary condition for optimal control is established.

Source-Sink Dynamics in a Two-Patch SI Epidemic Model with Life Stages and No Recovery from Infection.

This study presents a comprehensive analysis of a two-patch, two-life stage SI model without recovery from infection, focusing on the dynamics of disease spread and host population viability in natural populations. The model, inspired by real-world ecological crises like the decline of amphibian populations due to chytridiomycosis and sea star populations due to Sea Star Wasting Disease, aims to understand the conditions under which a sink host population can present ecological rescue from a healthier, source population. Mathematical and numerical analyses reveal the critical roles of the basic reproductive numbers of the source and sink populations, the maturation rate, and the dispersal rate of juveniles in determining population outcomes. The study identifies basic reproduction numbers R 0 for each of the patches, and conditions for the basic reproduction numbers to produce a receiving patch under which its population. These findings provide insights into managing natural populations affected by disease, with implications for conservation strategies, such as the importance of maintaining reproductively viable refuge populations and considering the effects of dispersal and maturation rates on population recovery. The research underscores the complexity of host-pathogen dynamics in spatially structured environments and highlights the need for multi-faceted approaches to biodiversity conservation in the face of emerging diseases.

Memory impacts in hepatitis C: A global analysis of a fractional-order model with an effective treatment.

BACKGROUND AND OBJECTIVE: Hepatitis virus infections are affecting millions of people worldwide, causing death, disability, and considerable expenditure. Chronic infection with hepatitis C virus (HCV) can cause severe public health problems because of their high prevalence and poor long-term clinical outcomes. Thus a fractional-order epidemic model of the hepatitis C virus involving partial immunity under the influence of memory effect to know the transmission patterns and prevalence of HCV infection is studied. Investigating the transmission dynamics of HCV makes the issue more interesting. The HCV epidemic model and worldwide dynamics are examined in this study. Calculate the basic reproduction number for the HCV model using the next-generation matrix technique. We determine the model's global dynamics using reproduction numbers, the Lyapunov functional approach, and the Routh-Hurwitz criterion. The model's reproduction number shows how the disease progresses. METHODS: A fractional differential equation model of HCV infection has been created. Maximum relevant parameters, such as fractional power, HCV transmission rate, reproduction number, etc., influencing the dynamic process, have been incorporated. The model's numerical solutions are obtained using the fractional Adams method. Finally, numerical simulations support the theoretical conclusions, showing the great agreement between the two. RESULTS: In the fractional-order HCV infection model, the memory effect, which is not seen in the classical model, was shown on graphs so that disease dynamics and vector compartments could be seen. We found that the fractional-order HCV infection model has more stages of freedom than regular derivatives. Fractional-order derivations, which may be the best and most reliable, explained bodily approaches better than classical order. CONCLUSION: The current study modeled and analyzed a fractional-order HCV infection model. The current approach results in a much better understanding of HCV transmission in a population, which leads to important insights into its spread and control, such as better treatment dosage for different age groups, identifying the best control measure, improving health, prolonging life, reducing the risk of HCV transmission, and effectively increasing the quality of life of HCV patients. The creation of a fractional-order HCV infection model, which provides a better understanding of HCV transmission dynamics and leads to significant insights for better treatment dosages, identification of optimal control measures, and ultimately improvement of the quality of life for HCV patients, is the study's major outcome.

Stochastic two-strain epidemic model with saturated incidence rates driven by Lévy noise.

In this paper, we introduce a stochastic two-strain epidemic model driven by Lévy noise describing the interaction between four compartments; susceptible, infected individuals by the first strain, infected ones by the second strain and the recovered individuals. The forces of infection, for both strains, are represented by saturated incidence rates. Our study begins with the investigation of unique global solution of the suggested mathematical model. Then, it moves to the determination of sufficient conditions of extinction and persistence in mean of the two-strain disease. In order to illustrate the theoretical findings, we give some numerical simulations.

Analysis of a COVID-19 model with media coverage and limited resources.

The novel coronavirus disease (COVID-19) pandemic has profoundly impacted the global economy and human health. The paper mainly proposed an improved susceptible-exposed-infected-recovered (SEIR) epidemic model with media coverage and limited medical resources to investigate the spread of COVID-19. We proved the positivity and boundedness of the solution. The existence and local asymptotically stability of equilibria were studied and a sufficient criterion was established for backward bifurcation. Further, we applied the proposed model to study the trend of COVID-19 in Shanghai, China, from March to April 2022. The results showed sensitivity analysis, bifurcation, and the effects of critical parameters in the COVID-19 model.

Threshold dynamics of a switching diffusion SIR model with logistic growth and healthcare resources.

In this article, we have constructed a stochastic SIR model with healthcare resources and logistic growth, aiming to explore the effect of random environment and healthcare resources on disease transmission dynamics. We have showed that under mild extra conditions, there exists a critical parameter, i.e., the basic reproduction number $ R_0/ $, which completely determines the dynamics of disease: when $ R_0/ < 1 $, the disease is eradicated; while when $ R_0/ > 1 $, the disease is persistent. To validate our theoretical findings, we conducted some numerical simulations using actual parameter values of COVID-19. Both our theoretical and simulation results indicated that (1) the white noise can significantly affect the dynamics of a disease, and importantly, it can shift the stability of the disease-free equilibrium; (2) infectious disease resurgence may be caused by random switching of the environment; and (3) it is vital to maintain adequate healthcare resources to control the spread of disease.

An SEIR network epidemic model with manual and digital contact tracing allowing delays.

We consider an SEIR epidemic model on a network also allowing random contacts, where recovered individuals could either recover naturally or be diagnosed. Upon diagnosis, manual contact tracing is triggered such that each infected network contact is reported, tested and isolated with some probability and after a random delay. Additionally, digital tracing (based on a tracing app) is triggered if the diagnosed individual is an app-user, and then all of its app-using infectees are immediately notified and isolated. The early phase of the epidemic with manual and/or digital tracing is approximated by different multi-type branching processes, and three respective reproduction numbers are derived. The effectiveness of both contact tracing mechanisms is numerically quantified through the reduction of the reproduction number. This shows that app-using fraction plays an essential role in the overall effectiveness of contact tracing. The relative effectiveness of manual tracing compared to digital tracing increases if: more of the transmission occurs on the network, when the tracing delay is shortened, and when the network degree distribution is heavy-tailed. For realistic values, the combined tracing case can reduce R0 by 20%-30%, so other preventive measures are needed to reduce the reproduction number down to 1.2-1.4 for contact tracing to make it successful in avoiding big outbreaks.

Optimal control strategies for toxoplasmosis disease transmission dynamics via harmonic mean-type incident rate.

Toxoplasma infection in humans is considered due to direct contact with infected cats. Toxoplasma infection (an endemic disease) has the potential to affect various organs and systems (brain, eyes, heart, lungs, liver, and lymph nodes). Bilinear incidence rate and constant population (birth rate is equal to death rate) are used in the literature to explain the dynamics of Toxoplasmosis disease transmission in humans and cats. The goal of this study is to consider the mathematical model of Toxoplasma disease with harmonic mean type incident rate and also consider that the population of humans and cats is not equal (birth rate and the death rate are not equal). In examining Toxoplasma transmission dynamics in humans and cats, harmonic mean incidence rates are better than bilinear incidence rates. The disease dynamics are first schematically illustrated, and then the law of mass action is applied to obtain nonlinear ordinary differential equations (ODEs). Analysis of the boundedness, positivity, and equilibrium points of the system has been analyzed. The reproduction number is calculated using the next-generation matrix technique. The stability of disease-free and endemic equilibrium are analyzed. Sensitivity analysis is also done for reproduction number. Numerical simulation shows that the infection is spread in the population when the contact rate ß h and ß c increases while the infection is reduced when the recovery rate δ h increases. This study investigates the impact of various optimal control strategies, such as vaccinations for the control of disease and the awareness of disease awareness, on the management of disease.

A risk-induced dispersal strategy of the infected population for a disease-free state in the SIS epidemic model.

This article proposes a dispersal strategy for infected individuals in a spatial susceptible-infected-susceptible (SIS) epidemic model. The presence of spatial heterogeneity and the movement of individuals play crucial roles in determining the persistence and eradication of infectious diseases. To capture these dynamics, we introduce a moving strategy called risk-induced dispersal (RID) for infected individuals in a continuous-time patch model of the SIS epidemic. First, we establish a continuous-time n-patch model and verify that the RID strategy is an effective approach for attaining a disease-free state. This is substantiated through simulations conducted on 7-patch models and analytical results derived from 2-patch models. Second, we extend our analysis by adapting the patch model into a diffusive epidemic model. This extension allows us to explore further the impact of the RID movement strategy on disease transmission and control. We validate our results through simulations, which provide the effects of the RID dispersal strategy.

Infection-induced increases to population size during cycles in a discrete-time epidemic model.

One-dimensional discrete-time population models, such as those that involve Logistic or Ricker growth, can exhibit periodic and chaotic dynamics. Expanding the system by one dimension to incorporate epidemiological interactions causes an interesting complexity of new behaviors. Here, we examine a discrete-time two-dimensional susceptible-infectious (SI) model with Ricker growth and show that the introduction of infection can not only produce a distinctly different bifurcation structure than that of the underlying disease-free system but also lead to counter-intuitive increases in population size. We use numerical bifurcation analysis to determine the influence of infection on the location and types of bifurcations. In addition, we examine the appearance and extent of a phenomenon known as the 'hydra effect,' i.e., increases in total population size when factors, such as mortality, that act negatively on a population, are increased. Previous work, primarily focused on dynamics at fixed points, showed that the introduction of infection that reduces fecundity to the SI model can lead to a so-called 'infection-induced hydra effect.' Our work shows that even in such a simple two-dimensional SI model, the introduction of infection that alters fecundity or mortality can produce dynamics can lead to the appearance of a hydra effect, particularly when the disease-free population is at a cycle.

Influence of waning immunity on vaccination decision-making: A multi-strain epidemic model with an evolutionary approach analyzing cost and efficacy.

In this research, we introduce a comprehensive epidemiological model that accounts for multiple strains of an infectious disease and two distinct vaccination options. Vaccination stands out as the most effective means to prevent and manage infectious diseases. However, when there are various vaccines available, each with its costs and effectiveness, the decision-making process for individuals becomes paramount. Furthermore, the factor of waning immunity following vaccination also plays a significant role in influencing these choices. To understand how individuals make decisions in the context of multiple strains and waning immunity, we employ a behavioral model, allowing an epidemiological model to be coupled with the dynamics of a decision-making process. Individuals base their choice of vaccination on factors such as the total number of infected individuals and the cost-effectiveness of the vaccine. Our findings indicate that as waning immunity increases, people tend to prioritize vaccines with higher costs and greater efficacy. Moreover, when more contagious strains are present, the equilibrium in vaccine adoption is reached more rapidly. Finally, we delve into the social dilemma inherent in our model by quantifying the social efficiency deficit (SED) under various parameter combinations.

Comprehensive analysis of a stochastic wireless sensor network motivated by Black-Karasinski process.

Wireless sensor networks (WSNs) encounter a significant challenge in ensuring network security due to their operational constraints. This challenge stems from the potential infiltration of malware into WSNs, where a single infected node can rapidly propagate worms to neighboring nodes. To address this issue, this research introduces a stochastic S E I R S model to characterize worm spread in WSNs. Initially, we established that our model possesses a globally positive solution. Subsequently, we determine a threshold value for our stochastic system and derive a set of sufficient conditions that dictate the persistence or extinction of worm spread in WSNs based on the mean behavior. Our study reveals that environmental randomness can impede the spread of malware in WSNs. Moreover, by utilizing various parameter sets, we obtain approximate solutions that showcase these precise findings and validate the effectiveness of the proposed S E I R S model, which surpasses existing models in mitigating worm transmission in WSNs.

Dynamical Analysis of an Improved Bidirectional Immunization SIR Model in Complex Network.

In order to investigate the impact of two immunization strategies-vaccination targeting susceptible individuals to reduce their infection rate and clinical medical interventions targeting infected individuals to enhance their recovery rate-on the spread of infectious diseases in complex networks, this study proposes a bilinear SIR infectious disease model that considers bidirectional immunization. By analyzing the conditions for the existence of endemic equilibrium points, we derive the basic reproduction numbers and outbreak thresholds for both homogeneous and heterogeneous networks. The epidemic model is then reconstructed and extensively analyzed using continuous-time Markov chain (CTMC) methods. This analysis includes the investigation of transition probabilities, transition rate matrices, steady-state distributions, and the transition probability matrix based on the embedded chain. In numerical simulations, a notable concordance exists between the outcomes of CTMC and mean-field (MF) simulations, thereby substantiating the efficacy of the CTMC model. Moreover, the CTMC-based model adeptly captures the inherent stochastic fluctuation in the disease transmission, which is consistent with the mathematical properties of Markov chains. We further analyze the relationship between the system's steady-state infection density and the immunization rate through MCS. The results suggest that the infection density decreases with an increase in the immunization rate among susceptible individuals. The current research results will enhance our understanding of infectious disease transmission patterns in real-world scenarios, providing valuable theoretical insights for the development of epidemic prevention and control strategies.

Analysis and comparative study of a deterministic mathematical model of SARS-COV-2 with fractal-fractional operators: a case study.

In this paper, we investigate a fractal-fractional-order mathematical model with the influence of hospitalized patients and the impact of vaccination with fractal-fractional operators. The respective derivatives are considered in the Caputo, Caputo Fabrizio, and Atangana-Baleanu senses of fractional order α and fractal dimension τ . For the proposed problem, some results regarding basic reproduction number and stability are given. Using the next-generation matrix approach, we have investigated the global and local stability of several types of equilibrium points. We provide a detailed analysis of the existence and uniqueness of the solution. Moreover, we fit the model with the real data of Pakistan from June 01, 2020, till March 24, 2021. Then, we use the fractal-fractional derivative to find a numerical solution for the model. MATLAB software is used for numerical illustration. Graphical presentations corresponding to different parameteric values are given as well.

Final epidemic size of a two-community SIR model with asymmetric coupling.

Communities are commonly not isolated but interact asymmetrically with each other, allowing the propagation of infectious diseases within the same community and between different communities. To reveal the impact of asymmetrical interactions and contact heterogeneity on disease transmission, we formulate a two-community SIR epidemic model, in which each community has its contact structure while communication between communities occurs through temporary commuters. We derive an explicit formula for the basic reproduction number R 0 , give an implicit equation for the final epidemic size z, and analyze the relationship between them. Unlike the typical positive correlation between R 0 and z in the classic SIR model, we find a negatively correlated relationship between counterparts of our model deviating from homogeneous populations. Moreover, we investigate the impact of asymmetric coupling mechanisms on R 0 . The results suggest that, in scenarios with restricted movement of susceptible individuals within a community, R 0 does not follow a simple monotonous relationship, indicating that an unbending decrease in the movement of susceptible individuals may increase R 0 . We further demonstrate that network contacts within communities have a greater effect on R 0 than casual contacts between communities. Finally, we develop an epidemic model without restriction on the movement of susceptible individuals, and the numerical simulations suggest that the increase in human flow between communities leads to a larger R 0 .

Delay epidemic models determined by latency, infection, and immunity duration.

We propose new single and two-strain epidemic models represented by systems of delay differential equations and based on the number of newly exposed individuals. Transitions between exposed, infectious, recovered, and back to susceptible compartments are determined by the corresponding time delays. Existence and positiveness of solutions are proved. Reduction of delay differential equations to integral equations allows the analysis of stationary solutions and their stability. In the case of two strains, they compete with each other, and the strain with a larger individual basic reproduction number dominates the other one. However, if the basic reproduction number exceeds some critical values, stationary solution loses its stability resulting in periodic time oscillations. In this case, both strains are present and their dynamics is not completely determined by the basic reproduction numbers but also by other parameters. The results of the work are illustrated by comparison with data on seasonal influenza.

On cognitive epidemic models: spatial segregation versus nonpharmaceutical interventions.

Fick's law and the Fokker-Planck law of diffusion are applied to manifest the cognitive dispersal of individuals in two reaction-diffusion SEIR epidemic models, where the disease transmission is illustrated by nonlocal infection mechanisms in heterogeneous environments. Building upon the well-posedness of solutions, threshold dynamics are discussed in terms of the basic reproduction numbers for the two cognitive epidemic models. The numerical investigation reveals that the Fokker-Planck law can better describe the diffusion of individuals by taking different dispersal strategies of exposed individuals in our cognitive epidemic models, and provides some insights on spatial segregation and nonpharmaceutical interventions: (i) spatial segregation occurs in the random diffusion model when the nonlocal infection radius is small, while it appears in the symmetric diffusion model when the radius is large; (ii) nonpharmaceutical interventions on restricting the dispersal of exposed and infected individuals do not contribute to reducing the infection proportion, but rather eliminate the disease in a region, which expands as the nonlocal infection radius increases. We additionally find that the final infection size in the random diffusion model is significantly smaller than that in the symmetric diffusion model and decreases as the nonlocal infection radius increases.

SubEpiPredict: A tutorial-based primer and toolbox for fitting and forecasting growth trajectories using the ensemble n-sub-epidemic modeling framework.

An ensemble n-sub-epidemic modeling framework that integrates sub-epidemics to capture complex temporal dynamics has demonstrated powerful forecasting capability in previous works. This modeling framework can characterize complex epidemic patterns, including plateaus, epidemic resurgences, and epidemic waves characterized by multiple peaks of different sizes. In this tutorial paper, we introduce and illustrate SubEpiPredict, a user-friendly MATLAB toolbox for fitting and forecasting time series data using an ensemble n-sub-epidemic modeling framework. The toolbox can be used for model fitting, forecasting, and evaluation of model performance of the calibration and forecasting periods using metrics such as the weighted interval score (WIS). We also provide a detailed description of these methods including the concept of the n-sub-epidemic model, constructing ensemble forecasts from the top-ranking models, etc. For the illustration of the toolbox, we utilize publicly available daily COVID-19 death data at the national level for the United States. The MATLAB toolbox introduced in this paper can be very useful for a wider group of audiences, including policymakers, and can be easily utilized by those without extensive coding and modeling backgrounds.

Stochastic analysis and control for a delayed Hepatitis B epidemic model.

Considering the time delay originating from a certain incubation period or asymptomatic state, we propose a delayed epidemic system within the noisy environment of the hepatitis B virus to analyze the mechanism of disease transmission and elucidate how to control it by applying the strategy of vaccinating and treatment. Applying stochastic Lyapunov functional theory, we first construct an integral Lyapunov function coupling the time delay and stochastic fluctuation to investigate whether there exists a unique global solution to the model. Next, we yield the threshold condition for controlling disease extinction, and persistence, as well as its stationary distribution. Governed by these sufficient conditions, we study the existence of optimal control solutions in deterministic and stochastic scenarios to uncover how to accelerate disease extinction through vaccination and treatment. The results indicate that the time delay will prolong the duration of the disease for the original system but suppress the peak value of HBV in the controlled system. Finally, we verify the versatility of theoretical results by numerical simulations. These results will effectively decipher the importance of the time delay in the control of hepatitis B.