Bayesian hierarchical statistical SIRS models

Taking cues from machine learning, compartmental and time series models for SARS-CoV-2 omicron infection in Indian provinces

SARS-CoV-2, the virus responsible for COVID-19, posed a significant threat to the world. We analyzed COVID-19 dissemination data in the top ten Indian provinces by infection incidences using the Susceptible-Infectious-Removed (SIR) model, an Autoregressive Integrated Moving Average (ARIMA) time series model, a machine learning model based on the Random Forest, and distribution fitting. Outbreaks are expected to continue if the Basic Reproduction Number (R0) > 1, and infection waves are anticipated to end if the R0 < 1, as determined by the SIR model. Different parametric probability distributions are also fitted. Data collected from December 12, 2021, to March 31, 2022, encompassing data from both before and during the implementation of strict control measures. Based on the estimates of the model parameters, health agencies and government policymakers can develop strategies to combat the spread of the disease in the future, and the most effective technique can be recommended for real-world application for other outbreaks of COVID-19. The best method out of these could be also implemented further on the epidemiological data of other similar infectious agents.

REAL-TIME MECHANISTIC BAYESIAN FORECASTS OF COVID-19 MORTALITY

The COVID-19 pandemic emerged in late December 2019. In the first six months of the global outbreak, the U.S. reported more cases and deaths than any other country in the world. Effective modeling of the course of the pandemic can help assist with public health resource planning, intervention efforts, and vaccine clinical trials. However, building applied forecasting models presents unique challenges during a pandemic. First, case data available to models in real time represent a nonstationary fraction of the true case incidence due to changes in available diagnostic tests and test-seeking behavior. Second, interventions varied across time and geography leading to large changes in transmissibility over the course of the pandemic. We propose a mechanistic Bayesian model that builds upon the classic compartmental susceptible-exposed-infected-recovered (SEIR) model to operationalize COVID-19 forecasting in real time. This framework includes nonparametric modeling of varying transmission rates, nonparametric modeling of case and death discrepancies due to testing and reporting issues, and a joint observation likelihood on new case counts and new deaths; it is implemented in a probabilistic programming language to automate the use of Bayesian reasoning for quantifying uncertainty in probabilistic forecasts. The model has been used to submit forecasts to the U.S. Centers for Disease Control through the COVID-19 Forecast Hub under the name MechBayes. We examine the performance relative to a baseline model as well as alternate models submitted to the forecast hub. Additionally, we include an ablation test of our extensions to the classic SEIR model. We demonstrate a significant gain in both point and probabilistic forecast scoring measures using MechBayes, when compared to a baseline model, and show that MechBayes ranks as one of the top two models out of nine which regularly submitted to the COVID-19 Forecast Hub for the duration of the pandemic, trailing only the COVID-19 Forecast Hub ensemble model of which which MechBayes is a part.