Propensity score analysis with missing data using a multi-task neural network

BACKGROUND

Propensity score analysis is increasingly used to control for confounding factors in observational studies. Unfortunately, unavoidable missing values make estimating propensity scores extremely challenging. We propose a new method for estimating propensity scores in data with missing values.

MATERIALS AND METHODS

Both simulated and real-world datasets are used in our experiments. The simulated datasets were constructed under 2 scenarios, the presence (T = 1) and the absence (T = 0) of the true effect. The real-world dataset comes from LaLonde's employment training program. We construct missing data with varying degrees of missing rates under three missing mechanisms: MAR, MCAR, and MNAR. Then we compare MTNN with 2 other traditional methods in different scenarios. The experiments in each scenario were repeated 20,000 times. Our code is publicly available at https://github.com/ljwa2323/MTNN .

RESULTS

Under the three missing mechanisms of MAR, MCAR and MNAR, the RMSE between the effect and the true effect estimated by our proposed method is the smallest in simulations and in real-world data. Furthermore, the standard deviation of the effect estimated by our method is the smallest. In situations where the missing rate is low, the estimation of our method is more accurate.

CONCLUSIONS

MTNN can perform propensity score estimation and missing value filling at the same time through shared hidden layers and joint learning, which solves the dilemma of traditional methods and is very suitable for estimating true effects in samples with missing values. The method is expected to be broadly generalized and applied to real-world observational studies.

Resurgence of different influenza types in China and the US in 2021

Various nonpharmaceutical interventions (NPIs) were implemented to alleviate the COVID-19 pandemic since its outbreak. The transmission dynamics of other respiratory infectious diseases, such as seasonal influenza, were also affected by these interventions. The drastic decline of seasonal influenza caused by such interventions would result in waning of population immunity and may trigger the seasonal influenza epidemic with the lift of restrictions during the post-pandemic era. We obtained weekly influenza laboratory confirmations from FluNet to analyse the resurgence patterns of seasonal influenza in China and the US. Our analysis showed that due to the impact of NPIs including travel restrictions between countries, the influenza resurgence was caused by influenza virus A in the US while by influenza virus B in China.

Twice evasions of Omicron variants explain the temporal patterns in six Asian and Oceanic countries

BACKGROUND

The ongoing coronavirus 2019 (COVID-19) pandemic has emerged and caused multiple pandemic waves in the following six countries: India, Indonesia, Nepal, Malaysia, Bangladesh and Myanmar. Some of the countries have been much less studied in this devastating pandemic. This study aims to assess the impact of the Omicron variant in these six countries and estimate the infection fatality rate (IFR) and the reproduction number [Formula: see text] in these six South Asia, Southeast Asia and Oceania countries.

METHODS

We propose a Susceptible-Vaccinated-Exposed-Infectious-Hospitalized-Death-Recovered model with a time-varying transmission rate [Formula: see text] to fit the multiple waves of the COVID-19 pandemic and to estimate the IFR and [Formula: see text] in the aforementioned six countries. The level of immune evasion and the intrinsic transmissibility advantage of the Omicron variant are also considered in this model.

RESULTS

We fit our model to the reported deaths well. We estimate the IFR (in the range of 0.016 to 0.136%) and the reproduction number [Formula: see text] (in the range of 0 to 9) in the six countries. Multiple pandemic waves in each country were observed in our simulation results.

CONCLUSIONS

The invasion of the Omicron variant caused the new pandemic waves in the six countries. The higher [Formula: see text] suggests the intrinsic transmissibility advantage of the Omicron variant. Our model simulation forecast implies that the Omicron pandemic wave may be mitigated due to the increasing immunized population and vaccine coverage.

Effects of vaccination on mitigating COVID-19 outbreaks: a conceptual modeling approach

