Culture And COVID-19 Related Mortality- A Cross Sectional Study of 53 Countries

DOI: https://doi.org/10.21203/rs.3.rs-942150/v1

Abstract

The aim of this paper is to investigate the influence of culture on COVID-19 related mortality relying on Hofstede’s six dimensional concept of culture. Retrieving data from 53 countries a multivariable regression model was fitted that controls for health related, economic and policy related variables that have been found to be associated with mortality. In particular we included the percentage of population above 65 years of age, the prevalence of relevant co-morbidities (i.e. diabetes, hypertension, cancer. obesity), and tobacco usage as health related variables. Economic variables were GDP per capita and the connectedness of a country as measured by the number of international arrivals. As policy variables we included the Oxford stringency index as well as stringency speed, and the Global Health Security index. The results suggests that individualistic societies can navigate well through the pandemic when they are also characterized by low power distance which reflects a low acceptance of inequality within a society. This finding contradicts with previous studies that supported the popular narrative that autocratic societies with an obedient population are better positioned to manage the pandemic. 

1. Introduction

Due to the relevance and omnipresence of COVID-19, a considerable steam of literature exists that analyses factors influencing the spread of infection and resulting mortality. So far COVID -19 has claimed more than 25 million lives globally [1]. Most studies take clinical and socioeconomic variables into account to explain mortality while culture as an influencing factor has received little attention in this context. Among some notable exceptions is the study by Gokmen et al [2] who build upon Hofstede’s work on national culture. Hofstede's pioneering model consists of six dimensions that constitute cultural norms [3,4]: Power distance, individualism, masculinity, uncertainty avoidance, long term orientation, and indulgence. In their empirical analysis Gokmen et al. reported significant and positive associations between COVID-19 incident rates in selected European countries’ and both individualism and indulgence while power distance was negatively related to COVID-19 incidence. Individualism in Hofstede’s sense is defined as a social preference for a “loosely-knit social framework in which individuals are expected to take care of only themselves and their immediate families” [5]. The link between individualism and COVID-19 cases was interpreted in the way that individualistic cultures feel less responsible for achieving collective goals such as the containment of a pandemic and, therefore, are less supportive towards policies that restrict individual rights even though the society at large might benefit from them. Moreover, in individualistic countries privacy is highly valued [6] which make some containment measures such as contact tracking difficult. Conversely, a more collectivistic society people are more willing to sacrifice personal freedom for the sake of collective benefits [7].

A similar argument was put forward for the impact of indulgence. Indulgence stands for  “a society that allows relatively free gratification of basic and natural human drives related to enjoying life and having fun” [3]. The antipode of indulgent cultures are restraint societies with strict social norms [8]. It was argued, that such societies more easily accept restrictions in daily life at the expense of quality of life. Lastly, power distance describes to which extent to less powerful members of societies accept an unequal distribution of power. It is argued that societies with a high power distance are less likely to follow orders by state authorities which makes containment measures less effective. 

Another study that has empirically linked cultural variables with COVID-19 outcomes is by Chen et al [9]. They found that individualistic cultures are slow to respond to the COVID-19 pandemic and that a higher speed of implementing containment measures leads to a lower COVID-19 related mortality. This result partly echoes the findings Cao et al.[10] who established a negative correlation between Hofstede’s individualism dimension and both incidence and mortality rates. 

Notwithstanding the substantial new insights that those studies provided to the reader, none of them controlled for clinical parameters that are known to be associated with COVID-19 related mortality. Clinical variables might be might be even more important in terms of effect size and their omittance would result in an underspecified model. Omitted variables can lead to biased and inconsistent coefficient estimates. We therefore want to have a fresh look on the relationship between culture and COVID-19 related mortality while taking set of variables into consideration. For instance, age, obesity, and co-morbidities such as cancer, cardiovascular diseases, and high blood pressure are found to be associated with both COVID-19 incidence and mortality [11,12,13,14]. Controlling for those factors allows for a better specification of the model and therefore might lead to different results. 

2. Methods

Variables and Data Sources

We assembled a cross sectional dataset consisting of 53 countries spanning from the beginning of the pandemic (Jan 1st 2020) until April 27th 2021 

