K-Track-Covid: interactive web-based dashboard for analyzing geographical and temporal spread of COVID-19 in South Korea.

The COVID-19 pandemic has necessitated the development of robust tools for tracking and modeling the spread of the virus. We present 'K-Track-Covid,' an interactive web-based dashboard developed using the R Shiny framework, to offer users an intuitive dashboard for analyzing the geographical and temporal spread of COVID-19 in South Korea. Our dashboard employs dynamic user interface elements, employs validated epidemiological models, and integrates regional data to offer tailored visual displays. The dashboard allows users to customize their data views by selecting specific time frames, geographic regions, and demographic groups. This customization enables the generation of charts and statistical summaries pertinent to both daily fluctuations and cumulative counts of COVID-19 cases, as well as mortality statistics. Additionally, the dashboard offers a simulation model based on mathematical models, enabling users to make predictions under various parameter settings. The dashboard is designed to assist researchers, policymakers, and the public in understanding the spread and impact of COVID-19, thereby facilitating informed decision-making. All data and resources related to this study are publicly available to ensure transparency and facilitate further research.

Threshold dynamics of a reaction-advection-diffusion schistosomiasis epidemic model with seasonality and spatial heterogeneity.

Most water-borne disease models ignore the advection of water flows in order to simplify the mathematical analysis and numerical computation. However, advection can play an important role in determining the disease transmission dynamics. In this paper, we investigate the long-term dynamics of a periodic reaction-advection-diffusion schistosomiasis model and explore the joint impact of advection, seasonality and spatial heterogeneity on the transmission of the disease. We derive the basic reproduction number R 0 and show that the disease-free periodic solution is globally attractive when R 0 < 1 whereas there is a positive endemic periodic solution and the system is uniformly persistent in a special case when R 0 > 1 . Moreover, we find that R 0 is a decreasing function of the advection coefficients which offers insights into why schistosomiasis is more serious in regions with slow water flows.

Markovian Approach for Exploring Competitive Diseases with Heterogeneity-Evidence from COVID-19 and Influenza in China.

Due to the complex interactions between multiple infectious diseases, the spreading of diseases in human bodies can vary when people are exposed to multiple sources of infection at the same time. Typically, there is heterogeneity in individuals' responses to diseases, and the transmission routes of different diseases also vary. Therefore, this paper proposes an SIS disease spreading model with individual heterogeneity and transmission route heterogeneity under the simultaneous action of two competitive infectious diseases. We derive the theoretical epidemic spreading threshold using quenched mean-field theory and perform numerical analysis under the Markovian method. Numerical results confirm the reliability of the theoretical threshold and show the inhibitory effect of the proportion of fully competitive individuals on epidemic spreading. The results also show that the diversity of disease transmission routes promotes disease spreading, and this effect gradually weakens when the epidemic spreading rate is high enough. Finally, we find a negative correlation between the theoretical spreading threshold and the average degree of the network. We demonstrate the practical application of the model by comparing simulation outputs to temporal trends of two competitive infectious diseases, COVID-19 and seasonal influenza in China.

Real-Time Use of a Dynamic Model To Measure the Impact of Public Health Interventions on Measles Outbreak Size and Duration - Chicago, Illinois, 2024.

Measles is a highly infectious, vaccine-preventable disease that can cause severe illness, hospitalization, and death. A measles outbreak associated with a migrant shelter in Chicago occurred during February-April 2024, in which a total of 57 confirmed cases were identified, including 52 among shelter residents, three among staff members, and two among community members with a known link to the shelter. CDC simulated a measles outbreak among shelter residents using a dynamic disease model, updated in real time as additional cases were identified, to produce outbreak forecasts and assess the impact of public health interventions. As of April 8, the model forecasted a median final outbreak size of 58 cases (IQR = 56-60 cases); model fit and prediction range improved as more case data became available. Counterfactual analysis of different intervention scenarios demonstrated the importance of early deployment of public health interventions in Chicago, with a 69% chance of an outbreak of 100 or more cases had there been no mass vaccination or active case-finding compared with only a 1% chance when those interventions were deployed. This analysis highlights the value of using real-time, dynamic models to aid public health response, set expectations about outbreak size and duration, and quantify the impact of interventions. The model shows that prompt mass vaccination and active case-finding likely substantially reduced the chance of a large (100 or more cases) outbreak in Chicago.

Modelling the impact of precaution on disease dynamics and its evolution.

