COVID-19 case-fatality rate and demographic and socioeconomic influencers: a worldwide spatial regression analysis based on country-level data
==============================================================================================================================================

* Yang Cao
* Ayako Hiyoshi
* Scott Montgomery

## Abstract

We used the COVID-19 dataset obtained from the Our World in Data website and investigated the associations between COVID-19 CFR and nine country-level indices of 209 countries/territories using the Matern correlation regression model. Spatial dependence among the data was controlled using the latitude and longitude of the centroid of the countries/territories. Stratified analyses were conducted by economic level and COVID-19 testing policy. The average of country/territory-specific COVID-19 CFR is about 2-3% worldwide, which is higher than previously reported at 0.7-1.3%. Statistically significant associations were observed between COVID-19 CFR and population size and proportion of female smokers. The open testing policies are associated with decreased CFR. Strictness of anti-COVID-19 measures was not statistically significantly associated with CFR overall, but the higher stringency index was associated with higher CFR in higher income countries with active testing policies. The statistically significant association between population size and COVID-19 CRF suggests the healthcare strain and lower treatment efficiency in countries with large populations. The observed association between smoking in females and COVID-19 CFR might be due to that the proportion of female smokers reflected broadly income level of a country. When testing is warranted and healthcare resources are sufficient, strict quarantine and/or lockdown measures might result in excess deaths in underprivileged populations.

Key words
*   COVID-19
*   case-fatality rate
*   spatial regression analysis
*   global health
*   policy
*   demographic factor
*   socioeconomic factor

## Introduction

The pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-COV2 Virus/COVID-19), which was initially reported in Wuhan, a city in the Hubei province in China in December 2019, has become a major global health concern (Margallo II et al. 2020). Poor outcomes in those with COVID□19 infections correlate with clinical and laboratory features of cytokine storm syndrome, an exaggerated systemic inflammatory phenomenon due to over-production of proinflammatory cytokines by the immune system that results in diffuse inflammatory lung disease and acute respiratory distress syndrome (Henderson et al. 2020). It may be complicated by sepsis, respiratory failure, acute respiratory distress syndrome (ARDS), and subsequent multi-organ failure (Zaim et al. 2020). Although COVID-19 related deaths are not clearly defined in the international reports available so far, many governments are warning people to be particularly stringent in following the recommended prevention measures because COVID-19 may results in severe conditions that need critical care including ventilation or death (Jordan et al. 2020). Till the end of June, 2020, the pandemic has resulted in over 10 million confirmed cases and 500 thousand deaths worldwide (Worldometer 2020). COVID-19 related fatality rates are difficult to assess with certainty, but according to the estimates based on the data from China, the United Kingdom, Italy, and the Diamond Princess cruise ship, the overall death rate from the confirmed COVID-19 cases is around 0.7-1.3%, with a sharply rising from less than 0.002% in children aged 9 years or younger to 8% in people aged over 80, which is much greater than seasonal influenza at about 0.1% (Baud et al. 2020; Jordan et al. 2020; Mahase 2020; Onder et al. 2020; Rajgor et al. 2020; Russell et al. 2020).

During the COVID-19 pandemic, numerous studies on the global public health emergency have covered a significant range of disciplines including medicine, mathematics and social sciences. The spatial spread is one of the most important properties of COVID-19, a characteristic which mainly depends on the epidemic mechanism, human mobility, and control strategy (Gross et al. 2020). Spatial statistical methods are frequently used to uncover relationships between spatiotemporal patterns of infectious diseases and host or environmental characteristics (Lawson 2013), generate detailed maps to visualize the distribution of the diseases’ morbidity or mortality (Palk and Blower 2018; Zulu et al. 2014), and identify hotspots, clusters, and potential risk factors (Anselin et al. 2006; Kulldorff and Nagarwalla 1995; Pfeiffer et al. 2008; Sasaki et al. 2008).

Clinical risk factors for mortality of adult patients with COVID-19 have been investigated in numerous studies, and the identified factors include older age (Brooke and Jackson 2020), male sex (Jin et al. 2020), higher sequential organ failure assessment score (Zhou et al. 2020), obesity (Klang et al. 2020; Zhang et al. 2020), preexisting concurrent diabetes (Guo et al. 2020), cardiovascular, cerebrovascular (Du et al. 2020), and kidney diseases (Cheng et al. 2020b), and macro-economic and environmental risk factors such as socioeconomic deprivation (Bibbins-Domingo 2020), air pollution (Ogen 2020; Wu et al. 2020) and diurnal temperature variation (Ma et al. 2020), etc. However, there is a lack of published studies on the effects of demographic and socioeconomic factors on COVID-10 case-fatality rates at the population level. It is an important issue for governments and regional or international non-governmental organizations to identify country characteristics that are associated with high case fatality rate (CRF) and help developing prevention and intervention measures to fight against this global public health crisis.

Using the publicly available data from the non-governmental organization Our World in Data (Our World in Data 2020a), we aimed to investigate the relationship of key country level demographic and socioeconomic indices and COVID-19 case-fatality, and explore factors associated with CFR, which may indicate treatment efficiency and strain in healthcare resources, while controlling for the spatial dependence of the data collected at different locations.

The data used in the study are freely available in the open source for research. There are no individual data available in the dataset. The owner of the data gives the permission to use, distribute, and reproduce the data in any medium. Therefore, no private and confidential information could be disclosed in the study, and ethical approval is not applicable. In the study, estimates of health indicators at the global level were reported according to the GATHER statement (Supplemental material 1) (Stevens et al. 2016).

## Results

### Descriptive characteristics of the variables

In total, 209 countries and territories and 17 variables (including latitude and longitude) were included in the study. The descriptive statistics of the variables are shown in Tables 1. The pairwise relationships of the variables are shown in Figure 1. High multicollinearity was found between number of confirmed cases and population size (pairwise Pearson’s r = 0.76, P < 0.001, Figure 1), and gross domestic product (GDP) per capita and life expectancy, median age, proportion of aged 65 years or older, proportion of extreme poverty, and hospital beds per 1000 people (all pairwise Pearson’s r greater than or approximate to 0.70, and P < 0.001, Figure 1; multiple correlation coefficient = 0.92, P < 0.001). Therefore, median age, proportion of aged 65 or older, and hospital beds per 1000 people were excluded in later regression analysis. Although proportion of extreme poverty is an interesting factor to be investigated, it is highly correlated with GDP per capita (r = −0.83) and the latter is available in more countries, therefore proportion of extreme poverty was also excluded from analysis.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T1)

