COVID-19 mortality: educational inequalities and socio-spatial context in two provinces of Argentina

Carlos M. Leveau Guillermo A. Velázquez About the authors

ABSTRACT

With the aim of describing the association between sociodemographic characteristics and contextual factors with COVID-19 mortality during 2020-2021 in the provinces of Mendoza and San Juan in Argentina, we conducted an ecological study, which included the sociodemographic factors: age, sex and educational level, and the contextual factors: poverty and urbanization at the departmental level. The analyses were estimated using negative binomial Bayesian hierarchical models. Educational inequalities existed regardless of socioeconomic context and level of urbanization. The exception was the age group 65 years and older during 2021, which, regardless of educational level, showed a higher risk of death by COVID-19 in departments with high levels of structural poverty. In conclusion, educational inequality is an indicator of social inequality that increases vulnerability to COVID-19 mortality.

Keywords:
Spatial analysis; socioeconomic disparities in health; mortality; SARS-CoV-2; medical geography

KEY MESSAGES

Motivation for the study. There are very few studies on the educational inequalities in COVID-19 mortality, taking into account social contextual factors.

Main findings. We found educational inequalities of COVID-19 mortality during both the 2020 and 2021 waves, regardless of the level of poverty and urbanization in the departments of Mendoza and San Juan provinces (Argentina).

Implications. Preventive policies should focus not only in areas with high levels of poverty, but also in areas with adults of low educational level.

Keywords:
Spatial analysis; socioeconomic disparities in health; mortality; SARS-CoV-2; medical geography

INTRODUCTION

There are ecological studies that have shown the relationship between social and geographical inequalities in mortality due to COVID-19 in Latin American countries 11. Bermudi PMM, Lorenz C, de Aguiar BS, Failla MA, Barrozo LV, Chiaravalloti-Neto F. Spatiotemporal ecological study of COVID-19 mortality in the city of São Paulo, Brazil: shifting of the high mortality risk from areas with the best to those with the worst socio-economic conditions. Travel Med Infect Dis. 2021;39:101945. doi: 10.1016/j.tmaid.2020.101945.
https://doi.org/10.1016/j.tmaid.2020.101...

2. Mena G, Martinez PP, Mahmud AS, Marquet PA, Buckee CO, Santillana M. Socioeconomic status determines COVID-19 incidence and related mortality in Santiago, Chile. Science. 2021;372(6545):eabg5298. doi: 10.1101/2021.01.12.21249682.
https://doi.org/10.1101/2021.01.12.21249...

3. Silva J, Ribeiro-Alves M. Social inequalities and the pandemic of COVID-19: the case of Rio de Janeiro. J Epidemiol Community Health. 2021;75(10):975-979. doi: 10.1136/jech-2020-214724.
https://doi.org/10.1136/jech-2020-214724...
-44. Leveau CM, Soares Bastos L. Desigualdades socio-espaciales de la mortalidad por COVID-19 en tres olas de propagación: un análisis intra-urbano en Argentina. Cad Saúde Pública. 2022;38(5):e00163921. doi: 10.1590/0102-311XES163921.
https://doi.org/10.1590/0102-311XES16392...
, but few studies have analyzed these inequalities considering socioeconomic characteristics of the deceased. The socioeconomic level of the area constitutes a contextual factor that may contribute to a greater spread of COVID-19 55. Albuquerque MV de, Ribeiro LHL. Desigualdade, situação geográfica e sentidos da ação na pandemia da COVID-19 no Brasil. Cad Saúde Pública. 2021;36(12): e00208720. doi: 10.1590/0102-311X00208720.
https://doi.org/10.1590/0102-311X0020872...
, associated with situations of structural poverty (overcrowding of households, low access to green spaces and recreation or low access to healthy food). Urban areas favor greater mobility and contact between different populations, increasing the probability of contagion of communicable diseases, besides, rural areas have very low geographic accessibility to health services, but there seems to be no clear consensus on the differential impact of COVID-19 along the urban-rural gradient 66. Anzalone AJ, Horswell R, Hendricks BM, Chu S, Hillegass WB, Beasley WH, et al. Higher hospitalization and mortality rates among SARS-CoV-2-infected persons in rural America. J Rural Health. 2023;39:39-54. doi: 10.1111/jrh.12689.
https://doi.org/10.1111/jrh.12689...

