# Explaining among-country variation in COVID-19 case fatality rate.

Gabriele Sorci, Bruno Faivre, Serge Morand

<figure><img src="image_1.png" /></figure>

## To cite this version:

Gabriele Sorci, Bruno Faivre, Serge Morand. Explaining among-country variation in COVID-19 case fatality rate.. Scientific Reports, 2020, 10 (1), pp.18909. 10.1038/s41598-020-75848-2. hal-02992811

HAL Id: hal-02992811

https://hal.science/hal-02992811

Submitted on 19 Nov 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# Explaining among-country variation in COVID-19 case fatality rate

Gabriele Sorci<fnref n="1" />✉, Bruno Faivre<fnref n="1" /> & Serge Morand²,³

While the epidemic of SARS-CoV-2 has spread worldwide, there is much concern over the mortality rate that the infection induces. Available data suggest that COVID-19 case fatality rate had varied temporally (as the epidemic has progressed) and spatially (among countries). Here, we attempted to identify key factors possibly explaining the variability in case fatality rate across countries. We used data on the temporal trajectory of case fatality rate provided by the European Center for Disease Prevention and Control, and country-specific data on different metrics describing the incidence of known comorbidity factors associated with an increased risk of COVID-19 mortality at the individual level. We also compiled data on demography, economy and political regimes for each country. We found that temporal trajectories of case fatality rate greatly vary among countries. We found several factors associated with temporal changes in case fatality rate both among variables describing comorbidity risk and demographic, economic and political variables. In particular, countries with the highest values of DALYs lost to cardiovascular, cancer and chronic respiratory diseases had the highest values of COVID-19 CFR. CFR was also positively associated with the death rate due to smoking in people over 70 years. Interestingly, CFR was negatively associated with share of death due to lower respiratory infections. Among the demographic, economic and political variables, CFR was positively associated with share of the population over 70, GDP per capita, and level of democracy, while it was negatively associated with number of hospital beds ×1000. Overall, these results emphasize the role of comorbidity and socio-economic factors as possible drivers of COVID-19 case fatality rate at the population level.

At the end of 2019, a novel coronavirus (SARS-CoV-2) emerged from an animal reservoir in the city of Wuhan, China<fnref n="1" />,². Having established a human-to-human transmission, the virus rapidly spread, first within China, and subsequently outside China, worldwide. On 11th March 2020 the WHO declared SARS-CoV-2 as a pandemic³. The coronavirus disease 2019 (COVID-19) produces a series of respiratory symptoms that vary in severity⁴⁻¹⁰. Although in most people the infection produces an asymptomatic disease or mild symptoms that do not require particular medical care, in a fraction of patients the disease develops into severe respiratory distress, due to an overreacting inflammatory response, requiring hospitalization in intensive care units. In the absence of specific treatments that prevent or block viral replication, a serious matter of concern is the proportion of infected patients that will eventually die. Current data indicate that, worldwide, case fatality rate (CFR, the ratio between number of deaths and number of confirmed cases) might be around 4%. However, at the country level, CFR ranges from 0 to more than 20%. There are many possible reasons for such a variation¹¹. First, the epidemic has spread in some countries earlier than in others and therefore the difference in CFR might reflect different stages during the spread of the disease. According to this hypothesis, one might expect an increase in CFR with time in countries where CFR is currently low. Second, at the individual level, clinical data have reported several risk factors associated with a poor prognosis. Age and comorbidities (cardiovascular diseases, cancers, diabetes mellitus, chronic lung diseases) seem to greatly increase the mortality risk¹²⁻¹⁵. Therefore, countries with a greater share of elderly in the population, or with higher incidence of recognized comorbidity factors might pay the highest toll to the infection. Third, as said before, in the absence of effective treatments, patients that develop severe symptoms require hospitalization in intensive care units for respiratory assistance. A major concern is that, as the number of infected people increases, the healthcare system will be overwhelmed¹⁶. Under this scenario, mortality rate might reflect the country-specific capacity to tackle a large number of patients requiring

<footnote n="1">¹Biogéosciences, CNRS UMR 6282, Université de Bourgogne Franche-Comté, 6 Boulevard Gabriel, 21000 Dijon, France. ²CNRS ISEM - CIRAD ASTRE - Montpellier Université, Montpellier, France. ³Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand. ✉email: gabriele.sorci@u-bourgogne.fr</footnote>
respiratory assistance and intensive care. Finally, given that CFR is defined as the ratio between number of deaths and number of confirmed cases, countries might simply differ in the accuracy with which they detect the infection. Indeed, CFR estimates are prone to error if the actual number of infected people is much higher than the number of PCR confirmed cases (producing an overestimated CFR), or if mortality occurs with a delay (producing an underestimated CFR). Therefore, variation in CRF might reflect among country variation (i) in population screening (i.e., the number of tests performed), which affects the denominator of the ratio between number of deaths and number of confirmed cases, (ii) in counting and communicating the actual number of patients that have succumbed from SARS-CoV-2 infection.

Here, we conducted an analysis of the factors that might account for the variability of COVID-19 CFR among countries. Using data updated to June 11th 2020, we first investigated if CFR significantly varies between countries, independently from (i) the stage of the epidemic wave, (ii) the testing strategy, and (iii) the social distancing policies adopted by each country. Second, we used country-specific variables assessing the occurrence of comorbidities, as well as demographic, economic and political variables to uncover any association pattern with COVID-19 CFR.

## Methods

We used data on daily number of confirmed cases and deaths for each country reported by the European Center for Disease Prevention and Control (ECDC). We computed the case fatality rate as the ratio between deaths and confirmed cases. We restricted the dataset to countries with at least 100 confirmed cases to avoid spurious results due to small numbers. For each country, we also counted the number of days between the 100th case and 11th June 2020, and the number of days from the occurrence of at least one death × 1,000,000, indicating the progression of the epidemic.

We used the online resource https://ourworldindata.org/ to retrieve data on the incidence of known comorbidity factors for each country, as well as information on demographics, economics and political regimes. In particular, we used different metrics to describe comorbidities: (1) disability-adjusted life years (DALYs); (2) share of total disease burden; (3) age-standardized death rates per 100,000; (4) share of deaths. We focused on known comorbidities such as cardiovascular diseases, cancers, chronic respiratory diseases, diabetes mellitus, chronic kidney diseases. We also included metrics related to factors that might impinge on the severity of respiratory syndromes, such as smoking and air pollution. We also used https://ourworldindata.org/ to retrieve information on demographic, economic and political indicators. Table S1, in the supplementary information, reports the list, description and source of the variables used here, according to the GATHER statement¹⁷.