In this paper, we introduce the notion of practically susceptible population, which is a fraction of the biologically susceptible population. Assuming that the fraction depends on the severity of the epidemic and the public's level of precaution (as a response of the public to the epidemic), we propose a general framework model with the response level evolving with the epidemic. We firstly verify the well-posedness and confirm the disease's eventual vanishing for the framework model under the assumption that the basic reproduction number R 0 < 1 . For R 0 > 1 , we study how the behavioural response evolves with epidemics and how such an evolution impacts the disease dynamics. More specifically, when the precaution level is taken to be the instantaneous best response function in literature, we show that the endemic dynamic is convergence to the endemic equilibrium; while when the precaution level is the delayed best response, the endemic dynamic can be either convergence to the endemic equilibrium, or convergence to a positive periodic solution. Our derivation offers a justification/explanation for the best response used in some literature. By replacing "adopting the best response" with "adapting toward the best response", we also explore the adaptive long-term dynamics.

Conventional and frugal methods of estimating COVID-19-related excess deaths and undercount factors.

Across the world, the officially reported number of COVID-19 deaths is likely an undercount. Establishing true mortality is key to improving data transparency and strengthening public health systems to tackle future disease outbreaks. In this study, we estimated excess deaths during the COVID-19 pandemic in the Pune region of India. Excess deaths are defined as the number of additional deaths relative to those expected from pre-COVID-19-pandemic trends. We integrated data from: (a) epidemiological modeling using pre-pandemic all-cause mortality data, (b) discrepancies between media-reported death compensation claims and official reported mortality, and (c) the "wisdom of crowds" public surveying. Our results point to an estimated 14,770 excess deaths [95% CI 9820-22,790] in Pune from March 2020 to December 2021, of which 9093 were officially counted as COVID-19 deaths. We further calculated the undercount factor-the ratio of excess deaths to officially reported COVID-19 deaths. Our results point to an estimated undercount factor of 1.6 [95% CI 1.1-2.5]. Besides providing similar conclusions about excess deaths estimates across different methods, our study demonstrates the utility of frugal methods such as the analysis of death compensation claims and the wisdom of crowds in estimating excess mortality.

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.

Importance of social inequalities to contact patterns, vaccine uptake, and epidemic dynamics.

Individuals' socio-demographic and economic characteristics crucially shape the spread of an epidemic by largely determining the exposure level to the virus and the severity of the disease for those who got infected. While the complex interplay between individual characteristics and epidemic dynamics is widely recognised, traditional mathematical models often overlook these factors. In this study, we examine two important aspects of human behaviour relevant to epidemics: contact patterns and vaccination uptake. Using data collected during the COVID-19 pandemic in Hungary, we first identify the dimensions along which individuals exhibit the greatest variation in their contact patterns and vaccination uptake. We find that generally higher socio-economic groups of the population have a higher number of contacts and a higher vaccination uptake with respect to disadvantaged groups. Subsequently, we propose a data-driven epidemiological model that incorporates these behavioural differences. Finally, we apply our model to analyse the fourth wave of COVID-19 in Hungary, providing valuable insights into real-world scenarios. By bridging the gap between individual characteristics and epidemic spread, our research contributes to a more comprehensive understanding of disease dynamics and informs effective public health strategies.

Advanced forecasting of COVID-19 epidemic: Leveraging ensemble models, advanced optimization, and decomposition techniques.

In the global effort to address the outbreak of the new coronavirus pneumonia (COVID-19) pandemic, accurate forecasting of epidemic patterns has become crucial for implementing successful interventions aimed at preventing and controlling the spread of the disease. The correct prediction of the course of COVID-19 outbreaks is a complex and challenging task, mainly because of the significant volatility in the data series related to COVID-19. Previous studies have been limited by the exclusive use of individual forecasting techniques in epidemic modeling, disregarding the integration of diverse prediction procedures. The lack of attention to detail in this situation can yield worse-than-ideal results. Consequently, this study introduces a novel ensemble framework that integrates three machine learning methods (kernel ridge regression (KRidge), Deep random vector functional link (dRVFL), and ridge regression) within a linear relationship (L-KRidge-dRVFL-Ridge). The optimization of this framework is accomplished through a distinctive approach, specifically adaptive differential evolution and particle swarm optimization (A-DEPSO). Moreover, an effective decomposition method, known as time-varying filter empirical mode decomposition (TVF-EMD), is employed to decompose the input variables. A feature selection technique, specifically using the light gradient boosting machine (LGBM), is also implemented to extract the most influential input variables. The daily datasets of COVID-19 collected from two countries, namely Italy and Poland, were used as the experimental examples. Additionally, all models are implemented to forecast COVID-19 at two-time horizons: 10- and 14-day ahead (t+10 and t+14). According to the results, the proposed model can yield higher correlation coefficient (R) for both case studies: Italy (t+10 = 0.965, t+14 = 0.961) and Poland (t+10 = 0.952, t+14 = 0.940) than the other models. The experimental results demonstrate that the model suggested in this paper has outstanding results in various kinds of complex epidemic prediction situations. The proposed ensemble model demonstrates exceptional accuracy and resilience, outperforming all similar models in terms of efficacy.