Table 1. 
Descriptive statistics of the variables

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F1)

Figure 1. 
Pairwise scatter plots and Pearson’s correlation coefficients of the variables \***|, P<0.001; **, P<0.01; *, P<0.05

### Worldwide COVID-19 CFR distribution

Distribution of COVID-19 case fatality rates of the 209 countries/territories is show in Figure 2. The mean and median CFR worldwide are 3.31% and 2.19%, respectively (Table 1), with the highest rates found in Yeman (27%), West and North Europe (14-19%), and North America (9-12%, Figure 2).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F2)

Figure 2. 
Global COVID-19 case-fatality rates (%) nan: no data available

### Spatial autocorrelation of the COVID-19 CFR

Statistically significant spatial autocorrelation was found among the countries/territories’ COVID-19 case-fatality rates. The residuals from the common (non-spatial) multivariate linear regression models show apparent spatial dependence around the countries/territories with high fatality (Figure 3). The P-value from the Moran’s I test for the spatial autocorrelation of the residuals is 2.32×10−5.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F3)

Figure 3. 
Residuals of the common (non-spatial) multivariate linear regression

The estimated spatial autocorrelation coefficient of COVID-19 case-fatality rates between two locations against their distance is shown in Figure 4, with a strength parameter *ν* of 2.48 and a decay parameter *ρ* of 0.12. Basically, locations more than 20 degrees (in longitude or latitude) away have an autocorrelation coefficient below 0.5 (Figure 4).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F4)

Figure 4. 
Strength and decay of the spatial autocorrelation between pair of locations

### Associations of demographic and socioeconomic variables with COVID-19 CFR

Overall, after controlling for the spatial dependence, the statistically significant variables associated with COVID-19 case-fatality rate are population size and proportion of female smokers in a country/region (Table 2). The multivariate adjusted results indicate that, approximately, a doubling in size of population is associated with 0.48 (95% confidence interval (CI): 0.25, 0.70) percent increase in COVID-19 case-fatality rate, and a doubling in proportion of female smokers is associated with 0.55 (95% CI: 0.09, 1.02) percent increase in in COVID-19 case-fatality rate. Open public testing policy is associated with decreased case-fatality rate (Beta = 2.23, 95% CI: −4.25, −0.21) percent compared to no testing policy. Although not statically significant, stringency index showed the smallest coefficient with COVID-19 case-fatality rate (Beta = −0.02, 95% CI: −0.08, 0.03) among all the investigated variables (Table 2).

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T2)

Table 2. 
Estimation for variables’ regression coefficients (Beta) for CFR

The estimated contour of COVID-19 case-fatality rates worldwide is shown in Figure 5. The areas with the higher risks are mainly around North America and West Europe.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F5)

Figure 5. 
Contour plot of estimated COVID-19 case-fatality rate (%)

The subgroup analysis by economic level indicates that population size, CVD death rate, diabetes prevalence, and open public testing policy; population size, population density, diabetes prevalence, and testing only symptomatic and specified policy and testing anyone symptomatic policy; and population size, and population density are statistically significantly associated with COVID-19 case-fatality rate in lower-middle, upper-middle, and high income countries, respectively (Table 2).

However, the subgroup analysis in upper-middle-income and high-income countries by testing policy indicates that, if testing was ensured (testing policy =2 or 3), stringency index is associated with increased case-fatality rate (Beta = 0.14, 95% CI: 0.01, 0.27) (Table 3). The find suggests an open society policy would be good when testing is ensured. Meanwhile, diabetes prevalence is inversely associated with case-fatality rate.

View this table:
[Table 3.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T3)

Table 3. 
Estimation for variables’ regression coefficients (Beta) for CFR by testing policy in upper-middle-income and high-income countries

## Discussion

### Geospatial analysis of the COVID-19 pandemic

The COVID-19 pandemic is still full of unknowns, and many of them have a spatial dimension that lead to understanding the emergency geographically. Analysis of the COVID-19 data requires an interdisciplinary approach, including spatial statistics that may provide solutions to address the spatial issues in the pandemic(Meade 2014). A recently published review summarized studies by May 1, 2020 on geospatial and spatial-statistical analysis of the COVID-19 pandemic. In total 63 studies were categorized into five subjects: spatiotemporal analysis, health and social geography, environmental variables, data mining, and web-based mapping (Franch-Pardo et al. 2020). Although 15 of the studies address the question globally, none of them investigate the association of COVID-19 related deaths with country level demographic and/or socioeconomic factors. From a global health perspective, there is apparent knowledge gap in the research field. Our geospatial analysis fills this gap and help to understand the consequences of COVID-19 and their related factors from a global level, and contribute to the predictive modelling and decision-making to combat the pandemic.

### Population size and density and COVID-19 CFR

Our results indicate that larger population size is the most consistently associated with higher COVID-19 CFR, but population density was not associated with the outcome, controlled for some demographic and socioeconomic variables and spatial dependence worldwide. The association of population size with CFR can be interpreted at least two ways. It could be because larger countries experienced greater number of death and conducted relatively less test compared with countries with smaller population. Alternatively, in larger countries, healthcare might have been strained and resulted in relatively higher number of deaths among confirmed cases than smaller countries (World Health Organization 2020a). We are unable to disentangle the mechanisms of the association. Therefore, it is recommended that the analysis is replicated by a study with more detailed healthcare information on both individual- and country-levels. However, in national studies, higher density has been shown to associate with higher Covid-19 prevalence in Japan (Bassino and Ladmiral 2020), Italy(Sjodin et al. 2020), and Iran (Ahmadi et al. 2020). The lack of the association of COVID-19 CFR with population density globally might be due to the confounding by testing strategies and economic levels. In countries like Germany and South Korea, which took more active testing strategy than, for example, the United Kingdom (UK) where polymerase chain reaction (PCR) test for COVID-19 was only performed among those who were with severe symptoms and hospital admitted at the beginning (Cheng et al. 2020a), CFR naturally show lower values. There are weak but statistically significant negative associations between population size and population density (r = −0.15, p = 0.042) and GDP per capita (r = −0.20, p = 0.006). In order to minimize the confounding, we conducted stratified analysis by economic level (Table 2) and testing policy (only within upper-middle-income and high-income countries, Table 3). The results indicate that in high-income countries, higher population density was associated with increased COVID-19 CFR (Table 2). Furthermore, we conducted a sensitivity analysis (Appendix) for CDR and the results were similar (Table S1 and S2). The results suggests that, globally, healthcare strain should be firstly relieved and treatment efficiency should be improved in countries with large populations.