This paper is devoted to investigating the impact of vaccination on mitigating COVID-19 outbreaks. In this work, we propose a compartmental epidemic ordinary differential equation model, which extends the previous so-called SEIRD model ^[1,2,3,4] by incorporating the birth and death of the population, disease-induced mortality and waning immunity, and adding a vaccinated compartment to account for vaccination. Firstly, we perform a mathematical analysis for this model in a special case where the disease transmission is homogeneous and vaccination program is periodic in time. In particular, we define the basic reproduction number $ \mathcal{R}_0 $ for this system and establish a threshold type of result on the global dynamics in terms of $ \mathcal{R}_0 $. Secondly, we fit our model into multiple COVID-19 waves in four locations including Hong Kong, Singapore, Japan, and South Korea and then forecast the trend of COVID-19 by the end of 2022. Finally, we study the effects of vaccination again the ongoing pandemic by numerically computing the basic reproduction number $ \mathcal{R}_0 $ under different vaccination programs. Our findings indicate that the fourth dose among the high-risk group is likely needed by the end of the year.

Estimation of the serial interval of monkeypox during the early outbreak in 2022

With increased transmissibility and novel transmission mode, monkeypox poses new threats to public health globally in the background of the ongoing COVID-19 pandemic. Estimates of the serial interval, a key epidemiological parameter of infectious disease transmission, could provide insights into the virus transmission risks. As of October 2022, little was known about the serial interval of monkeypox due to the lack of contact tracing data. In this study, public-available contact tracing data of global monkeypox cases were collected and 21 infector-infectee transmission pairs were identified. We proposed a statistical method applied to real-world observations to estimate the serial interval of the monkeypox. We estimated a mean serial interval of 5.6 days with the right truncation and sampling bias adjusted and calculated the reproduction number of 1.33 for the early monkeypox outbreaks at a global scale. Our findings provided a preliminary understanding of the transmission potentials of the current situation of monkeypox outbreaks. We highlighted the need for continuous surveillance of monkeypox for transmission risk assessment.

The mechanism shaping the logistic growth of mutation proportion in epidemics at population scale

Virus evolution is a common process of pathogen adaption to host population and environment. Frequently, a small but important fraction of virus mutations are reported to contribute to higher risks of host infection, which is one of the major determinants of infectious diseases outbreaks at population scale. The key mutations contributing to transmission advantage of a genetic variant often grow and reach fixation rapidly. Based on classic epidemiology theories of disease transmission, we proposed a mechanistic explanation of the process that between-host transmission advantage may shape the observed logistic curve of the mutation proportion in population. The logistic growth of mutation is further generalized by incorporating time-varying selective pressure to account for impacts of external factors on pathogen adaptiveness. The proposed model is implemented in real-world data of COVID-19 to capture the emerging trends and changing dynamics of the B.1.1.7 strains of SARS-CoV-2 in England. The model characterizes and establishes the underlying theoretical mechanism that shapes the logistic growth of mutation in population.

Differences in the superspreading potentials of COVID-19 across contact settings

BACKGROUND

Superspreading events (SSEs) played a critical role in fueling the COVID-19 outbreaks. Although it is well-known that COVID-19 epidemics exhibited substantial superspreading potential, little is known about the risk of observing SSEs in different contact settings. In this study, we aimed to assess the potential of superspreading in different contact settings in Japan.

METHOD

Transmission cluster data from Japan was collected between January and July 2020. Infector-infectee transmission pairs were constructed based on the contact tracing history. We fitted the data to negative binomial models to estimate the effective reproduction number (R) and dispersion parameter (k). Other epidemiological issues relating to the superspreading potential were also calculated.

RESULTS