Transmission dynamics of a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods.

In this paper, we propose a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods. Firstly, we establish the well-posedness of the model. Secondly, we define the basic reproduction number ℜ 0 for this model and show that ℜ 0 is a threshold parameter: if ℜ 0 < 1 , then the disease-free periodic solution is globally attractive; if ℜ 0 > 1 , the system is uniformly persistent. Thirdly, we study the global attractivity of the positive steady state when the spatial environment is homogeneous and the advection of mosquitoes is ignored. As an example, we use the model to investigate the dengue fever transmission case in Guangdong Province, China, and explore the impact of model parameters on ℜ 0 . Our findings indicate that ignoring seasonality may underestimate ℜ 0 . Additionally, the spatial heterogeneity of transmission may increase the risk of disease transmission, while the increase of seasonal developmental durations, intrinsic incubation periods and advection rates can all reduce the risk of disease transmission.

Informing policy via dynamic models: Cholera in Haiti.

Public health decisions must be made about when and how to implement interventions to control an infectious disease epidemic. These decisions should be informed by data on the epidemic as well as current understanding about the transmission dynamics. Such decisions can be posed as statistical questions about scientifically motivated dynamic models. Thus, we encounter the methodological task of building credible, data-informed decisions based on stochastic, partially observed, nonlinear dynamic models. This necessitates addressing the tradeoff between biological fidelity and model simplicity, and the reality of misspecification for models at all levels of complexity. We assess current methodological approaches to these issues via a case study of the 2010-2019 cholera epidemic in Haiti. We consider three dynamic models developed by expert teams to advise on vaccination policies. We evaluate previous methods used for fitting these models, and we demonstrate modified data analysis strategies leading to improved statistical fit. Specifically, we present approaches for diagnosing model misspecification and the consequent development of improved models. Additionally, we demonstrate the utility of recent advances in likelihood maximization for high-dimensional nonlinear dynamic models, enabling likelihood-based inference for spatiotemporal incidence data using this class of models. Our workflow is reproducible and extendable, facilitating future investigations of this disease system.

An immuno-epidemiological model with waning immunity after infection or vaccination.

In epidemics, waning immunity is common after infection or vaccination of individuals. Immunity levels are highly heterogeneous and dynamic. This work presents an immuno-epidemiological model that captures the fundamental dynamic features of immunity acquisition and wane after infection or vaccination and analyzes mathematically its dynamical properties. The model consists of a system of first order partial differential equations, involving nonlinear integral terms and different transfer velocities. Structurally, the equation may be interpreted as a Fokker-Planck equation for a piecewise deterministic process. However, unlike the usual models, our equation involves nonlocal effects, representing the infectivity of the whole environment. This, together with the presence of different transfer velocities, makes the proved existence of a solution novel and nontrivial. In addition, the asymptotic behavior of the model is analyzed based on the obtained qualitative properties of the solution. An optimal control problem with objective function including the total number of deaths and costs of vaccination is explored. Numerical results describe the dynamic relationship between contact rates and optimal solutions. The approach can contribute to the understanding of the dynamics of immune responses at population level and may guide public health policies.

Impact of immune evasion, waning and boosting on dynamics of population mixing between a vaccinated majority and unvaccinated minority.