The outcome variable in our model is the COVID-19 related mortality rate defined as the cumulative number of confirmed deaths per million population and are taken from a publicly accessible data set assembled by the University of Oxford [15]. We prefer the mortality rate over infection counts because a large proportion of infections with COVID-19 from non-risk individuals are mild and therefore, measurement errors might occur with regard to the incidence rate. In addition, there are large differences between countries in the availability and quality of testing. In contrast, data on death rates can be more easily and reliably collected. 

The independent variables of interest are the six cultural dimensions of Hofstede. The current six cultural dimensions according to Hofstede are the power distance index (PDI), individualism versus collectivism (IDV), masculinity versus feminism (MAS), the uncertainty avoidance index (UAI), long term versus short term orientation (LTO) and indulgence versus restraint (IVR). PDI describes how a society deals with inequality, with a higher value associated with greater acceptance of inequality. High IDV scores indicate a preference for a loose society with a focus on the individual. The collective self-image here is defined more as "me" as opposed to "we." The MAS dimension attempts to enact a society between competitive and consensual values. UAI is concerned with a society's attitude toward an uncertain future. A high value in this dimension is associated with intolerance of deviant behavior, while a low value suggests a more relaxed attitude. High LTO scores describe a society with a high focus on the future. Low values describe more traditional attitudes, with aversion to major change. The sixth dimension, IVR, places nations between free enjoyment and restraint [16]. The original research consisted of the results of a value study of over 110,000 IBM employees from 40 countries, between 1967 and 1973. Hofstede used factor analysis to define four cultural dimensions based on this data set, which were later expanded to six dimensions after replication studies in other countries. Although the model and its simplification of culture to a few cultural dimensions is not considered uncontroversial, especially for multiethnic societies [17], it still represents an important tool for understanding culture in an international comparison as well as in a professional context [18]. 

The following covariates were used in our analysis: The prevalence of co-morbidities that might be associated with a severe course of the diseases are from the World Health Organization. This included the percentage of the population with obesity (i.e. body mass index [BMI] above 30) [19], and hypertension [20], that was reported to be a risk factor for severe disease aggravation [21]. From the same database we also include the percentage of the population that regularly uses tobacco products [22], because patients with a smoking history have a higher likelihood of developing a more severe disease course [23]. The prevalence of cancer as another relevant co-morbidity is from the Global Health Data Exchange [24]. We also include the proportion of the total population over age 65 from the Worldbank database [25] because age is major risk factor for COVID-19 related mortality [26,27]. Finally, we use the proportion of the population between 20 and 79 with diabetes [28] because a recent meta-analysis suggested that diabetes is one major cause of COVID-19 related mortality [29]. 

Besides the clinical covariates we use the GDP per capita as an economic performance indicator. People in richer countries have better access to health care and modern treatments. Values are in constant 2017 USD using purchasing power parities as exchange rates. Those data are taken from the World Bank as well [30]. Moreover, we take the Global Health Security Index that represents an attempt to quantify, how prepared a country is for a pandemic outbreak [31]. Other variables that were used to explain COVID-19 mortality include the connectedness of a country measured as the number of international arrivals [32]. We use this as well and take the data from the Worldbank [33].

Finally, we include two variables that aim to assess government measures of pandemic containment. First, we take the “COVID-19 government response tracker” stringency index [15]. This index consists of several subscores on different restriction measures, as well as an overall index, which is used in our analysis It is intended to quantify the severity of national pandemic containment measures. Here, only the overall index is included as the mean of the daily total index since the beginning of the pandemic as a control variable. Second, we include the stringency speed index that relates to the time it takes for the government to impose containment measures. Reaction speed is considered to be more effective than the level of stringency. For example, a study suggests the elapsed time between the first national death and the ban on public events as a significant factor [34]. Stringency Speed is calculated as the as the marginal rate of change of stringency index and is taken from Chen et al. [9].

Statistical Analysis

We employ a multivariate regression model that examines the association of Hofstede's six cultural dimensions with the cumulative number of COVID-19 deaths in a country. The value a country achieves on each dimension is included as a single variable in the model. Furthermore, the model includes a number of control variables that were outlined   above. Absolute data are divided by the number of million inhabitants to make them comparable, relative data and indices are kept unchanged. Correlations between individual variables are calculated and displayed in a correlation matrix. Combinations of variables with the greatest correlation coefficients are then examined by means of an interaction plot to examine whether interaction effects between variables are present.  

3. Results

