Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Face mask wearing rate predicts country’s COVID-19 death rates: with supplementary state-by-state data in the United States

View ORCID ProfileDaisuke Miyazawa, View ORCID ProfileGen Kaneko
doi: https://doi.org/10.1101/2020.06.22.20137745
Daisuke Miyazawa
1Miyazawa Clinic, 1-6-5 Akuraminami Takarazuka, Hyogo 665-0823, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daisuke Miyazawa
  • For correspondence: kusami1@ybb.ne.jp
Gen Kaneko
2School of Arts & Sciences, University of Houston-Victoria. 3007 Ben Wilson, Victoria, TX 77901, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Gen Kaneko
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Identification of biomedical and socioeconomic predictors for the number of deaths by COVID-19 among countries will lead to the development of effective intervention. While previous multiple regression studies have identified several predictors for the number of COVID-19-related deaths, little is known for the association with mask non-wearing rate possibly because the data is available for limited number of countries, which constricts the application of traditional multiple regression approach to screen a large number of potential predictors. In this study, we used the hypothesis-driven regression approach to test the association with limited number of predictors based on the hypothesis that the mask non-wearing rate can predict the number of deaths to a large extent together with age and BMI, other relatively independent risk factors for hospitalized patients of COVID-19. The mask non-wearing rate, percentage of age ≥ 80 (male), and male BMI showed Spearman’s correlations up to about 0.8, 0.7, and 0.6, respectively, with the numbers of deaths per million in 22 countries from mid-March to mid-June, 2020. The observed numbers of deaths per million were significantly correlated with those predicted by the lasso regression model including four predictors, age ≥ 80 (male), male BMI, and mask non-wearing rates from mid-March and late April to early May (Pearson’s coefficient = 0.919). The multiple linear regression models including the mask non-wearing rates, age, and obesity-related predictors explained up to 75% variation of the number of deaths per million in the 22 countries with little concerns about multicollinearity. Furthermore, linear regressions using the mask non-wearing rate in mid-March as the sole predictor still explained up to 72% of the variation of the numbers of deaths from March to mid-June, emphasizing the importance of the strongest predictor. Although further verification is needed to identify causes of the national differences in COVID-19 mortality rates, these findings highlight the importance of the mask, age, and BMI in predicting the COVID-19-related deaths, providing a useful strategy for future regression analyses that attempt to contribute to the mechanistic understanding of COVID-19.

Introduction

There have been considerable differences in the numbers of deaths by COVID-19 among countries. Western countries (Europe and the United States), in particular, have been suffered from the high number of deaths and the high death rate compared to East Asian countries. Several regression studies have identified predictors, or candidates for true risk factors, for the difference among countries such as age, obesity, and previous BCG vaccination (Bellali, Chtioui, & Chahed, 2020; Leung, Bulterys, & Bulterys, 2020; Squalli, 2020; Knittel & Ozaltun, 2020). However, it is important to continue characterizing novel and known predictors under various biomedical and socioeconomic contexts for potential future interventions to counteract the ongoing COVID-19 pandemic.

Face masks have been the subject of increasing attention throughout the pandemic period. Leffler et al. (2020) used issued dates for each government’s mandates or recommendations for face mask wearing as a predictor in their regression analyses, demonstrating that the duration of infection in the country, duration of wearing masks, and percentage of the population over age 60 significantly predicted each country’s COVID-19 mortality rate. Also, a systematic review and meta-analysis identifying 172 observational studies indicated that face mask use could result in a large reduction in risk of COVID-19 (Chu et al., 2020). These findings, along with the big difference in face mask wearing rates between Western and other East Asian countries, suggest that the importance of mask wearing rate as a predictor for the number of deaths by COVID-19. However, to our knowledge, there has been no analytical studies including the face mask wearing rate in regression models. This is possibly because the common regression approach uses a large number of samples and predictors to achieve the high prediction rates. The mask wearing rate is available only for limited number of countries.

Other studies identified age, BMI, and male sex as relatively independent risk factors for hospitalized patients of COVID-19. For example, in a cohort of hospitalized patients in a minority-predominant population, severe obesity, increasing age, and male sex were associated with higher in-hospital mortality (Palaiodimos et al., 2020). The critical illness after admission was strongly associated with age, heart failure, BMI > 40, and male sex in 5279 patients with laboratory confirmed COVID-19 patients in New York City (Petrilli et al., 2020). Docherty et al. (2020) reported using data from 20133 patients admitted to hospital with COVID-19 in UK that independent risk factors for mortality are increasing age, male sex, and chronic comorbidity including obesity. Furthermore, several mechanistic implications have been explored for these factors (Miyazawa, 2020; Mueller, McNamara, & Sinclair, 2020; Grandi, Facchinetti, & Bitzer, 2020), which has supported their independent nature.

In this study, we employed the hypothesis-driven regression approach to identify predictors for the number of deaths by COVID-19 across countries based on the hypothesis that the mask non-wearing rate can predict the number of deaths to a large extent together with age, BMI, and sex, other relatively independent risk factors for hospitalized patients of COVID-19. This approach uses limited number of a priori determined predictors based on a certain hypothesis. The hypothesis-driven regression is especially useful in providing biological insights related to the association with small number of samples, while the implication is restricted to the hypothesis, and has been used in various fields of research (Klein et al., 2013; Sabuncu et al., 2016; Carriedo et al., 2020). Based on the results we discussed the cultural factors for the observed difference among countries, which may have lowered the number of COVID-19-related deaths in East Asian countries.