### Economic level and COVID-19 CFR

In our analysis, high COVID-19 CFR was found mainly around North America and West Europe (Figure 1). One of the possible reasons might be that these countries counted COVID-19 deaths by including those who died with it, not only from it (Onder et al. 2020; Schellekens and Sourrouille 2020). Determination of COVID-19 deaths also differed by country. Some countries recorded a COVID-19 death as any death once the patient became a confirmed case, even the death happened after 2 months possibly by other reasons (such as an accident), while in some other countries, a COVID-19 death was recorded as the death occurring within a certain time period (ranging from 2 to 8 weeks) after COVID-19 symptom onset (Our World in Data 2020b). Furthermore, the extent that the counting covered home, institutions, and hospitals in high-income countries is different with that in low-income countries (Schellekens and Sourrouille 2020).

It has been reported in earlier studies that CFR was more favorable in low-income countries (Cash and Patel 2020; Schellekens and Sourrouille 2020). There are three possibilities to explain this unusual pattern: it may be because of younger population, poor data quality, or it was still the early stage (till we wrote this paper) of COVID-19 infection (Schellekens and Sourrouille 2020). There is a tight relationship between income level of a country and demographic structure. For example, many of African countries were classified as low-income, with median age of 20 years, and 61% of population of them is 24 years or younger, and merely 3% is equal to or older than 65 years in 2015 (United Nations 2019). It has been shown that younger age is associated with lower likelihood of severe COVID-19 (Jordan et al. 2020; Li et al. 2020; Team 2020). However, age-related explanation may not be satisfactory given other factors that are usually associated with higher spread and severity of COVID-19, lower treatment efficiency, and higher healthcare strain. In terms of absolute count, older population in low income countries is larger than that in high income countries, and the prevalence of risk factors such as lack of hygiene facilities, handwashing soap, and water are greater (Schellekens and Sourrouille 2020). Higher viral load has been suggested to be linked to more severe disease (Zheng et al. 2020). Healthcare resources are usually low in low-income countries (Fryatt et al. 2010). Even before the pandemic, developing countries particularly had challenges to collect, verify, and aggregate reliable data in timely manner due to lack of resource, communication and technological development (Bram et al. 2015). And the pandemic might have accentuated the pre-existing challenges (Schellekens and Sourrouille 2020). The extent of bias is difficult to know, including whether it is still in early stage of infection in developing countries. Therefore, other possibilities, such as data quality and the stage of infection spread may need to be considered in interpreting the results and need further investigation in the future.

### Proportion of smokers and COVID-19 CFR

The proportion of female smokers was positively associated with COVID-19 CFR in the overall analysis, but the association diminished when the analysis was stratified by economical level of the countries/territories. This is counterintuitive given that severer COVID-19 was associated with male sex due to possibly immune system and hormone levels (Channappanavar et al. 2017) and smoking (Cai 2020; Vardavas and Nikitara 2020). The observed association between smoking in females and COVID-19 CFR might be due to that the proportion of female smoker reflected broadly income level of a country (Figure 6a). Linking to the theory of diffusion of innovation, the spread of smoking habit has been illustrated to take several stages of rise, leveling and decline, from rich to poor, men to women, and young to old (Huisman et al. 2005; Khlat et al. 2016; Pampel 2006). In the early phrase, the prevalence of smoking increase in men, and women take up smoking about a few decades later. Subsequently, male smoking start to decline and female follows later on. This pattern was found to spread from rich to poor countries. In general, Asian and African countries tended to have low female smoking but high male smoking, while in European countries prevalence was similar between men and women (Gallus et al. 2020; Hitchman and Fong 2011). Therefore, female smoking in the overall analysis was a marker of the development of a country and it diminished when the analysis was stratified by it.

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F6/graphic-11.medium.gif)

[](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F6/graphic-11)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F6/graphic-12.medium.gif)

[](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F6/graphic-12)

Figure 6. 
Contour plot of projected (a) proportion of female smoker, (b) population, (c) GDP per capita, and (d) stringency index

### Stringency of measures against the COVID-19 epidemic and COVID-19 CFR

According to the currently available data from the Oxford COVID-19 Government Response Tracker (Hale et al. 2020), South America and Asian took the strictest measures (Table 4 and Figure 6d), and they also had relatively lower COVID-19 CFR (Table 4). However, in our multivariable analysis controlled for other variables and spatial dependence, we did not observe a statistically significant association between stringency index and COVID-19 CRF. In contrast, stricter measures were even found being associated with higher CFR in upper-middle-income and high-income countries with active testing policies (Table 3), which seems supporting the current argument that lockdown measures might result in excess deaths in underprivileged populations and those in need are hit harder by the crisis (Melnick and Ioannidis 2020). So far, the evidence that stricter response reduced healthcare strain or treatment efficiency reflected by COVID-19 CFR globally is lacking. However, the findings need to be further examined by comparing the all-cause mortalities in previous years. Meanwhile, the reliability of the stringency index also needs to be further investigated. The relationship between socioeconomic measures against the pandemic and COVID-19 CFR is a complicated issue, which needs deeper spatiotemporal analysis with more detailed and reliable information in the future.

View this table:
[Table 4.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T4)

Table 4. 
Stringency of measures against COVID-19 epidemic and COVID-19 CRF by continent

Noticeably, we also observed negative associations between COVID-19 CFR and CVD death rate and diabetes prevalence in some analysis, which might be partially explained by the competing risk between the deaths and/or comorbidities (Berry et al. 2010), because most of COVID-19 deaths are elderly and have one or more comorbidities (Garnier-Crussard et al. 2020; Richardson et al. 2020; Yang et al. 2020). Therefore, the COVID-19 CFR worldwide deserves deeper investigation with more detailed comorbidity information.

