A Hypothesis-Free Bridging of Disease Dynamics and Non-pharmaceutical Policies

Based on epidemiological parameter data, probe into a stochastically perturbed dominant variant of the COVID-19 pandemic model

During the COVID-19 pandemic, the SARS-CoV-2 gene mutation has been rapidly emerging and spreading all over the world. Experts worldwide regularly monitor genetic mutations and variants through genome-sequence-based surveillance, laboratory testing, outbreak investigation, and epidemiological probing. Clinical pathologists and medical laboratory scientists prefer developing or endorsing COVID-19 vaccines with a broader immune response involving various antibodies and cells to protect against mutations or new variants. Randomness plays an enormous role in pathology and epidemiology. Hence, based on epidemiological parameter data, we construct and probe a stochastically perturbed dominant variant of the coronavirus epidemic model with three nonlinear saturated incidence rates. We reveal the existence of a unique global positive solution to the constructed stochastic COVID-19 model. The Lyapunov function method is used to determine the presence of a stationary distribution of positive solutions. We derive sufficient conditions for the coronavirus to be eradicated. Eventually, numerical simulations validate the effectiveness of our theoretical outcomes.

Spatial heterogeneity analysis for the transmission of syphilis disease in China via a data-validated reaction-diffusion model

Based on the distinctive spatial diffusion characteristics observed in syphilis transmission patterns, this paper introduces a novel reaction-diffusion model for syphilis disease dynamics, incorporating general incidence functions within a heterogeneous environment. We derive the basic reproduction number essential for threshold dynamics and investigate the uniform persistence of the model. We validate the model and estimate its parameters by employing the multi-objective Markov Chain Monte Carlo (MCMC) method, using real syphilis data from the years 2004 to 2018 in China. Furthermore, we explore the impact of spatial heterogeneity and intervention measures on syphilis transmission. Our findings reveal several key insights: (1) In addition to the original high-incidence areas of syphilis, Xinjiang, Guizhou, Hunan and Northeast China have also emerged as high-incidence regions for syphilis in China. (2) The latent syphilis cases represent the highest proportion of newly reported cases, highlighting the critical importance of considering their role in transmission dynamics to avoid underestimation of syphilis outbreaks. (3) Neglecting spatial heterogeneity results in an underestimation of disease prevalence and the number of syphilis-infected individuals, undermining effective disease prevention and control strategies. (4) The initial conditions have minimal impact on the long-term spatial distribution of syphilis-infected individuals in scenarios of varying diffusion rates. This study underscores the significance of spatial dynamics and intervention measures in assessing and managing syphilis transmission, which offers insights for public health policymakers.

From Policy to Prediction: Assessing Forecasting Accuracy in an Integrated Framework with Machine Learning and Disease Models

To improve the forecasting accuracy of the spread of infectious diseases, a hybrid model was recently introduced where the commonly assumed constant disease transmission rate was actively estimated from enforced mitigating policy data by a machine learning (ML) model and then fed to an extended susceptible-infected-recovered model to forecast the number of infected cases. Testing only one ML model, that is, gradient boosting model (GBM), the work left open whether other ML models would perform better. Here, we compared GBMs, linear regressions, k-nearest neighbors, and Bayesian networks (BNs) in forecasting the number of COVID-19-infected cases in the United States and Canadian provinces based on policy indices of future 35 days. There was no significant difference in the mean absolute percentage errors of these ML models over the combined dataset [H ( 3 ) = 3.10 , p = 0.38 ]. In two provinces, a significant difference was observed [H ( 3 ) = 8.77 , H ( 3 ) = 8.07 , p < 0.05 ], yet posthoc tests revealed no significant difference in pairwise comparisons. Nevertheless, BNs significantly outperformed the other models in most of the training datasets. The results put forward that the ML models have equal forecasting power overall, and BNs are best for data-fitting applications.

An early warning indicator trained on stochastic disease-spreading models with different noises

The timely detection of disease outbreaks through reliable early warning signals (EWSs) is indispensable for effective public health mitigation strategies. Nevertheless, the intricate dynamics of real-world disease spread, often influenced by diverse sources of noise and limited data in the early stages of outbreaks, pose a significant challenge in developing reliable EWSs, as the performance of existing indicators varies with extrinsic and intrinsic noises. Here, we address the challenge of modelling disease when the measurements are corrupted by additive white noise, multiplicative environmental noise and demographic noise into a standard epidemic mathematical model. To navigate the complexities introduced by these noise sources, we employ a deep learning algorithm that provides EWS in infectious disease outbreaks by training on noise-induced disease-spreading models. The indicator's effectiveness is demonstrated through its application to real-world COVID-19 cases in Edmonton and simulated time series derived from diverse disease spread models affected by noise. Notably, the indicator captures an impending transition in a time series of disease outbreaks and outperforms existing indicators. This study contributes to advancing early warning capabilities by addressing the intricate dynamics inherent in real-world disease spread, presenting a promising avenue for enhancing public health preparedness and response efforts.

A novel indicator in epidemic monitoring through a case study of Ebola in West Africa (2014-2016)

The E/S (exposed/susceptible) ratio is analyzed in the SEIR model. The ratio plays a key role in understanding epidemic dynamics during the 2014-2016 Ebola outbreak in Sierra Leone and Guinea. The maximum value of the ratio occurs immediately before or after the time-dependent reproduction number (Rt) equals 1, depending on the initial susceptible population (S(0)). It is demonstrated that transmission rate curves corresponding to various incubation periods intersect at a single point referred to as the Cross Point (CP). At this point, the E/S ratio reaches an extremum, signifying a critical shift in transmission dynamics and aligning with the time when Rt approaches 1. By plotting transmission rate curves, β(t), for any two arbitrary incubation periods and tracking their intersections, we can trace CP over time. CP serves as an indicator of epidemic status, especially when Rt is close to 1. It provides a practical means of monitoring epidemics without prior knowledge of the incubation period. Through a case study, we estimate the transmission rate and reproduction number, identifying CP and Rt = 1 while examining the E/S ratio across various values of S(0).