Methods

Data collection World

All data was collected from publicly available secondary sources. Analyses in the present study included 13 Western (UK, France, Italy, USA, Spain, Mexico, Germany, Canada, Sweden, Norway, Finland, Denmark, and Australia) and 9 Asian (Malaysia, China, Saudi Arabia, India, Indonesia, Philippines, Japan, Singapore, and Thailand) countries. These countries were chosen because of the availability for mid-March mask-wearing data. The face mask wearing rates in mid-March (3/9-18) 2020 and late April to early May (4/26-5/1) 2020 across countries were derived from “% of people in each country who answered they are: Wearing a face mask when in public spaces” in Smith (2020). This database (YouGov) has partnered with the Institute of Global Health Innovation at Imperial College London and summarizes interviews with nationally representative sample sizes (150– 2000/week depending on countries). Total COVID deaths per million were obtained from Ritchie (2020). BMI data was obtained from Global Status Report on Noncommunicable Diseases 2014. 5 Oct. 2015, www.who.int/nmh/publications/ncd-status-report-2014/en/. Population percent by age data was obtained from https://data.worldbank.org on June 15, 2020.

United States

The face mask wearing rates in late March to late April (3/26-4/29) 2020 across states in the United States were derived from “State-by-state: Face mask adoption across the US; Which, if any, of the following measures have you taken in the past 2 weeks to protect yourself from the Coronavirus (COVID-19)? (% of US adults in each state who say they wear a face mask when in public)” in Nguyen (2020). Total COVID deaths per million by state in the United States on July 3 were obtained from CDC COVID Data Tracker July 3, 2020, https://www.cdc.gov/covid-data-tracker/index.html. (Accessed July 3, 2020).

Data analysis

R version 3.6.2 was used for all statistical analyses. R packages used in this study include fields, ggplot2, glmnet, tidyverse, broom, car, and ggcorrplot. Pearson’s and Spearman’s correlations were calculated by cor and cor.test functions. Lasso regression was used to find potential association of predictors with the number of deaths per million (Tibshirani, 2011). Multiple and single linear regression models were built using the lm function.

Results

Spearman’s correlations

We first calculated Spearman’s correlations between the number of deaths per million from various dates and predictors related to the mask non-wearing rate, age, and obesity (Fig. 1). The correlation coefficients between the number of deaths per million and mask non-wearing rates (both mid-March and late April to early May) showed time-dependent increases, reaching the apparent plateau in early May. The mask non-wearing rate in mid-March generally showed higher correlation coefficients than those that in late April to early May. The highest correlation of 0.79 was observed between the mid-March mask non-wearing rate and the number of deaths per million on June 6, 2020 (P = 6.865e-06).

Fig. 1.
  • Download figure
  • Open in new tab
Fig. 1.

Spearman’s correlation (over time) among parameters related to the number of COVID19-related deaths per million, obesity, age, and mask non-wearing rate. Age-related predictors indicate the percentage of the population of the specified range.

Age-related predictors were also highly correlated with the number of deaths per million with Spearman’s correlation coefficients of up to ∼0.68. In both sexes, correlation coefficients were lower in younger age groups (age 65–69) compared to those of higher age groups, and the highest correlations were found in age ≥ 80 groups. In contrast to the mask non-wearing rate, correlations between the number of deaths and age-related predictors were highest in mid-May and decreased thereafter.

In general BMI showed lower correlations with the number of deaths per million with a clear sex-dependent difference. Male BMI was more closely associated with the number of deaths per million than female BMI. The highest correlation of 0.59 was observed between male BMI and the number of deaths per million on May 26.

Time course profiles for Spearman’s correlation coefficients of several highly correlated parameters demonstrated the unique feature of age ≥ 80 (male) among parameters used (Fig. 2). Spearman’s correlation coefficients for this parameter remained almost constant throughout the study period, whereas those for other three parameters showed time-dependent increase and reached the apparent plateau around day 50 (mid-May).

Fig. 2.
  • Download figure
  • Open in new tab
Fig. 2.

Time-dependent changes in Spearman’s correlations between the number of deaths per million on May 13 and mask non-wearing rates (mid-March and late April to early May), age ≥ 80, and male BMI.

Scatter plots

We subsequently made scatter plots showing the relationship between the number of deaths per million and predictors that showed high Spearman’s correlations (Fig. 3). The number of deaths per million on May 13 was used because it showed the highest correlations with most important predictors including the mask non-wearing rates (mid-March), age ≥ 80 (male), and male BMI. The number of deaths per million showed exponential type increases in most plots, in the plot with the mask non-wearing rate (mid-March) in particular, with clear separations between Western and Asian countries. As expected, the number of deaths per million did not show a clear positive correlation with the mask non-wearing rate (late April to early May). Rather, in the Western countries there was a tendency toward a negative association between the number of deaths per million and mask non-wearing rate in April - May (Spearman’s correlation; r = −0.4126551, P = 0.1611). This is attributed to the marked reduction in the mask non-wearing rate in some countries with high number of deaths, possibly caused by the fear of the disease. Likewise, within the United States, the mask non-wearing rate in late March to late April across states also showed a negative correlation with the total number of deaths per million on July 3 (Spearman’s correlation; r = −0.477556, P = 0.0003946) (Fig. S1).