The country selection includes 53 culturally and geographically diverse nations, as shown in figure 1. Regional differences in deaths per million population are represented by different color shadings. 

Figure 2 displays the correlation matrix of the variables. It is evident that some of the independent variables are related, mostly strongly power distance and individualism with correlation coefficient of -0.74. 

The results of the best fitted multivariable model is reported in table 2. 

The result of the complete regression model reaches an adjusted R-squared of 0.53. The model as a whole is highly significant according to the F-test result. After performing a Breusch-Pagan test, the null hypothesis of homoscedasticity of the data cannot be rejected. Furthermore, scale-location, normal Q-Q and residual vs. fitted plots suggest a good model fit.

According to our results prevalence of obesity and cancer, as well as the share of the population older than 65 years is significantly associated with the COVID-19 related death rate. The stringency index is positively associated with mortality as well while stringency speed exhibit a significant negative correlation. 

Masculinity and long term orientation is positively associated with the mortality rate. The effect of Individualism is moderated by power distance in that its associated with a higher mortality rate in countries that are characterized by a high power distance index. Conversely, individualism is negatively related to the mortality rate when power distance is low. However, the significance level of this interaction effect is only 10%. The interaction plot in figure 4 further illustrates the relation between individualism and mortality contingent on power distance. Other interaction effects were not detected and other variables were insignificant.

4. Discussion

The results partly mirror those of previous studies. In particular, the influence of obesity and age on death rates is supported by the model. The speed of implementation of containment measures is associated with lower mortality, which was reported in previous studies as well [35;36]. Somewhat surprising, the average level of the stringency index was positively correlated with mortality, which is at odds with most of the literature [29,37,38]. Therefore, our results could be due to reverse causality that would lead to a biased parameter estimate. Controlling for this so called simultaneity bias would require the utilization of instruments which are usually hard to find. On the other hand, a recent study did not find evidence that COVID-19 related death trigger interventions, which would cause reverse causality issues. Instead it was shown that are more likely to follow the policies of nearby countries leading to badly timed lookdowns that are implemented either to early or too late [39]. This finding suggests that reverse causality might be less of a problem. Using instrumental variables and other methods to control for endogeneity Bjørnskov (2021) found no evidence of reversed causality as well and concluded that stringency measures had no or even a positive effect on mortality [40]. He used the discretionary power of governments as instrument variable because policy responses are stronger in countries in which the constitution allows for strong political responses irrespective of the pandemic situation [41]. For this reason, the instrument is correlated with the stringency index but not with mortality.  Appositive relation between stringency and mortality could be observed when people subtitle gatherings in restaurants with meetings at home where they pay less attention to hygiene and distance keeping. 

Turning to the culture variables, we find that individualism is related with mortality and moderated by power distance. Specifically, higher scores on the individualism dimension are related to higher death counts in presence of high values in power distance and vice versa. However, this interaction effect is only weakly significant. 

The results suggests that individualistic societies can navigate well through the pandemic when they are also characterized by low power distance. Power distance on the other hand reflects decision making of authorities in a persuasive vis a vis an autocratic style and also reflects to what degree “critical thinking” is valued by a society. Therefore, our finding contradicts with the popular narrative that autocratic societies such as China are better positioned to manage  the COVID 19 pandemic. For instance the UK based Daily Mail newspaper recently titled that “China was able to get a quicker grip on the virus, because people are more obedient and follow the rules” [42]. With regard to China’s containment policies the term ‘Authoritarian Advantage’ was coined [43]. Our results rather indicate that individualism that goes along with personal responsibility and critical attitude towards authority does not necessarily constitute an obstacle to the management of a pandemic.

We also observe that long term orientation shows a significant positive influence on mortality. According to Hofstede, pragmatic societies score high in this dimension in their strive to prepare for the future. Traditional societies on the other hand that value time-honored   norms have relatively low scores. Moreover, in low scoring countries there is more focus on short-term and quick results. This mind set might have helped to handle the pandemic. 

Finally, we observe a positive coefficient for the masculinity variable. The counterpart of masculinity is feminity that stands for a preference for cooperation, and caring for the weak and the society at large is more consensus oriented in those countries. Apparently, those characteristics are well suited to navigate through the pandemic. 

