The total impact of the COVID-19 pandemic on mortality should be the least controversial outcome to measure. However, its analysis is complicated by the lack of real time cause specific data; a potential additional issue concerns the quality of coding on the death certificates, particularly in the earlier stages of the pandemic. Furthermore, there are important differences in the recording systems, both across and within countries (e.g. in the UK, England and Wales have consistently reported data on daily deaths based on different time of recording in comparison to Scotland and Northern Ireland). In this context, estimating excess deaths for all causes at national level, with respect to past year trends has been used in several countries as an effective way to evaluate the total burden of the COVID-19 pandemic, including direct COVID-19-related, as well as indirect effects (e.g. people not being able to access health-care). At the same time, all-cause mortality is not affected by mis-coding on the death certificates.

Despite these positive features, this approach can only present global pictures of the total burden of the first wave of the infection. However, to understand the dynamics of the pandemic, we need to analyse data at sub-national level; this would allow to account for geographical differences due to the infectious nature of the disease, as well as those in the population characteristics and health system provision. Additionally, time trends can vary substantially, rendering comparisons even more complex.

The modelling is generally based on the following structure.

$y_{j t s k} =$ number of all-cause deaths in week $j$ of year $t$ for NUTS3 area $s$ and age-sex group $k$

$k = 1, \dots, K = 10 =$ age-sex group (male/female and $<$ 40, 40-59, 60-69, 70-79, 80)
$P_{j t s k} =$ population at risk in week $j$ of year $t$ for NUTS3 area $s$ and age-sex group $k$
$ρ_{j t s k} =$ the risk of death (mortality rate) in week $j$ of year $t$ for NUTS3 area $s$ and age-sex group $k$
$z_{j} =$ dummy variable for public holiday
$x_{j t s} =$ average weekly temperature in each area

The observed counts are modelled as

$y_{j t s k} \sim Poisson (ρ_{j t s k} P_{j t s k}) \log (ρ_{j t s k}) = β_{0 t} + β_{1} z_{j} + f (x_{j t s}) + b_{s} + w_{j}$

and

“Unstructured” effects $β$
- $β_{0 t} = β_{0} + ε_{t}$ : year-specific intercept

$β_{0} \sim Normal (0, 10^{3})$ : global intercept
$ε_{t} \sim Normal (0, τ_{ε}^{- 1})$ : unstructured random effect
$β_{1} \sim Normal (0, 10^{3})$ : effect of public holidays

Non-linear effect of average weekly temperature $f (x_{j t s})$ - RW2 model: $x_{j t s} ∣ x_{(j - 1) t s}, x_{(j - 2) t s}, τ_{x} \sim Normal (2 x_{(j - 1) t s} + x_{(j - 2) t s}, τ_{x}^{- 1})$

Spatial component - $b_{s} = \frac{1}{\sqrt{τ_{b}}} (\sqrt{1 - ϕ} τ_{v}^{0.5} v_{s} + \sqrt{ϕ} τ_{u}^{0.5} u_{s})$ : Besag-York-Mollié (BYM)-type model

$v_{s} \sim Normal (0, τ_{v}^{- 1})$ : unstructured random effect
$u_{s} ∣ u_{- s} \sim Normal (\frac{\sum_{r = 1}^{R} n_{r s} u_{r}}{\sum_{r = 1}^{R} n_{r s}}, \frac{1}{τ_{u} \sum_{r = 1}^{R} n_{r s}})$ : spatially structured random effect
$ϕ \in [0, 1]$ : mixing parameter (measures proportion of variance explained by the structured effect)

Temporal component (non-linear weekly effect) - RW1 model (accounts for seasonality): $w_{j} ∣ w_{j - 1}, τ_{w} \sim Normal (w_{j - 1}, τ_{w}^{- 1})$

Last updated: 20 November 2023

Relevant publications

Filters Active - 0


20202
20221
20233


Journal article4
Preprint/Working paper1
Uncategorised1

Search:


2023
G Konstantinoudis, V Gomez Rubio, M Cameletti, M Pirani, G Baio, M Blangiardo. (2023). "A Workflow for Estimating and Visualising Excess Mortality During the COVID-19 Pandemic". The R Journal Publication GitHub repo DOI
F Hassard, M Vu, S Rahimzadeh, V Castro-Gutierrez, I Stanton, B Burczynska, D Wildeboer, G Baio, M Brown, H Garelick, J Hofman, B Kasprzyk-Hordern, R Kwiatkowska, A Majeed, S Priest, H Denise, M Khalifa, I Bassano, M Wade, J Grimsley, L Lundy, A Singer, M Di Cesare. (2023). "Wastewater monitoring for detection of public health markers during the COVID-19 pandemic: Near-source monitoring of schools in England over an academic year". PLoS One Publication DOI
G Konstantinoudis, V Gomez Rubio, M Cameletti, M Pirani, G Baio, M Blangiardo. (2023). "A framework for estimating and visualising excess mortality during the COVID-19 pandemic". Arxiv
2022
G Konstantinoudis, M Cameletti, V Gomez Rubio, I Leon Gomez, M Pirani, G Baio, A Larrauri, J Riou, M Egger, P Vineis, M Blangiardo. (2022). "Regional excess mortality during the 2020 COVID-19 pandemic: a study of five European countries". Nature Communications Publication GitHub repo DOI
2020
M Blangiardo, M Cameletti, M Pirani, G Corsetti, M Battaglini, G Baio. (2020). "Estimating weekly excess mortality at sub-national level in Italy during the COVID-19 pandemic". PLOS One 15(10): e0240286 DOI
G Baio, M Blangiardo. (2020). "Why counting coronavirus deaths is not an exact science". The Guardian, 19 April 2020 Publication