Combining and comparing regional epidemic dynamics in Italy: Bayesian meta-analysis of compartmental models and model assessment via Global Sensitivity Analysis

During autumn 2020, Italy faced a second important SARS-CoV-2 epidemic wave. We explored the time pattern of the instantaneous reproductive number, R0(t), and estimated the prevalence of infections by region from August to December calibrating SIRD models on COVID19-related deaths, fixing at values from literature Infection Fatality Rate (IFR) and infection duration. A Global Sensitivity Analysis (GSA) was performed on the regional SIRD models. Then, we used Bayesian meta-analysis and meta-regression to combine and compare the regional results and investigate their heterogeneity. The meta-analytic R0(t) curves were similar in the Northern and Central regions, while a less peaked curve was estimated for the South. The maximum R0(t) ranged from 2.61 (North) to 2.15 (South) with an increase following school reopening and a decline at the end of October. Average temperature, urbanization, characteristics of family medicine and health care system, economic dynamism, and use of public transport could partly explain the regional heterogeneity. The GSA indicated the robustness of the regional R0(t) curves to different assumptions on IFR. The infectious period turned out to have a key role in determining the model results, but without compromising between-region comparisons.