### Strengths and Limitations

To our knowledge, this is the first study that investigated relationship between COVID-19 CFR and demographic and socioeconomic factors globally. The study included in total 10,445,656 confirmed COVID-19 cases and 511,030 deaths worldwide. Although numerous studies have investigated the aforementioned factors related to the COVID-19 CFR, either they used much smaller sample size and investigated the question locally, or they did not approach this issue from a geospatial perspective. Our study may inspire new reflections from the healthcare workers to work together against the COVID-19 pandemic geographically and globally. International comparison of CFR may be challenging when the ascertainment of COVID-19 cases differed by country. To tackle this, we performed a sensitivity analysis using CDR. Although some risk factors, such as CVD and diabetes, showed different pattern of association, population showed consistent and positive association (see Appendix).

There are some limitations in our study. Firstly, the case-fatality analyzed here was based on the reported COVID-19 cases and deaths by countries/territories. According to the recent estimations, asymptomatic carriers and victims of COVID-19 could be as high as 10-80% in a population (Anastassopoulou et al. 2020; Day 2020; Kimball et al. 2020; Mizumoto et al. 2020; Nishiura et al. 2020; World Health Organization 2020b). However, this fraction was not taken into account in our analysis. Therefore, the case-fatality rates presented in the study might be significantly higher than the real ones. Secondly, age structure of population influences both prevalence and mortality of COVID-19, although we adjusted our analysis using the proportion of age over 60 years in populations, residual confounding largely remains. Thirdly, no detailed comorbidity information from the COVID-19 cases available in current study, which might bias the association towards an unknown direction. Finally, during an ongoing pandemic, delayed reporting occurs for both the number of cases and deaths, and strategies against the crisis also change time by time. Although we did the analysis two times using the data obtained on 17 June 17 and on July 2, and produced the same results, which suggests the bias due to delayed report might be negligible, the dynamic of the problem need to be addressed incorporating with temporal statistics methods.

## Conclusions

The average of country/territory-specific COVID-19 CFR is about 2-3% worldwide, which is higher than previously reported 0.7-1.3% and possibly due to the unreported asymptomatic cases. The COVID-19 CFR is statistically significantly associated with population size, especially in middle-income and high-income countries, which may indicate the healthcare strain and/or lower treatment efficiency in the countries with large populations. No statistically significant findings were found in low-income countries, which might be due to the challenges in data collection, communication, and verification in the countries and need to be further investigated in follow-up studies. To make global joint strategy and/or policy against the COVID-19 pandemic, spatial dependence and temporal trends must be consider in data analysis and decision making.

## Materials and methods

### Data on COVID-19 by Our World in Data

The COVID-19 dataset used in the study was downloaded from the Our World in Data website on July 2, 2020, which is a collection of the COVID-19 data maintained by the organization Our World in Data and updated daily. The dataset includes country level daily data on confirmed cases, deaths, and testing, as well as other variables of potential interest (Our World in Data 2020a; Roser et al. 2020). The data sources of the dataset including the European Centre for Disease Prevention and Control, the International Organization for Standardization (ISO), national government reports, the Department of Economic and Social Affairs of the United Nations (UN), UN Population Division, UN Statistics Division, Oxford COVID-19 Government Response Tracker, the World Bank, the Global Burden of Disease Collaborative Network, and Eurostat of the Organization for Economic Cooperation and Development (Our World in Data 2020a). There are in total 34 indices from 209 countries and territories in the dataset by July 2, 2020. The dataset was linked to the global geospatial vector database using the ISO 3166-1 alpha-2 codes for the spatial modelling (International Organization for Standardization (ISO) 2006).

### Case-fatality rate (CFR)

CFR of COVID-19 was calculated as the number of total deaths due to COVID-19 divided by the number of total confirmed COVID-19 cases by July 2, 2020, multiplied by 100. CFR was investigated in our study because it may reflect disease severity as well as efficiency of treatment and healthcare response and strain. CRF is not constant. It can vary between populations and over time, depending on the interplay between the causative agent of disease, the host, and the environment as well as available treatments and quality of patient care. For example, it can increase if the health care system is overwhelmed by the sudden increase of cases (Harrington 2020).

We also calculated crude cause-specific death rate (CDR) of COVID-19 in a sensitivity analysis and compared it with CFR. The CDR was calculated as the number of total deaths due to Covid-19 divided by the production of population and months of the data collected, multiplied by 1,000,000 (Centers for Disease Control Prevention 2012).

### Statistical analysis

The numbers of confirmed cases, tests, and tests per thousand people were not included in the analysis because they are much depending on the population, detection policy, and quarantine and isolation policy of the countries and territories. Instead, we included the stringency index in the analysis, which is a composite measure based on 9 response indicators including school closures, workplace closures, testing policy, and travel bans etc., rescaled to a value from 0 to 100 (100 = strictest response) (Haleet al. 2020). The stringency index data were obtained from the World Intellectual Property Organization website on July 1, 2020 (World Intellectual Property Organization 2020). Because the variable “proportion of the population with basic handwashing facilities on premises” has missing values in more than 50% of the countries/territories, it was also excluded from the analysis.

A subcomponent of the stringency index is the government policy on the access to COVID-19 test by 4 groups: 0, no policy; 1, only those who were both symptomatic and met specific criteria; 2, anyone symptomatic; 3, and open public testing, such as drive through testing (Haleet al. 2020). The testing policy indicator was used for stratification of the analysis.

Collinearity and multiple collinearity between the variables were examined using the Pearson’s correlation coefficient and multiple correlation coefficient, respectively (Mason and Perreault 1991). Spatial autocorrelation (or spatial dependence) is defined as the relationship between spatial proximity among some observational units and similarity among their values; positive spatial autocorrelation refers to situations in which the nearer the observational units, the more similar their values (and vice versa for its negative counterpart) (Lee 2017). This feature violates the assumption of independence among observations upon which many regression analyses are applied. Spatial autocorrelation among the fatality rates of the countries/territories was examined using a multivariate linear regression model and the Moran’s I test (Ren et al. 2014). The autocorrelation was visualized using the Matern correlation coefficient (Matérn 2013).