In conclusion, our results suggests a significant influence of some of the culture dimensions on COVID-19 related mortality. To us, the most striking finding is that after controlling for major clinical variables, we observe that low power distance is associated with lower death rates in individualistic societies while it is associated with higher death rates more collectivistic societies. We therefore conclude that collectivism and obedience is not necessarily a good prerequisite for handling this pandemic as previous studies suggested. We also conclude that there is no universal blue print for COVID-19 containment measures and an effective policy needs to find country-specific approaches that build on local cultural circumstances. 

Declarations

Conflict of interest: On behalf of all authors, the corresponding author states that there is no conflict of interest.

Competing interests: The authors declare no competing interests

References

1   Pifarré i Arolas, H., Acosta, E., López-Casasnovas, G. et al. Years of life lost to COVID-19 in 81 countries. Scientific Reports 2021;11:3504.

2   Gokmen Y, Baskici C, Ercil Y. The impact of national culture on the increase of COVID-19: A cross-country analysis of European countries. International Journal of Intercultural Relations 2021;81:1-8.

3   Hofstede G. Culture’s consequences: Comparing values, behaviors, institutions

and organizations across nations. Sage 2001

4   Hofstede G, McCrae R. Personality and culture revisited: Linking traits and dimensions of culture. Cross-Cultural Research 2004;38(1):52–88.

5   Hofstede G. Nationl cultures, available at : https://hi.hofstede-insights.com/national-culture, accessed on Sep 7th, 2021.

6   Markus H, Kitayama S. Culture and the self: Implications for cognition, emotion, and motivation. Psychol Review 1991;98(2):224–253.

7   Porcher S. Culture and the Quality of Government. Public Administration Review 2021;81:333-343.

8   DeBode J, Haggard D, Haggard K. Economic freedom and Hofstede’s cultural dimensions. International 

Journal of Organization Theory and Behavior 2002;23(1):65–84.

9   Chen D, Peng D, Rieger M, Wang M. Institutional and cultural determinants of speed of government responses during COVID-19 pandemic. Humanities and Social Sciences Communications 2021;8:171.

10  Cao C, Li N, Liu L. Do national cultures matter in the containment of COVID-19? International Journal of Sociology and Social Policy 2020;40(9/10):939–961.

11  Chaudhry R, Dranitsaris G, Mubashir T, Bartoszko J, Riazi S. A country level analysis measuring the impact of government actions, country preparedness and socioeconomic factors on COVID-19 mortality and related health outcomes. EClinicalMedicine 2020;25:100464.

12  Sá F. Socioeconomic Determinants of Covid-19 Infections and Mortality: Evidence from England and Wales. CEPR Discussion Paper 2020:DP14781

13  Toraih E, Elshazli R, Hussein M, Elgaml A, Amin M, El‐Mowafy M, El‐Mesery M, Ellythy A, Duchesne J, Killackey M, Ferdinand K, Kandil E, Fawzy M. Association of cardiac biomarkers and comorbidities with increased mortality, severity, and cardiac injury in COVID‐19 patients: A meta‐regression and decision tree analysis. Journal of Medical Virology 2020;92(11):2473-2488.

14  Sannigrahi S, Pilla F, Basu B, Basu A, Molter A. Examining the association between socio-demographic composition and COVID-19 fatalities in the European region using spatial regression approach. Sustainable Cities and Society 2020;62:102418.

15  Hale T, Angrist N, Goldszmidt R, Kira B, Petherick A, Phillips T, Webster S, Cameron-Blake E, Hallas L, Majumdar S, Tatlow H. A global panel database of pandemic policies (Oxford COVID-19 Government Response Tracker). Nature Human Behaviour 2021;5:529–538.

16  Hofstede G. Culture's consequences. 2nd ed. Newbury Park: Sage 2001

17  Baskerville R. Hofstede never studied culture. Accounting, Organizations and Society 2003;28(1):1-14.

18  de Mooij M, Hofstede G. The Hofstede model. International Journal of Advertising 2010;29(1):85-110.

19  WHO, 2021. GHO | By indicator | Prevalence of obesity among adults, BMI & GreaterEqual; 30 (crude estimate) (%) (Noncommunicable diseases). [online] Available at: <https://apps.who.int/gho/data/node.imr.NCD_BMI_30C?lang=en> [Accessed 16 September 2021].

20  WHO, 2021. GHO | By indicator | Raised blood pressure (SBP>=140 OR DBP>=90) (age-standardized estimate) (Global Dementia Observatory (GDO)). [online] Available at: <https://apps.who.int/gho/data/node.imr.BP_04?lang=en> [Accessed 16 September 2021].