The overall estimated R and k are 0.561 (95% CrI: 0.496, 0.640) and 0.221 (95% CrI: 0.186, 0.262), respectively. The transmission in community, healthcare facilities and school manifest relatively higher superspreading potentials, compared to other contact settings. We inferred that 13.14% (95% CrI: 11.55%, 14.87%) of the most infectious cases generated 80% of the total transmission events. The probabilities of observing superspreading events for entire population and community, household, health care facilities, school, workplace contact settings are 1.75% (95% CrI: 1.57%, 1.99%), 0.49% (95% CrI: 0.22%, 1.18%), 0.07% (95% CrI: 0.06%, 0.08%), 0.67% (95% CrI: 0.31%, 1.21%), 0.33% (95% CrI: 0.13%, 0.94%), 0.32% (95% CrI: 0.21%, 0.60%), respectively.

CONCLUSION

The different potentials of superspreading in contact settings highlighted the need to continuously monitoring the transmissibility accompanied with the dispersion parameter, to timely identify high risk settings favoring the occurrence of SSEs.

Global stability and optimal control for a COVID-19 model with vaccination and isolation delays

COVID-19 pandemic remains serious around the world and causes huge deaths and economic losses. To investigate the effect of vaccination and isolation delays on the transmission of COVID-19, we propose a mathematical model of COVID-19 transmission with vaccination and isolation delays. The basic reproduction number is computed, and the global dynamics of the model are proved. When R 0 < 1 , the disease-free equilibrium is globally asymptotically stable. The unique endemic equilibrium is globally asymptotically stable if R 0 > 1 . Based on the public information, parameter values are estimated, and sensitivity analysis is carried out by the partial rank correlation coefficients (PRCCs) and the extended version of the Fourier amplitude sensitivity test (eFAST). Our results suggest that the isolation rates of asymptomatic and symptomatic infectious individuals have a significant impact on the transmission of COVID-19. When the COVID-19 is epidemic, the optimal control strategies of our model with vaccination and isolation delays are analyzed. Under the limited resource with constant and time-varying isolation rates, we find that the optimal isolation rates may minimize the cumulative number of infected individuals and the cost of disease control, and effectively contain the transmission of COVID-19. Our study may help public health to prevent and control the COVID-19 spread.

Unravelling the dynamics of Lassa fever transmission with differential infectivity: Modeling analysis and control strategies

Epidemic models have been broadly used to comprehend the dynamic behaviour of emerging and re-emerging infectious diseases, predict future trends, and assess intervention strategies. The symptomatic and asymptomatic features and environmental factors for Lassa fever (LF) transmission illustrate the need for sophisticated epidemic models to capture more vital dynamics and forecast trends of LF outbreaks within countries or sub-regions on various geographic scales. This study proposes a dynamic model to examine the transmission of LF infection, a deadly disease transmitted mainly by rodents through environment. We extend prior LF models by including an infectious stage to mild and severe as well as incorporating environmental contributions from infected humans and rodents. For model calibration and prediction, we show that the model fits well with the LF scenario in Nigeria and yields remarkable prediction results. Rigorous mathematical computation divulges that the model comprises two equilibria. That is disease-free equilibrium, which is locally-asymptotically stable (LAS) when the basic reproduction number, $ {\mathcal{R}}_{0} $, is $ < 1 $; and endemic equilibrium, which is globally-asymptotically stable (GAS) when $ {\mathcal{R}}_{0} $ is $ > 1 $. We use time-dependent control strategy by employing Pontryagin's Maximum Principle to derive conditions for optimal LF control. Furthermore, a partial rank correlation coefficient is adopted for the sensitivity analysis to obtain the model's top rank parameters requiring precise attention for efficacious LF prevention and control.

Evaluation of Effectiveness of Global COVID-19 Vaccination Campaign

To model estimated deaths averted by COVID-19 vaccines, we used state-of-the-art mathematical modeling, likelihood-based inference, and reported COVID-19 death and vaccination data. We estimated that >1.5 million deaths were averted in 12 countries. Our model can help assess effectiveness of the vaccination program, which is crucial for curbing the COVID-19 pandemic.

Estimating the incubation period of SARS-CoV-2 Omicron BA.1 variant in comparison with that during the Delta variant dominance in South Korea