The Matern correlation model, a commonly used model for spatial correlated data, was fitted for our data to investigate the relationship between COVID-19 case-fatality and the demographic and socioeconomic variables. The latitude and longitude of the centroid of the countries/territories were used as random effects in the Matern correlation model (Ward and Gleditsch 2018).

Variables with skewed distribution were log transformed before entering the regression models. Multiple imputation method was used to handle the missing values in the data, and 10 imputed datasets were used for the regression analyses (Royston 2004).

Subgroup analysis was conducted by economic levels of the countries/territories according to the World Bank’s newest classification (World Bank 2020).

The associations of the studied variables with COVID-19 CDR (per 1,000,000 person months) of the countries/territories for December 31, 2019 to July 1, 2020 were also evaluated using the same methodology but using Poisson regression model, and the results were presented as appendix.

All the analysis were conducted in R 4.02 (the R Foundation for Statistical Computing, Vienna, Austria) using the package *spaMM* (Rousset and Ferdy 2014) and in Python 3.6 (Python Software Foundation) (van Rossum 1995) using the packages *geopandas* and *geoplot* (Jordahl et al. 2019).

## Data Availability

All data used in this study are publicly available and are referenced in the manuscript.

[https://ourworldindata.org/](https://ourworldindata.org/) 

## Competing interests

The authors declare no competing interests.

## Author contributions

YC initially designed the study. All authors refined the study design and developed the concept. YC collected and analyzed the data, and drafted the manuscript with AH. All authors interpreted data and critically revised the manuscript. All authors have approved the final version for publication.

## Data sharing

All data used in this study are publicly available and are referenced in the manuscript. The R code used to run the analyses is available on request.

## Supplementary files

*   GATHER statement checklist

## Appendix

![Figure S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F7.medium.gif)

[Figure S1.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F7)

Figure S1. 
Global COVID-19 deaths (per 1,000,000 person months)

![Figure S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F8.medium.gif)

[Figure S2.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F8)

Figure S2. 
Residuals of COVID-19 deaths per 1,000,000 person months based on the common (non-spatial) multivariate Poisson regression The P-value from the Moran’s I test for the spatial autocorrelation of the residuals < 0.001.

![Figure S3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F9.medium.gif)

[Figure S3.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F9)

Figure S3. 
Strength and decay of the spatial autocorrelation between pair of locations (COVID-19 deaths per 1,000,000 person months) The estimated spatial autocorrelation coefficient of COVID-19 case-fatality rates between two locations against their distance is shown in Figure S3, with a strength parameter *ν* of 0.11 and a decay parameter *ρ* of 0.01.

View this table:
[Table S1.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T5)

Table S1. 
Estimation for variables’ incidence rate ratio (IRR) for COVID-19 CDR

![Figure S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/31/2020.07.31.20165811/F10.medium.gif)

[Figure S4.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/F10)

Figure S4. 
Contour plot of estimated COVID-19 deaths (per 1,000,000 person months)

View this table:
[Table S2.](http://medrxiv.org/content/early/2020/07/31/2020.07.31.20165811/T6)

Table S2. 
Estimation for variables’ incidence rate ratio (IRR) for COVID-19 CDR by testing policy in upper-middle-income and high-income countries

*   Received July 31, 2020.
*   Revision received July 31, 2020.
*   Accepted July 31, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  Ahmadi, M.; Sharifi, A.; Dorosti, S.; Ghoushchi, S.J.; Ghanbari, N. Investigation of effective climatology parameters on COVID-19 outbreak in Iran. Science of the Total Environment 2020;729
    
    

2.  Anastassopoulou, C.; Russo, L.; Tsakris, A.; Siettos, C. Data-based analysis, modelling and forecasting of the COVID-19 outbreak. PLoS One 2020;15:e0230405
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0230405&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

3.  Anselin, L.; Syabri, I.; Kho, Y. GeoDa: An introduction to spatial data analysis. Geogr Anal 2006;38:5–22
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.0016-7363.2005.00671.x&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234927600002&link_type=ISI) 

4.  Bassino, J.-P.; Ladmiral, G. Socio-Economic Factors Influencing COVID-19 Spread in Japan; Evidence from the Second Wave (March-May 2020). Evidence from the Second Wave (March-May 2020)(June 30, 2020) 2020;[http://dx.doi.org/10.2139/ssrn.3637994](http://dx.doi.org/10.2139/ssrn.3637994)
    
    

5.  Baud, D.; Qi, X.; Nielsen-Saines, K.; Musso, D.; Pomar, L.; Favre, G. Real estimates of mortality following COVID-19 infection. Lancet Infect Dis 2020;
    
    

6.  Berry, S.D.; Ngo, L.; Samelson, E.J.; Kiel, D.P. Competing risk of death: an important consideration in studies of older adults. J Am Geriatr Soc 2010;58:783–787
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1532-5415.2010.02767.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20345862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

7.  Bibbins-Domingo, K. This Time Must Be Different: Disparities During the COVID-19 Pandemic. Ann Intern Med 2020;
    
    

8.  Bram, J.T.; Warwick-Clark, B.; Obeysekare, E.; Mehta, K. Utilization and Monetization of Healthcare Data in Developing Countries. Big Data 2015;3:59–66
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/big.2014.0053&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26487984&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

9.  Brooke, J.; Jackson, D. Older people and COVID-19: Isolation, risk and ageism. J Clin Nurs 2020;29:2044–2046
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

10. Cai, H. Sex difference and smoking predisposition in patients with COVID-19. Lancet Respir Med 2020;8:e20
    
    

11. Cash, R.; Patel, V. Has COVID-19 subverted global health? The Lancet 2020;395:1687–1688
    
    

12. Centers for Disease Control Prevention. Principles of epidemiology in public health practice ed^eds. Atlanda, US: US Department of Health and Human Services; 2012
    
    

13. Channappanavar, R.; Fett, C.; Mack, M.; Ten Eyck, P.P.; Meyerholz, D.K.; Perlman, S. Sex-Based Differences in Susceptibility to Severe Acute Respiratory Syndrome Coronavirus Infection. J Immunol 2017;198:4046–4053
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE5OC8xMC80MDQ2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMzEvMjAyMC4wNy4zMS4yMDE2NTgxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. Cheng, H.Y.; Li, S.Y.; Yang, C.H. Initial rapid and proactive response for the COVID-19 outbreak - Taiwan’s experience. J Formos Med Assoc 2020a;119:771–773
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

15. Cheng, Y.; Luo, R.; Wang, K.; Zhang, M.; Wang, Z.; Dong, L.; Li, J.; Yao, Y.; Ge, S.; Xu, G. Kidney disease is associated with in-hospital death of patients with COVID-19. Kidney Int 2020b;97:829–838
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.kint.2020.03.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

16. Day, M. Covid-19: four fifths of cases are asymptomatic, China figures indicate. BMJ 2020;369:m1375
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMDJfMTEvbTEzNzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8zMS8yMDIwLjA3LjMxLjIwMTY1ODExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

17. Du, R.H.; Liang, L.R.; Yang, C.Q.; Wang, W.; Cao, T.Z.; Li, M.; Guo, G.Y.; Du, J.; Zheng, C.L.; Zhu, Q.; Hu, M.; Li, X.Y.; Peng, P.; Shi, H.Z. Predictors of mortality for patients with COVID-19 pneumonia caused by SARS-CoV-2: a prospective cohort study. Eur Respir J 2020;55
    
    

18. Franch-Pardo, I.; Napoletano, B.M.; Rosete-Verges, F.; Billa, L. Spatial analysis and GIS in the study of COVID-19. A review. Sci Total Environ 2020;739:140033
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2020.140033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

19. Fryatt, R.; Mills, A.; Nordstrom, A. Financing of health systems to achieve the health Millennium Development Goals in low-income countries. Lancet 2010;375:419–426
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(09)61833-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20113826&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274305500029&link_type=ISI) 

20. Gallus, S.; Lugo, A.; Liu, X.; Behrakis, P.; Boffi, R.; Bosetti, C.; Carreras, G.; Chatenoud, L.; Clancy, L.; Continente, X.; Dobson, R.; Effertz, T.; Filippidis, F.T.; Fu, M.; Geshanova, G.; Gorini, G.; Keogan, S.; Ivanov, H.; Lopez, M.J.; Lopez-Nicolas, A.; Precioso, J.; Przewozniak, K.; Radu-Loghin, C.; Ruprecht, A.; Semple, S.; Soriano, J.B.; Starchenko, P.; Trapero-Bertran, M.; Tigova, O.; Tzortzi, A.S.; Vardavas, C.; Vyzikidou, V.K.; Colombo, P.; Fernandez, E.; Tack, S.H.S.P.I. Who smokes in Europe? Data from 12 European countries in the TackSHS survey (2017-2018). J Epidemiol 2020;
    
    

21. Garnier-Crussard, A.; Forestier, E.; Gilbert, T.; Krolak-Salmon, P. Novel Coronavirus (COVID-19) Epidemic: What Are the Risks for Older Patients? J Am Geriatr Soc 2020;68:939–940
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

22. Gross, B.; Zheng, Z.; Liu, S.; Chen, X.; Sela, A.; Li, J.; Li, D.; Havlin, S. Spatio-temporal propagation of COVID-19 pandemics. medRxiv 2020;2020.03.23.20041517
    
    

23. Guo, W.N.; Li, M.Y.; Dong, Y.L.; Zhou, H.F.; Zhang, Z.L.; Tian, C.X.; Qin, R.J.; Wang, H.J.; Shen, Y.; Du, K.Y.; Zhao, L.; Fan, H.; Luo, S.S.; Hu, D.S. Diabetes is a risk factor for the progression and prognosis of COVID-19. Diabetes-Metab Res 2020;
    
    

24. Hale, T.; Webster, S.; Petherick, A.; Phillips, T.; Kira, B. Oxford covid-19 government response tracker. Blavatnik School of Government 2020;25
    
    

25. Harrington, R.A. Case fatality rate in Encyclopedia Britannica. 2020
    
    

26. Henderson, L.A.; Canna, S.W.; Schulert, G.S.; Volpi, S.; Lee, P.Y.; Kernan, K.F.; Caricchio, R.; Mahmud, S.; Hazen, M.M.; Halyabar, O.; Hoyt, K.J.; Han, J.; Grom, A.A.; Gattorno, M.; Ravelli, A.; De Benedetti, F.; Behrens, E.M.; Cron, R.Q.; Nigrovic, P.A. On the Alert for Cytokine Storm: Immunopathology in COVID-19. Arthritis Rheumatol 2020;
    
    

27. Hitchman, S.C.; Fong, G.T. Gender empowerment and female-to-male smoking prevalence ratios. Bull World Health Organ 2011;89:195–202
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.10.079905&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21379415&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000288406200010&link_type=ISI) 