Supplementary Fig. 1:
  • Download figure
  • Open in new tab
Supplementary Fig. 1:

Scatter plot showing cross-state correlation between ln(total number of deaths per million on July 3) and mask non-wearing rate (late March to late April) by state in the United States. Spearman’s correlation; r = −0.477556, P = 0.0003946.

Fig. 3.
  • Download figure
  • Open in new tab
Fig. 3.

Scatter plots showing the correlations among parameters related to COVID19 deaths, obesity, age, and mask non-wearing rate. The number of deaths per million on May 13 was used with (bottom panels) and without (top panels) logarithmic transformation.

Based on the cross-country scatter plots, we applied the logarithmic transformation to the number of deaths per million. The transformed number of deaths per million on May 13 showed significant linear correlation with the mask non-wearing rate in mid-March (r = 0.842, P = 8.754e-07), mask non-wearing rate in late April to early May (r = 0.496, P = 0.01897), age ≥ 80 (male) (r = 0.658, P = 0.0008717), and BMI (male) (r = 0.682, P = 0.0004762). The relationship between the number of deaths per million (May 13) and age ≥ 80 (female) appeared to be similar in the scatter plot (data not shown), but their linear correlation did not reach the statistical significance with a lower Pearson correlation coefficient compared to that of age ≥ 80 of males (r = 0.416, P = 0.05381).

Lasso regression

Next, we used the least absolute shrinkage and selection operator (lasso) regression to find the potential association of these predictors with the log-transformed number of deaths per million. The lasso regression is a machine learning method that estimates the regression coefficients with the lasso penalty terms, retaining important explanatory variables, while less important variables are removed from the regression (Tibshirani, 2011). When all parameters in Fig. 2 were included together with age ≥ 80 (female), the lasso regression selected BMI (male) (beta = 0.21830178), age ≥ 80 (male) (0.35009932), the mask non-wearing rate (mid-March) (0.06088824), and the mask non-wearing rate in late April to early May (−0.02080248) with an intercept of −6.73379877. Age ≥ 80 (female) was omitted from the regression. The observed number of deaths per million showed high correlation with the number predicted by the lasso regression model (Fig. 4; Pearson’s coefficient = 0.919, P = 1.601e-09).

Fig. 4.
  • Download figure
  • Open in new tab
Fig. 4.

Pearson’s correlation between observed and predicted values of log-transformed number of deaths per million. The minimum lambda of 0.01995262 was used for the lasso regression. Pearson’s coefficient = 0.919, P = 1.601e-09.

Multiple linear regression

We subsequently attempted to predict the number of deaths per million by the traditional multiple linear regression approach. In contrast to the machine learning-based lasso regression that is designed primarily for prediction, the multiple linear regression provides a statistical interference, with which the role of each predictor can be discussed in detail. The multiple linear regression, however, has two major constraints: 1) The number of predictors should be conservative, not exceeding ∼10% of the sample size (Steyerberg & Harrell, 2004); 2) The multiple linear regressions is vulnerable to multicollinearity, a severe distortion of the model caused by high correlations among predictors. The definition of high correlation depends on the field and model, but in general |r| > 0.7 (Dormann et al., 2013) or |r| > 0.8–0.9 (Franke, 2010) are considered to be highly correlated. In order to deal with the first constraint, we employed the hypothesis-driven regression using two or three predictors based on the priori determined hypothesis. Predictors selected by the lasso regression were used. For the second constraint, we calculated the correlation between the predictors prior to conducting the multiple regression.

The mask non-wearing rate in mid-March showed significant correlations with that in late April to early May, age ≥ 80 (male), and male BMI with Spearman’s correlation coefficients of 0.87, 0.52, and 0.54, respectively. Compared to the mask non-wearing rate in mid-March, that in late April to early May showed a weak association with male BMI, and its correlation with age ≥ 80 (male) was not statistically significant. Age ≥ 80 (male) and male BMI were not significantly correlated. Therefore, we concluded that our regression model may have a problem of multicollinearity when it contains the two mask non-wearing rates, but other combinations of the predictors would not severely affect the accuracy of the regression. For further verification, we calculated the variance inflation factor (VIF) in all regression models.

We first selected BMI (male), age ≥ 80 (male), and the mask non-wearing rate in mid-March for this regression (Table 1, model 1). This multiple regression model explained 76.0% variation of the number of deaths per million on May 13 (P = 2.075e-06), identifying the significant association with mask non-wearing rate in mid-March and age ≥ 80 (male). BMI (male) also showed a trend for association, but the association was not statistically significant. This regression model may have a weak multicollinearity problem because the VIF of the mask non-wearing rate exceeded 2.5, a strict threshold of VIF in evaluating the presence of multicollinearity.

View this table:
  • View inline
  • View popup
Table 1.

Multiple linear regressions for predicting the number of deaths per million (May 13) using three or two predictors