Climate-dependent effectiveness of nonpharmaceutical interventions on COVID-19 mitigation

Environmental factors have a significant impact on the transmission of infectious diseases. Existing results show that the novel coronavirus can persist outside the host. We propose a susceptible-exposed-presymptomatic-infectious-asymptomatic-recovered-susceptible (SEPIARS) model with a vaccination compartment and indirect incidence to explore the effect of environmental conditions, temperature and humidity, on the transmission of the SARS-CoV-2 virus. Using climate data and daily confirmed cases data in two Canadian cities with different atmospheric conditions, we evaluate the mortality rates of the SARS-CoV-2 virus and further estimate the transmission rates by the inverse method, respectively. The numerical results show that high temperature or humidity can be helpful in mitigating the spread of COVID-19 during the warm summer months. Our findings verify that nonpharmaceutical interventions are less effective if the virus can persist for a long time on surfaces. Based on climate data, we can forecast the transmission rate and the infection cases up to four weeks in the future by a generalized boosting machine learning model.

Assessing potential insights of an imperfect testing strategy: Parameter estimation and practical identifiability using early COVID-19 data in India

A deterministic model with testing of infected individuals has been proposed to investigate the potential consequences of the impact of testing strategy. The model exhibits global dynamics concerning the disease-free and a unique endemic equilibrium depending on the basic reproduction number when the recruitment of infected individuals is zero; otherwise, the model does not have a disease-free equilibrium, and disease never dies out in the community. Model parameters have been estimated using the maximum likelihood method with respect to the data of early COVID-19 outbreak in India. The practical identifiability analysis shows that the model parameters are estimated uniquely. The consequences of the testing rate for the weekly new cases of early COVID-19 data in India tell that if the testing rate is increased by 20% and 30% from its baseline value, the weekly new cases at the peak are decreased by 37.63% and 52.90%; and it also delayed the peak time by four and fourteen weeks, respectively. Similar findings are obtained for the testing efficacy that if it is increased by 12.67% from its baseline value, the weekly new cases at the peak are decreased by 59.05% and delayed the peak by 15 weeks. Therefore, a higher testing rate and efficacy reduce the disease burden by tumbling the new cases, representing a real scenario. It is also obtained that the testing rate and efficacy reduce the epidemic's severity by increasing the final size of the susceptible population. The testing rate is found more significant if testing efficacy is high. Global sensitivity analysis using partial rank correlation coefficients (PRCCs) and Latin hypercube sampling (LHS) determine the key parameters that must be targeted to worsen/contain the epidemic.

Estimating the time-dependent effective reproduction number and vaccination rate for COVID-19 in the USA and India

The effective reproduction number, $ R_t $, is a vital epidemic parameter utilized to judge whether an epidemic is shrinking, growing, or holding steady. The main goal of this paper is to estimate the combined $ R_t $ and time-dependent vaccination rate for COVID-19 in the USA and India after the vaccination campaign started. Accounting for the impact of vaccination into a discrete-time stochastic augmented SVEIR (Susceptible-Vaccinated-Exposed-Infectious-Recovered) model, we estimate the time-dependent effective reproduction number $ (R_t) $ and vaccination rate $ (\xi_t) $ for COVID-19 by using a low pass filter and the Extended Kalman Filter (EKF) approach for the period February 15, 2021 to August 22, 2022 in India and December 13, 2020 to August 16, 2022 in the USA. The estimated $ R_t $ and $ \xi_t $ show spikes and serrations with the data. Our forecasting scenario represents the situation by December 31, 2022 that the new daily cases and deaths are decreasing for the USA and India. We also noticed that for the current vaccination rate, $ R_t $ would remain greater than one by December 31, 2022. Our results are beneficial for the policymakers to track the status of the effective reproduction number, whether it is greater or less than one. As restrictions in these countries ease, it is still important to maintain safety and preventive measures.

From Policy to Prediction: Forecasting COVID-19 Dynamics Under Imperfect Vaccination

Understanding the joint impact of vaccination and non-pharmaceutical interventions on COVID-19 development is important for making public health decisions that control the pandemic. Recently, we created a method in forecasting the daily number of confirmed cases of infectious diseases by combining a mechanistic ordinary differential equation (ODE) model for infectious classes and a generalized boosting machine learning model (GBM) for predicting how public health policies and mobility data affect the transmission rate in the ODE model (Wang et al. in Bull Math Biol 84:57, 2022). In this paper, we extend the method to the post-vaccination period, accordingly obtain a retrospective forecast of COVID-19 daily confirmed cases in the US, and identify the relative influence of the policies used as the predictor variables. In particular, our ODE model contains both partially and fully vaccinated compartments and accounts for the breakthrough cases, that is, vaccinated individuals can still get infected. Our results indicate that the inclusion of data on non-pharmaceutical interventions can significantly improve the accuracy of the predictions. With the use of policy data, the model predicts the number of daily infected cases up to 35 days in the future, with an average mean absolute percentage error of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20.15\%$$\end{document}20.15%, which is further improved to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$14.88\%$$\end{document}14.88% if combined with human mobility data. Moreover, the most influential predictor variables are the policies of restrictions on gatherings, testing and school closing. The modeling approach used in this work can help policymakers design control measures as variant strains threaten public health in the future.