28. Huisman, M.; Kunst, A.E.; Mackenbach, J.P. Inequalities in the prevalence of smoking in the European Union: comparing education and income. Prev Med 2005;40:756–764
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ypmed.2004.09.022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15850876&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229006700021&link_type=ISI) 

29. International Organization for Standardization (ISO). ISO 3166-1: Codes for the representation of names of countries and their subdivisions–Part 1: Country codes. ISO Geneva, Switzerland; 2006
    
    

30. Jin, J.M.; Bai, P.; He, W.; Wu, F.; Liu, X.F.; Han, D.M.; Liu, S.; Yang, J.K. Gender Differences in Patients With COVID-19: Focus on Severity and Mortality. Front Public Health 2020;8:152
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

31. Jordahl, K.; Bossche, J.d.; Wasserman, J.; McBride, J.; Gerard, J.; Fleischmann, M.; Tratner, J. Geopandas/geopandas: V0. 6.1. Zenodo. doi; 2019
    
    

32. Jordan, R.E.; Adab, P.; Cheng, K.K. Covid-19: risk factors for severe disease and death. BMJ 2020;368:m1198
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjgvbWFyMjZfMTEvbTExOTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8zMS8yMDIwLjA3LjMxLjIwMTY1ODExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

33. Khlat, M.; Pampel, F.; Bricard, D.; Legleye, S. Disadvantaged Social Groups and the Cigarette Epidemic: Limits of the Diffusion of Innovations Vision. Int J Environ Res Public Health 2016;13
    
    