7. Cuadros DF, Branscum AJ, Mukandavire Z, Miller FD, MacKinnon N. Dynamics of the COVID-19 epidemic in urban and rural areas in the United States. Ann Epidemiol. 2021;59:16-20. doi: 10.1016/j.annepidem.2021.04.007.
https://doi.org/10.1016/j.annepidem.2021...
-88. Petrelli A, Ventura M, Di Napoli A, Mateo-Urdiales A, Pezzotti P, Fabiani M. Geographic heterogeneity of the epidemiological impact of the COVID-19 pandemic in Italy using a socioeconomic proxy-based classification of the national territory. Front Public Health. 2023;11:1143189. doi: 10.3389/fpubh.2023.1143189.
https://doi.org/10.3389/fpubh.2023.11431...
. Regarding demographic factors, there is a higher risk of death in older people and in men 99. Angelici L, Sorge C, Di Martino M, Cappai G, Stafoggia M, Agabiti N, et al. Incidence of SARS-CoV-2 Infection and Related Mortality by Education Level during Three Phases of the 2020 Pandemic: A Population-Based Cohort Study in Rome. J Clin Med. 2022;11(3):877. doi: 10.3390/jcm11030877.
https://doi.org/10.3390/jcm11030877...

10. Chiaravalloti Neto F, Bermudi PMM, Aguiar BS de, Failla MA, Barrozo LV, Toporcov TN. Covid-19 hospital mortality using spatial hierarchical models: cohort design with 74,994 registers. Rev Saúde Pública. 2023; (suppl 1):2s. doi: 10.11606/s1518-8787.2023057004652.
https://doi.org/10.11606/s1518-8787.2023...
-1111. Feldman JM, Bassett MT. Variation in COVID-19 mortality in the US by race and ethnicity and educational attainment. JAMA Netw Open. 2021;4(11):e2135967. doi: 10.1001/jamanetworkopen.2021.35967.
https://doi.org/10.1001/jamanetworkopen....
. Although the social position of individuals, measured by educational level, is associated with a higher risk of COVID-19 mortality 1212. Spijker JJA, Trias-Llimós S. Cause-specific mortality in Spain during the pandemic. Educational differences and its impact on life expectancy. Eur J Public Health 2023;33(3):543-9. doi: 10.1093/eurpub/ckad036.
https://doi.org/10.1093/eurpub/ckad036...
, few studies have studied these inequalities considering both the social position of patients who died due to COVID-19 and the socioeconomic characteristics of the areas where they resided 1313. Li SL, Pereira RH, Prete Jr CA, Zarebski AE, Emanuel L, Alves PJ, et al. Higher risk of death from COVID-19 in low-income and non-White populations of São Paulo, Brazil. BMJ Glob Health. 2021;6(4):e004959. doi: 10.1136/bmjgh-2021-004959.
https://doi.org/10.1136/bmjgh-2021-00495...
.

Our study area were the Argentine provinces of Mendoza and San Juan, both of which have a high proportion of information on educational level 1414. Leveau CM, Hussein M, Tapia-Granados JA, Velázquez GA. Economic fluctuations and educational inequalities in premature ischemic heart disease mortality in Argentina. Cad Saúde Pública. 2023;39(5):e00181222. doi: 10.1590/0102-311xen181222.
https://doi.org/10.1590/0102-311xen18122...
. These provinces had, during 2020-2021, two waves of COVID-19 deaths (Figure 1). The first wave had its mortality peaks in October (Mendoza) and November (San Juan) 2020, and was characterized by measures of social distancing and restriction of population mobility 1515. Boletín Oficial de la República Argentina. Boletín Oficial República Argentina - Aislamiento Social Preventivo Y Obligatorio - Decreto 297/2020. 2020 [consultado el 21 de agosto de 2022]. Disponible en: https://www.boletinoficial.gob.ar/detalleAviso/primera/227042.
https://www.boletinoficial.gob.ar/detall...
. The second wave had mortality peaks in May (Mendoza) and June (San Juan) 2021, characterized by the appearance of the gamma and lambda variants 1616. Freitas ARR, Beckedorff OA, Cavalcanti LP de G, Siqueira AM, Castro DB de, Costa CF da, et al. The emergence of novel SARS-CoV-2 variant P.1 in Amazonas (Brazil) was temporally associated with a change in the age and sex profile of COVID-19 mortality: A population based ecological study. Lancet Reg Health Am. 2021;1:100021. doi: 10.1016/j.lana.2021.100021.
https://doi.org/10.1016/j.lana.2021.1000...
, the former associated with an increase in mortality in the young population, and the beginning of mass vaccination against COVID-19.

Figure 1
Crude mortality rate by COVID-19 (per 100,000 inhabitants) in residents of the provinces of Men-doza (blue line) and San Juan (red line) in Argentina 2020-2021.

Therefore, the aim of this study was to analyze the association between socio-demographic characteristics and contextual factors with COVID-19 mortality, during 2020-2021 in the Argentine provinces of Mendoza and San Juan.