Based on exposure history and symptom onset of 22 Omicron BA.1 cases in South Korea from November to December 2021, we estimated mean incubation period of 3.5 days (95% CI: 2.5, 3.8), and then compared to that of 6.5 days (95% CI: 5.3, 7.7) for 64 cases during Delta variants' dominance in June 2021. For Omicron BA.1 variants, we found that 95% of symptomatic cases developed clinical conditions within 6.0 days (95% CI: 4.3, 6.6) after exposure. Thus, a shorter quarantine period may be considered based on symptoms, or similarly laboratory testing, when Omicron BA.1 variants are circulating.

A simple model to estimate the transmissibility of the Beta, Delta, and Omicron variants of SARS-COV-2 in South Africa

The COVID-19 pandemic caused multiple waves of mortality in South Africa, where three genetic variants of SARS-COV-2 and their ancestral strain dominated consecutively. State-of-the-art mathematical modeling approach was used to estimate the time-varying transmissibility of SARS-COV-2 and the relative transmissibility of Beta, Delta, and Omicron variants. The transmissibility of the three variants were about 73%, 87%, and 276% higher than their preceding variants. To the best of our knowledge, our model is the first simple model that can simulate multiple mortality waves and three variants' replacements in South Africa. The transmissibility of the Omicron variant is substantially higher than that of previous variants.

Reduction in the infection fatality rate of Omicron variant compared with previous variants in South Africa

OBJECTIVE

The SARS-CoV-2 Omicron (B.1.1.529) variant has caused global concern. Previous studies have shown that the variant has enhanced immune evasion ability and transmissibility and reduced severity.

METHODS

In this study, we developed a mathematical model with time-varying transmission rate, vaccination, and immune evasion. We fit the model to reported case and death data up to February 6, 2022 to estimate the transmissibility and infection fatality ratio of the Omicron variant in South Africa.

RESULTS

We found that the high relative transmissibility of the Omicron variant was mainly due to its immune evasion ability, whereas its infection fatality rate substantially decreased by approximately 78.7% (95% confidence interval: 66.9%, 85.0%) with respect to previous variants.

CONCLUSION

On the basis of data from South Africa and mathematical modeling, we found that the Omicron variant is highly transmissible but with significantly lower infection fatality rates than those of previous variants of SARS-CoV-2.

Transmission dynamics of COVID-19 pandemic with combined effects of relapse, reinfection and environmental contribution: A modeling analysis

Reinfection and reactivation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have recently raised public health pressing concerns in the fight against the current pandemic globally. In this study, we propose a new dynamic model to study the transmission of the coronavirus disease 2019 (COVID-19) pandemic. The model incorporates possible relapse, reinfection and environmental contribution to assess the combined effects on the overall transmission dynamics of SARS-CoV-2. The model's local asymptotic stability is analyzed qualitatively. We derive the formula for the basic reproduction number ( R 0 ) and final size epidemic relation, which are vital epidemiological quantities that are used to reveal disease transmission status and guide control strategies. Furthermore, the model is validated using the COVID-19 reported situations in Saudi Arabia. Moreover, sensitivity analysis is examined by implementing a partial rank correlation coefficient technique to obtain the ultimate rank model parameters to control or mitigate the pandemic effectively. Finally, we employ a standard Euler technique for numerical simulations of the model to elucidate the influence of some crucial parameters on the overall transmission dynamics. Our results highlight that contact rate, hospitalization rate, and reactivation rate are the fundamental parameters that need particular emphasis for the prevention, mitigation and control.

Regional heterogeneity of in-hospital mortality of COVID-19 in Brazil

Background

The ongoing Coronavirus disease of 2019 (COVID-19) pandemic has hit Brazil hard in period of different dominant variants. Different COIVD-19 variants have swept through the region, resulting that the total number of cases in Brazil is the third highest in the world. This study is aimed at investigating the regional heterogeneity of in-hospital mortality of COVID-19 in Brazil and the effects of vaccination and social inequality.

