Background: Vaccination was the single most effective measure in mitigating the impact of the COVID-19 pandemic. Our study aims to quantify the impact of vaccination programmes during the initial year of vaccination (2021) by estimating the number of case fatalities avoided, using Sweden as a case study.

Methods: Using Swedish data on age-specific reported incidence and vaccination uptake, along with vaccine efficacies, age-specific contact patterns and under-reporting from the literature, we fit a Bayesian SEIR epidemic model with time-varying community contact rate for COVID-19 incidence. Age-specific fatality rates from the literature are adjusted proportionally to fit the observed number of case fatalities in the factual analysis, resulting in 5,510 (95% PI: 5,370-5,650) matching the observed number 5,430. The estimated time-varying community contact rate is then used in a counterfactual analysis where the population is unvaccinated, leading to more infections and fatalities. A sensitivity analysis is performed to identify which parameters influence our conclusions.

Findings: The counterfactual analysis result in a severe epidemic outbreak during the early autumn of 2021, resulting in about 37,100 (36,700–37,500) number of case fatalities. Consequently, the number of lives saved by the vaccination programme is estimated to be about 31,600 (31,300–32,000), out of which 5,170 are directly saved and 26,400 are indirectly saved, mainly by drastically reducing the severe outbreak in the early autumn of 2021, which would have occurred without vaccination and unchanged community contact rate.

Interpretation: Our mathematical model is used to analyse the impact of COVID-19 vaccination on lives saved in Sweden during 2021, but the same methodology can be applied to other countries. The counterfactual analysis offers insights into an alternative trajectory of the pandemic without vaccination. The results show the direct impact of vaccination on reducing deaths for infected individuals and shed light on the indirect effects of reduced transmission dynamics.

For evaluating the quality of care provided by hospitals, special interest lies in the identification of performance outliers. The classification of healthcare providers as outliers or non-outliers is a decision under uncertainty, because the true quality is unknown and can only be inferred from an observed result of a quality indicator. We propose to embed the classification of healthcare providers into a Bayesian decision theoretical framework that enables the derivation of optimal decision rules with respect to the expected decision consequences. We propose paradigmatic utility functions for two typical purposes of hospital profiling: the external reporting of healthcare quality and the initiation of change in care delivery. We make use of funnel plots to illustrate and compare the resulting optimal decision rules and argue that sensitivity and specificity of the resulting decision rules should be analyzed. We then apply the proposed methodology to the area of hip replacement surgeries by analyzing data from 1,277 hospitals in Germany which performed over 180,000 such procedures in 2017. Our setting illustrates that the classification of outliers can be highly dependent upon the underlying utilities. We conclude that analyzing the classification of hospitals as a decision theoretic problem helps to derive transparent and justifiable decision rules. The methodology for classifying quality indicator results is implemented in an R package (iqtigbdt) and is available on GitHub. 

Background Globally, since 1 January 2020 and as of 24 January 2023, there have been over 664 million cases of COVID-19 and over 6.7 million deaths reported to WHO. WHO developed an evidence-based alert system, assessing public health risk on a weekly basis in 237 countries, territories and areas from May 2021 to June 2022. This aimed to facilitate the early identification of situations where healthcare capacity may become overstretched. Methods The process involved a three-stage mixed methods approach. In the first stage, future deaths were predicted from the time series of reported cases and deaths to produce an initial alert level. In the second stage, this alert level was adjusted by incorporating a range of contextual indicators and accounting for the quality of information available using a Bayes classifier. In the third stage, countries with an alert level of 'High' or above were added to an operational watchlist and assistance was deployed as needed. Results Since June 2021, the system has supported the release of more than US$27 million from WHO emergency funding, over 450 000 rapid antigen diagnostic testing kits and over 6000 oxygen concentrators. Retrospective evaluation indicated that the first two stages were needed to maximise sensitivity, where 44% (IQR 29%-67%) of weekly watchlist alerts would not have been identified using only reported cases and deaths. The alerts were timely and valid in most cases; however, this could only be assessed on a non-representative sample of countries with hospitalisation data available. Conclusions The system provided a standardised approach to monitor the pandemic at the country level by incorporating all available data on epidemiological analytics and contextual assessments. While this system was developed for COVID-19, a similar system could be used for future outbreaks and emergencies, with necessary adjustments to parameters and indicators.

The real-time analysis of infectious disease surveillance data is essential in obtaining situational awareness about the current dynamics of a major public health event such as the COVID-19 pandemic. This analysis of e.g., time-series of reported cases or fatalities is complicated by reporting delays that lead to under-reporting of the complete number of events for the most recent time points. This can lead to misconceptions by the interpreter, for instance the media or the public, as was the case with the time-series of reported fatalities during the COVID-19 pandemic in Sweden. Nowcasting methods provide real-time estimates of the complete number of events using the incomplete time-series of currently reported events and information about the reporting delays from the past. In this paper we propose a novel Bayesian nowcasting approach applied to COVID-19-related fatalities in Sweden. We incorporate additional information in the form of time-series of number of reported cases and ICU admissions as leading signals. We demonstrate with a retrospective evaluation that the inclusion of ICU admissions as a leading signal improved the nowcasting performance of case fatalities for COVID-19 in Sweden compared to existing methods.