Other combinations of three predictors also resulted in significant associations between the number of deaths per million on May 13 and the selected predictors with ∼0.7 adjusted R square values (Table 1). We found several important features in these multiple regressions. First, the association with mask non-wearing rate in mid-March was significant in all regressions tested with low P values, and exclusion of this predictor markedly decreased the adjusted R square (model 2). Second, BMI (male) showed a relatively low association with the number of deaths on May 13, being significant only when the strongest predictor, mask non-wearing rate in mid-Match, was not included the regression (model 3). Third, the highest adjusted R square of 0.7934 was obtained when we include predictors age ≥ 80 (male) and mask non-wearing rates of mid-March and late April to early May (model 4). All regression coefficients were statistically significant in this model. Forth, VIF values were generally low unless the two mask non-wearing rates were together included in the model.

Regressions including two predictors also underscored the importance of the mask non-wearing rate in mid-March as a predictor for the number of deaths per million on May 13. Namely, the mask non-wearing rate in mid-March was always a significant predictor (models 5–7) with adjusted R square of up to ∼0.75. While BMI (male) and age ≥ 80 (male) were also significantly associated with the number of deaths per million on May 13, exclusion of the mask non-wearing rate in mid-March dropped the adjusted R square to ∼0.65. It is also noted that the regression coefficient of the mask non-wearing late from April to early May became negative when this predictor was included in the model with that of mid-March (models 4 and 7). This is probably related to the weak association of this predictor (Fig. 3) and its correlation with the mask non-wearing rate from April to May (Fig. 5; see Discussion for details). Altogether, we concluded that the model 5 would be the best regression model that predict about 75% of the variation of the number of deaths per million by satisfying all constrains in the linear regression.

Fig. 5.
  • Download figure
  • Open in new tab
Fig. 5.

Spearman’s correlation matrix between mask non-wearing rates, age ≥ 80 (male), and male BMI. Only significant correlations are shown (P < 0.05).

Single regression

Lastly, we built several single linear regression models for predicting the number of deaths from other dates using this predictor (Table 2) since the mask non-wearing rate in mid-March showed the strong association with the number of deaths on May 13. Overall, changes in adjusted R square showed the similar tendency with those of the Spearman’s correlation coefficients (Fig. 1), in which the value gradually increased during March and April, reaching the plateau in mid-May. Overall, the mask non-wearing rate in mid-March explained up to 72% of variations of the number of deaths per million.

View this table:
  • View inline
  • View popup
Table 2.

Linear regressions for predicting the log-transformed number of deaths per million from various dates using the mask non-wearing date in mid-March as the single predictor

Discussion

The number of deaths per million on a certain date highly depends on when the epidemic began in the country, in particular in the early stages of the pandemic. In this study, we have tentatively determined March 20 as the start date of a global pandemic. It is likely that this assumption is a reasonable approximation for our hypothesis because all employed predictors (mask non-wearing rates, age and BMI) tended to show weak correlations with the number of deaths per million at the onset of the pandemic followed by marked increases as the progress of the pandemic. These features enabled us to select a specific date (May 13), at which most parameters showed the high correlations with the number of deaths per million, and to build regression models that well explain the variation of the response variable with limited number of predictors. On the other hand, it is noteworthy that the age-related predictors were highly correlated to the number of deaths per million from the very beginning of the pandemic and then remained constant or gradually declined. Below we discuss the potential biomedical and socioeconomic mechanisms by which each predictor affects the number of deaths by COVID-19.

Obesity

We have proposed possible mechanisms by which obesity could be a risk factor for the severity of COVID-19 and H1N1 influenza infections (Miyazawa, 2020). Namely, respiratory failure is the most important pathology that contributes to the severity of both COVID-19 and H1N1 influenza infections. Since obese patients generally show a restrictive breathing pattern and reduced lung volumes, the obesity-hypoventilation syndrome can lead to respiratory failure in COVID-19 patients, being a risk factor especially for patients with severe symptoms. Also, obesity has been reported to be a risk factor for the development of acute respiratory distress syndrome (Zhi et al, 2016), which is a serious clinical manifestation of COVID-19 (Simonnet et al., 2020). According to Moriconi et al. (2020), in patients of COVID-19, inflammatory markers are higher in obese group than non-obese group at admission, and obese group show a worse pulmonary clinical picture with lower PaO2. Thus, obesity could be a risk factor for already hospitalized COVID-19 patients or even for those with serious symptoms.

These potential mechanisms may explain the reason why the association of BMI (male) and the number of deaths increased May to June 2020. Typical COVID-19 patients are hospitalized a few days to a week after infection and in some cases the symptom become severe a few weeks after hospitalization. If the obesity affects the last phase of the transition of the disease, there must be some delay in the increase in correlations between BMI and the number of deaths per million compared to the spread of COVID-19. However, we cannot rule out the possibility that obese people have high chance to be infected because obese adults are known to inhale air average 50% more per day than non-obese adults (Brochu, 2014), which may lead to increase the chance to inhale the virus.

Age

In contrast to BMI, age ≥ 80 (male) was correlated with the number of deaths from the onset of the global pandemic. This may reflect the common etiological feature that old people are more susceptible to infection and have higher risk of death after infection compared to young people because of the dysfunction of immunity.