Methods

We fitted a multivariate mixed-effects Cox model to a national database of inpatient data in Brazil who were admitted for COVID-19 from February 27, 2020 to March 15, 2022. The in-hospital mortality risks of vaccinated and unvaccinated patients were compared, with adjustment for age, state, ethnicity, education and comorbidities. And the effects of variables to in-hospital mortality were also compared. Stratified analysis was conducted across different age groups and vaccine types.

Results

By fitting the multivariate mixed-effects Cox model, we concluded that age was the most important risk factor for death. With regards to educational level, illiterate patients (hazard ratio: 1.63, 95% CI: 1.56–1.70) had a higher risk than those with a university or college degree. Some common comorbidities were more dangerous for hospitalized patients, such as liver disease (HR: 1.46, 95% CI: 1.34–1.59) and immunosuppression (HR:1.32, 95% CI: 1.26–1.40). In addition, the states involving Sergipe (HR: 1.75, 95% CI: 1.46–2.11), Roraima (HR: 1.65, 95% CI: 1.43–1.92), Maranhão (HR: 1.57, 95% CI: 1.38–1.79), Acre (HR: 1.44, 95% CI: 1.12–1.86), and Rondônia (HR: 1.26, 95% CI: 1.10–1.44) in the north and the northeast region tended to have higher hazard ratios than other area. In terms of vaccine protection, vaccination did not significantly reduce mortality among hospitalized patients. Sinovac and AstraZeneca offered different protection in different regions, and no vaccine provided high protection in all regions.

Conclusion

The study revealed the regional heterogeneity of in-hospital mortality of Covid-19 in Brazil and the effects of vaccination and social inequality. We found that ethnic concentrations were consistent with higher proportion of death cases relative to population size. White Brazilians had more frequent international travel opportunities. As race revealed the intersection of social connections, we speculated that uneven interactions with residential communities partially contribute to the spread of the epidemic. Additionally, the vaccine showed different protection in different regions. In the northern and northeastern regions, AstraZeneca was much more protective than Sinovac, while Sinovac was more protective for hospitalized patients with varying numbers of comorbidities in the Central-west, Southeast and South regions.

Characterizing superspreading potential of infectious disease: Decomposition of individual transmissibility

In the context of infectious disease transmission, high heterogeneity in individual infectiousness indicates that a few index cases can generate large numbers of secondary cases, a phenomenon commonly known as superspreading. The potential of disease superspreading can be characterized by describing the distribution of secondary cases (of each seed case) as a negative binomial (NB) distribution with the dispersion parameter, k. Based on the feature of NB distribution, there must be a proportion of individuals with individual reproduction number of almost 0, which appears restricted and unrealistic. To overcome this limitation, we generalized the compound structure of a Poisson rate and included an additional parameter, and divided the reproduction number into independent and additive fixed and variable components. Then, the secondary cases followed a Delaporte distribution. We demonstrated that the Delaporte distribution was important for understanding the characteristics of disease transmission, which generated new insights distinct from the NB model. By using real-world dataset, the Delaporte distribution provides improvements in describing the distributions of COVID-19 and SARS cases compared to the NB distribution. The model selection yielded increasing statistical power with larger sample sizes as well as conservative type I error in detecting the improvement in fitting with the likelihood ratio (LR) test. Numerical simulation revealed that the control strategy-making process may benefit from monitoring the transmission characteristics under the Delaporte framework. Our findings highlighted that for the COVID-19 pandemic, population-wide interventions may control disease transmission on a general scale before recommending the high-risk-specific control strategies.

