We adapted the BSTI Model for infectious diseases to the current Covid19 pandemic and fit county specific models in germany based on a 3 week window of data. The model is calculated on the HPC Systems at the Juelich Supercomputing Center. The results are can be explored on https://covid19-bayesian.fz-juelich.de/.
The BSTI model is a Bayesian spatio-temporal interaction model, a probabilistic generalized linear model, that predicts aggregated case counts within spatial regions (counties) and time intervals (calendar weeks) using a history of reported cases, temporal features (seasonality and trend) and region-specific as well as demographic information.
Key features of the model:
It is recommended to install a local virtual environment for python3, a requirements.txt is provided.
virtualenv --no-site-packages --python=python3 .venv followed by
pip install requirements.txt to get going.
The experiment is then run by first sampling interactions effects and then the posterior, currently implemented via a SGE scripts
./src/gridjob_sample_posterior.sge . Adapt to different runtimes as needed.
The data is provided by the Robert Koch Institute, and consists of weekly reports of case counts for three diseases, campylobacteriosis, rotavirus infections and Lyme borreliosis. They are aggregated by county and collected over a time period spanning from the 1st of January 2011 (2013 for borreliosis) to the 31st of December 2017 via the SurvNet surveillance system. Aggregated case counts of diseases with mandatory reporting in Germany is available online.
Covid-19 data is provided by the Robert Koch Institute via the publically accessiable this link. The downloadable CSV table can be found in ./data/raw and is preprocessed to fit the BSTI Model implementation via the notebook ./notebooks/preprocess_covid19_data.ipynb.
Information about the shape of counties within Germany is publicly provided by the German federal agency for cartography and geodesy (Bundesamt für Kartographie und Geodäsie) (© GeoBasis-DE / BKG 2018) under the dl-de/by-2-0 license.