BACKGROUND: We previously demonstrated that when vaccines prevent infection, the dynamics of mixing between vaccinated and unvaccinated sub-populations is such that use of imperfect vaccines markedly decreases risk for vaccinated people, and for the population overall. Risks to vaccinated people accrue disproportionately from contact with unvaccinated people. In the context of the emergence of Omicron SARS-CoV-2 and evolving understanding of SARS-CoV-2 epidemiology, we updated our analysis to evaluate whether our earlier conclusions remained valid. METHODS: We modified a previously published Susceptible-Infectious-Recovered (SIR) compartmental model of SARS-CoV-2 with two connected sub-populations: vaccinated and unvaccinated, with non-random mixing between groups. Our expanded model incorporates diminished vaccine efficacy for preventing infection with the emergence of Omicron SARS-CoV-2 variants, waning immunity, the impact of prior immune experience on infectivity, "hybrid" effects of infection in previously vaccinated individuals, and booster vaccination. We evaluated the dynamics of an epidemic within each subgroup and in the overall population over a 10-year time horizon. RESULTS: Even with vaccine efficacy as low as 20%, and in the presence of waning immunity, the incidence of COVID-19 in the vaccinated subpopulation was lower than that among the unvaccinated population across the full 10-year time horizon. The cumulative risk of infection was 3-4 fold higher among unvaccinated people than among vaccinated people, and unvaccinated people contributed to infection risk among vaccinated individuals at twice the rate that would have been expected based on the frequency of contacts. These findings were robust across a range of assumptions around the rate of waning immunity, the impact of "hybrid immunity", frequency of boosting, and the impact of prior infection on infectivity in unvaccinated people. INTERPRETATION: Although the emergence of the Omicron variants of SARS-CoV-2 has diminished the protective effects of vaccination against infection with SARS-CoV-2, updating our earlier model to incorporate loss of immunity, diminished vaccine efficacy and a longer time horizon, does not qualitatively change our earlier conclusions. Vaccination against SARS-CoV-2 continues to diminish the risk of infection among vaccinated people and in the population as a whole. By contrast, the risk of infection among vaccinated people accrues disproportionately from contact with unvaccinated people.

A New Auto-Regressive Multi-Variable Modified Auto-Encoder for Multivariate Time-Series Prediction: A Case Study with Application to COVID-19 Pandemics.

The SARS-CoV-2 global pandemic prompted governments, institutions, and researchers to investigate its impact, developing strategies based on general indicators to make the most precise predictions possible. Approaches based on epidemiological models were used but the outcomes demonstrated forecasting with uncertainty due to insufficient or missing data. Besides the lack of data, machine-learning models including random forest, support vector regression, LSTM, Auto-encoders, and traditional time-series models such as Prophet and ARIMA were employed in the task, achieving remarkable results with limited effectiveness. Some of these methodologies have precision constraints in dealing with multi-variable inputs, which are important for problems like pandemics that require short and long-term forecasting. Given the under-supply in this scenario, we propose a novel approach for time-series prediction based on stacking auto-encoder structures using three variations of the same model for the training step and weight adjustment to evaluate its forecasting performance. We conducted comparison experiments with previously published data on COVID-19 cases, deaths, temperature, humidity, and air quality index (AQI) in São Paulo City, Brazil. Additionally, we used the percentage of COVID-19 cases from the top ten affected countries worldwide until May 4th, 2020. The results show 80.7% and 10.3% decrease in RMSE to entire and test data over the distribution of 50 trial-trained models, respectively, compared to the first experiment comparison. Also, model type#3 achieved 4th better overall ranking performance, overcoming the NBEATS, Prophet, and Glounts time-series models in the second experiment comparison. This model shows promising forecast capacity and versatility across different input dataset lengths, making it a prominent forecasting model for time-series tasks.

Long-term effects of non-pharmaceutical interventions on total disease burden in parsimonious epidemiological models.

The recent global COVID-19 pandemic resulted in governments enacting non-pharmaceutical interventions (NPIs) targeted at reducing transmission of SARS-CoV-2. But the NPIs also affected the transmission of viruses causing non-target seasonal respiratory diseases, including influenza and respiratory syncytial virus (RSV). In many countries, the NPIs were found to reduce cases of such seasonal respiratory diseases, but there is also evidence that subsequent relaxation of NPIs led to outbreaks of these diseases that were larger than pre-pandemic ones, due to the accumulation of susceptible individuals prior to relaxation. Therefore, the net long-term effects of NPIs on the total disease burden of non-target diseases remain unclear. Knowledge of this is important for infectious disease management and maintenance of public health. In this study, we shed light on this issue for the simplified scenario of a set of NPIs that prevent or reduce transmission of a seasonal respiratory disease for about a year and are then removed, using mathematical analyses and numerical simulations of a suite of four epidemiological models with varying complexity and generality. The model parameters were estimated using empirical data pertaining to seasonal respiratory diseases and covered a wide range. Our results showed that NPIs reduced the total disease burden of a non-target seasonal respiratory disease in the long-term. Expressed as a percentage of population size, the reduction was greater for larger values of the basic reproduction number and the immunity loss rate, reflecting larger outbreaks and hence more infections averted by imposition of NPIs. Our study provides a foundation for exploring the effects of NPIs on total disease burden in more-complex scenarios.

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 .