Superspreading is one of the key transmission features of many infectious diseases and is considered a consequence of the heterogeneity in infectiousness of individual cases. To characterize the superspreading potential, we divided individual infectiousness into two independent and additive components, including a fixed baseline and a variable part. Such decomposition produced an improvement in the fit of the model explaining the distribution of real-world datasets of COVID-19 and SARS that can be captured by the classic statistical tests. Disease control strategies may be developed by monitoring the characteristics of superspreading. For the COVID-19 pandemic, population-wide interventions are suggested first to limit the transmission at a scale of general population, and then high-risk-specific control strategies are recommended subsequently to lower the risk of superspreading.

Modelling COVID-19 outbreak on the Diamond Princess ship using the public surveillance data

The novel coronavirus disease 2019 (COVID-19) outbreak on the Diamond Princess (DP) ship has caused over 634 cases as of February 20, 2020. We model the transmission process on DP ship as a stochastic branching process, and estimate the reproduction number at the innitial phase of 2.9 (95%CrI: 1.7–7.7). The epidemic doubling time is 3.4 days, and thus timely actions on COVID-19 control were crucial. We estimate the COVID-19 transmissibility reduced 34% after the quarantine program on the DP ship which was implemented on February 5. According to the model simulation, relocating the population at risk may sustainably decrease the epidemic size, postpone the timing of epidemic peak, and thus relieve the tensive demands in the healthcare. The lesson learnt on the ship should be considered in other similar settings.

Heterogeneous epidemic modelling within an enclosed space and corresponding Bayesian estimation

Since March 11th, 2020, COVID-19 has been a global pandemic for more than one years due to a long and infectious incubation period. This paper establishes a heterogeneous epidemic model that divides the incubation period into infectious and non-infectious and employs the Bayesian framework to model the ‘Diamond Princess’ enclosed space incident. The heterogeneity includes two different identities, two transmission methods, two different-size rooms, and six transmission stages. This model is also applicable to similar mixed structures, including closed schools, hospitals, and communities. As the COVID-19 pandemic continues, our mathematical modeling can provide management insights to the governments and policymakers on how the COVID-19 disease has spread and what prevention strategies still need to be taken.

•

A completely data-driven linear heterogeneous epidemic model for a relatively enclosed space composed of subspaces up to three people is developed.

•

The heterogeneity of the model includes the following five aspects: two different identities, two infection sources, two transmission methods, two different-size rooms, and six transmission stages.

•

A probabilistic algorithm to calculate the change in the subspace distribution of people is constructed.

•

A simple method is introduced to calculate the basic reproduction number R0 for different infection sources.

•

An improved approximate Bayesian updating computation method based on entirely non-informative priors that simultaneously infers any number of unknown parameters is proposed.

Superspreading potential of SARS-CoV-2 Delta variants under intensive disease control measures in China

Given the heterogeneity in individual transmissibility, we estimated the superspreading potential of SARS-CoV-2 Delta variants. Using case series of Delta variants in Guangdong, China, we found 15% (95%CrI: 12, 19) of cases seeded 80% of offspring cases.

Impact of low vaccine coverage on the resurgence of COVID-19 in Central and Eastern Europe

The coronavirus disease 2019 (COVID-19) pandemic has caused a tremendous global impact both socially and economically. The mechanisms behind the disparity in the severity, vaccine coverage, and variant replacement patterns across European countries are unclear. In this work, we aim to reveal the possible reasons via data visualization and model fitting. We developed a model with a vaccination component to simulate the mortality waves in these countries. Deaths averted by the vaccination campaign were estimated. Finally, we discuss the potential reasons behind the differences in vaccine coverage across European countries. Contemporary transportation and global trade bring significant convenience to our daily life but also facilitate the spread of the novel virus COVID-19 to anywhere globally within a short time. The observations and results in this work highlight the importance of the global campaign to mitigate the COVID-19 pandemic and future pandemics under the One Health approach.

•

We reveal disparity in COVID-19 vaccine coverage across European counties.

•

We reveal different patterns of COVID-19 variants across European countries.

•

Using a mathematical model, we calculate deaths averted by the vaccine in Europe.

•

We discuss the reasons behind the disparity in vaccine coverage in Europe.