Socioeconomic, demographic and healthcare determinants of the COVID-19 pandemic: an ecological study of Spain
BMC Public Health volume 21, Article number: 606 (2021)
The coronavirus disease 2019 (COVID-19) pandemic has posed a major challenge to health, economic and political systems around the world. Understanding the socioeconomic, demographic and health determinants affecting the pandemic is of interest to stakeholders. The purpose of this ecological study is to analyse the effect of the different socioeconomic, demographic and healthcare determinants on the mortality rate and estimated cumulative incidence of COVID-19 first wave in the Spanish regions.
From the available data of the 17 Spanish regions (Autonomous Communities), we have carried out an ecological study through multivariate linear regression using ordinary least squares. To do this, we conducted an analysis using two distinct dependent variables: the logarithm of mortality rate per 1,000,000 inhabitants and the estimated cumulative incidence. The study has 12 explanatory variables.
After applying the backward stepwise multivariate analysis, we obtained a model with nine significant variables at different levels for mortality rate and a model with seven significant variables for estimated cumulative incidence. Among them, six variables are statistically significant and of the same sign in both models: “Nursing homes beds”, “Proportion of care homes over 100 beds”, “Log GDP per capita”, “Aeroplane passengers”, “Proportion of urban people”, and the dummy variable “Island region”.
The different socioeconomic, demographic and healthcare determinants of each region have a significant effect on the mortality rate and estimated cumulative incidence of COVID-19 in territories where the measures initially adopted to control the pandemic have been identical.
The emergence of the coronavirus disease 2019 (COVID-19) pandemic has created a global health emergency  that threatens to bring about profound social, health and economic changes in all countries.
The analysis of the effect of the different social, economic and demographic determinants on the incidence of COVID-19 is a subject of interest to researchers, as it can help to understand how the particular conditions of each geographical area have affected the spread and the impact of the disease, and can help policymakers to make decisions that can help to deal with future pandemics safely [2, 3]. Different ecological studies on the COVID-19 pandemic have shown significant associations between different socioeconomic, demographic and healthcare determinants, and the outcomes of the pandemic, both at the regional level in Europe [4,5,6], China , Iran , and at the country level [9, 10]. It can also be noted that the results in terms of the number of infected and dead in each country or region, as well as the effect that the pandemic is having on the economies, present enormous variations. Therefore, it is considered necessary to analyse the effect that the different social, demographic and economic determinants of an exogenous nature may have had on the effect of the pandemic in the different territories.
Spain has been one of the developed countries most affected by the first wave of the COVID-19 pandemic, with one of the highest mortality rates per inhabitant in the world during that period. The first death from COVID-19 in Spain occurred on 13 February 2020 and the country decreed a very restrictive, total lockdown on 15 March 2020, when the country already had 121 deaths. In any case, the impact of the pandemic in Spain has not been similar in all its regions, with significant differences in mortality rates and cumulative incidence depending on the geographical location. Spain has its own territorial and political characteristics that make it suitable for this analysis since, due to territorial management issues, two administrations with direct responsibilities for the management of the pandemic coexist. On the one hand, the national government is responsible for defining both public health policy and legal framework for all its territories in a comprehensive manner, so that the measures adopted to contain and mitigate the effects of the pandemic have been applied uniformly in all Spanish regions; it also corresponds to the central government to take general economic and social measures for the country as a whole. On the other hand, the organisation of the Spanish territory is characterised by the existence of different Spanish regions (called Autonomous Communities) that have a wide level of decentralisation, being responsible for, among other things, healthcare in their territories and long-term care policies. Besides, every Autonomous Community has its own social and economic determinants, which are the result of its degree of economic and social development, geographical location, natural resources and local history.
In this way, given the adoption of the same measures of containment of the effects of the pandemic for all the territories, the present ecological study allows us to deepen in the effect that the own different determinants of each territory have had on mortality and estimated cumulative incidence of the COVID-19 pandemic. The ecological study is a type of aggregate study that allows us to make inferences about the effect of different variables on the dependent variables, using groups for analysis instead of individuals, that has been shown to be effective in the analysis of similar situations.
Administratively, Spain is divided into 17 regions (Autonomous Communities) with a high degree of decentralisation and autonomy in the management of their resources, where each is directly responsible for both the expenditure and management of the healthcare systems in its region, as well as for other aspects of economic and social nature.
This study uses two distinct dependent variables, used separately or together in studies similar to the present [4,5,6,7,8,9,10]: first, the mortality rate of COVID-19 per 1,000,000 inhabitants in each of the Spanish regions, calculated from the official figures of cumulative mortality rate as of 23 May 2020 in each of the regions, and provided by the Spanish Ministry of Health. The population data corresponds to 2019 and have been taken from the Spanish National Institute of Statistics. The variable is presented in logarithmic form.
The second dependent variable used is the estimated cumulative incidence, extracted from the second wave (May 2020) of the national seroprevalence studies carried out by the Spanish government . This variable is expressed as a decimal. We chose this variable due to the lack of reliability of infection data confirmed by tests during the period of the study, which significantly underestimated the actual data on infected people (235,290 cumulative infections on 23 May 2020 confirmed by tests compared to the 2.44 M infected people estimated by the seroprevalence study; which is 5.2% of the total Spanish population).
As explanatory variables, 12 variables have been finally chosen, and classified into three groups: (a) healthcare and long-term care variables, (b) economic variables and (c) environmental, sociodemographic and geographical variables. The definition of each variable is given in Table 1.
A study of collinearity of the variables (see next section) was carried out through the variance inflation factor (VIF), from which it was concluded that there was no collinearity among the selected variables. Previously, two independent variables initially incorporated into the model were eliminated due to high collinearity: population density and research and development (R&D) expenditure per capita.
Likewise, in order to deepen the analysis of the social determinants, the dependent variable “Proportion of immigrant population in each region” was also included in the statistical analyses of both models, being in all cases not significant and not having any effect on the parameters of the two models; therefore, it has not been included in the final models to simplify their presentation.
We have used multivariate ordinary least squares (OLS) regression model to study the relationship between each of the dependent variables and the set of explanatory variables.
The model is specified as follows:
where y represents the dependent variable, xi are the explanatory variables, β0 is the intercept, βi are the regression coefficients and ε is a random component.
For the construction of the model, regressions analysis was performed through backward stepwise multiple regression analysis. Starting with the fully saturated model, the least significant independent variable was eliminated in every step, until a model with all significant variables was achieved.
An important assumption in multiple linear regressions is the independence between independent variables, so that there is no relationship between each of them, to ensure the validity of the estimation using OLS, avoiding multicollinearity problems. As indicated, in this study the evaluation of multicollinearity was performed with VIF, defined as:
where p is the number of explanatory variables, and R2 is the square of the coefficient of determination of the regression of Xj on other covariates (p - 1), so that:
if 0 < VIF < 5, there is no evidence of a multicollinearity problem;
if 5 < VIF < 10, there is a moderate multicollinearity problem; and.
if VIF > 10, there is a serious multicollinearity problem of variables.
The significance tests are two-sided. We consider three levels of significance for the p-value: < 0.1, < 0.05 and < 0.01. All statistical analyses were performed using R software.
We have used the standardised residues to perform the analysis of the assumptions of the multiple linear regression.
In this study, all the 17 Spanish regions (Autonomous Communities) have been included. The analysis of the descriptive statistics can be seen in Table 2. Descriptive analysis of non-log transformed variables can be seen in Supplementary material, Table 1.
For the analysis of the relationship between variables, multivariate linear regressions using OLS have been carried out, using a backward stepwise approach. For the dependent variable “Log mortality per 1,000,000 inh”, the application of the selection of variables using the backward method eliminated three variables, generating the most parsimonious model with nine variables, all of them significant in different degrees. The four resulting models are summarised in Table 3. All models are statistically significant. In the final model (model 4), the VIF of the independent variables are below 5, so there are no multicollinearity problems.
With regard to the parameters obtained, in the three models generated, the logarithm of average public health expenditure over the last 10 years is of a positive sign and statistically significant at 10% in model 4 (1.774; p-value = 0.091).
Also, the variables “Number of physicians per 100,000 inhabitants” (negative sign), “Nursing home beds per 100 elderly people” (positive sign), “Aeroplane passengers” (positive sign), “Island region” (negative sign) and “Log GDP per capita” (positive sign) have statistical significance to different degrees in the three models generated.
Within model 4, the variables “Proportion of nursing homes >100 beds” and “Proportion of urban population” reach statistical significance at 10%, so that model 4 retains 9 of the 12 independent variables used, with statistical significance to different degrees.
For the dependent variable “Estimated cumulative incidence”, the backward selection of the stepwise multivariate regression analysis eliminates 5 variables, resulting in the model with 7 independent variables (model 6); all of them statistically significant in different degrees and with absence of collinearity problems. The result of the backward stepwise regression is summarised in Table 4.
In the aforementioned model 6, the variables “Nursing home beds per 100 elderly people” (positive sign), “Aeroplane passengers” (positive sign), “Island region” (negative sign) “Log GDP per capita” (positive sign), “Proportion of nursing homes >100 beds” (positive sign), “Proportion of urban population” (negative sign) and “Proportion of elderly people” (negative sign) are statistically significant at different levels.
In this ecological study, we have evaluated the impact of the different determinants at regional level that have influenced the degree of mortality rate and estimated cumulative incidence in the first wave of COVID-19 in Spain, assuming that the systems of surveillance and control of the pandemic, carried out by the national government, have been identical in all regions.
The difficulty of having individual data for analysis in this phase of the pandemic makes ecological studies a useful tool for understanding the factors that have affected the development of the pandemic. In any case, it cannot be ignored that this type of study is susceptible to bias, where the so-called ecological bias—the temptation to extrapolate the results obtained through this type of studies with aggregate data to an individual level— is the most important.
Knowing this limitation, all necessary precautions have been adopted to minimise such ecological bias, to the extent that the available data have allowed it. These precautions include the homogeneity of the analysis groups  by using the smallest regional unit (the Autonomous Community) for which data are available, and the mitigation of the existence of confounding factors when working with sociodemographic factors , especially the multicollinearity bias. In any case, the results obtained in this study should be taken with due caution, given the limitations indicated.
Having pointed this out, the present analysis carried out in Spain allows us to compare regions whose response to the initial control of the pandemic has been identical, since the measures were taken by a higher decision-making body, which allows us to isolate this aspect in the analysis of the effect of the different determinants that may have influenced the expansion and degree of affectation of the pandemic.
The statistical analysis carried out for each of the two dependent variables presents six variables with statistical significance, which are associated with both mortality and cumulative incidence.
The first finding in this study is the strong association between long-term care resources, both with mortality rate and the estimated cumulative prevalence. Thus, we can observe that the presence of a greater number of care homes is the variable most strongly associated with the mortality rate of COVID-19 in Spain, where the increase of 1 care home per 100 elderly inhabitants is significantly associated with an increase in mortality rate of 43%, keeping all other variables constant. This variable is also positively associated with a greater cumulative incidence of the disease. These findings confirm previous findings during the current pandemic [14,15,16]. In a similar manner, we also found a positive and statistically significant association between mortality rate and estimated cumulative incidence with the presence in the region of a greater proportion of nursing homes which size is greater than 100 places; this reinforces the idea of the relevance that the population of nursing homes has had in all the outcome figures of this pandemic in Spain.
Concerning the above, it should be noted that the proportion of the population over 64 years of age has no statistically significant effect on the mortality rate of COVID-19, showing a negative and statistically significant relationship with the estimated cumulative incidence, similar to what has been described in other regional studies, such as in the case of Germany  and northern Italy . It seems that the effects of the pandemic in Spain have been associated not so much to the high ageing of its population as to the institutionalised elderly people. This may be explained by the fact that the non-institutionalised elderly population in Spain has followed stricter measures of self-isolation and restriction of social relations than the rest of the population, as advised by the health authorities (in Spain, during the lockdown, specific walking hours were decreed exclusively for the elderly). These measures have allowed greater control of the spread of the disease in this age group. But those institutionalised in care homes are concentrated, of advanced age and with multiple pathologies, which facilitates the spread of the virus.
Another important finding in the present study is related to the economic variables, where we find a positive and statistically significant association between the GDP per capita and both mortality and estimated cumulative incidence. Bearing in mind that the variable is presented in its logarithmic form, as is the mortality rate, we can interpret that a 1% increase in the GDP per capita is associated with a 3.21% increase in mortality rate, keeping all other variables constant. This positive association between the health outcomes of the pandemic and the income level has been observed in similar studies [10, 17]. One possible explanation is that a higher GDP per capita is an indicator of greater economic and trade dynamism, which would be indicative of greater mobility, both within and between regions, as well as with foreign countries, such that the likelihood of contagion and death would increase. In addition, we find that the number of air passengers received by each region during February 2020 (the month before the great expansion of the pandemic in Spain) presents a positive sign and is statistically significant for both dependent variables used so that an association is observed between the increased mobility of people and both incidence and mortality of the disease; this result is in agreement with other studies conducted during this pandemic [9, 18].
From the territorial and demographic point of view, it should be noted that the insularity of the territories has meant, keeping all other variables constant, a protective and statistically significant element against the estimated cumulative incidence and mortality rate of COVID-19, which seems to indicate that the geographical barriers have been, in the case of Spain, important aspects in the control of the pandemic. Likewise, our analysis shows a negative and statistically significant association between the proportion of inhabitants living in municipalities of more than 100,000 inhabitants in each region and the mortality rate of COVID-19.
The variables that consider health resources only present a significant association with mortality, not with incidence. This finding is plausible due to the fact that, in Spain, specialised care has a high weight within the Spanish public healthcare system in public health expenditure (47.47% of the total for the year 2017), which influence on the incidence of the disease is less decisive than on mortality; more decisive could be the relatively speaking, substantially lower expenditure in primary care (21.9% of the total) and public health (0.7% of the total).
Considering the above, public health expenditure, measured as the logarithm of average per capita expenditure over the last 10 years, shows a positive and statistically significant association with the mortality rate of COVID-19 in Spain. The average of the last 10 years has been chosen to incorporate into the analysis the lagging effect of public expenditure on the management of the 2020 pandemic . Given that the variable is presented in its logarithmic form, as is the dependent variable, it can be interpreted as the semi-elasticity of the average public health expenditure of the last 10 years over the mortality rate, such that an increase in the average public health expenditure per capita of 1% is associated with an increase in mortality of 1.77%, all other factors remaining constant (model 4).
Although it is common to find in the literature a positive relationship between public health expenditure and the improvement of different health outcomes, such as different mortality rates or life expectancy, the evidence does not always offer conclusive results. Even in developing countries, where it is common to find evidences of a strong positive association between public health expenditure and a reduction in mortality rates and an increase in life expectancy [20, 21]. In some cases, this relationship is not so clear , and this is essentially related to the quality of governance in each analysed country. When the analysis is carried out in developed countries, as is our case, this relationship also appears sometimes rather blurred. In a systematic review of the topic , it is concluded from the analysis of 17 papers that, although the complexity of the delivery system of public healthcare services makes it difficult to demonstrate definitive associations between expenditure and outcomes, financial investments in health often have the potential to improve the health of the community.
In relation to the effect of health expenditure on the mortality rate during a pandemic, there are evidences of a significant relationship between higher health expenditure and lower mortality rates during the 2009 H1N1 influenza pandemic . A study of the current pandemic finds, as does the present study, a positive and significant association between public health expenditure and mortality rate of COVID-19, in an analysis of 96 countries using Bayesian model averaging techniques .
There is evidence [25,26,27] of the existence of decreasing returns on public health expenditure in relation to health outcomes. The findings of our study seem to be compatible with the idea of decreasing returns on public health expenditure, such that, once a certain level of public health expenditure is reached, the additional gains from an additional increase of expenditure disappear, in this case in relation to the mortality rate of COVID-19.
Our study also shows a negative and statistically significant association between the number of practitioners, as well as the number of hospital beds, and the mortality rate of COVID-19, results that confirm previous findings within the present pandemic outcomes ; in the case of practitioners, an increase of 1 unit in the ratio of practitioners per 100,000 inhabitants is associated with a reduction in the mortality rate of COVID-19 by 0.5%. In the case of hospital beds, an increase of 1 unit in the ratio of beds per 100,000 inhabitants reduces the mortality rate of COVID-19 by 0.3%, all other factors remaining constant. It should be kept in mind that in this study both variables represent the total of both resources in both the public and private health sectors, the latter representing approximately 30% of total Spanish healthcare, measured in terms of expenditure.
With regard to environmental aspects, our study does not find a significant relationship between PM2.5 concentration levels and increased mortality rate of COVID-19, as has been shown by other studies of this pandemic [29, 30].
Our study does not either find significant associations between the variables of social structure and inequity, and the outcomes of the pandemic.
The ecological analysis carried out in the present study indicates that the effect of the estimated cumulative incidence and mortality rate of COVID-19 in the different regions of the Spanish state is notably influenced by the social, economic, healthcare and demographic characteristics of each region. In particular, after the analysis of two distinct dependent variables for the same group of independent variables, we found a total of six variables with statistical significance and of the same sign in both models, which leads us to greater robustness of the results obtained in the present study. The detailed analysis of these determinants can help policymakers to make decisions to mitigate the effects of future hypothetic pandemic scenarios.
Availability of data and materials
The datasets supporting the conclusions of this article are public and are available in the repositories listed in Table 1.
Spanish National Aeronautics and Air-Traffic Administration
Gross Domestic Product
Spanish National Institute of Statistics
Particulate Matter below 2.5 μm
Variance Inflation Factor
Ordinary Least Squares
Research and Development
World Health Organization. Coronavirus disease 2019 (COVID-19) situation report—124. Available at: https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200523-covid-19-sitrep-124.pdf?sfvrsn=9626d639_2. Accessed 24 May2020.
Stojkoski V, Utkovski Z, Jolakoski P, Tevdovski D, Kocarev L. The socioeconomic determinants of the coronavirus disease (COVID-19) pandemic; 2020. https://doi.org/10.1101/2020.04.15.20066068.
Khalatbari-Soltani S, Cumming R, Delpierre C, Kelly-Irving M. Importance of collecting data on socioeconomic determinants from the early stage of the COVID-19 outbreak onwards. Epidemiol Community Health. 2020. https://doi.org/10.1136/jech-2020-214297.
Ehlert A. The socio-economic determinants of COVID-19: a spatial analysis of German county level data. medRxiv. 2020. https://doi.org/10.1101/2020.06.25.20140459.
Ginsburgh V, Magerman G. Natali I. COVID-19 and the role of economic conditions in French regional departments. (May), 2020. Working papers ECARES 2020-17, ULB -- Universite Libre de Bruxelles.
Buja A, Paganini S, Cocchio S, Scioni M, Rebba V, Baldo V. Demographic and socio-economic factors, and healthcare resource indicators associated with the rapid spread of COVID-19 in northern Italy: an ecological Study. SSRN Electron J. 2020. https://doi.org/10.2139/ssrn.3582855.
Qiu Y, Chen X, Shi W. Impacts of social and economic factors on the transmission of coronavirus disease 2019 (COVID-19) in China. J Popul Econ. 2020;33(4):1127–72. https://doi.org/10.1007/s00148-020-00778-2.
Ramirez-Aldana R, Gomez-Verjan J, Bello-Chavolla O. Spatial analysis of COVID-19 spread in Iran: Insights into geographical and structural transmission determinants at a province level. medRxiv, page 2020.04.19.20071605, 2020.
Azarpazhooh MR, Morovatdar N, Avan A, Phan TG, Divani AA, Yassi N, et al. COVID-19 pandemic and burden of non-communicable diseases: an ecological study on data of 185 countries. J Stroke Cerebrovasc Dis. 2020;29(9):105089. https://doi.org/10.1016/j.jstrokecerebrovasdis.2020.105089.
Chaudhry R, Dranitsaris G, Mubashirc 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. https://doi.org/10.1016/j.eclinm.2020.100464.
Ene-Covid Study. Final report. National Study os Sero-epidemioligy of infection by SARS-COV-2 in Spain. 2020. Available in: https://www.mscbs.gob.es/ciudadanos/ene-covid/docs/ESTUDIO_ENE-COVID19_INFORME_FINAL.pdf (spanish).
Gili M, Alvárez M, Cortés M, García C, Nieto A, Najera E. Factores determinantes de la Falacia Ecológica y procedimientos Para su control. Rev San Hig Pub. 1986;60:297–307.
Morgenstern H. Ecologic studies in epidemiology: concepts, principles, and methods. Annu Rev Public Health. 1995;16:611.
McMichael TM, Currie DW, Clark S, Pogosjans S, Kay M, Schwartz NG, et al. Epidemiology of COVID-19 in a long-term care facility in king county, Washington. N Engl J Med. 2020;382(21):2005–11. https://doi.org/10.1056/NEJMoa2005412.
Comas-Herrera A, Zalakain J. Mortality associated with COVID-19 outbreaks in care homes: early international evidence, April 12 2020. Article in LTCcovid.org, international long-term care policy network, CPEC-LSE.
Kimball A, Hatfield KM, Arons M, James A, Taylor J, Spicer K, et al. 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(13):377–81. https://doi.org/10.15585/mmwr.mm6913e1.
Gangemi S, Billeci L, Tonacci A. Rich at risk: socio-economic drivers of COVID-19 pandemic spread. Clin Mol Allergy. 2020;18(1):12. https://doi.org/10.1186/s12948-020-00127-4.
Oztig L, Askin O. Human mobility and COVID-19: a negative binomial regression analysis. Public Health. 2020. https://doi.org/10.1016/j.puhe.2020.07.002.
Rivera B. The effects of public health spending on self-assessed health status: an orderer probit model. Appl Econ. 2001;33(10):1313–9. https://doi.org/10.1080/00036840010007146.
Kojo R, Emecheta C, Ogechi M. Public health expenditure and health outcomes in Nigeria. Am J Biomed Life Sci. 2017;5(5):96–102. https://doi.org/10.11648/j.ajbls.20170505.13.
Novignon J, Lawanson A. Health expenditure and child health outcomes in Sub-Saharan Africa. Afr Rev Econ Finance. 2017;9(1):96–121.
Sunil A, Swaroop V. Public spending and outcomes: does governance matter? J Dev Econ. 2008;86(1):96–111.
Singh S. Public health spending and population health. A systematic review. 2014. Am J Prev Med. 2014;47(5):634–40.
Nikolopoulos G, Bagos P, Lytras T, Bonovas S. An ecological study of the determinants of differences in 2009 pandemic influenza mortality rates between countries in Europe. PLoS One. 2011;6(5):e19432. https://doi.org/10.1371/journal.pone.0019432.
Nixon J, Ulmann P. The relationship between health care expenditure and health outcomes. Eur J Health Econom. 2006;7(1):7–18. https://doi.org/10.1007/s10198-005-0336-8.
Murphy KM, Topel R. Diminishing returns?: The costs and benefits of improving health. Perspect Biol Med. 2003;46(3):S108–28. Project MUSE. https://doi.org/10.1353/pbm.2003.0048.
Rodríguez L, Sobrino C. Public HEs, diminishing returns, and life expectancy in Latin American and Caribbean countries: a panel data análisis. J Bras Econ Saúde. 2015;7(2):91–8.
Szulczuk K, Cheema M. COVID-19: the impact of socioeconomic characteristics on the fatality rate. 2020. https://doi.org/10.2139/ssrn.3591551.
Wu X, Nethery RC, Sabath MB, Braun D, Dominici F. Exposure to air pollution and COVID-19 mortality in the United States: a nationwide cross-sectional study. medRxiv. 2020. https://doi.org/10.1101/2020.04.05.20054502.
Pozzer A, Dominici F, Haines A, Witt C, Munzel T, Lelieveld J. Regional and global contributions of air pollution to risk of death from COVID-19. Cardiovasc Res. 2020;116:cvaa288. https://doi.org/10.1093/cvr/cvaa288.
I would like to thank Dr. Antonio Sarria for his support in revising this paper.
This research was fully supported by Alfonso X University and Alfonso X University Foundation. The funding body didn’t have any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Not applicable. Any administrative permissions weren’t required to access the raw data/samples mentioned in the methods. The data used in this study were anonymised before its use.
Consent for publication
The author declares that he has no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
García, C.N. Socioeconomic, demographic and healthcare determinants of the COVID-19 pandemic: an ecological study of Spain. BMC Public Health 21, 606 (2021). https://doi.org/10.1186/s12889-021-10658-3