Introduction: Excess mortality (EM) is a valid indicator of COVID-19's impact on public health. Several studies regarding the estimation of EM have been conducted in Italy, and some of them have shown conflicting values. We focused on three estimation models and compared their results with respect to the same target population, which allowed us to highlight their strengths and limitations.
Methods: We selected three estimation models: model 1 (Maruotti et al.) is a Negative-Binomial GLMM with seasonal patterns; model 2 (Dorrucci et al.) is a Negative Binomial GLM epidemiological approach; and model 3 (Scortichini et al.) is a quasi-Poisson GLM time-series approach with temperature distributions. We extended the time windows of the original models until December 2021, computing various EM estimates to allow for comparisons.
Results: We compared the results with our benchmark, the ISS-ISTAT official estimates. Model 1 was the most consistent, model 2 was almost identical, and model 3 differed from the two. Model 1 was the most stable towards changes in the baseline years, while model 2 had a lower cross-validation RMSE.
Discussion: Presently, an unambiguous explanation of EM in Italy is not possible. We provide a range that we consider sound, given the high variability associated with the use of different models. However, all three models accurately represented the spatiotemporal trends of the pandemic waves in Italy.
Keywords: COVID-19; all-cause mortality; coronavirus; excess deaths; statistical models.