21  Tang, Q., Liu, Y., Fu, Y. et al. A comprehensive evaluation of early potential risk factors for disease aggravation in patients with COVID-19. Scientific Reports 2021;11:8062

22  WHO, 2021. GHO | By indicator | Estimate of current tobacco use prevalence (%) (age-standardized rate) (Tobacco control). [online] Available at: <https://apps.who.int/gho/data/node.imr.M_Est_tob_curr_std?lang=en> [Accessed 16 September 2021].

23  Shastri MD, Shukla SD, Chong WC, Kc R, Dua K, Patel RP, Peterson GM, O'Toole RF. Smoking and COVID-19: What we know so far. Respiratory Medicine 2021;176:106237. 

24  Global Health Data Exchange GHDx healthdata, available at: http://ghdx.healthdata.org/, accessed on 26th. July 2021.

25  Worldbank. Population ages 65 and above in 2019 (% of total population), available at: https://data.worldbank.org/indicator/SP.POP.65UP.TO.ZS, accessed on July 31rd 2021.

26  Kar S, Chawla R, Haranath S . Multivariable mortality risk prediction using machine learning for COVID-19 patients at admission (AICOVID). Scientific Reports 2021;11:12801.

27  Undurraga EA, Chowell G, Mizumoto K. COVID-19 case fatality risk by age and gender in a high testing setting in Latin America: Chile, March-August 2020. Infectious Diseases and Poverty 2021;10(1):11

28  World Bank. Diabetes prevalence (% of population ages 20 to 79) | Data. [online] Available at: <https://data.worldbank.org/indicator/SH.STA.DIAB.ZS> [Accessed 15 September 2021]. 

29  Corona G, Pizzocaro A, Vena W, Rastrelli G, Semeraro F, Isidori AM, Pivonello R, Salonia A, Sforza A, Maggi M. Diabetes is most important cause for mortality in COVID-19 hospitalized patients: Systematic review and meta-analysis. Reviews in Endocrine and Metabolic Disorders 2021;22(2):275-296. 

30  Worldbank. GDP per capita, PPP (constant 2017 international $) in 2019, available at https://data.worldbank.org/indicator/NY.GDP.PCAP.PP.KD, accessed on 11 September 2021.

31  GHS Index. 2021. The Global Health Security Index. <https://www.ghsindex.org/> [Accessed 26. July 2021].

32  Pana T, Bhattacharya S, Gamble D, Pasdar Z, Szlachetka W, Perdomo-Lampignano J, Ewers K, McLernon D, Myint P. Country-level determinants of the severity of the first global wave of the COVID-19 pandemic: an ecological study. BMJ Open 2021;11(2):e042034.

33  Worldbank. International tourism, number of arrivals, 2019, available at https://data.worldbank.org/indicator/ST.INT.ARVL, accessed on 11 September 2021.

34  Fountoulakis K, Fountoulakis N, Koupidis S, Prezerakos P. Factors determining different death rates because of the COVID-19 outbreak among countries. Journal of Public Health 2020;42(4):681-687.

35  Deb P, Furceri D, Ostry J, Tawk N. The Effect of Containment Measures on the COVID-19 Pandemic. IMF Working Paper 2020;WP/20/159.

36  Fuller JA, Hakim A, Victory KR, Date K, Lynch M, Dahl B, Henao O; CDC COVID-19 Response Team. Mitigation Policies and COVID-19-Associated Mortality - 37 European Countries, January 23-June 30, 2020. MMWR Morbidity and Mortality Weekly Reports 2021;70(2):58-62.

37  Bilgili F, Dundar M, Kuskaya S, Lorente D, Ünlü F, Gençoglu P, Mugaloglu E. The Age Structure, Stringency Policy, Income, and Spread of Coronavirus Disease 2019: Evidence From 209 Countries. Front. Psychol 2021 11:632192.

38  Hale T, Angrist N, Hale AJ, Kira B, Majumdar S, et al. Government responses and COVID-19 deaths: Global 

evidence across multiple pandemic waves. PLOS ONE 2021;16(7): e0253116

39  Sebhatu A, Wennberg K, Arora-Jonsson S, Lindberg S. Explaining the homogeneous diffusion of COVID-19 nonpharmaceutical interventions across heterogeneous countries. Proceedings of the National Academy of Sciences 2020;117(35):21201-21208.