34. Kimball, A.; Hatfield, K.M.; Arons, M.; James, A.; Taylor, J.; Spicer, K.; Bardossy, A.C.; Oakley, L.P.; Tanwar, S.; Chisty, Z.; Bell, J.M.; Methner, M.; Harney, J.; Jacobs, J.R.; Carlson, C.M.; McLaughlin, H.P.; Stone, N.; Clark, S.; Brostrom-Smith, C.; Page, L.C.; Kay, M.; Lewis, J.; Russell, D.; Hiatt, B.; Gant, J.; Duchin, J.S.; Clark, T.A.; Honein, M.A.; Reddy, S.C.; Jernigan, J.A.; Public Health, S.; King, C.; Team, C.C.-I. Asymptomatic and Presymptomatic SARS-CoV-2 Infections in Residents of a Long-Term Care Skilled Nursing Facility - King County, Washington, March 2020. MMWR Morb Mortal Wkly Rep 2020;69:377–381
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm6913e1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32240128&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

35. Klang, E.; Kassim, G.; Soffer, S.; Freeman, R.; Levin, M.A.; Reich, D.L. Morbid Obesity as an Independent Risk Factor for COVID-19 Mortality in Hospitalized Patients Younger than 50. Obesity (Silver Spring) 2020;
    
    

36. Kulldorff, M.; Nagarwalla, N. Spatial disease clusters: detection and inference. Stat Med 1995;14:799–810
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.4780140809&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7644860&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QY45400007&link_type=ISI) 

37. Lawson, A.B. Bayesian disease mapping: hierarchical modeling in spatial epidemiology ed^eds: CRC press; 2013
    
    

38. 1.  Shekhar S.; 
    2.  Xiong H.; 
    3.  Zhou X.
    
    Lee, S.-I. Correlation and Spatial Autocorrelation. in: Shekhar S.; Xiong H.; Zhou X., eds. Encyclopedia of GIS. Heidelberg, Germany: Springer Nature; 2017
    
    

39. Li, X.; Xu, S.; Yu, M.; Wang, K.; Tao, Y.; Zhou, Y.; Shi, J.; Zhou, M.; Wu, B.; Yang, Z.; Zhang, C.; Yue, J.; Zhang, Z.; Renz, H.; Liu, X.; Xie, J.; Xie, M.; Zhao, J. Risk factors for severity and mortality in adult COVID-19 inpatients in Wuhan. J Allergy Clin Immunol 2020;146:110–118
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

40. Ma, Y.; Zhao, Y.; Liu, J.; He, X.; Wang, B.; Fu, S.; Yan, J.; Niu, J.; Zhou, J.; Luo, B. Effects of temperature variation and humidity on the death of COVID-19 in Wuhan, China. Sci Total Environ 2020;724:138226
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2020.138226&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32408453&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

41. Mahase, E. Covid-19: death rate is 0.66% and increases with age, study estimates. Bmj-Brit Med J 2020;369
    
    

42. Margallo II, L.N.; Diaz, M.; Lim, P.P. 2019 Novel Coronavirus Pandemic: What Do We Know? South Dakota Medicine 2020;73
    
    

43. Mason, C.H.; Perreault, W.D. Collinearity, Power, and Interpretation of Multiple-Regression Analysis. J Marketing Res 1991;28:268–280
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/3172863&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991FW93700002&link_type=ISI) 

44. Matérn, B. Spatial variation ed^eds: Springer Science & Business Media; 2013
    
    

45. Meade, M.S. Medical geography. The Wiley Blackwell Encyclopedia of Health, Illness, Behavior, and Society. Hoboken, NJ: John Wiley & Sons, Inc., 2014
    
    

46. Melnick, E.R.; Ioannidis, J.P. Should governments continue lockdown to slow the spread of covid-19? BMJ 2020;369:m1924
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNjkvanVuMDNfNS9tMTkyNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzMxLzIwMjAuMDcuMzEuMjAxNjU4MTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

47. Mizumoto, K.; Kagaya, K.; Zarebski, A.; Chowell, G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill 2020;25
    
    

48. Nishiura, H.; Kobayashi, T.; Miyama, T.; Suzuki, A.; Jung, S.M.; Hayashi, K.; Kinoshita, R.; Yang, Y.; Yuan, B.; Akhmetzhanov, A.R.; Linton, N.M. Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19). Int J Infect Dis 2020;94:154–155
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.03.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

49. Ogen, Y. Assessing nitrogen dioxide (NO2) levels as a contributing factor to coronavirus (COVID-19) fatality. Science of the Total Environment 2020;726
    
    

50. Onder, G.; Rezza, G.; Brusaferro, S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. JAMA 2020;
    
    

51. Our World in Data. Our World in Data; 2020a
    
    

52. Our World in Data. Coronavirus (COVID-19) Deaths. 2020b
    
    

53. Palk, L.; Blower, S. Geographic variation in sexual behavior can explain geospatial heterogeneity in the severity of the HIV epidemic in Malawi. Bmc Med 2018;16
    
    

54. Pampel, F.C. Global Patterns and Determinants of Sex Differences in Smoking. Int J Comp Sociol 2006;47:466–487
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0020715206070267&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21874066&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

55. Pfeiffer, D.; Robinson, T.P.; Stevenson, M.; Stevens, K.B.; Rogers, D.J.; Clements, A.C. Spatial analysis in epidemiology ed^eds: Oxford University Press Oxford; 2008
    
    

56. 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 2020;20:776–777
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

57. Ren, T.X.; Long, Z.H.; Zhang, R.G.; Chen, Q.Q. Moran’s I test of spatial panel data model - Based on bootstrap method. Econ Model 2014;41:9–14
    
    

58. Richardson, S.; Hirsch, J.S.; Narasimhan, M.; Crawford, J.M.; McGinn, T.; Davidson, K.W.; and the Northwell, C.-R.C., Barnaby, D.P.; Becker, L.B.; Chelico, J.D.; Cohen, S.L.; Cookingham, J.; Coppa, K.; Diefenbach, M.A.; Dominello, A.J.; Duer-Hefele, J.; Falzon, L.; Gitlin, J.; Hajizadeh, N.; Harvin, T.G.; Hirschwerk, D.A.; Kim, E.J.; Kozel, Z.M.; Marrast, L.M.; Mogavero, J.N.; Osorio, G.A.; Qiu, M.; Zanos, T.P. Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area. JAMA 2020;
    
    