THE STUDY

Study design

This is a spatial ecological study. The study area is comprised of the provinces of Mendoza and San Juan, two jurisdictions that, together with 22 others, make up the territory of the Argentine Republic. Both provinces are located in the west-central part of the country and are subdivided into departments, 18 in Mendoza and 19 in San Juan. According to the 2022 Census, the provinces of Mendoza and San Juan are home to 2,014,533 and 818,234 inhabitants, respectively 1717. Instituto Nacional de Estadística y Censos. Datos provisionales del CENSO 2022. Censo Nac Poblac Hogares Viviendas 2023. [consultado el 31 de mayo de 2023]. Disponible en: https://censo.gob.ar/index.php/datos_provisionales/.
https://censo.gob.ar/index.php/datos_pro...
.

Data and study variables

COVID-19 deaths were identified using the following 10th International Classification of Diseases (ICD-10) codes: U071, U072, U109. Mortality data by age, sex, province of residence, and educational level were obtained from the Ministry of Health 1818. Ministerio de Salud de la Nación. Datos Abiertos del Ministerio de Salud - Defunciones ocurridas y registradas en la República Argentina 2023. [consultado el 16 de diciembre de 2023]. Disponible en: http://datos.salud.gob.ar/.
http://datos.salud.gob.ar/...
of Argentina for the years 2020 and 2021.

Age was divided into four categories: 25 to 44, 45 to 64, 65 to 74, and 75 years and older, whereas educational level included those who never attended education and those who completed or not any of the primary, secondary (including basic general education and polymodal), and higher-university levels. Sex and age were included due to a higher risk of mortality in older persons and in men 1919. Scruzzi GF, Aballay LR, Carreño P, Díaz Rousseau GA, Franchini CG, Cecchetto E, et al. Vacunación contra SARS-CoV-2 y su relación con enfermedad y muerte por COVID-19 en Argentina. Rev Panam Salud Pública. 2023;46:e39. doi: 10.26633/RPSP.2022.39.
https://doi.org/10.26633/RPSP.2022.39...
.

Regarding age, young adults (25 to 64 years) were analyzed separately from older people (65 and older), due to the increase in deaths in the first age group during the second wave. Educational level was grouped into two categories: low and medium-high, as an indicator of socioeconomic level commonly used in mortality studies 2020. Huisman M, Kunst AE, Bopp M, Borgan J-K, Borrell C, Costa G, et al. Educational inequalities in cause-specific mortality in middle-aged and older men and women in eight western European populations. Lancet. 2005;365(9458):493-500. doi: 10.1016/S0140-6736(05)17867-2.
https://doi.org/10.1016/S0140-6736(05)17...
,2121. Brønnum-Hansen H, Baadsgaard M. Increasing social inequality in life expectancy in Denmark. Eur J Public Health. 2007;17:585-6. doi: 10.1093/eurpub/ckm045.
https://doi.org/10.1093/eurpub/ckm045...
. During 2020, with the exception of the department of Angaco in San Juan province (1 death with unknown educational level out of a total of 5 COVID-19 deaths), all remaining departments had 80% or more COVID-19 deaths with data on the educational level of the deceased). During 2021 all departments had 80% or more COVID-19 deaths with this data.

The low educational level category included participants with incomplete secondary education. The medium-high educational level considered people who completed secondary school (6 or more years of elementary school, 5 or more years of secondary school) or had incomplete or complete tertiary education (1 or more years of tertiary education or no tertiary education after secondary education). The departments of the provinces of Mendoza and San Juan constituted the spatial units. Two variables were used to measure the socioeconomic and urbanization level in each department. The percentage of households with unsatisfied basic needs (UBN) was used as an indicator of socioeconomic level, commonly used in Argentina as a measure of structural poverty 2222. Instituto Nacional de Estadística y Censos de la República Argentina. Necesidades básicas insatisfechas 2024. [consultado el 30 de noviembre de 2023]. Disponible en: https://www.indec.gob.ar/indec/web/Nivel4-Tema-4-47-156.
https://www.indec.gob.ar/indec/web/Nivel...
, while population density (individuals per km2) was used as an indicator of the level of urbanization, since the classic definition implies population concentration 2323. Tisdale H. The Process of Urbanization. Soc Forces. 1942;20:311. doi: 10.2307/3005615.
https://doi.org/10.2307/3005615...
. To interpret the results, the population density and UBN values were transformed into z-scores. Data for both variables was obtained from the 2010 Census. Linear projections of these variables were calculated using the 2001 and 2010 censuses because currently there is no data available on population by age structure and educational level at the departmental level during 2020-2021.