The precise reason why older age contributes to the number of COVID-19-related deaths remains to be elucidated, while the immune dysfunction has been proposed as a potential mechanism (Mueller, McNamara, & Sinclair, 2020). Meanwhile, it is controversial whether the reduced immunity simply contributes to the higher mortality after infection. Excessive immunity could cause the cytokine storm and ARDS in COVID-19 in some cases (Ye, Wang, & Mao, 2020). Furthermore, two patients with X-linked agammaglobulinemia have recovered from COVID-19, suggesting that the B-cell response may not be necessary for the recovery from this disease (Soresina et al., 2020). While other dysfunction of immunity, such as age-related T-cell dysfunction, may be involved in the severity of COVID-19 (Minato, Hattori, & Hamazaki, 2020), the mechanism should be carefully investigated.

Along with immune dysfunction, old individuals tend to provoke unwanted inflammation, which may contribute to the severity of COVID-19 (Mueller, McNamara, & Sinclair, 2020). This tendency has been observed in SARS-CoV-infected old nonhuman primate (Smits et al., 2010). Also, To et al. (2020) reported that the older age is correlated with the higher salivary viral load, which becomes the highest during the first week after symptom onset. Furthermore, as in the case of obese patients, decreased respiratory function of elderly could also be the reason. In elder individuals, both forced expiratory volume in one second (FEV1) and forced vital capacity (FVC) decrease dramatically (Falaschetti et al., 2004) to almost half of their lifetime maximum values, especially in males (Leem et al., 2019).

Sex differences

We identified a clear sex-dependent difference in the Spearman’s correlation coefficients between the number of deaths per million and obesity-related parameters. This is in line with the fact that males are more susceptible to severe outcomes from COVID-19 than females (Montopoli et al., 2020; Grandi, Facchinetti, & Bitzer, 2020). Likewise, it has been proposed that estrogen may be protective against COVID-19 in females (Grandi, Facchinetti, & Bitzer, 2020), while androgens may worsen COVID-19 in males (Montopoli et al., 2020). It is also possible that smoking, which is more common in males globally, may play a role in worsening COVID-19 (Grandi, Facchinetti, & Bitzer, 2020). Also, females generally have more fear to COVID-19 and tend to take more preventive behaviors than males (Yildirim, 2020)

Face mask

Face masks are considered to be effective in preventing transmission of COVID-19 to others, and also have a protective effect in preventing transmission from others. A recent simulation study showed that universal masking at 80% adoption suppresses COVID-19 deaths significantly more than maintaining a lockdown. Comparison between the validated modeling results and empirical data from Asian regions showed an almost perfect correlation between early universal masking and successful suppression of COVID-19 outbreaks (Kai et al., 2020). Also, Zhang et al. (2020) reported that face mask alone could have reduced the number of infections by over 78,000 in Italy from April 6 to May 9 and over 66,000 in New York City from April 17 to May 9, claiming that the failure in containing the propagation of COVID-19 pandemic worldwide is largely attributed to the unrecognized importance of airborne virus transmission.

In line with these previous studies, the mask non-wearing rate in mid-March was found to be the strongest predictor for the number of deaths per million in this study, where even the single regression explained more than 70% of the variation of the response variable. We also visualized that there is a big difference in face mask wearing rates between Western countries and Asian countries. The county’s policy for wearing face mask alone cannot explain this big difference because, for example, face mask has never been mandated in Japan, despite its high face mask wearing rate. We speculate that the cultural may be the major reason for the difference. Jack, Caldara, & Schyns (2012) mentioned that “whereas Western Caucasian internal representations predominantly featured the eyebrows and mouth, East Asian internal representations showed a preference for expressive information in the eye region.” This tendency may be the major reason why it is more considered rude wearing sunglasses among East Asia (Gesteland, 2020), whereas it is considered more suspicious to wear face masks in Western countries.

Correlation between variables

There was a significant positive correlation between male BMI and mask non-wearing rate in mid-March, although the degree of correlation was far below to introduce severe multicollinearity in our regression. While this does not simply indicate that obese people tend not to wear masks, it is possible to speculate that when people become more obese, they feel more uncomfortable to wear masks since obese adult inhale air average 50% more per day than non-obese adults (Brochu, 2014). Also, body temperature is positively associated with obesity (Bastardot, 2019), and face masks could rise body temperature (Yip et al., 2005; Hayashi & Tokura, 2004). Therefore, obese people may feel more heat and discomfort when wearing face masks depending on the temperature and humidity (Li et al., 2005). The small size of universal face masks may be a simple reason for the correlation.

There was a significant positive correlation between age ≥ 80 (male) and mask non-wearing rate in mid-March. Due to the lack of data on mask wearing rates by age, this cannot be fully discussed either. It has been reported that age is correlated with fear and preventive behaviors (Yildirim, 2020). Face mask wearing in public correlates with the proportion of people who are afraid of COVID-19, except for Scandinavian countries and United Kingdom (Smith, 2020). These facts are inconsistent with the correlation found in this study. Further data is required for detailed investigation.

Conclusion

The mask non-wearing rate in mid-March alone explained up to 72% of variations of the number of deaths per million, whereas together with age ≥ 80 (male) it explained 75% of the variation. The fact that the face mask non-wearing rate in mid-March, the early phase of the pandemic, more strongly predicted the total number of deaths per population over the entire three-month course of the pandemic than the face mask non-wearing rate in late April to early May may suggest that face mask wearing from the early phase of the pandemic may be very important for suppressing COVID-19-related deaths, if mask wearing is assumed to be an independent risk factor.