40  Bjørnskov C. Did Lockdown Work? An Economist’s Cross-Country Comparison. CESifo Economic Studies 2021:ifab003.

41  Bjørnskov C, Voigt S. This Time is Different? On the Use of Emergency Measures during the Corona Pandemic, ILE Working Paper 2020;36, University of Hamburg.

42  Daily Mail. In China if there is a rule, we follow the rule': Chinese London resident who was trapped in Wuhan during peak of pandemic says country was able to control virus quicker than UK because people are more obedient. Oct 18, 2020, available at https://www.dailymail.co.uk/news/article-8852629/Chinese-London-resident-says-country-control-virus-people-follow-rules.html, accessed on Sept 13, 2021.

43  Schwartz J. Coronavirus and China’s ‘Authoritarian Advantage. The Diplomat 2020, available at: https://thediplomat.com/2020/02/coronavirus-and-chinas-authoritarian-advantage/

Tables

Table 1 Descriptive Statistics

Variable

Designation in model

Min

Max

Median

Std

Cumulative deaths per million population

DeathsMio

1.5

2 577.7

962.9

713.8

% of population with elevated blood pressure (SBP>=140 or DBP>=90)

Blood_Pressure

11.0

32.4

20.6

5.1

% of population with BMI > 30

Obesity

3.4

37.3

24.5

8.4

% of the population with regular tobacco use

Tobacco_Use

7.9

44.7

24.3

8.3

GDP per capita (PPP) in constant 2017 USD

GdpCapita

4 753.7

114 323.4

33 221.5

22 833.5

Annual foreign arrivals per million inhabitants

ArrivalsMio

1 981.0

14 756 238.0

763 516.3

2 473 137.0

% of population over 65

Population65

5.2

28.0

17.3

5.7

% of the population between 20 and 79 with diabetes

Diabetes

3.2

16.7

6.6

2.5

Cancer cases per 100 000 inhabitants

Cancer

100.2

1 278.5

319.4

196.1

Average daily Oxford Stringency Index since the beginning of the pandemic

Ox

33.3

72.4

55.3

8.8

Stringency Speed

speed

0.09

1.524

0.313

0.2521

Global Health Security overall index

Ghs

35.0

83.5

57.6

11.3

Hofstede – power distance

pdi

11.0

104.0

60.0

21.0

Hofstede – individualism

idv

13.0

91.0

46.0

23.3

Hofstede – uncertainty avoidance

uai

8.0

112.0

70.0

23.6

Hofstede – masculinity

mas

5.0

95.0

48.0

19.8

Hofstede – long term orientation

lto

13.0

100.0

47.0

22.3

Hofstede – indulgence

ivr

13.0

97.0

48.0

20.9

Interaction term – power distance * individualism

pdi:idv

605.0

4 875.0

2 378.0

899.9

 Table 2 Regression Results

 

Estimate

Standard Error

P-Value

Significance

Intercept

-2294

1768

0.203899

 

Blood_Pressure

43.50

27.32

0.121443

 

Obesity

40.69

16.59

0.019980

*

GdpCapita

-0.002069

0.005079

0.686620

 

Tobacco_Use

-18.85

12.06

0.128266

 

ArrivalsMio

0.00001913

0.00003957

0.632150

 

Population65

69.02

31.15

0.034174

*

Diabetes

-21.47

37.21

0.568062

 

Cancer

1.436

0.7699

0.071593

.

Ox

51.80

11.83

0.000126

***

Speed

-1346

505.9

0.012274

*

Ghs

0.6626

11.22

0.953279

 

PDI

-30.57

20.05

0.137378

 

UAI

-11.46

8.357

0.180112

 

IDV

-60.24

27.18

0.034142

*

MAS

8.156

4.774

0.097521

.

LTO

16.09

5.818

0.009479

**

LVR

8.510

6.271

0.184582

 

Pdi:idv

0.8129

0.4107

0.056697

.

Residual Standard Error: 504.2 on 31 degrees of freedom

Multiple R-squared: 0.7002

Adjusted R-squared: 0.5262

F-statistic: 4.023 on 18 and 31 degrees of freedom, p-value 0.0003376

Significance level: 0.001         (***); 0.01(**), 0.05 (*); 0.1 (.1)