We comment the paper by Jahn et al. (On the role of data, statistics and decisions in a pandemic, 2022).

Background: The European Centre for Disease Prevention and Control (ECDC) systematically collates information from sources to rapidly detect early public health threats. The lack of a freely available, customisable and automated early warning tool using data from Twitter prompted the ECDC to develop epitweetr, which collects, geotocates and aggregates tweets generating signals and email alerts. Aim: This study aims to compare the performance of epitweetr to manually monitoring tweets for the purpose of early detecting public health threats. Methods: We calculated the general and specific positive predictive value (PPV) of signals generated by epitweetr between 19 October and 30 November 2020. Sensitivity, specificity, timeliness and accuracy and performance of tweet geolocation and signal detection algorithms obtained from epitweetr and the manual monitoring of 1,200 tweets were compared. Results: The epitweetr geolocation algorithm had an accuracy of 30.1% at national, and 25.9% at subnational levels. The signal detection algorithm had 3.0% general PPV and 74.6% specific PPV. Compared to manual monitoring, epitweetr had greater sensitivity (47.9% and 78.6%, respectively), and reduced PPV (97.9% and 74.6%, respectively). Median validation time difference between 16 common events detected by epitweetr and manual monitoring was -48.6 hours (IQR: -102.8 to -23.7). Conclusion: Epitweetr has shown sufficient performance as an early warning toot for public health threats using Twitter data. Since epitweetr is a free, open-source tool with configurable settings and a strong automated component, it is expected to increase in usability and usefulness to public health experts.

We determined secondary attack rates (SAR) among close contacts of 59 asymptomatic and symptomatic coronavirus disease case-patients by presymptomatic and symptomatic exposure. We observed no transmission from asymptomatic case-patients and highest SAR through presymptomatic exposure. Rapid quarantine of close contacts with or without symptoms is needed to prevent presymptomatic transmission.

We analysed the coronavirus disease 2019 epidemic curve from March to the end of April 2020 in Germany. We use statistical models to estimate the number of cases with disease onset on a given day and use back-projection techniques to obtain the number of new infections per day. The respective time series are analysed by a trend regression model with change points. The change points are estimated directly from the data. We carry out the analysis for the whole of Germany and the federal state of Bavaria, where we have more detailed data. Both analyses show a major change between 9 and 13 March for the time series of infections: from a strong increase to a decrease. Another change was found between 25 March and 29 March, where the decline intensified. Furthermore, we perform an analysis stratified by age. A main result is a delayed course of the pandemic for the age group 80 + resulting in a turning point at the end of March. Our results differ from those by other authors as we take into account the reporting delay, which turned out to be time dependent and therefore changes the structure of the epidemic curve compared to the curve of newly reported cases.

To assess the current dynamics of an epidemic, it is central to collect information on the daily number of newly diseased cases. This is especially important in real-time surveillance, where the aim is to gain situational awareness, for example, if cases are currently increasing or decreasing. Reporting delays between disease onset and case reporting hamper our ability to understand the dynamics of an epidemic close to now when looking at the number of daily reported cases only. Nowcasting can be used to adjust daily case counts for occurred-but-not-yet-reported events. Here, we present a novel application of nowcasting to data on the current COVID-19 pandemic in Bavaria. It is based on a hierarchical Bayesian model that considers changes in the reporting delay distribution over time and associated with the weekday of reporting. Furthermore, we present a way to estimate the effective time-varying case reproduction number Re(t) based on predictions of the nowcast. The approaches are based on previously published work, that we considerably extended and adapted to the current task of nowcasting COVID-19 cases. We provide methodological details of the developed approach, illustrate results based on data of the current pandemic, and evaluate the model based on synthetic and retrospective data on COVID-19 in Bavaria. Results of our nowcasting are reported to the Bavarian health authority and published on a webpage on a daily basis (https://corona.stat.uni-muenchen.de/). Code and synthetic data for the analysis are available from https://github.com/FelixGuenther/nc_covid19_bavaria and can be used for adaption of our approach to different data.

As several countries gradually release social distancing measures, rapid detection of new localized COVID-19 hotspots and subsequent intervention will be key to avoiding large-scale resurgence of transmission. We introduce ASMODEE (automatic selection of models and outlier detection for epidemics), a new tool for detecting sudden changes in COVID-19 incidence. Our approach relies on automatically selecting the best (fitting or predicting) model from a range of user-defined time series models, excluding the most recent data points, to characterize the main trend in an incidence. We then derive prediction intervals and classify data points outside this interval as outliers, which provides an objective criterion for identifying departures from previous trends. We also provide a method for selecting the optimal breakpoints, used to define how many recent data points are to be excluded from the trend fitting procedure. The analysis of simulated COVID-19 outbreaks suggests ASMODEE compares favourably with a state-of-art outbreak-detection algorithm while being simpler and more flexible. As such, our method could be of wider use for infectious disease surveillance. We illustrate ASMODEE using publicly available data of National Health Service (NHS) Pathways reporting potential COVID-19 cases in England at a fine spatial scale, showing that the method would have enabled the early detection of the flare-ups in Leicester and Blackburn with Darwen, two to three weeks before their respective lockdown. ASMODEE is implemented in the free R package trendbreaker.