59. Roser, M.; Ritchie, H.; Ortiz-Ospina, E.; Hasell, J. Coronavirus disease (COVID-19)–Statistics and research. Our World in data 2020;
    
    

60. Rousset, F.; Ferdy, J.B. Testing environmental and genetic effects in the presence of spatial autocorrelation. Ecography 2014;37:781–790
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ecog.00566&link_type=DOI) 

61. Royston, P. Multiple imputation of missing values. The Stata Journal 2004;4:227–241
    
    

62. Russell, T.W.; Hellewell, J.; Jarvis, C.I.; van Zandvoort, K.; Abbott, S.; Ratnayake, R.; Cmmid Covid-Working, G.; Flasche, S.; Eggo, R.M.; Edmunds, W.J.; Kucharski, A.J. Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship, February 2020. Euro Surveill 2020;25
    
    

63. Sasaki, S.; Suzuki, H.; Igarashi, K.; Tambatamba, B.; Mulenga, P. Spatial analysis of risk factor of cholera outbreak for 2003-2004 in a peri-urban area of Lusaka, Zambia. Am J Trop Med Hyg 2008;79:414–421
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI3OS8zLzQxNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzMxLzIwMjAuMDcuMzEuMjAxNjU4MTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

64. Schellekens, P.; Sourrouille, D.M. COVID-19 Mortality in Rich and Poor Countries: A Tale of Two Pandemics? World Bank Policy Research Working Paper 2020;[https://ssrn.com/abstract=3614141](https://ssrn.com/abstract=3614141)
    
    

65. Sjodin, H.; Wilder-Smith, A.; Osman, S.; Farooq, Z.; Rocklov, J. Only strict quarantine measures can curb the coronavirus disease (COVID-19) outbreak in Italy, 2020. Eurosurveillance 2020;25:7–12
    
    

66. Stevens, G.A.; Alkema, L.; Black, R.E.; Boerma, J.T.; Collins, G.S.; Ezzati, M.; Grove, J.T.; Hogan, D.R.; Hogan, M.C.; Horton, R.; Lawn, J.E.; Marusic, A.; Mathers, C.D.; Murray, C.J.; Rudan, I.; Salomon, J.A.; Simpson, P.J.; Vos, T.; Welch, V.; Group, G.W. Guidelines for Accurate and Transparent Health Estimates Reporting: the GATHER statement. PLoS Med 2016;13:e1002056
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27351744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

67. Team, C.C.-R. Severe Outcomes Among Patients with Coronavirus Disease 2019 (COVID-19) - United States, February 12-March 16, 2020. Mmwr-Morbid Mortal W 2020;69:343–346
    
    

68. United Nations. World Population Prospects 2019: Highlights. in: United Nations D.o.E.a.S.A., ed. New York city: United Nations; 2019
    
    

69. van Rossum, G. Python tutorial, Technical Report CS-R9526, Centrum voor Wiskunde en Informatica (CWI), Amsterdam.”. 1995;
    
    

70. Ward, M.D.; Gleditsch, K.S. Spatial regression models ed^eds: Sage Publications; 2018
    
    

71. Vardavas, C.I.; Nikitara, K. COVID-19 and smoking: A systematic review of the evidence. Tob Induc Dis 2020;18:20
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18332/tid/119324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

72. World Bank. World Bank Country and Lending Groups. World Bank; 2020
    
    

73. World Health Organization. COVID-19 strategy update. Geneva, Switzerland; 2020a
    
    

74. World Health Organization. Q&A: Similarities and differences – COVID-19 and influenza. 2020b
    
    

75. World Intellectual Property Organization. WIPO’s COVID-19 Response. in: Organization W.I.P., ed. New York: World Intellectual Property Organization; 2020
    
    

76. Worldometer. COVID-19 Coronavirus Pandemic. Worldometer; 2020
    
    

77. Wu, X.; Nethery, R.C.; Sabath, B.M.; Braun, D.; Dominici, F. Exposure to air pollution and COVID-19 mortality in the United States: A nationwide cross-sectional study. medRxiv 2020;
    
    

78. Yang, J.; Zheng, Y.; Gou, X.; Pu, K.; Chen, Z.; Guo, Q.; Ji, R.; Wang, H.; Wang, Y.; Zhou, Y. Prevalence of comorbidities and its effects in patients infected with SARS-CoV-2: a systematic review and meta-analysis. Int J Infect Dis 2020;94:91–95
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.03.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

79. Zaim, S.; Chong, J.H.; Sankaranarayanan, V.; Harky, A. COVID-19 and Multiorgan Response. Curr Probl Cardiol 2020;45:100618
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cpcardiol.2020.100618&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

80. Zhang, F.; Xiong, Y.; Wei, Y.; Hu, Y.; Wang, F.; Li, G.; Liu, K.; Du, R.; Wang, C.Y.; Zhu, W. Obesity predisposes to the risk of higher mortality in young COVID-19 patients. J Med Virol 2020;
    
    

81. Zheng, S.; Fan, J.; Yu, F.; Feng, B.; Lou, B.; Zou, Q.; Xie, G.; Lin, S.; Wang, R.; Yang, X.; Chen, W.; Wang, Q.; Zhang, D.; Liu, Y.; Gong, R.; Ma, Z.; Lu, S.; Xiao, Y.; Gu, Y.; Zhang, J.; Yao, H.; Xu, K.; Lu, X.; Wei, G.; Zhou, J.; Fang, Q.; Cai, H.; Qiu, Y.; Sheng, J.; Chen, Y.; Liang, T. Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January-March 2020: retrospective cohort study. BMJ 2020;369:m1443
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvYXByMjFfMjAvbTE0NDMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8zMS8yMDIwLjA3LjMxLjIwMTY1ODExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

82. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; Guan, L.; Wei, Y.; Li, H.; Wu, X.; Xu, J.; Tu, S.; Zhang, Y.; Chen, H.; Cao, B. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet 2020;395:1054–1062
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30566-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F31%2F2020.07.31.20165811.atom) 

83. Zulu, L.C.; Kalipeni, E.; Johannes, E. Analyzing spatial clustering and the spatiotemporal nature and trends of HIV/AIDS prevalence using GIS: the case of Malawi, 1994-2010. Bmc Infect Dis 2014;14