In the lasso regression, four factors, BMI (male), age ≥ 80 (male), the mask non-wearing rate (mid-March), and the mask non-wearing rate in late April to early May were chosen, and the observed and predicted model of number of deaths per million showed high correlation (Pearson’s coefficient = 0.919, P = 1.635e-09). Although this has not been proven at this time to be the cause of the differences in COVID-19 mortality by country, it may be a useful method for predicting COVID-19 mortality.

The death rate of the country in the early stages of a pandemic could be further predicted if factors such as the timing and number of infected people entering the country from abroad are taken into account.

Data Availability

All data was collected from publicly available secondary sources. % of people in each country who say they are: Wearing a face mask when in an public spaces as face mask wearing rate in mid-March(3/9-18) 2020 and late April to early May(4/26-5/1) 2020 of countries listed below were obtained from: Smith M. International COVID-19 Tracker Update: 2 May. 2 May 2020, yougov.co.uk/topics/international/articles-reports/2020/05/01/international-covid-19-tracker-update-2-may Smith M. International COVID-19 Tracker Update: 18 May. 18 May 2020, yougov.co.uk/topics/international/articles-reports/2020/05/18/international-covid-19-tracker-update-18-may Total COVID death per million of countries listed below except for Taiwan were obtained from: Ritchie, Hannah. Coronavirus Source Data. 16 June 2020, ourworldindata.org/coronavirus-source-data Data of BMIs of countries listed below except for Taiwan were obtained from: Global Status Report on Noncommunicable Diseases 2014. 5 Oct. 2015, www.who.int/nmh/publications/ncd-status-report-2014/en/. Population% by age of countries listed below except for Taiwan were obtained from: Kose, Ayhan, et al. World Bank Open Data. Data, 15 June 2020, data.worldbank.org/.

Acknowledgements

The authors express sincere thanks to a journalist Mr. Junya Iwai for his help in data collection and constructive discussion.

