Skip to content

A stochastic epidemiological model that supplements the conventional reported cases with pooled samples from wastewater for assessing the overall SARS-CoV-2 burden at the community level.

Notifications You must be signed in to change notification settings

mehrdadfazli/COVID-SEIR-1

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

43 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

IEEE medRxiv mendeley contributions-welcome researchgate

Referenced paper : Wastewater-Based Epidemiological Modeling for Continuous Surveillance of COVID-19 Outbreak

COVID-19 Wastewater Surveillance: An Epidemiological Model 😷

A Simple SEIR

Using wastewater surveillance as a continuous pooled sampling technique has been in place in many countries since the early stages of the outbreak of COVID-19. Since the beginning of the outbreak, many research works have emerged, studying different aspects of viral SARS-CoV-2 DNA concentrations (viral load) in wastewater and its potential as an early warning method. However, one of the questions that has remained unanswered is the quantitative relation between viral load and clinical indicators such as daily cases, deaths, and hospitalizations. Few studies have tried to couple viral load data with an epidemiological model to relate the number of infections in the community to the viral burden. We propose a stochastic wastewater-based SEIR model to showcase the importance of viral load in the early detection and prediction of an outbreak in a community. We built three models based on whether or not they use the case count and viral load data and compared their simulations and forecasting quality.

We consider a stochastic wastewater-based epidemiological model with four compartments (hidden states) of susceptible (S), exposed (E), infectious (I), and recovered/removed (R).


Model Specifications

We define the force of infection as below:

ι is the imported infection, α is the mixing parameter, and w is the Gamma white noise.

The rest of the parameters are summarized in table below:

We used iterated filtering in "pomp" package in R to estimate the parameters for the three models.

Models to compare:

  • SEIR-VY: Model that uses both viral load data and case count data
  • SEIR-V: Model that uses only viral load data
  • SEIR-Y: Model that uses only case count data

Maximum Likelihood Estimation (MLE) parameters of the three models:


Results

Estimated Basic Reproduction Number R0

  • SEIR-VY model: R0=1.46
  • SEIR-V model: R0=1.08
  • SEIR-Y model: R0=1.92

Simulation Comparison

Forecasting Comparison

In addition to the three models we also fitted an Autoregressive Integrated Moving Average (ARIMA) model.


Conclusion

In this study, we implemented a SEIR model for three cases of using merely case count data, viral load data, and both. We allowed for stochasticity in the transmission rates and extra-demographic stochasticity accounting for the unforeseen events. We fitted our models using a simulation-based inference called Iterated Filtering. We compared the models from three different aspects, parameter estimation, simulation, and prediction. Our results suggest that the viral load data is an informative data source for monitoring the spread of COVID-19 cases on a community level. The viral load has enough information, which enables it to approximate the number of infected cases when employed with a proper epidemiological model. Additionally, viral load data is more consistent and less uncertain, making it a critical part of COVID-19 surveillance.


For more detailed information of our modeling and references please refer to our paper


Cite us:

@INPROCEEDINGS{9671543,
  author={Fazli, Mehrdad and Sklar, Samuel and Porter, Michael D. and French, Brent A. and Shakeri, Heman},
  booktitle={2021 IEEE International Conference on Big Data (Big Data)}, 
  title={Wastewater-Based Epidemiological Modeling for Continuous Surveillance of COVID-19 Outbreak}, 
  year={2021},
  volume={},
  number={},
  pages={4342-4349},
  doi={10.1109/BigData52589.2021.9671543}}

About

A stochastic epidemiological model that supplements the conventional reported cases with pooled samples from wastewater for assessing the overall SARS-CoV-2 burden at the community level.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • R 99.9%
  • Jupyter Notebook 0.1%