Modeling county level COVID-19 transmission in the greater St. Louis area: Challenges of uncertainty and identifiability when fitting mechanistic models to time-varying processes.

We use a compartmental model with a time-varying transmission parameter to describe county level COVID-19 transmission in the greater St. Louis area of Missouri and investigate the challenges in fitting such a model to time-varying processes. We fit this model to synthetic and real confirmed case and hospital discharge data from May to December 2020 and calculate uncertainties in the resulting parameter estimates. We also explore non-identifiability within the estimated parameter set. We find that the death rate of infectious non-hospitalized individuals, the testing parameter and the initial number of exposed individuals are not identifiable based on an investigation of correlation coefficients between pairs of parameter estimates. We also explore how this non-identifiability ties back into uncertainties in the estimated parameters and find that it inflates uncertainty in the estimates of our time-varying transmission parameter. However, we do find that R0 is not highly affected by non-identifiability of its constituent components and the uncertainties associated with the quantity are smaller than those of the estimated parameters. Parameter values estimated from data will always be associated with some uncertainty and our work highlights the importance of conducting these analyses when fitting such models to real data. Exploring identifiability and uncertainty is crucial in revealing how much we can trust the parameter estimates.

Estimating household contact matrices structure from easily collectable metadata.

Contact matrices are a commonly adopted data representation, used to develop compartmental models for epidemic spreading, accounting for the contact heterogeneities across age groups. Their estimation, however, is generally time and effort consuming and model-driven strategies to quantify the contacts are often needed. In this article we focus on household contact matrices, describing the contacts among the members of a family and develop a parametric model to describe them. This model combines demographic and easily quantifiable survey-based data and is tested on high resolution proximity data collected in two sites in South Africa. Given its simplicity and interpretability, we expect our method to be easily applied to other contexts as well and we identify relevant questions that need to be addressed during the data collection procedure.

R 0 May Not Tell Us Everything: Transient Disease Dynamics of Some SIR Models Over Patchy Environments.

This paper examines the short-term or transient dynamics of SIR infectious disease models in patch environments. We employ reactivity of an equilibrium and amplification rates, concepts from ecology, to analyze how dispersals/travels between patches, spatial heterogeneity, and other disease-related parameters impact short-term dynamics. Our findings reveal that in certain scenarios, due to the impact of spatial heterogeneity and the dispersals, the short-term disease dynamics over a patch environment may disagree with the long-term disease dynamics that is typically reflected by the basic reproduction number. Such an inconsistence can mislead the public, public healthy agencies and governments when making public health policy and decisions, and hence, these findings are of practical importance.

Predicting Drug-Drug Interactions Involving Rifampicin Using a Semi-mechanistic Hepatic Compartmental Model.

AIMS: To develop a semi-mechanistic hepatic compartmental model to predict the effects of rifampicin, a known inducer of CYP3A4 enzyme, on the metabolism of five drugs, in the hope of informing dose adjustments to avoid potential drug-drug interactions. METHODS: A search was conducted for DDI studies on the interactions between rifampicin and CYP substrates that met specific criteria, including the availability of plasma concentration-time profiles, physical and absorption parameters, pharmacokinetic parameters, and the use of healthy subjects at therapeutic doses. The semi-mechanistic model utilized in this study was improved from its predecessors, incorporating additional parameters such as population data (specifically for Chinese and Caucasians), virtual individuals, gender distribution, age range, dosing time points, and coefficients of variation. RESULTS: Optimal parameters were identified for our semi-mechanistic model by validating it with clinical data, resulting in a maximum difference of approximately 2-fold between simulated and observed values. PK data of healthy subjects were used for most CYP3A4 substrates, except for gilteritinib, which showed no significant difference between patients and healthy subjects. Dose adjustment of gilteritinib co-administered with rifampicin required a 3-fold increase of the initial dose, while other substrates were further tuned to achieve the desired drug exposure. CONCLUSIONS: The pharmacokinetic parameters AUCR and CmaxR of drugs metabolized by CYP3A4, when influenced by Rifampicin, were predicted by the semi-mechanistic model to be approximately twice the empirically observed values, which suggests that the semi-mechanistic model was able to reasonably simulate the effect. The doses of four drugs adjusted via simulation to reduce rifampicin interaction.