References

  1. Bastardot F, Marques-Vidal P, Vollenweider P. Association of body temperature with obesity. The CoLaus study. Int J Obes (Lond). 2019;43(5):1026–1033. doi:10.1038/s41366-018-0218-7
    OpenUrlCrossRef
  2. ↵
    Bellali H, Chtioui N, Chahed M. Factors associated with country-variation in COVID-19 morbidity and mortality worldwide: an observational geographic study. medRxiv 2020.05.27.20114280; doi: https://doi.org/10.1101/2020.05.27.20114280
  3. ↵
    Brochu P, Bouchard M, Haddad S. Physiological daily inhalation rates for health risk assessment in overweight/obese children, adults, and elderly. Risk Anal. 2014;34(3):567–582. doi:10.1111/risa.12125
    OpenUrlCrossRef
  4. Carriedo A, Cecchini JA, Gonzalez C. Soccer spectators’ moral functioning and aggressive tendencies in life and when watching soccer matches. Ethics & Behav, published online. doi:10.17605/OSF.IO/TPA6U.
    OpenUrlCrossRef
  5. ↵
    Chu DK, Akl EA, Duda S et al. Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis. Lancet 2020;395:1973–1987. doi:10.1016/S0140-6736(20)31142-9
    OpenUrlCrossRefPubMed
  6. ↵
    Docherty AB, Harrison EM, Green CA et al. Features of 20?133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ. 2020;369:m1985. doi:10.1136/bmj.m1985
    OpenUrlAbstract/FREE Full Text
  7. ↵
    Dormann CF, Elith J, Bacher S et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 2013:36;027–046.
    OpenUrlCrossRef
  8. ↵
    Falaschetti E, Laiho J, Primatesta P, Purdon S. Prediction equations for normal and low lung function from the Health Survey for England. Eur Respir J 2004;23(3):456–463. doi: 10.1183/09031936.04.00055204
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Sheth, J.,
    2. Malhotra, N.
    Franke GR. Multicollinearity. In Wiley International Encyclopedia of Marketing; Sheth, J., Malhotra, N., Eds.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2010.
  10. Gesteland R, Cross-Cultural Business Behavior: A Guide for Global Management (Fifth Edition), Copenhagen Business School Press, Norway, 2012.
  11. ↵
    Grandi G. Facchinetti F & Bitzer J (2020) The gendered impact of coronavirus disease (COVID-19): do estrogens play a role? Eur J Contracept Reprod Health Care. doi: 10.1080/13625187.2020.1766017
    OpenUrlCrossRef
  12. Hartman WR, Hess AS, Connor J. Persistent viral RNA shedding after COVID-19 symptom resolution in older convalescent plasma donors [published online ahead of print, 2020 Jun 13]. Transfusion. 2020;10.1111/trf.15927. doi:10.1111/trf.15927
    OpenUrlCrossRef
  13. ↵
    Hayashi C, Tokura H. The effects of two kinds of mask (with or without exhaust valve) on clothing microclimates inside the mask in participants wearing protective clothing for spraying pesticides. Int Arch Occup Environ Health. 2004;77(1):73–78. doi:10.1007/s00420-003-0472-3
    OpenUrlCrossRefPubMedWeb of Science
  14. ↵
    Jack RE, Caldara R, Schyns PG. Internal representations reveal cultural diversity in expectations of facial expressions of emotion. J Exp Psychol Gen. 2012;141(1):19–25. doi:10.1037/a0023463
    OpenUrlCrossRef
  15. Kai D, Goldstein G, Morgunov A et al. Universal masking is urgent in the COVID-19 pandemic: SEIR and agent based models, empirical validation, policy recommendations. arxiv:2004.13553 doi: 10.13140/RG.2.2.21662.08001
    OpenUrlCrossRef
  16. ↵
    Klein NJ, Beck AJ, Hutchins DA, Sanudo-Wilhelmy SA. Regression modeling of the North East Atlantic Spring Bloom suggests previously unrecognized biological roles for V and Mo. Front Microbiol 2020;4:45. doi: 10.3389/fmicb.2013.00045
    OpenUrlCrossRef
  17. ↵
    Knittel C, Ozaltun B. What does and does not correlate with COVID-19 death rates. medRxiv 2020.06.09.20126805; doi: https://doi.org/10.1101/2020.06.09.20126805
  18. ↵
    Leem, AY, Park B, Kim, YS et al. Longitudinal decline in lung function: a community-based cohort study in Korea. Sci Rep 2019;9:13614. doi: doi.org/10.1038/s41598-019-49598-9
    OpenUrlCrossRef
  19. Leffler CT, Ing EB, Lykins VJD et al. Association of country-wide coronavirus mortality with demographics, testing, lockdowns, and public wearing of masks. Update June 2, 2020. medRxiv 2020.05.22.20109231; doi: https://doi.org/10.1101/2020.05.22.20109231
  20. ↵
    Leung NY, Bulterys MA, Bulterys PL. Predictors of COVID-19 incidence, mortality, and epidemic growth rate at the country level. medRxiv 2020.05.15.20101097; doi: https://doi.org/10.1101/2020.05.15.20101097
  21. ↵
    Li Y, Tokura H, Guo YP et al. Effects of wearing N95 and surgical facemasks on heart rate, thermal stress and subjective sensations. Int Arch Occup Environ Health 2005;78(6):501-509. doi:10.1007/s00420-004-0584-4 May 16]. Metabolism. 2020;108:154262. doi:10.1016/j.metabol.2020.154262
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    Minato N, Hattori M, Hamazaki Y. Physiology and pathology of T-cell aging. Int Immunol. 2020;32(4):223–231. doi:10.1093/intimm/dxaa006
    OpenUrlCrossRef
  23. Miyazawa D. Why obesity, hypertension, diabetes, and ethnicities are common risk factors for COVID-19 and H1N1 influenza infections. J Med Virol, published online. doi: doi.org/10.1002/jmv.26220
    OpenUrlCrossRef
  24. ↵
    Montopoli M, Zumerle S, Vettor R, et al. Androgen-deprivation therapies for prostate cancer and risk of infection by SARS-CoV-2: a population-based study (n=4532). Ann Oncol, published online. doi:10.1016/j.annonc.2020.04.479
    OpenUrlCrossRef
  25. ↵
    Moriconi D, Masi S, Rebelos E et al. Obesity prolongs the hospital stay in patients affected by COVID-19, and may impact on SARS-COV-2 shedding. Obes Res Clin Pract. 2020;14:205–209. doi:10.1016/j.orcp.2020.05.009
    OpenUrlCrossRef
  26. ↵
    Mueller AL, McNamara MS, Sinclair DA. Why does COVID-19 disproportionately affect older people? Aging (Albany NY). 2020;12(10):9959–9981. doi:10.18632/aging.103344
    OpenUrlCrossRefPubMed
  27. Nguyen H. The states that are more and less likely to adopt face masks. YouGov. https://today.yougov.com/topics/politics/articles-reports/2020/05/08/states-are-more-and-less-likely-adopt-face-masks. published May 8, 2020. accessed July 3, 2020.
  28. ↵
    Palaiodimos L, Kokkinidis DG, Li W et al. Severe obesity, increasing age and male sex are independently associated with worse in-hospital outcomes, and higher in-hospital mortality, in a cohort of patients with COVID-19 in the Bronx, New York. Metabolism. 2020;108:154262. doi: 10.1016/j.metabol.2020.154262
    OpenUrlCrossRefPubMed
  29. ↵
    Petrilli CM, Jones SA, Yang J 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. 2020;369:m1966. doi:10.1136/bmj.m1966
    OpenUrlAbstract/FREE Full Text
  30. ↵
    Ritchie H. Coronavirus Source Data. 16 June 2020. https://ourworldindata.org/coronavirus-source-data
  31. ↵
    Sabuncu MR, Ge T, Holmes AJ et al. Morphometricity as a measure of the neuroanatomical signature of a trait. Proc Natl Acad Sci USA 2016:113;E5749–5756. doi: doi.org/10.1073/pnas.1604378113
    OpenUrlAbstract/FREE Full Text
  32. Simonnet A, Chetboun M, Poissy J et al. High prevalence of obesity in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) requiring invasive mechanical ventilation. Obesity, published online. doi:10.1002/oby.22831
    OpenUrlCrossRefPubMed
  33. ↵
    Smith M. International COVID-19 Tracker Update: 2 May. 2 May 2020, yougov.co.uk/topics/international/articles-reports/2020/05/01/international-covid-19-tracker-update-2-may
  34. ↵
    Smits SL, de Lang A, van den Brand JMA, Leijten LM, van IJcken WF, Eijkemans MJC, et al. (2010) Exacerbated innate host response to SARS-CoV in aged non-human primates. PLoS Pathog 2010;6(2):e1000756. doi: doi.org/10.1371/journal.ppat.1000756
    OpenUrlCrossRefPubMed
  35. Soresina A, Moratto D, Chiarini M et al. Two X-linked agammaglobulinemia patients develop pneumonia as COVID-19 manifestation but recover. Pediatr Allergy Immunol, published online. doi: 10.1111/pai.13263. doi:10.1111/pai.13263
    OpenUrlCrossRef
  36. ↵
    Squalli J. Evaluating the determinants of COVID-19 mortality: A cross-country study. medRxiv 2020.05.12.20099093v1; https://doi.org/10.1101/2020.05.12.20099093
  37. Steyerberg E & Harrel FE Jr.. Chapter 8: Statistical Models for Prediction. In Max MD & Lynn J. Symptom Research: Methods and Opportunities. National Institutes of Health, Bethesda, MD, USA, 2004.
  38. ↵
    Tibshirani R. Regression shrinkage and selection via the lasso: a retrospective. J R Stat Soc 2011;73(3):273–282. doi: doi.org/10.1111/j.1467-9868.2011.00771.x
    OpenUrlCrossRef
  39. ↵
    To KK, Tsang OT, Leung WS et al. Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study. Lancet Infect Dis 2020;20(5):565–574. doi:10.1016/S1473-3099(20)30196-1
    OpenUrlCrossRefPubMed
  40. ↵
    Ye Q, Wang B, Mao J. The pathogenesis and treatment of the ‘cytokine storm’ in COVID-19. J Infect 2020;80(6):607–613. doi: doi.org/10.1016/j.jinf.2020.03.037
    OpenUrlCrossRefPubMed
  41. ↵
    Yip WL, Leung LP, Lau PF, & Tong HK. The Effect of wearing a face mask on body temperature. Hong Kong J Emerg Med 2005;12:23–27.
    OpenUrl
  42. ↵
    Yildirim M, Geçer E, Akgül Ö. The impacts of vulnerability, perceived risk, and fear on preventive behaviours against COVID-19. Psychol Health Med. 2020;1–9. doi: 10.1080/13548506.2020.1776891
    OpenUrlCrossRef
  43. Zhang R, Li Y, Zhang AL, Wang Y, Molina MJ. Identifying airborne transmission as the dominant route for the spread of COVID-19. Proc Natl Acad Sci USA, published online. doi:10.1073/pnas.2009637117
    OpenUrlAbstract/FREE Full Text
  44. ↵
    Zhi G, Xin W, Ying W, Guohong X, Shuying L. “Obesity Paradox” in acute respiratory distress syndrome: asystematic review and meta-analysis. PLoS One 2016;11(9):e0163677. doi: doi.org/10.1371/journal.pone.0163677
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted July 14, 2020.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Face mask wearing rate predicts country’s COVID-19 death rates: with supplementary state-by-state data in the United States
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Face mask wearing rate predicts country’s COVID-19 death rates: with supplementary state-by-state data in the United States
Daisuke Miyazawa, Gen Kaneko
medRxiv 2020.06.22.20137745; doi: https://doi.org/10.1101/2020.06.22.20137745
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Face mask wearing rate predicts country’s COVID-19 death rates: with supplementary state-by-state data in the United States
Daisuke Miyazawa, Gen Kaneko
medRxiv 2020.06.22.20137745; doi: https://doi.org/10.1101/2020.06.22.20137745

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Infectious Diseases (except HIV/AIDS)
Subject Areas
All Articles
  • Addiction Medicine (215)
  • Allergy and Immunology (495)
  • Anesthesia (106)
  • Cardiovascular Medicine (1095)
  • Dentistry and Oral Medicine (196)
  • Dermatology (141)
  • Emergency Medicine (274)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (500)
  • Epidemiology (9765)
  • Forensic Medicine (5)
  • Gastroenterology (480)
  • Genetic and Genomic Medicine (2308)
  • Geriatric Medicine (222)
  • Health Economics (462)
  • Health Informatics (1557)
  • Health Policy (734)
  • Health Systems and Quality Improvement (603)
  • Hematology (236)
  • HIV/AIDS (503)
  • Infectious Diseases (except HIV/AIDS) (11638)
  • Intensive Care and Critical Care Medicine (617)
  • Medical Education (237)
  • Medical Ethics (67)
  • Nephrology (257)
  • Neurology (2141)
  • Nursing (134)
  • Nutrition (336)
  • Obstetrics and Gynecology (426)
  • Occupational and Environmental Health (517)
  • Oncology (1175)
  • Ophthalmology (364)
  • Orthopedics (128)
  • Otolaryngology (220)
  • Pain Medicine (146)
  • Palliative Medicine (50)
  • Pathology (311)
  • Pediatrics (695)
  • Pharmacology and Therapeutics (300)
  • Primary Care Research (266)
  • Psychiatry and Clinical Psychology (2180)
  • Public and Global Health (4653)
  • Radiology and Imaging (776)
  • Rehabilitation Medicine and Physical Therapy (457)
  • Respiratory Medicine (623)
  • Rheumatology (274)
  • Sexual and Reproductive Health (225)
  • Sports Medicine (210)
  • Surgery (251)
  • Toxicology (43)
  • Transplantation (120)
  • Urology (94)