**Statistical analyses.** We first aimed at exploring whether CFR varied among countries while controlling for the differences in the epidemic progression. To this purpose, we run a linear mixed model (LMM) where CFR was the dependent variable, time since 100th case (in days) and squared time since 100th case (as to model non-linear variation), country and the two-way interactions were the fixed factors. Country was also included as a random effect. For computational reasons, in this model the covariance structure of the R matrix was modeled using variance components, and degrees of freedom were computed using the between-within method. Only countries for which at least 10 days had elapsed between the record of the 100th case and 11th June 2020, and for which number of deaths was higher than 1 per 1,000,000 inhabitants were included in this model, to allow a better estimate of the variation between CFR and time. This model included 143 countries and 8441 daily values of CFR. In a second model, we also included the number of tests performed × 1000 as to control for differences in testing strategies among countries, and a stringency index describing the severity of the social distancing rules adopted by each country. This reduced the number of countries included in the model to 72 and the number of total observations to 3778.

To assess the pattern of association between country-specific CFR and comorbidities, we ran four LMMs that included different metrics as fixed factors. The first model included DALYs lost to cardiovascular, cancer and chronic respiratory diseases (due to the strong correlations among variables, we summed the three DALYs and used the sum in the model), and DALYs lost × 100,000 for people older than 70 years. The second model included share of disease burden (due to cardiovascular, cancer and chronic respiratory diseases). The third model included age-standardized death rates (due to cardiovascular diseases, cancer, air pollution, ambient particulate matter pollution, and smoking for people older than 70 years). The last model included share of deaths (due to cardiovascular, cancer, chronic respiratory, kidney diseases, lower respiratory infections, diabetes mellitus, outdoor air pollution). The four models always included a set of covariates that described demographic, economic and political variables, namely population size, share of the population over 70 years, gross domestic product (GDP) per capita, total health care expenditure as share of GDP, number of hospital beds (×1000 inhabitants), political regime, stringency index, and the number of tests performed × 1000. Political regime is scored according to the level of democracy, between – 10 (full autocracy) to + 10 (full democracy) (Polity IV as reported in https://ourworldindata.org/). The stringency index describes the severity of the policies implemented by each country to limit the spread of the virus, according to the following categories: school closure, workplace closure, public events cancelled, restrictions on gatherings, public transport closure, public information campaigns, stay at home, restriction on internal movements, international travel controls, testing policy, contact tracing. Finally, all models also included time since 1 death × 1,000,000 was reached (in days), and squared time as to take into account the progression of the epidemic in each country. In each model, we tested the interaction between the different comorbidity, demographic, economic factors and time. This allowed us to ascertain whether the temporal changes in CFR differed as a function of country-specific comorbidity, demographic, and socio-economic factors. All variables were standardized with mean = 0 and standard deviation = 1 to make parameter estimates directly comparable (variables with asymmetrical distribution were previously log-transformed to reduce skewness). To
<table><thead><tr><th>Fixed effects</th><th>F</th><th>df</th><th>p</th></tr></thead><tbody><tr><td>Country</td><td>2484.95</td><td>142, 8012</td><td>&lt;0.0001</td></tr><tr><td>Time since 100th case (days)</td><td>0.10</td><td>1, 8012</td><td>0.7518</td></tr><tr><td>Squared time since 100th case</td><td>0.01</td><td>1, 8012</td><td>0.9293</td></tr><tr><td>Country × time since 100th case</td><td>193.48</td><td>142, 8012</td><td>&lt;0.0001</td></tr><tr><td>Country × squared time since 100th case</td><td>109.93</td><td>142, 8012</td><td>&lt;0.0001</td></tr></tbody></table>

**Table 1.** Linear mixed model exploring variation of COVID-19 case fatality rate (CFR) as a function of time since the 100th case for each country. The model also included squared time since 100th case and the interactions between country and time. Country was also declared as a random effect in the model. The model was restricted to countries that had 10 or more days elapsed between the occurrence of the 100th case and 11th June 2020 and for which number of deaths was higher than 1 per 1,000,000 inhabitants. The analysis is based on 143 countries and 8441 observations. Significant p-values are in bold.

<figure><img src="image_4.png" /></figure>

**Figure 1.** Time-dependent variation in COVID-19 case fatality rate (CFR) among countries. Time refers to the period between 30 and 90 days post 100th case. For illustrative reasons, only 20 countries are reported here.

take into account the covariation between observations at different geographical scale, the LMMs also included three nested effects as random variables: continent, geographic region¹⁸ nested within continent, country nested within geographic region nested within continent. The covariance structure of the R matrix was modeled using a first-order autoregressive structure to model the temporal autocorrelation between daily CFR values. Degrees of freedom were computed using the Satterthwaite approximation. In order to minimize the risk of false discovery rate, we set the value of α to 0.01. Therefore, only p values lower than 0.01 were considered as indicative of statistically significant associations.

All the analyses were conducted with SAS 14.3 (PROC MIXED).

## Results

We found strong evidence for among-country differences in COVID-19 CFR. The LMM showed a highly significant interaction between time since the 100th case and country, indicating that CFR trajectories (both linear and quadratic components), did differ among countries as the epidemic progressed (Table 1, Fig. 1). CFR is defined as the ratio between number of deaths and number of confirmed cases. Since many cases of asymptomatic infection (or infections with mild symptoms) might get unnoticed, CFR is certainly an overestimate of the actual risk of dying when infected with SARS-Cov-2. In the light of this argument, a strategy of massive screening of the population (not only patients who are admitted to the hospitals, but also those with no or mild symptoms) might provide a more realistic estimate of the denominator of the CFR. Although the number of tests performed over the course of the epidemic is available for a smaller number of countries, it nevertheless varies tremendously among them, indicating different screening strategies adopted by each country. Focusing, for instance, on day 80 since the 100th case, and including only countries with more than 1 death per 1,000,000 inhabitants (n = 50 countries), number of tests performed per 1000 inhabitants varied from 0.9 (Indonesia) to 179 (Iceland). Interestingly,
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Changes in COVID-19 CFR as a function of the number of tests performed (×1000). For illustrative reasons, we report some representative countries showing how the relationship between CFR and number of tests can vary from negative to positive.</figcaption></figure>

however, countries that better screened their population did not always suffer from the lowest CFR. At day 80 since the 100th case, there was no correlation between number of tests per 1000 inhabitants and CFR (Pearson's r = -0.101, *p* = 0.4873, n = 50). Using the whole dataset, and including the number of tests ×1000, showed that the relationship between CFR and number of tests was country specific, with some countries where screening resulted in a decline in CFR, whereas in other the relationship was positive (interaction between country and number of tests ×1000: F<sub>71,3782</sub> = 23.17, *p* < 0.0001; Fig. 2). Importantly, the interaction between time and country remained highly significant even when adding both the number of tests ×1000 and the stringency index in the model (time × country: F<sub>69,3429</sub> = 43.50, *p* < 0.0001; squared time × country: F<sub>69,3429</sub> = 45.88, *p* < 0.0001). This result, therefore, suggests that variation in CFR among countries does not merely depend on the screening effort provided by each country, nor on the severity of social distancing and isolation rules.

In order to uncover the factors that might account for such among-country variation in CFR, we ran four LMMs that included different fixed factors related to comorbidities, demographic, economic and political variables. These models provided evidence for several associations between country-specific risk factors and the temporal dynamics of CFR (Table 2). In particular, countries with high values of DALYs lost to cardiovascular, cancer and chronic respiratory diseases had high CFR (model 1, Table 2, Fig. 3A). Countries with high share of disease burden due to chronic respiratory diseases had high CFR (model 2, Table 2). Countries with the lowest death rate due to cardiovascular diseases had the lowest CFR and countries with the highest death rate due to smoking in people over 70 year-old had the highest CFR (model 3, Table 2, Fig. 3B). Finally, model 4 showed that CFR was positively associated with high share of death due to chronic respiratory, cardiovascular, kidney diseases and outdoor air pollution (Table 2, Fig. 3C). Interestingly, this model also showed a negative correlation between CFR and share of death due to lower respiratory infections (Table 2, Fig. 3D).

Among the demographic, economic and political variables, the four models consistently provided evidence for positive associations between the temporal dynamic of CFR and population size, GDP per capita, total health expenditure as share of GDP, share of the population over 70 years, stringency index (Table 3, Fig. 4A–C). As mentioned above, number of tests performed ×1000 showed a more complex pattern of association with CFR (Table 3, Fig. 4D). Two models (out of four) provided evidence for a negative association between number of hospital beds ×1000 and CFR (Table 3, Fig. 4E), while only one model showed a positive association between CFR and political regime, with democracies having the highest CFR (Table 3, Fig. 4F).

## Discussion

As the SARS-CoV-2 epidemic continues to spread worldwide, the mortality induced by the disease has been a serious matter of concern, with some countries paying a high toll to the infection. In this light, it is important to understand why some countries seem to experience lower mortality rate than others and possibly to uncover the associated factors. Here, we showed that CFR greatly differs among countries (even after controlling for the number of tests and the severity of social distancing rules), and several comorbidity (DALYs lost to cardiovascular, cancer and chronic respiratory diseases; death rate due to smoking in people older than 70; share of deaths due to cardiovascular, chronic respiratory, kidney diseases), and socio-economic factors (population size, GDP
<table><thead><tr><th>Fixed effects</th><th>Estimate</th><th>SE</th><th>t</th><th>p</th><th>95% CI</th></tr></thead><tbody><tr><td colspan="6"><strong>Model 1 (DALYs)</strong></td></tr><tr><td>Intercept</td><td>1.243</td><td>2.000</td><td></td><td></td><td></td></tr><tr><td>Time</td><td>-1.607</td><td>0.257</td><td>-6.26</td><td><strong>&lt;0.0001</strong></td><td>-2.110/-1.104</td></tr><tr><td>Time<sup>2</sup></td><td>1.267</td><td>0.261</td><td>4.85</td><td><strong>&lt;0.0001</strong></td><td>0.754/1.779</td></tr><tr><td>DALYs lost to cardiovascular, cancer and chronic respiratory diseases</td><td>-0.569</td><td>5.381</td><td>-0.11</td><td>0.9161</td><td>-11.344/10.206</td></tr><tr><td>DALYs lost × 100,000 for people older than 70 year-old</td><td>0.161</td><td>0.943</td><td>0.17</td><td>0.8647</td><td>-1.728/2.051</td></tr><tr><td>Time × DALYs lost to cardiovascular, cancer and chronic respiratory diseases</td><td>1.691</td><td>0.876</td><td>1.93</td><td>0.0537</td><td>-0.027/3.408</td></tr><tr><td>Time × DALYs lost × 100,000 for people older than 70 year-old</td><td>0.164</td><td>0.190</td><td>0.87</td><td>0.3866</td><td>-0.208/0.537</td></tr><tr><td>Time<sup>2</sup> × DALYs lost to cardiovascular, cancer and chronic respiratory diseases</td><td>-3.022</td><td>1.002</td><td>-3.02</td><td><strong>0.0026</strong></td><td>-4.986/-1.058</td></tr><tr><td>Time<sup>2</sup> × DALYs lost × 100,000 for people older than 70 year-old</td><td>0.169</td><td>0.228</td><td>0.74</td><td>0.4604</td><td>-0.279/0.617</td></tr><tr><td><strong>Covariance parameters</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Variance</td><td>7.212</td><td>1.363</td><td>5.29</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td>First-order autoregression</td><td>-0.107</td><td>0.146</td><td>-0.73</td><td>0.4633</td><td></td></tr><tr><td>Residual</td><td>0.756</td><td>0.018</td><td>41.86</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td colspan="6"><strong>Model 2 (share of disease burden)</strong></td></tr><tr><td>Intercept</td><td>1.498</td><td>1.417</td><td></td><td></td><td></td></tr><tr><td>Time</td><td>-2.014</td><td>0.159</td><td>-12.71</td><td><strong>&lt;0.0001</strong></td><td>-2.324/-1.703</td></tr><tr><td>Time<sup>2</sup></td><td>1.861</td><td>0.142</td><td>13.07</td><td><strong>&lt;0.0001</strong></td><td>1.582/2.141</td></tr><tr><td>Cardiovascular</td><td>0.080</td><td>0.836</td><td>0.10</td><td>0.9242</td><td>-1.597/1.757</td></tr><tr><td>Cancer</td><td>2.209</td><td>1.327</td><td>1.66</td><td>0.1017</td><td>0.451/4.869</td></tr><tr><td>Chronic respiratory</td><td>0.314</td><td>0.602</td><td>0.52</td><td>0.6037</td><td>-0.892/1.520</td></tr><tr><td>Time × cardiovascular</td><td>-0.158</td><td>0.100</td><td>-1.58</td><td>0.1151</td><td>-0.354/0.039</td></tr><tr><td>Time × cancer</td><td>0.457</td><td>0.197</td><td>2.33</td><td>0.0201</td><td>0.072/0.843</td></tr><tr><td>Time × chronic respiratory</td><td>0.311</td><td>0.093</td><td>3.35</td><td><strong>0.0008</strong></td><td>0.129/0.493</td></tr><tr><td>Time<sup>2</sup> × cardiovascular</td><td>0.183</td><td>0.110</td><td>1.67</td><td>0.0960</td><td>-0.032/0.398</td></tr><tr><td>Time<sup>2</sup> × cancer</td><td>0.048</td><td>0.215</td><td>0.22</td><td>0.8234</td><td>-0.374/0.470</td></tr><tr><td>Time<sup>2</sup> × chronic respiratory</td><td>-0.630</td><td>0.098</td><td>-6.43</td><td><strong>&lt;0.0001</strong></td><td>-0.822/-0.438</td></tr><tr><td><strong>Covariance parameters</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Variance</td><td>6.978</td><td>1.330</td><td>5.25</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td>First-order autoregression</td><td>-0.104</td><td>0.148</td><td>-0.70</td><td>0.4835</td><td></td></tr><tr><td>Residual</td><td>0.738</td><td>0.018</td><td>41.84</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td colspan="6"><strong>Model 3 (age-standardized death rate)</strong></td></tr><tr><td>Intercept</td><td>1.192</td><td>1.250</td><td></td><td></td><td></td></tr><tr><td>Time</td><td>-1.561</td><td>0.171</td><td>-9.13</td><td><strong>&lt;0.0001</strong></td><td>-1.896/-1.226</td></tr><tr><td>Time<sup>2</sup></td><td>1.665</td><td>0.158</td><td>10.57</td><td><strong>&lt;0.0001</strong></td><td>1.356/1.974</td></tr><tr><td>Cardiovascular</td><td>-0.941</td><td>0.784</td><td>-1.20</td><td>0.2357</td><td>-2.513/0.632</td></tr><tr><td>Cancer</td><td>0.441</td><td>0.616</td><td>0.72</td><td>0.4775</td><td>-0.796/1.677</td></tr><tr><td>Air pollution</td><td>-0.866</td><td>1.842</td><td>-0.47</td><td>0.6402</td><td>-4.559/2.827</td></tr><tr><td>Ambient particulate matter pollution</td><td>0.178</td><td>1.252</td><td>0.14</td><td>0.8877</td><td>-2.333/2.689</td></tr><tr><td>Smoking in people older than 70 years</td><td>1.446</td><td>0.616</td><td>2.35</td><td>0.0226</td><td>0.211/2.681</td></tr><tr><td>Time × cardiovascular</td><td>-0.349</td><td>0.108</td><td>-3.24</td><td><strong>0.0012</strong></td><td>-0.560/-0.138</td></tr><tr><td>Time × cancer</td><td>0.120</td><td>0.099</td><td>1.22</td><td>0.2244</td><td>-0.074/0.315</td></tr><tr><td>Time × air pollution</td><td>0.611</td><td>0.363</td><td>1.69</td><td>0.0920</td><td>-0.100/1.323</td></tr><tr><td>Time × ambient particulate matter pollution</td><td>0.089</td><td>0.253</td><td>0.35</td><td>0.7252</td><td>-0.407/0.585</td></tr><tr><td>Time × smoking in people older than 70 years</td><td>0.730</td><td>0.115</td><td>6.35</td><td><strong>&lt;0.0001</strong></td><td>0.505/0.955</td></tr><tr><td>Time<sup>2</sup> × cardiovascular</td><td>0.206</td><td>0.116</td><td>1.77</td><td>0.0769</td><td>-0.022/0.434</td></tr><tr><td>Time<sup>2</sup> × cancer</td><td>0.055</td><td>0.114</td><td>0.48</td><td>0.6299</td><td>-0.169/0.279</td></tr><tr><td>Time<sup>2</sup> × air pollution</td><td>-0.285</td><td>0.439</td><td>-0.65</td><td>0.5164</td><td>-1.145/0.575</td></tr><tr><td>Time<sup>2</sup> × ambient particulate matter pollution</td><td>-0.125</td><td>0.316</td><td>-0.40</td><td>0.6915</td><td>-0.744/0.494</td></tr><tr><td>Time<sup>2</sup> × smoking in people older than 70 years</td><td>-0.860</td><td>0.127</td><td>-6.77</td><td><strong>&lt;0.0001</strong></td><td>-1.110/-0.611</td></tr><tr><td><strong>Covariance parameters</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Variance</td><td>6.181</td><td>1.197</td><td>5.16</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td>First-order autoregression</td><td>-0.198</td><td>0.150</td><td>-1.32</td><td>0.1877</td><td></td></tr><tr><td>Residual</td><td>0.732</td><td>0.018</td><td>41.82</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td colspan="6"><strong>Model 4 (share of death)</strong></td></tr><tr><td>Intercept</td><td>1.533</td><td>1.314</td><td></td><td></td><td></td></tr><tr><td>Continued</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 2.** Linear mixed models investigating the association between COVID-19 case fatality rate (CFR) and several descriptors of comorbidities, demographics, economics and political regime for each country. Each model included the same demographic, economic and political regime variables (GDP per capita, population size, total health care expenditure as share of GDP, number of hospital beds ×1000 inhabitants, share of the population over 70 years, political regime, stringency index and number of tests performed ×1000). In addition, model 1 included DALYs lost to cardiovascular, cancer and chronic respiratory diseases, and DALYs lost ×100,000 for people older than 70 years. Model 2 included share of disease burden (cardiovascular, cancer and chronic respiratory diseases). Model 3 included age-standardized death rates ×100,000 due to cardiovascular diseases, cancer, air pollution, ambient particulate air pollution, and smoking over 70 years. Model 4 included share of deaths for cardiovascular diseases, cancer, chronic respiratory diseases, lower respiratory diseases, diabetes, and outdoor air pollution. All models included time since number of deaths ×1,000,000 higher than 1 and squared time. Three nested factors were also included as random factors (continent, region within continent and country within region within continent). The table reports parameter estimates (with SE and 95% CI), t and p values (in bold significant p-values at the 0.01 threshold) for the comorbidity factors. Sample size is five continents, 17 geographical regions, 67 countries and 3596 observations.</caption><thead><tr><th>Fixed effects</th><th>Estimate</th><th>SE</th><th>t</th><th>p</th><th>95% CI</th></tr></thead><tbody><tr><td>Time</td><td>-1.994</td><td>0.168</td><td>-11.88</td><td><strong>&lt;0.0001</strong></td><td>-2.324/-1.665</td></tr><tr><td>Time<sup>2</sup></td><td>2.017</td><td>0.153</td><td>13.22</td><td><strong>&lt;0.0001</strong></td><td>1.718/2.316</td></tr><tr><td>Cardiovascular</td><td>0.499</td><td>1.078</td><td>0.46</td><td>0.6456</td><td>-1.667/2.664</td></tr><tr><td>Cancer</td><td>1.366</td><td>1.139</td><td>1.20</td><td>0.2357</td><td>-0.920/3.652</td></tr><tr><td>Chronic respiratory</td><td>0.348</td><td>0.648</td><td>0.54</td><td>0.5937</td><td>-0.953/1.648</td></tr><tr><td>Lower respiratory infections</td><td>-0.085</td><td>0.596</td><td>-0.14</td><td>0.8869</td><td>-1.285/1.115</td></tr><tr><td>Kidney</td><td>0.263</td><td>0.723</td><td>0.36</td><td>0.7182</td><td>-1.190/1.715</td></tr><tr><td>Diabetes mellitus</td><td>1.169</td><td>0.648</td><td>1.80</td><td>0.0771</td><td>-0.132/2.469</td></tr><tr><td>Outdoor air pollution</td><td>-0.236</td><td>0.765</td><td>-0.31</td><td>0.7589</td><td>-1.771/1.299</td></tr><tr><td>Time × cardiovascular</td><td>-0.847</td><td>0.157</td><td>-5.41</td><td><strong>&lt;0.0001</strong></td><td>-1.153/-0.540</td></tr><tr><td>Time × cancer</td><td>-0.402</td><td>0.171</td><td>-2.34</td><td>0.0192</td><td>-0.738/-0.065</td></tr><tr><td>Time × chronic respiratory</td><td>0.720</td><td>0.106</td><td>6.81</td><td><strong>&lt;0.0001</strong></td><td>0.512/0.927</td></tr><tr><td>Time × lower respiratory infections</td><td>-0.601</td><td>0.101</td><td>-5.94</td><td><strong>&lt;0.0001</strong></td><td>-0.800/-0.403</td></tr><tr><td>Time × kidney</td><td>-0.403</td><td>0.100</td><td>-4.04</td><td><strong>&lt;0.0001</strong></td><td>-0.599/-0.208</td></tr><tr><td>Time × diabetes mellitus</td><td>-0.098</td><td>0.118</td><td>-0.83</td><td>0.4066</td><td>-0.330/0.134</td></tr><tr><td>Time × outdoor air pollution</td><td>0.374</td><td>0.111</td><td>3.35</td><td><strong>0.0008</strong></td><td>0.155/0.592</td></tr><tr><td>Time<sup>2</sup> × cardiovascular</td><td>0.241</td><td>0.183</td><td>1.32</td><td>0.1870</td><td>-0.117/0.600</td></tr><tr><td>Time<sup>2</sup> × cancer</td><td>0.294</td><td>0.197</td><td>1.49</td><td>0.1355</td><td>-0.092/0.679</td></tr><tr><td>Time<sup>2</sup> × chronic respiratory</td><td>-0.919</td><td>0.121</td><td>-7.59</td><td><strong>&lt;0.0001</strong></td><td>-1.156/-0.681</td></tr><tr><td>Time<sup>2</sup> × lower respiratory infections</td><td>0.332</td><td>0.117</td><td>2.84</td><td><strong>0.0045</strong></td><td>0.103/0.561</td></tr><tr><td>Time<sup>2</sup> × kidney</td><td>0.333</td><td>0.113</td><td>2.96</td><td><strong>0.0031</strong></td><td>0.112/0.554</td></tr><tr><td>Time<sup>2</sup> × diabetes mellitus</td><td>-0.068</td><td>0.134</td><td>-0.50</td><td>0.6137</td><td>-0.331/0.195</td></tr><tr><td>Time<sup>2</sup> × outdoor air pollution</td><td>-0.221</td><td>0.132</td><td>-1.68</td><td>0.0937</td><td>-0.479/0.037</td></tr><tr><td>Covariance parameters</td><td>Estimate</td><td>SE</td><td>z</td><td>p</td><td></td></tr><tr><td>Variance</td><td>6.878</td><td>1.358</td><td>5.07</td><td><strong>&lt;0.0001</strong></td><td></td></tr><tr><td>First-order autoregression</td><td>-0.206</td><td>0.145</td><td>-1.43</td><td>0.1535</td><td></td></tr><tr><td>Residual</td><td>0.711</td><td>0.017</td><td>41.80</td><td><strong>&lt;0.0001</strong></td><td></td></tr></tbody></table>

per capita, share of the population over 70, number of hospital beds ×1000, political regime) were positively associated with COVID-19 CFR.

Assessing case fatality rate during an ongoing outbreak is particularly difficult for, at least two reasons¹¹,¹⁹. CFR at a given time might be an underestimate of the final mortality rate because cumulative number of deaths will eventually keep increasing as some patients are still hosted in intensive care units (i.e., right censoring). Conversely, CFR might overestimate actual mortality if a large fraction of infected people do not develop disease symptoms (or develop mild symptoms), get unnoticed and are not included into the population of confirmed cases. Several countries have implemented a strategy of massive screening, aiming at identifying positives and isolating them, as to avoid further spreading of the disease. Data on the number of tests performed per 1000 people, therefore, allow investigating whether countries with a massive screening policy were also those with the lowest CFR. The results showed no simple association between CFR and testing, with countries showing the expected negative relationship, while others showing a positive correlation. This suggests that differences in population screening are not enough to explain the tremendous variation in CFR among countries. As mentioned above, the other possible bias when computing CFR is that some of the patients who are still in intensive care units may eventually die, and therefore if the total number of confirmed cases remains unchanged, CFR might
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Figure 3.** Time-dependent variation in COVID-19 case fatality rate (CFR) according to comorbidity factors. Time refers to the number of days between the date of occurrence of 1 death ×1,000,000 and June 11th 2020. (A) DALYs lost to cardiovascular, cancer and chronic respiratory diseases; (B) death rate (×100,000) due to smoking in people over 70 years; (C) share of death due to chronic respiratory diseases; (D) share of death due to lower respiratory infections. The surfaces were generated using a smoothed spline interpolation on the predicted values of the LMMs described in the text. Darker colors indicated higher values of CFR. X- and Y-axis are standardized values, allowing to have similarly scaled axis.

still increase. However, visual inspection of the temporal trajectories of CFR shows that values have reached an asymptote for the vast majority of countries (Fig. 1). This suggests that, unless a second wave hits in the following weeks/months, we should not expect CFR to substantially vary as a consequence of residual mortality.

Despite the uncertainty associated with the estimation of CFR, its comparison between countries can provide useful insights into the heterogeneity in the burden paid to the disease. We showed that the trajectories of time dependent variation in CFR greatly differed among countries, while controlling for differences in testing strategies and stringency index. This quantitatively corroborates and statistically validates the intuition that some countries better dealt with the disease than others. The following step was to try to understand whether such heterogeneity arises as the consequence of predictable factors. Previous reports on the clinical outcome of the disease have identified two major factors associated with poor prognosis: age and the presence of comorbidities. Elderly people have been shown to suffer the highest mortality rate following infection with SARS-CoV-2<sup>11,19</sup>. Similarly, previous history of cardiovascular disorders, cancer and diabetes have been reported to substantially increase COVID-19 mortality risk<sup>20,21</sup>. We therefore predicted that countries with a higher share of elderly people and a higher incidence of known comorbidity factors might suffer from the highest CFR.

Our integrated modelling approach provided some evidence in support to these predictions. In particular, several metrics of comorbidity factors were positively associated with the temporal dynamics of CFR (DALYs lost to cardiovascular, cancer and chronic respiratory diseases; share of burden due to chronic respiratory diseases; death rate due to cardiovascular diseases; death rate due to smoking in people older than 70; share of death due to chronic respiratory and kidney diseases). Interestingly, our model also showed that share of death due to lower respiratory infections was negatively associated with COVID-19 CFR. This negative association is
<table><caption>**Table 3.** Linear mixed models investigating the association between COVID-19 case fatality rate (CFR) and several descriptors of comorbidities, demographics, economics and political regime for each country. Each model included the same demographic, economic and political regime variables (GDP per capita, population size, total health care expenditure as share of GDP, number of hospital beds × 1000 inhabitants, share of the population over 70 years, political regime, stringency index and number of tests performed × 1000). In addition, model 1 included DALYs lost to cardiovascular, cancer and chronic respiratory diseases, and DALYs lost × 100,000 for people older than 70 years. Model 2 included share of disease burden (cardiovascular, cancer and chronic respiratory diseases). Model 3 included age-standardized death rates × 100,000 due to cardiovascular diseases, cancer, air pollution, ambient particulate air pollution, and smoking over 70 years. Model 4 included share of deaths for cardiovascular diseases, cancer, chronic respiratory diseases, lower respiratory diseases, diabetes, and outdoor air pollution. All models included time since number of deaths × 1,000,000 higher than 1 and squared time. Three nested factors were also included as random factors (continent, region within continent and country within region within continent). The table reports parameter estimates (SE), and p values for socio-economic factors in each model (in bold significant p-values at the 0.01 threshold). Sample size is five continents, 17 geographical regions, 67 countries and 3596 observations.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Model 1</th><th colspan="2">Model 2</th><th colspan="2">Model 3</th><th colspan="2">Model 4</th></tr><tr><th>Estimate (SE)</th><th>p</th><th>Estimate (SE)</th><th>p</th><th>Estimate (SE)</th><th>p</th><th>Estimate (SE)</th><th>p</th></tr></thead><tbody><tr><td>GDP per capita</td><td>0.427 (0.994)</td><td>0.6689</td><td>-0.525 (1.003)</td><td>0.6028</td><td>-0.090 (1.114)</td><td>0.9361</td><td>-0.712 (1.050)</td><td>0.5010</td></tr><tr><td>Population size</td><td>3.308 (6.469)</td><td>0.6110</td><td>2.840 (0.807)</td><td><strong>0.0009</strong></td><td>2.815 (0.726)</td><td><strong>0.0003</strong></td><td>3.039 (0.894)</td><td><strong>0.0013</strong></td></tr><tr><td>Total health care expenditure as share of GDP</td><td>0.155 (0.602)</td><td>0.7977</td><td>0.089 (0.592)</td><td>0.8809</td><td>-0.201 (0.609)</td><td>0.7432</td><td>-0.138 (0.605)</td><td>0.8203</td></tr><tr><td>Number of hospital beds×1000</td><td>-1.592 (0.913)</td><td>0.0867</td><td>-1.525 (0.942)</td><td>0.1111</td><td>-1.712 (0.821)</td><td>0.0417</td><td>-0.978 (0.924)</td><td>0.2950</td></tr><tr><td>Share of the population over 70 years</td><td>1.709 (1.274)</td><td>0.1852</td><td>0.421 (1.208)</td><td>0.7287</td><td>1.330 (1.008)</td><td>0.1929</td><td>1.465 (1.113)</td><td>0.1938</td></tr><tr><td>Political regime</td><td>-0.183 (0.599)</td><td>0.761</td><td>-0.952 (0.699)</td><td>0.1786</td><td>-1.069 (0.735)</td><td>0.1518</td><td>-0.989 (0.778)</td><td>0.2091</td></tr><tr><td>Stringency index</td><td>0.137 (0.047)</td><td><strong>0.0033</strong></td><td>0.127 (0.046)</td><td><strong>0.0057</strong></td><td>0.143 (0.046)</td><td><strong>0.0019</strong></td><td>0.111 (0.045)</td><td>0.0145</td></tr><tr><td>Number of tests performed×1000</td><td>0.894 (0.119)</td><td><strong>&lt;0.0001</strong></td><td>1.001 (0.116)</td><td><strong>&lt;0.0001</strong></td><td>0.628 (0.123)</td><td><strong>&lt;0.0001</strong></td><td>0.948 (0.118)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time×GDP per capita</td><td>2.768 (0.168)</td><td><strong>&lt;0.0001</strong></td><td>1.874 (0.174)</td><td><strong>&lt;0.0001</strong></td><td>2.518 (0.200)</td><td><strong>&lt;0.0001</strong></td><td>2.111 (0.178)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time×population size</td><td>-0.969 (1.046)</td><td>0.3541</td><td>1.030 (0.120)</td><td><strong>&lt;0.0001</strong></td><td>0.851 (0.121)</td><td><strong>&lt;0.0001</strong></td><td>0.466 (0.132)</td><td><strong>0.0004</strong></td></tr><tr><td>Time×total health care expenditure as share of GDP</td><td>0.658 (0.101)</td><td><strong>&lt;0.0001</strong></td><td>0.322 (0.097)</td><td><strong>0.0009</strong></td><td>0.753 (0.112)</td><td><strong>&lt;0.0001</strong></td><td>0.331 (0.098)</td><td><strong>0.0007</strong></td></tr><tr><td>Time×number of hospital beds×1000</td><td>-0.268 (0.101)</td><td><strong>0.0082</strong></td><td>-0.148 (0.099)</td><td>0.1351</td><td>-0.252 (0.090)</td><td><strong>0.0054</strong></td><td>0.265 (0.105)</td><td>0.0116</td></tr><tr><td>Time×share of the population over 70 years</td><td>0.451 (0.173)</td><td><strong>0.0091</strong></td><td>0.706 (0.162)</td><td><strong>&lt;0.0001</strong></td><td>0.888 (0.125)</td><td><strong>&lt;0.0001</strong></td><td>0.959 (0.143)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time×political regime</td><td>0.164 (0.102)</td><td>0.1072</td><td>-0.117 (0.110)</td><td>0.2904</td><td>0.009 (0.133)</td><td>0.9487</td><td>0.214 (0.129)</td><td>0.0990</td></tr><tr><td>Time×stringency index</td><td>0.188 (0.073)</td><td>0.0106</td><td>0.107 (0.073)</td><td>0.1420</td><td>0.189 (0.075)</td><td>0.0115</td><td>0.342 (0.075)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time×number of tests performed×1000</td><td>-1.016 (0.096)</td><td><strong>&lt;0.0001</strong></td><td>-1.038 (0.094)</td><td><strong>&lt;0.0001</strong></td><td>-0.954 (0.097)</td><td><strong>&lt;0.0001</strong></td><td>-1.106 (0.097)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time²×GDP per capita</td><td>-2.090 (0.190)</td><td><strong>&lt;0.0001</strong></td><td>-1.135 (0.201)</td><td><strong>&lt;0.0001</strong></td><td>-1.804 (0.229)</td><td><strong>&lt;0.0001</strong></td><td>-1.589 (0.209)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time²×population size</td><td>2.397 (1.198)</td><td>0.0454</td><td>-1.131 (0.125)</td><td><strong>&lt;0.0001</strong></td><td>-0.985 (0.129)</td><td><strong>&lt;0.0001</strong></td><td>-0.0646 (0.136)</td><td><strong>&lt;0.0001</strong></td></tr><tr><td>Time²×total health care expenditure as share of GDP</td><td>-0.258 (0.105)</td><td>0.0146</td><td>0.079 (0.102)</td><td>0.4352</td><td>-0.341 (0.116)</td><td><strong>0.0033</strong></td><td>-0.037 (0.104)</td><td>0.7207</td></tr><tr><td>Time²×number of hospital beds×1000</td><td>0.398 (0.101)</td><td><strong>&lt;0.0001</strong></td><td>0.190 (0.096)</td><td>0.0479</td><td>0.273 (0.092)</td><td><strong>0.0030</strong></td><td>-0.212 (0.109)</td><td>0.0513</td></tr><tr><td>Time²×share of the population over 70 years</td><td>0.382 (0.188)</td><td>0.0429</td><td>-0.368 (0.182)</td><td>0.0432</td><td>-0.193 (0.131)</td><td>0.1399</td><td>-0.300 (0.157)</td><td>0.0557</td></tr><tr><td>Time²×political regime</td><td>-0.328 (0.108)</td><td><strong>0.0024</strong></td><td>-0.157 (0.117)</td><td>0.1787</td><td>-0.264 (0.151)</td><td>0.0809</td><td>-0.256 (0.151)</td><td>0.0901</td></tr><tr><td>Time²×stringency index</td><td>-0.026 (0.075)</td><td>0.7331</td><td>0.109 (0.075)</td><td>0.1481</td><td>-0.040 (0.077)</td><td>0.6042</td><td>-0.120 (0.077)</td><td>0.1233</td></tr><tr><td>Time²×number of tests performed×1000</td><td>0.340 (0.110)</td><td><strong>0.0020</strong></td><td>0.386 (0.108)</td><td><strong>0.0004</strong></td><td>0.342 (0.112)</td><td><strong>0.0023</strong></td><td>0.487 (0.114)</td><td><strong>&lt;0.0001</strong></td></tr></tbody></table>

intriguing, especially in the light of recent reports of possible cross-immunity that might confer partial protection to SARS-CoV-2²².

All the above-mentioned associations between comorbidities and CFR hold while controlling for several potential confounding factors describing the socio-economic context of different countries. As such, positive associations between comorbidities and CFR at the country level, do not merely reflect differences in the structure of the age-pyramid, or the amount of resources allocated to the health care system. Actually, focusing on such demographic and socio-economic factors allowed us to identify several other variables that contributed to explain among-country variation in CFR. As predicted, countries with the highest share of elderly people (over 70) also had the highest CFR.

Economic parameters might equally well contribute to shape COVID-19 mortality. As the number of severe cases increases during the epidemic, the health care system can get overwhelmed and might be unable to receive and treat all those who need intensive care. Mortality might therefore results from health care systems that are inadequate to deal with large number of cases requiring simultaneous admittance in intensive care units. We used several proxies describing the investment of each country into the health care system and found a negative association between the number of hospital beds per 1000 inhabitants and CFR. However, seemingly in contradiction
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 4.** Time-dependent variation in COVID-19 case fatality rate (CFR) according to socio-economic factors. Time refers to the number of days between the date of occurrence of 1 death × 1,000,000 and June 11th 2020. (A) share of the population over 70 years; (B) GDP per capita; (C) stringency index; (D) number of tests × 1000; (E) number of hospital beds × 1000; (F) political regime. The surfaces were generated using a smoothed spline interpolation on the predicted values of the LMMs described in the text. Darker colors indicated higher values of CFR. X- and Y-axis are standardized values, allowing to have similarly scaled axis.</figcaption></figure>

to this view, we also found that CFR was highest in countries with high GDP per capita and high total health expenditure as share of GDP. While odd, this result corroborates the impression that wealthy countries in Europe and North America have paid a severe toll to the infection. Overall, the relationship between investment into health care system and CFR appear to be more complex than one might expect.

A final source of variation accounting for differences in CFR might be due to differential reports of number of deaths and/or confirmed cases between countries. This might reflect different counting/reporting methodology (e.g., testing strategy, deciding whether or not a given patient died because of COVID-19). In addition, most of the countries have been implementing social distancing protocols that differed in the severity of the restrictions imposed and on the timing of policy execution. Moreover, different populations might follow governmental instructions more or less loosely, due to the perceived risk/benefit of applying such instructions. For instance, social distancing and isolation might be more easily applied in countries with autocratic regimes that exert a more stringent control over the population. The role of social and cultural traits in the emergence of zoonotic diseases has already been discussed in the past, including the idea that collectivistic societies might have built as a way to better control epidemic waves<sup>23,24</sup>. We explored how political regime and the severity of isolation policies were associated with CFR. We found moderate evidence (1 out of 4 models) suggesting that countries with a democratic regime were those with the highest CFR. The analysis of the stringency index, describing the severity of the restrictions implemented by each country, showed that highest values of CFR were reached for intermediate values of the stringency index. This might reflect different processes. First, countries where the epidemic wave was relatively low (perhaps because of the factors described above) could have implemented relatively mild restriction policies, compared to countries where the epidemic got out of control and that required imposing more stringent social distancing rules.

Although we report here some associations between comorbidities, demographic, socio-economic variables and COVID-19 CFR, we fully acknowledge that these factors do not perfectly explain the variation in CFR among countries. This might come from the coarse grain of the analyses (country level), the error associated with the metrics used in our study, the role played by other factors not taken into account in our study, or the uncertainties associated with the estimation of CFR while the epidemic is still ongoing. If serological tests will be used on a very large scale to assess the proportion of the population that has been infected by the virus, we will have a better estimate of the mortality rate and the possible factors explaining the among-country heterogeneity. With this in mind, we nevertheless believe that our results stress the role of comorbidities, socio-economic and political factors as potential drivers affecting how a country deals with globally threatening epidemics.

## Data availability

The full dataset is available in the online appendix.
Received: 22 April 2020; Accepted: 8 October 2020

Published online: 03 November 2020

## References

1. Zhu, N. et al. A novel coronavirus from patients with pneumonia in China, 2019. *N. Engl. J. Med.* **382**, 727–733. https://doi.org/10.1056/NEJMoa2001017 (2020).

2. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. *Nature* https://doi.org/10.1038/s41586-020-2012-7 (2020).

3. WHO Director-General's opening remarks at the media briefing on COVID-19—11 March 2020. https://www.who.int/dg/speec hes/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19--11-march-2020.

4. Wu, Z. & McGoogan, J. M. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China. Summary of a report of 72314 cases from the Chinese center for disease control and prevention. *JAMA* **323**, 1239–1242. https://doi.org/10.1001/jama.2020.2648 (2020).

5. Wang, D. et al. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China. *JAMA* **323**, 1061–1069. https://doi.org/10.1001/jama.2020.1585 (2020).

6. Petrilli, C. M. et al. Factors associated with hospital admission and critical illness among 5279 people with coronavirus disease 2019 in New York City: prospective cohort study. *BMJ* **369**, m1966 (2020).

7. Zheng, S. et al. Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January–March 2020: retrospective cohort study. *BMJ* **369**, m1443 (2020).

8. Yang, X. et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China. *Lancet Respir. Med.* **8**, 475–481 (2020).

9. Lewnard, J. A. et al. Incidence, clinical outcomes, and transmission dynamics of severe coronavirus disease 2019 in California and Washington: prospective cohort study. *BMJ* **369**, m1923 (2020).

10. Zhou, F. et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. *Lancet* **395**, 1054–1062. https://doi.org/10.1016/S0140-6736(20)30566-3 (2020).

11. Verity, R. et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. *Lancet Infect. Dis.* **20**, 669–677 (2020).

12. CDC COVID-19 Response Team. Preliminary estimates of the prevalence of selected underlying health conditions among patients with coronavirus disease—United States, February 12–March 28, 2020. *MMWR Morb. Mortal. Wkly. Rep.* **2020**(69), 382–386. https://doi.org/10.15585/mmwr.mm6913e2 (2019).

13. Chen, N. et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. *Lancet* **395**, 507–513. https://doi.org/10.1016/S0140-6736(20)30211-7 (2020).

14. Huang, C. et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. *Lancet* **395**, 497–506. https://doi.org/10.1016/S0140-6736(20)30183-5 (2020).

15. Lescure, F.-X., Bouadma, L. & Nguyen, D. Clinical and virological data of the first cases of COVID-19 in Europe: a case series. *Lancet Infect. Dis.* **20**, 697–706. https://doi.org/10.1016/S1473-3099(20)30200-0 (2020).

16. Ji, Y., Ma, Z., Peppelenbosch, M. P. & Pan, Q. Potential association between COVID-19 mortality and health-care resource availability. *Lancet Glob. Health* **8**, E480. https://doi.org/10.1016/S2214-109X(20)30068-1 (2020).

17. Stevens, G. A. et al. Guidelines for accurate and transparent health estimates reporting: the GATHER statement. *Lancet* **388**, e19-23. https://doi.org/10.1016/S0140-6736(16)30388-9 (2016).

18. United Nations. Standard country or area codes for statistical use (M49). 2017. https://unstats.un.org/unsd/methodology/m49/.

19. Rajgor, D. D., Lee, M. H., Archuleta, S., Bagdasarian, N. & Quek, S. C. The many estimates of the COVID-19 case fatality rate. *Lancet Infect. Dis.* **20**, 776–777. https://doi.org/10.1016/S1473-3099(20)30244-9 (2020).

20. Ruan, Q., Yang, K., Wang, W., Jiang, L. & Song, J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. *Intensive Care Med.* **46**, 1294–1297. https://doi.org/10.1007/s00134-020-05991-x (2020).

21. Yang, X. et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. *Lancet Respir. Med.* **8**, 475–481. https://doi.org/10.1016/S2213-2600(20)30079-5 (2020).

22. Grifoni, A. et al. Targets of T cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals. *Cell* **181**, 1489–1501 (2020).

23. Fincher, C. L., Thornhill, R., Murray, D. R. & Schaller, M. Pathogen prevalence predicts human cross-cultural variability in individualism/collectivism. *Proc. R. Soc. Lond. B* **275**, 1279–1285 (2008).

24. Morand, S. & Walther, B. Individualistic values are related to an increase in the outbreaks of infectious diseases and zoonotic diseases. *Sci. Rep.* **8**, 3866. https://doi.org/10.1038/s41598-018-22014-4 (2018).

## Author contributions

G.S., B.F. and S.M. conceived the study; G.S. compiled the data and analyzed them; G.S., B.F., and S.M. interpreted the results; G.S. wrote the first draft of the manuscript; G.S., B.F. and S.M. revised the manuscript and approved the final version.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-75848-2.

Correspondence and requests for materials should be addressed to G.S.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_16.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020