Statistical analysis

To test for associations between COVID-19 mortality and sociodemographic and contextual factors, we used a negative binomial Bayesian hierarchical model, with cells at level 1, consisting of deaths in population numerators and denominators (model offset) cross-tabulated by age group, sex and educational level, which were nested within the 37 departments of Mendoza and San Juan at level 2. Each department had data of the population density and UBN households. Models under negative binomial distributions showed a better fit regarding models under Poisson distributions, when considering the Watanabe-Akaike information criterion (Supplementary Material, Table A1). For each year and three large age groups (25 to older, 25 to 64, and 65 to older), a spatial mixed model 2424. Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol. 2013;4:33-49. doi: 10.1016/j.sste.2012.12.001.
https://doi.org/10.1016/j.sste.2012.12.0...
with spatially structured random effects (“BYM2” model) 2525. Riebler A, Sørbye SH, Simpson D, Rue H. An intuitive Bayesian spatial model for disease mapping that accounts for scaling. Stat Methods Med Res. 2016;25:1145-65. doi: 10.1177/0962280216660421.
https://doi.org/10.1177/0962280216660421...
was used to account for spatial dependence between departments (i.e., there is a tendency for similar mortality rates between neighboring departments). Thus, we included a parameter that takes into account spatially structured random effects plus a parameter of unstructured residuals 2626. Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991;43:1-20. doi: 10.1007/BF00116466.
https://doi.org/10.1007/BF00116466...
. Spatially structured random effects were calculated considering a spatial contiguity matrix, where the neighborhood criterion was determined if a department shared a boundary with another department. Relative risks (RR) were estimated as a measure of association between the risk of COVID-19 mortality and the independent variables. Then, hyperparameters, the precision (the inverse of the variance) of the random effects, and a parameter controlling for the importance of spatial structure were estimated 2525. Riebler A, Sørbye SH, Simpson D, Rue H. An intuitive Bayesian spatial model for disease mapping that accounts for scaling. Stat Methods Med Res. 2016;25:1145-65. doi: 10.1177/0962280216660421.
https://doi.org/10.1177/0962280216660421...
. Finally, residual relative risks were calculated after controlling for the independent variables included in the models, and posterior probabilities of risks greater than 1 2424. Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol. 2013;4:33-49. doi: 10.1016/j.sste.2012.12.001.
https://doi.org/10.1016/j.sste.2012.12.0...
. Statistical analysis was performed with the INLA package in the R program version R.3.2.5 (http://www.r-project.org), while the maps were made with the QGIS program version 2.14.3 (https://qgis.org/en/site/).

Ethical aspects

This study used secondary data available under Law 27.275 on the Right of Access to Public Information. Therefore, no approval from a research ethics committee was required since we worked with anonymous statistical data.

FINDINGS

During 2020 and 2021, 2715 and 4844 deaths due to COVID-19 were reported, respectively, in the population aged 25 years and older residing in the provinces of Mendoza and San Juan. We found that of all the deaths, 78% corresponded to Mendoza, with a higher frequency of men (57%), population aged 65 years and older (72%), and individuals with low educational level (73%). The highest values of population density were concentrated in the small departments comprising the capital cities of San Juan and Mendoza (Figure 2). The highest values of percentage of households with UBN were concentrated mainly in the eastern half of the province of San Juan.

Figure 2
Geographic distribution of the residual relative risk of COVID-19 mortality (adjusted for independent variables) and posterior probability of risk > 1, by year and age group and geographic distribution of population density, households with unmet basic needs in the Provinces of Mendoza and San Juan, Argentina, 2020-2021.

Considering all age groups, we detected a higher risk of COVID-19 mortality (95% confidence intervals - 95% CI - greater than 1) associated with increasing age, male sex, and low educational level (RR 2020: 1.81, 95% CI: 1.60-2.06; RR 2021: 1.64, 95% CI: 1.49-1.81; Table 1) during 2020-2021. Relative inequalities in mortality by age group were lower in 2021 compared to 2020. In general, the associations with the variables were maintained in the models that considered separately the 25-64 and 65+ age groups. A high percentage of households with UBN was associated with higher mortality due to COVID-19 (RR: 1.15, 95% CI: 1.02-1.31; Table 1) in population aged 65 and older during 2021.

Table 1
Sociodemographic and contextual variables associated with COVID-19 mortality in the provinces of Mendoza and San Juan, Argentina, 2020-2021.

We found a higher mortality risk mainly in the northern half of Mendoza province during 2020, after considering sex, age, and educational level of people who died due to COVID-19, in addition to population density and the percentage of households with UBN at the departmental level (Figure 2). During 2021, the high mortality risk appeared to be concentrated in the southern half of Mendoza province. When comparing the posterior probabilities of risk greater than 1 between the 25-64 and 65+ age groups, we observed greater geographic similarity in 2020, with higher mortality risks in the northern half of Mendoza. Then, in 2021, there was greater geographic differentiation in the 65 to older age group, with respect to 2020, reflected in a spread of mortality in the southern half of Mendoza province (Figure 2).

DISCUSSION

Educational inequalities related to COVID-19 mortality existed in the provinces of Mendoza and San Juan in Argentina during 2020-2021 regardless of socioeconomic background and urbanization level. The exception was the population aged 65 years and older, which, during 2021, showed a higher risk of death from COVID-19 in departments with high levels of structural poverty.

The risk of COVID-19 mortality was higher mostly in the province of Mendoza during 2020-2021. In contrast to the province of San Juan, a strategy of balance between health and economy was implemented in the province of Mendoza (2727. La Nación. Coronavirus: murió un hombre de 74 años en Mendoza y hay 110 víctimas en el país 2020. [consultado el 21 de agosto de 2020]. Disponible en: https://www.lanacion.com.ar/sociedad/coronavirus-argentina-murio-hombre-74-anos-mendoza-nid2354467.
https://www.lanacion.com.ar/sociedad/cor...
) that could have led to more population circulation in the province, higher number of infections and, therefore, higher number of deaths. This can be evidenced by comparing the variation of mobility in transport stations between the two provinces, which shows a more pronounced drop in mobility during 2020 compared to the pre-pandemic reference period in the province of San Juan 2828. Data Driven Argentina. Reporte de movilidad de Google. Data Driven Argent 2022. [consultado el 1 de noviembre de 2022]. Disponible en: https://datadriven.com.ar/movilidad-google-argentina/.
https://datadriven.com.ar/movilidad-goog...
. While the drop in mobility in transport stations in Mendoza province mostly did not reach 40% (considering 14-day moving averages), the drop in mobility in San Juan province was more pronounced, mostly with values between -60% and -70% 2828. Data Driven Argentina. Reporte de movilidad de Google. Data Driven Argent 2022. [consultado el 1 de noviembre de 2022]. Disponible en: https://datadriven.com.ar/movilidad-google-argentina/.
https://datadriven.com.ar/movilidad-goog...
.

Our results showed similar educational inequalities for COVID-19 mortality between the two waves, after controlling for the level of structural poverty and urbanization. Educational level would not only be an indicator of a person’s socioeconomic situation, in terms of salary, type of labor insertion (precarious-stable), but also of his or her health situation, with worse indicators in populations with low educational level 2929. Rodríguez López S, Bilal U, Ortigoza AF, Diez-Roux AV. Educational inequalities, urbanicity and levels of non-communicable diseases risk factors: evaluating trends in Argentina (2005-2013). BMC Public Health. 2021;21:1-12. doi:10.1186/s12889-021-11617-8.
https://doi.org/10.1186/s12889-021-11617...
. In the State of São Paulo (Brazil) 1313. Li SL, Pereira RH, Prete Jr CA, Zarebski AE, Emanuel L, Alves PJ, et al. Higher risk of death from COVID-19 in low-income and non-White populations of São Paulo, Brazil. BMJ Glob Health. 2021;6(4):e004959. doi: 10.1136/bmjgh-2021-004959.
https://doi.org/10.1136/bmjgh-2021-00495...
, the city of Rome (Italy) 99. Angelici L, Sorge C, Di Martino M, Cappai G, Stafoggia M, Agabiti N, et al. Incidence of SARS-CoV-2 Infection and Related Mortality by Education Level during Three Phases of the 2020 Pandemic: A Population-Based Cohort Study in Rome. J Clin Med. 2022;11(3):877. doi: 10.3390/jcm11030877.
https://doi.org/10.3390/jcm11030877...
, and in the United States 1111. Feldman JM, Bassett MT. Variation in COVID-19 mortality in the US by race and ethnicity and educational attainment. JAMA Netw Open. 2021;4(11):e2135967. doi: 10.1001/jamanetworkopen.2021.35967.
https://doi.org/10.1001/jamanetworkopen....
, higher risk of mortality from COVID-19 were found in populations with low level of education. In the case of São Paulo, populations with lower educational level had higher prevalence of comorbidities and worse working conditions (face-to-face work, less possibility of teleworking, unpaid leave) compared to populations with medium and high educational level 1313. Li SL, Pereira RH, Prete Jr CA, Zarebski AE, Emanuel L, Alves PJ, et al. Higher risk of death from COVID-19 in low-income and non-White populations of São Paulo, Brazil. BMJ Glob Health. 2021;6(4):e004959. doi: 10.1136/bmjgh-2021-004959.
https://doi.org/10.1136/bmjgh-2021-00495...
.

The risk of COVID-19 mortality was not higher in urban areas during 2020, despite restrictions on population mobility that may have protected against SARS-CoV-2 transmission to peripheral areas with low population density. We also did not find an increase in COVID-19 mortality in departments with lower population density (more rural) during 2021. A possible explanation for the absence of differences in mortality between rural and urban areas may be due to the possible absence of geographic inequalities in access to intensive care beds. A study conducted in the city of São Paulo (Brazil) showed that delayed access to hospitalization was not associated with an increased risk of mortality due to COVID-19, after considering the presence of comorbidities 1010. Chiaravalloti Neto F, Bermudi PMM, Aguiar BS de, Failla MA, Barrozo LV, Toporcov TN. Covid-19 hospital mortality using spatial hierarchical models: cohort design with 74,994 registers. Rev Saúde Pública. 2023; (suppl 1):2s. doi: 10.11606/s1518-8787.2023057004652.
https://doi.org/10.11606/s1518-8787.2023...
. Therefore, it is possible that the impossibility of protecting oneself from infection, together with the presence of comorbidities, may be more predominant factors linked to mortality from COVID-19 than access to hospitalization.

Higher levels of structural poverty were associated with higher mortality from COVID-19 during 2021 in the population aged 65 years and older. This could indicate a greater transmission of SARS-CoV-2 to peripheral departments of lower socioeconomic status, particularly in the eastern part of the province of San Juan, in a context of greater population mobility compared to 2020. In parallel to an increasing population mobility since the end of 2020, at the beginning of 2021, people aged 70 years and older started to be vaccinated, who were considered to be a high-priority population group. One possible explanation is the lower proportion of vaccination against COVID-19 in low socioeconomic level departments during the first months of 2021, when cases began to spread in the provinces of Mendoza and San Juan during the second wave.

This study has several limitations. One of the main limitations was that there is no current data on populations by sex, age, educational level, and structural poverty at the departmental level, so we made linear projections using census data from 2001 and 2010. Therefore, the use of data from the 2022 Census, not yet available, could modify our findings. Then, the spatial units used in this study may reflect a broad level of generalization that masks significant socioeconomic variations within spatial units. Another limitation is that mortality data may be affected by problems related to the registration of the cause of death due to COVID-19 (under- and over-registration), although we do not have information on geographic variations regarding this data limitation. Finally, we do not have information on the prevalence of comorbidities at the departmental level, which is why we did not include these variables in the statistical models.

In conclusion, our results would indicate the higher risk of mortality due to COVID-19 in populations with low educational level compared to populations with medium-high educational level. This would indicate the need to focus not only on areas with high levels of poverty, but also on adults with low education levels.

References

  • 1
    Bermudi PMM, Lorenz C, de Aguiar BS, Failla MA, Barrozo LV, Chiaravalloti-Neto F. Spatiotemporal ecological study of COVID-19 mortality in the city of São Paulo, Brazil: shifting of the high mortality risk from areas with the best to those with the worst socio-economic conditions. Travel Med Infect Dis. 2021;39:101945. doi: 10.1016/j.tmaid.2020.101945.
    » https://doi.org/10.1016/j.tmaid.2020.101945
  • 2
    Mena G, Martinez PP, Mahmud AS, Marquet PA, Buckee CO, Santillana M. Socioeconomic status determines COVID-19 incidence and related mortality in Santiago, Chile. Science. 2021;372(6545):eabg5298. doi: 10.1101/2021.01.12.21249682.
    » https://doi.org/10.1101/2021.01.12.21249682
  • 3
    Silva J, Ribeiro-Alves M. Social inequalities and the pandemic of COVID-19: the case of Rio de Janeiro. J Epidemiol Community Health. 2021;75(10):975-979. doi: 10.1136/jech-2020-214724.
    » https://doi.org/10.1136/jech-2020-214724
  • 4
    Leveau CM, Soares Bastos L. Desigualdades socio-espaciales de la mortalidad por COVID-19 en tres olas de propagación: un análisis intra-urbano en Argentina. Cad Saúde Pública. 2022;38(5):e00163921. doi: 10.1590/0102-311XES163921.
    » https://doi.org/10.1590/0102-311XES163921
  • 5
    Albuquerque MV de, Ribeiro LHL. Desigualdade, situação geográfica e sentidos da ação na pandemia da COVID-19 no Brasil. Cad Saúde Pública. 2021;36(12): e00208720. doi: 10.1590/0102-311X00208720.
    » https://doi.org/10.1590/0102-311X00208720
  • 6
    Anzalone AJ, Horswell R, Hendricks BM, Chu S, Hillegass WB, Beasley WH, et al. Higher hospitalization and mortality rates among SARS-CoV-2-infected persons in rural America. J Rural Health. 2023;39:39-54. doi: 10.1111/jrh.12689.
    » https://doi.org/10.1111/jrh.12689
  • 7
    Cuadros DF, Branscum AJ, Mukandavire Z, Miller FD, MacKinnon N. Dynamics of the COVID-19 epidemic in urban and rural areas in the United States. Ann Epidemiol. 2021;59:16-20. doi: 10.1016/j.annepidem.2021.04.007.
    » https://doi.org/10.1016/j.annepidem.2021.04.007
  • 8
    Petrelli A, Ventura M, Di Napoli A, Mateo-Urdiales A, Pezzotti P, Fabiani M. Geographic heterogeneity of the epidemiological impact of the COVID-19 pandemic in Italy using a socioeconomic proxy-based classification of the national territory. Front Public Health. 2023;11:1143189. doi: 10.3389/fpubh.2023.1143189.
    » https://doi.org/10.3389/fpubh.2023.1143189
  • 9
    Angelici L, Sorge C, Di Martino M, Cappai G, Stafoggia M, Agabiti N, et al. Incidence of SARS-CoV-2 Infection and Related Mortality by Education Level during Three Phases of the 2020 Pandemic: A Population-Based Cohort Study in Rome. J Clin Med. 2022;11(3):877. doi: 10.3390/jcm11030877.
    » https://doi.org/10.3390/jcm11030877
  • 10
    Chiaravalloti Neto F, Bermudi PMM, Aguiar BS de, Failla MA, Barrozo LV, Toporcov TN. Covid-19 hospital mortality using spatial hierarchical models: cohort design with 74,994 registers. Rev Saúde Pública. 2023; (suppl 1):2s. doi: 10.11606/s1518-8787.2023057004652.
    » https://doi.org/10.11606/s1518-8787.2023057004652
  • 11
    Feldman JM, Bassett MT. Variation in COVID-19 mortality in the US by race and ethnicity and educational attainment. JAMA Netw Open. 2021;4(11):e2135967. doi: 10.1001/jamanetworkopen.2021.35967.
    » https://doi.org/10.1001/jamanetworkopen.2021.35967
  • 12
    Spijker JJA, Trias-Llimós S. Cause-specific mortality in Spain during the pandemic. Educational differences and its impact on life expectancy. Eur J Public Health 2023;33(3):543-9. doi: 10.1093/eurpub/ckad036.
    » https://doi.org/10.1093/eurpub/ckad036
  • 13
    Li SL, Pereira RH, Prete Jr CA, Zarebski AE, Emanuel L, Alves PJ, et al. Higher risk of death from COVID-19 in low-income and non-White populations of São Paulo, Brazil. BMJ Glob Health. 2021;6(4):e004959. doi: 10.1136/bmjgh-2021-004959.
    » https://doi.org/10.1136/bmjgh-2021-004959
  • 14
    Leveau CM, Hussein M, Tapia-Granados JA, Velázquez GA. Economic fluctuations and educational inequalities in premature ischemic heart disease mortality in Argentina. Cad Saúde Pública. 2023;39(5):e00181222. doi: 10.1590/0102-311xen181222.
    » https://doi.org/10.1590/0102-311xen181222
  • 15
    Boletín Oficial de la República Argentina. Boletín Oficial República Argentina - Aislamiento Social Preventivo Y Obligatorio - Decreto 297/2020. 2020 [consultado el 21 de agosto de 2022]. Disponible en: https://www.boletinoficial.gob.ar/detalleAviso/primera/227042
    » https://www.boletinoficial.gob.ar/detalleAviso/primera/227042
  • 16
    Freitas ARR, Beckedorff OA, Cavalcanti LP de G, Siqueira AM, Castro DB de, Costa CF da, et al. The emergence of novel SARS-CoV-2 variant P.1 in Amazonas (Brazil) was temporally associated with a change in the age and sex profile of COVID-19 mortality: A population based ecological study. Lancet Reg Health Am. 2021;1:100021. doi: 10.1016/j.lana.2021.100021.
    » https://doi.org/10.1016/j.lana.2021.100021
  • 17
    Instituto Nacional de Estadística y Censos. Datos provisionales del CENSO 2022. Censo Nac Poblac Hogares Viviendas 2023. [consultado el 31 de mayo de 2023]. Disponible en: https://censo.gob.ar/index.php/datos_provisionales/
    » https://censo.gob.ar/index.php/datos_provisionales/
  • 18
    Ministerio de Salud de la Nación. Datos Abiertos del Ministerio de Salud - Defunciones ocurridas y registradas en la República Argentina 2023. [consultado el 16 de diciembre de 2023]. Disponible en: http://datos.salud.gob.ar/
    » http://datos.salud.gob.ar/
  • 19
    Scruzzi GF, Aballay LR, Carreño P, Díaz Rousseau GA, Franchini CG, Cecchetto E, et al. Vacunación contra SARS-CoV-2 y su relación con enfermedad y muerte por COVID-19 en Argentina. Rev Panam Salud Pública. 2023;46:e39. doi: 10.26633/RPSP.2022.39.
    » https://doi.org/10.26633/RPSP.2022.39
  • 20
    Huisman M, Kunst AE, Bopp M, Borgan J-K, Borrell C, Costa G, et al. Educational inequalities in cause-specific mortality in middle-aged and older men and women in eight western European populations. Lancet. 2005;365(9458):493-500. doi: 10.1016/S0140-6736(05)17867-2.
    » https://doi.org/10.1016/S0140-6736(05)17867-2
  • 21
    Brønnum-Hansen H, Baadsgaard M. Increasing social inequality in life expectancy in Denmark. Eur J Public Health. 2007;17:585-6. doi: 10.1093/eurpub/ckm045.
    » https://doi.org/10.1093/eurpub/ckm045
  • 22
    Instituto Nacional de Estadística y Censos de la República Argentina. Necesidades básicas insatisfechas 2024. [consultado el 30 de noviembre de 2023]. Disponible en: https://www.indec.gob.ar/indec/web/Nivel4-Tema-4-47-156
    » https://www.indec.gob.ar/indec/web/Nivel4-Tema-4-47-156
  • 23
    Tisdale H. The Process of Urbanization. Soc Forces. 1942;20:311. doi: 10.2307/3005615.
    » https://doi.org/10.2307/3005615
  • 24
    Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol. 2013;4:33-49. doi: 10.1016/j.sste.2012.12.001.
    » https://doi.org/10.1016/j.sste.2012.12.001
  • 25
    Riebler A, Sørbye SH, Simpson D, Rue H. An intuitive Bayesian spatial model for disease mapping that accounts for scaling. Stat Methods Med Res. 2016;25:1145-65. doi: 10.1177/0962280216660421.
    » https://doi.org/10.1177/0962280216660421
  • 26
    Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991;43:1-20. doi: 10.1007/BF00116466.
    » https://doi.org/10.1007/BF00116466
  • 27
    La Nación. Coronavirus: murió un hombre de 74 años en Mendoza y hay 110 víctimas en el país 2020. [consultado el 21 de agosto de 2020]. Disponible en: https://www.lanacion.com.ar/sociedad/coronavirus-argentina-murio-hombre-74-anos-mendoza-nid2354467
    » https://www.lanacion.com.ar/sociedad/coronavirus-argentina-murio-hombre-74-anos-mendoza-nid2354467
  • 28
    Data Driven Argentina. Reporte de movilidad de Google. Data Driven Argent 2022. [consultado el 1 de noviembre de 2022]. Disponible en: https://datadriven.com.ar/movilidad-google-argentina/
    » https://datadriven.com.ar/movilidad-google-argentina/
  • 29
    Rodríguez López S, Bilal U, Ortigoza AF, Diez-Roux AV. Educational inequalities, urbanicity and levels of non-communicable diseases risk factors: evaluating trends in Argentina (2005-2013). BMC Public Health. 2021;21:1-12. doi:10.1186/s12889-021-11617-8.
    » https://doi.org/10.1186/s12889-021-11617-8

  • Funding.

    CML was funded by the Agencia Nacional de Promoción de la Investigación, el Desarrollo Tecnológico y la Innovación, PICT 2021-I-INVI-00683.

  • Cite as:

    Leveau CM, Velázquez GA. COVID-19 mortality: educational inequalities and socio-spatial context in two provinces of Argentina. Rev Peru Med Exp Salud Publica. 2024;41(2): doi: 10.17843/rpmesp.2024.412.13201.

Publication Dates

  • Publication in this collection
    19 Aug 2024
  • Date of issue
    Apr-Jun 2024

History

  • Received
    07 Jan 2024
  • Accepted
    27 Mar 2024
Instituto Nacional de Salud Lima - Lima - Peru
E-mail: revmedex@ins.gob.pe