Next Article in Journal
Validation of the Spanish Version of the Double Knowledge Expectations and Received Knowledge Significant Other Scale for Informal Caregivers of People with Dementia (KESO-DEM/RKSO-DEM)
Next Article in Special Issue
Mitigation Planning and Policies Informed by COVID-19 Modeling: A Framework and Case Study of the State of Hawaii
Previous Article in Journal
Redrawing Cities with Children and Adolescents: Development of a Framework and Opportunity Index for Wellbeing—The REDibuja Study Protocol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Human Development Index Is Associated with COVID-19 Case Fatality Rate in Brazil: An Ecological Study

by
Camila Vantini Capasso Palamim
1,2,3,
Matheus Negri Boschiero
1,2,
Felipe Eduardo Valencise
1,2 and
Fernando Augusto Lima Marson
1,2,3,*
1
Laboratory of Cell and Molecular Tumor Biology and Bioactive Compounds, São Francisco University, Avenida São Francisco de Assis, 218, Jardim São José, Bragança Paulista 12916-900, SP, Brazil
2
Laboratory of Human and Medical Genetics, São Francisco University, Avenida São Francisco de Assis, 218, Jardim São José, Bragança Paulista 12916-900, SP, Brazil
3
Postgraduate Program in Health Science, São Francisco University, Avenida São Francisco de Assis, 218, Jardim São José, Bragança Paulista 12916-900, SP, Brazil
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(9), 5306; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph19095306
Submission received: 26 March 2022 / Revised: 22 April 2022 / Accepted: 24 April 2022 / Published: 27 April 2022

Abstract

:
The Human Development Index measures a region’s development and is a step for development debate beyond the traditional, economic perspective. It can also determine the success of a country’s response to the COVID-19 pandemic, mainly affecting the case fatality rate among severe cases of SARS-CoV-2 infection. We aimed to associate the Human Development Index with the case fatality rate due to COVID-19 in each Brazilian state and the Federal District, taking into account comorbidities and the need for invasive mechanical ventilation. We also evaluated the influence of the GINI index, number of intensive care unit beds, and occupied households in subnormal clusters on the case fatality rate. We performed an ecological study including two populations: COVID-19 individuals that did not require the mechanical ventilation protocol; and COVID-19 individuals under invasive mechanical ventilation. We performed a Pearson correlation test and a univariate linear regression analysis on the relationship between Human Development Index, Human Development Index—Education Level, Human Development Index—Life Expectancy, and Human Development Index—Gross National Income per capita and COVID-19 deaths. The same analyses were performed using the other markers. We grouped the patients with COVID-19 according to comorbidities and the need for invasive mechanical ventilation. Alpha = 0.05. We included 848,501 COVID-19 individuals, out of which 153,710 needed invasive mechanical ventilation and 314,164 died, and 280,533 COVID-19 individuals without comorbidity, out of which 33,312 needed invasive mechanical ventilation and 73,723 died. We observed a low negative Pearson correlation between the Human Development Index and death and a moderate negative Pearson correlation between the Human Development Index and deaths of individuals on invasive mechanical ventilation, with or without comorbidity. The univariate linear analysis showed the case fatality rate depends on at least 20–40% of the Human Development Index. In Brazil, regions with a low Human Development Index demonstrated a higher case fatality rate due to COVID-19, mainly in individuals who needed invasive mechanical ventilation, than regions with a higher Human Development Index. Although other indexes studied, such as intensive care unit beds and GINI, were also associated with the COVID-19 case fatality rate, they were not as relevant as the Human Development Index. Brazil is a vast territory comprising cultural, social, and economic diversity, which mirrors the diversity of the Human Development Index. Brazil is a model nation for the study of the Human Development Index’s influence on aspects of the COVID-19 pandemic, such as its impact on the case fatality rate.

1. Introduction

The coronavirus disease (COVID)-19 can spread through person-to-person contact and through the air, affecting people throughout whole countries. Individuals with comorbidities, such as diabetes mellitus, systemic arterial hypertension, and obesity, are more likely to develop a more severe course and progression of the disease [1,2,3]. In addition, factors such as social distancing compliance, personal hygiene methods, and governmental response to the pandemic seem to attenuate the impact of the COVID-19 [4,5]. In this context, markers such as the Human Development Index (HDI) can be valuable for measuring a region’s development. It can be considered a step for development debate beyond the traditional, economic perspective and can also determine the success of a country’s response to the COVID-19 pandemic, mainly affecting the case fatality rate among severe cases of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection [6,7].
Although one of the most important sociodemographic indexes, the HDI, solely, cannot deeply analyze a country’s response to COVID-19 due to its simplicity. Thus, jointly analyzing other socio, economic, and demographic indexes, such as the GINI index and the number of intensive care units (ICU) beds, which were already linked to COVID-19 in previous studies [8,9], is essential. Living conditions can also play a key role as a risk factor for COVID-19 since poor infrastructure and overcrowding, such as that observed in slums, can enhance SARS-CoV-2 dissemination [10,11]. Importantly, these indicators can help one to understand the dynamics of the COVID-19 pandemic and allow one to compare the general needs of certain areas, thus, allowing one to make decisions based on reliable indicators [12,13,14].
The HDI can be a valuable tool for better understanding the impact of the COVID-19 pandemic in certain areas and can also guide better public health measures since the factors used to calculate the HDI (life expectancy, education, and gross national income per capita) are associated with COVID-19 mortality and infection [15]. For instance, regions with high life expectancy perhaps suffer a higher impact from COVID-19 since older age is associated with the worst outcomes in COVID-19 [16]. Additionally, the higher the level of education, the more individuals are aware of the pandemic risks of not adhering to preventive measures, such as the use of masks or social distancing, thus, decreasing the spread of SARS-CoV-2 [15,17]. Finally, regions with higher gross national income per capita are perhaps more likely to have international travelers and can purchase more COVID-19 tests, which can enhance SARS-CoV-2 dissemination and decrease underreporting, respectively [15].
HDI scores can range from 0 (lower development; e.g., Niger (0.394), Central African Republic (0.397), and Chad (0.398) in 2019) to 1 (higher development; e.g., Norway (0.957), Ireland (0.955), and Switzerland (0.955) in 2019) [18,19]. Brazil, in 2020, was at position 84 worldwide, with an HDI of 0.765, behind other South American countries such as Argentina (0.845), Uruguay (0.817), and Chile (0.851) [19].
The literature described a correlation between SARS-CoV-2 infection and death prevalence due to COVID-19 and the HDI, as well as a country’s socioeconomic position [20,21,22], including two Brazilian studies performed in the Alagoas and Pernambuco states. These studies demonstrated higher mortality in regions with a low HDI score [23,24]. However, the data from these studies are conflicting. Regarding Brazil, there is a wide range of HDI scores among the Brazilian states and the Federal District (0.850 in the Federal District to 0.683 in the Alagoas state) [25]. To date, the main findings for the association between HDI and the case fatality rate in Brazil present only local perspectives alongside a scarce representation of Brazil’s national scenario.
Brazil’s vast territory explains the wide disparity among states and the Federal District in HDI and other markers, such as the GINI index and access to health support, which can compromise access to complex treatments, including mechanical ventilatory support. The Brazilian diversity can be a model to determine the influence of a country’s development on its ability to deal with a pandemic situation, such as COVID-19. In this context, we aimed to associate socio, economic, and demographic indexes, such as the HDI, ICU beds, the GINI index, and occupied households in subnormal clusters, with the case fatality rate due to COVID-19 in each Brazilian state and the Federal District taking into account comorbidities and the need for invasive mechanical ventilation (IMV).

2. Materials and Methods

We performed an ecological study using epidemiologic data (death due to COVID-19 and death due to COVID-19 among individuals who needed IMV) available at OpenDataSUS [26]. We computed the data from the Brazilian Ministry of Health according to the surveillance data of severe acute respiratory infection and data from the Information System platform for Epidemiological Surveillance of Influenza (SIVEP-Flu). The data were recorded during the first year of COVID-19 in Brazil after the first report in our country. In such a context, we retrieved the information for severe acute respiratory infection from December 2019 to April 2020 and for patients with COVID-19 (SARS-CoV-2 infection) from February 2020 to April 2021. Only four days of April’s first week were included in our study. Inclusion criteria: patients with a positive SARS-CoV-2 real-time polymerase chain reaction (RT-PCR) test and complete information on the need for ventilatory support, outcomes, and place of residence. Exclusion criteria: negative SARS-CoV-2 RT-PCR test or absence of classification of severe acute respiratory infection, absence of a description for the place of residence, or patients who lived in a country other than Brazil. We also excluded the patients who presented comorbidities or pregnancy in a second analysis. In addition, we excluded the individuals without gender data or outcomes information for both analyses. The complete flowchart of the patients included and excluded is shown in Figure 1.
We present the data as percentages calculated using the ratio between individuals who died due to COVID-19 and the sum of COVID-19 individuals who died or recovered in Brazilian states and the Federal District. Two analyses were performed for: (i) all COVID-19 individuals without mechanical ventilation protocol; and (ii) COVID-19 individuals who needed IMV. Additionally, for subgroups analysis, we considered the presence of confounders, such as comorbidities.
The HDI has three principles to classify a region as developed: life expectancy (HDI-LE); capacity to acquire knowledge, that is, mean and expected years of schooling (HDI-E); and access to resources for a decent standard of living, that is, gross national income per capita (HDI-GNI) [6,7,15,27]. We retrieved the HDI data of each Brazilian state and the Federal District from the AtlasBR [25], a private company licensed by the federal government. Unfortunately, the latest HDI score for each Brazilian state and the Federal District was from 2017. The Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) published the latest HDI in 2010 [28].
The GINI index and the number of occupied households in subnormal clusters were obtained from the IBGE website [29]. The GINI index measures the income inequality of a certain area. It varies from 0 to 1, with regions with values close to zero presenting lower inequality in contrast to regions with values close to 1, which present higher inequality [30]. Regarding the number of occupied households in subnormal clusters (slums, or favelas in Portuguese), the IBGE defines such clusters as “forms of irregular occupation of land owned by others for housing purpose, characterized by an irregular urban pattern, lack of essential public services and location in that have restrictions on occupancy” [29]. We retrieved the data from 2019 since it was the last update before the COVID-19 pandemic. We decided to present the absolute number of occupied households in subnormal clusters and the relative number (%) of occupied households in subnormal clusters, representing the proportion of this type of occupation and the total number of households [29].
Regarding ICU beds, the following markers were included in our study: the number of ICU beds in Brazil, ICU beds in the Brazilian public health system (SUS; in Portuguese, Sistema Único de Saúde), ICU beds in the Brazilian private health system, ICU beds per 10,000 inhabitants, ICU beds per 10,000 inhabitants in the Brazilian public health system, and ICU beds per 10,000 inhabitants in the Brazilian private health system (beneficiaries only). We retrieved the data from the Federal Council of Medicine, 2018 [9,31].
We performed a Pearson correlation test between HDI, HDI-E, HDI-LE, and HDI-GNI and the number of deaths due to COVID-19, considering the comorbidities and the need for IMV. The same statistical test was used to correlate the GINI index, the number/percentage of occupied households in subnormal clusters, and the number of ICU beds per 10,000 inhabitants (total, public health system, and private health system) with the number of deaths due to COVID-19, considering comorbidities and the need for IMV. We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; and negligible correlation, 0.00 to 0.30. We present the correlation coefficient (CC) and the 95% confidence interval (95%CI) for the Pearson correlation test.
We performed a univariate regression analysis between HDI, GINI index, the number/percentage of occupied households in subnormal clusters, and the number of ICU beds per 10,000 inhabitants (total, public health system, and private health system) and deaths due to COVID-19. The analysis considered the presence of comorbidities and the need for IMV. For the univariate regression analysis, we described the R-squared (goodness-of-fit measure for linear regression models that indicates the percentage of the variance in the dependent variable that the independent variables explain collectively). This equation represents how an independent variable X (HDI, GINI index, the number/percentage of occupied households in subnormal clusters, and the number of ICU beds per 10,000 inhabitants (total, public health system, and private health system)) is related to a dependent variable Y (case fatality rate) and the 95% confidence bands of the best-fit line. In brief, the primary marker evaluated in our study was the HDI and the others were the GINI index, the number/percentage of occupied households in subnormal clusters, and the number of ICU beds per 10,000 inhabitants (total, public health system, and private health system).
We performed the statistical analysis using the Statistical Package for the Social Sciences software (IBM SPSS Statistics for Macintosh, Version 27.0) and the GraphPad Prism version 8.00 for Apple Mac, GraphPad Software, San Diego, CA, USA, www.graphpad.com (accessed on 21 March 2022). We used an alpha of 0.05 in all statistical analyses.
The data used in our study were publicly available. As the data did not contain personal data of patients, the study was consent-free since it does not present risks to the research participants.

3. Results

3.1. Inclusion of the COVID-19 Individuals and the HDI of the Brazilian States and the Federal District

We included 848,501 individuals with positive RT-PCR tests for SARS-CoV-2 in the first analysis from the original cohort. Additionally, we performed a second analysis in which we enrolled individuals with SARS-CoV-2 infection and with no comorbidities screened; the second analysis comprised 280,533 individuals. For the first and second analyses, respectively, 153,710 and 33,312 individuals were under IMV. We present the description of inclusion and exclusion for the participants and the distribution of ventilatory support and outcome (case fatality rate and clinical recovery) in Figure 1.
We observed the highest HDI in the Federal District (0.850), followed by São Paulo (0.826) and Santa Catarina (0.808) states. In contrast, the lowest HDI occurred in Alagoas (0.683), Maranhão (0.687), and Piauí (0.697) states. Regarding the components of the HDI, we observed the highest HDI-E in São Paulo (0.828) state, followed by the Federal District (0.804) and Santa Catarina (0.779) state. The highest HDI-LE occurred in the Federal District (0.859), followed by São Paulo (0.796) and the Rio Grande do Sul (0.787) states. Finally, we observed the highest HDI-GNI in the Federal District (0.890), followed by Minas Gerais (0.875) and Santa Catarina (0.866) states (Figure 2 and Figure 3).
The Brazilian state which accounted for highest overall case fatality rate was Sergipe (3643/5806; 62.7%) state, followed by Espírito Santo (4618/7668; 60.2%) and Roraima (1121/1957; 57.3%) states, whereas the state with highest death rate in individuals without comorbidities was Sergipe (1010/1847; 54.7%), followed by Espírito Santo (1146/2191; 52.3%) and Rondônia (1475/3177; 46.4%) states (Table 1; Figure 2). Regarding IMV, the states which presented a higher case fatality rate among all individuals were Roraima (702/747; 94.0%), Pará (3990/4357; 91.5%), and Paraíba (2950/3233; 91.2%) states. However, a higher difference between the case fatality rate of all individuals and individuals who needed IMV occurred in Mato Grosso (63.8%), Piauí (52.3%), and the Federal District (52.1%). Regarding the individuals who needed IMV and did not have a comorbidity, the highest case fatality rate occurred in Roraima (226/248; 91.1%), Alagoas (394/436; 90.3%), and Paraíba (585/661; 88.5%) states. In contrast, the higher difference between the case fatality rate of individuals without comorbidities and individuals without comorbidities who needed IMV occurred in Mato Grosso (64.3%), Alagoas (60.3%), and Piauí (58.9%) states (Table 1, Figure 3, and Supplementary Material S2).

3.2. Association between Case Fatality Rate Due to COVID-19 and the HDI

We observed a negative correlation (p-value < 0.05) between HDI and case fatality rate (CC = −0.449; 95% CI = −0.708 to −0.083) and deaths without comorbidities (CC = −0.467; 95% CI = −0.720 to −0.106), between HDI-GNI and case fatality rate (CC = −0.397; 95% CI = −0.675 to −0.020) and deaths without comorbidities (CC = −0.446; 95% CI = −0.706 to −0.079), between HDI-E and case fatality rate (CC = −0.390; 95% CI = −0.670 to −0.011), and, finally, between HDI-LE and case fatality rate (CC = −0.444; 95% CI = −0.705 to −0.077) and deaths without comorbidities (CC = −0.481; 95% CI = −0.728 to −0.123) (Figure 2).

3.3. Association between Case Fatality Rate Due to COVID-19 in Individuals Who Needed IMV and the HDI

We observed a negative correlation (p-value < 0.05) between HDI and case fatality rate (CC = −0.629; 95% CI = −0.815 to −0.328) and deaths without comorbidities (CC = −0.619; 95% CI = −0.809 to −0.313), between HDI-GNI and case fatality rate (CC = −0.618; 95% CI = −0.808 to −0.311) and deaths without comorbidities (CC = −0.654; 95% CI = −0.828 to −0.365), between HDI-E and case fatality rate (CC = −0.569; 95% CI = −0.780 to −0.241) and deaths without comorbidities (CC = −0.509; 95% CI = −0.745 to −0.160), and, finally, between HDI-LE and case fatality rate (CC = −0.569; 95% CI = −0.780 to −0.241) and deaths without comorbidities (CC = −0.582; 95% CI = −0.788 to −0.260) (Figure 3). In brief, all these correlation indexes were moderate.

3.4. Univariate Regression Analysis between the Case Fatality Rate Due to COVID-19 and the HDI

We included patients with COVID-19, considering comorbidity and the need for IMV (Figure 4). We described the best fit R2 for case fatality rate in individuals who needed IMV (R2 = 0.396), followed by individuals who needed IMV and had no comorbidities (R2 = 0.384); overall, for individuals without comorbidities R2 = 0.218 and, overall, for COVID-19 individuals R2 = 0.202. In addition, the equations for each analysis were: overall case fatality rate, case fatality rate = −103.7 × (HDI) + 119.2 (Figure 4A); overall case fatality rate in COVID-19 individuals without comorbidities, case fatality rate = −130.4 × (HDI) + 128.9 (Figure 4B); case fatality rate in COVID-19 individuals who received IMV, case fatality rate = −69.94 × (HDI) + 137.2 (Figure 4C); and case fatality rate in COVID-19 individuals who needed IMV and did not have any comorbidities, case fatality rate = −120.1 × (HDI) + 168.1 (Figure 4D). The case fatality rate was the dependent marker, and the HDI was the independent marker.

3.5. Association between the Case Fatality Rate Due to COVID-19 and Other Populational Features

In our data, we described three other populational features (GINI coefficient, ICU bed distribution, and the number of occupied households in subnormal clusters) which contribute to the case fatality rate diversity in Brazil. We present the features in Table 2, Table S1, and Supplementary Material S2.
We observed the highest GINI index of household income per capita, at average prices for the year, in the Federal District (0.583), followed by Pernambuco (0.578) and Acre (0.575) states, which represented the states with the most significant inequalities. In contrast, the lowest GINI indexes of household income per capita, at average prices for the year, occurred in Santa Catarina (0.429), Mato Grosso (0.457), and Góias (0.474) states (Table 2). Regarding the component of the GINI index of the average real monthly income of people aged 14 and over received in the reference month, for all jobs, at average prices for the year, we observed the highest index in Piauí (0.592) state, followed by Sergipe (0.560) state and the Federal District (0.553). Additionally, the lowest GINI indexes occurred in Santa Catarina (0.417), Mato Grosso (0.440), and Goiás (0.441) states (Table 2).
The Brazilian state which accounted for the highest number of total ICU beds per 10,000 inhabitants was Rio de Janeiro (3.79) state, followed by the Federal District (3.39) and São Paulo (2.72) state. In contrast, the state with the lowest number of total ICU beds per 10,000 inhabitants was Acre (0.90) state, followed by Roraima (0.92) and Amapá (1.03) states (Table 2). Additionally, the Brazilian state which accounted for the highest number of total ICU beds per 10,000 inhabitants in the public health system was Paraná (1.54) state, followed by the Rio Grande do Sul (1.33) and Minas Gerais (1.30) states, whereas the state with the lowest number of total ICU beds per 10,000 inhabitants in the public health system was Amapá (0.33) state, followed by Piauí (0.56) and Roraima (0.57) states (Table 2). Finally, the Brazilian state which accounted for the highest number of total ICU beds per 10,000 inhabitants in the private health system was Mato Grosso (10.63) state, followed by the Federal District (8.78) and Rio de Janeiro (8.70) state. In contrast, the state with the lowest number of total ICU beds per 10,000 inhabitants in the private health system was Santa Catarina (2.61) state, followed by the Minas Gerais (3.14) and the Rio Grande do Sul (3.31) states (Table 2).
Additionally, the Brazilian state which accounted for the highest number of occupied households in subnormal clusters (%) was Amazonas (34.59%) state, followed by the Espiríto Santo (26.10%) and Amapá (21.58%) states. In contrast, the state with the lowest number of occupied households in subnormal clusters (%) was Mato Grosso do Sul (0.74) state, followed by Santa Catarina (1.42) and Goiás (1.55) states (Table 2).

3.5.1. Association between the Case Fatality Rate Due to COVID-19 and GINI Coefficient

We observed a positive correlation (p-value = 0.034) between the GINI of household income per capita, at average prices for the year, and the case fatality rate in patients without comorbidities who received IMV (CC = 0.409; 95% CI = 0.035 to 0.683) (Figure 5). The same index (GINI coefficient of household income per capita, at average prices for the year) had the best fit R2 for case fatality rate in individuals without comorbidities who needed IMV (R2 = 0.168) with the following equation: case fatality rate = 84.31 × (GINI) + 33.67 (Figure 6A). The other study groups did not present any significative correlation (Figure 5) or equation in the regression analysis (Figures S1 and S2).

3.5.2. Association between the Case Fatality Rate Due to COVID-19 and ICU

We observed a negative correlation (p-value = 0.040) between the number of ICU beds in the public health system per 10,000 inhabitants and the case fatality rate in patients who received IMV (CC = −0.397; 95% CI = −0.675 to −0.020) (Figure 7). No significant correlation occurred between the case fatality rate (all patients, patients who received IMV, all patients without comorbidities, and patients without comorbidities who received IMV) and total number of ICU beds per 10,000 inhabitants and the ICU beds in the private health system per 10,000 inhabitants (Figure 7). In addition, we observed a significative best fit R2 for case fatality rate in individuals who needed IMV (R2 = 0.158) and the distribution of the ICU beds in the public health system per 10,000 inhabitants with the following equation: case fatality rate = −7.00 × (ICU beds in the public health system per 10,000 inhabitants) + 91.14 (Figure 6B). The other study groups did not present any significative equation in the regression analysis (Figures S3–S5).

3.5.3. Association between the Case Fatality Rate Due to COVID-19 and Occupied Households in Subnormal Clusters

The number of occupied households in subnormal clusters and the percentage of occupied households in subnormal clusters did not present any significative correlation (Figure S6) with the case fatality rate due to COVID-19 or any significative equation in the regression analysis (Figures S7 and S8) that could predict death due to COVID-19.

4. Discussion

The COVID-19 pandemic became the most urgent public health crisis worldwide [4] as it caused thousands of deaths daily. In Brazil, COVID-19 was responsible for more than 660,000 deaths and millions of confirmed cases and is increasing daily, despite vaccinations [32]. Due to its large dimensions, 8,510,345.538 km2, climate differences, and governmental public health management differences [33], the HDI is significantly different for each Brazilian state and the Federal District. The HDI index varies between Brazilian states and the Federal District, ranging from 0.683 for the Alagoas state to 0.850 for the Federal District and 0.826 for the São Paulo state [25]. Indirectly, the index is associated with different levels of impact from the disease in Brazil, as shown in our data for case fatality rate due to COVID-19, mainly in cases who needed IMV and required hospitalization support, which is directly associated with the development of each state and the Federal District.
Although we studied the COVID-19 pandemic from February 2020 to April 2021, reports showed at least two different waves during this period in Brazil. To face the pandemic and its first wave, the Brazilian federal government transferred around 30% more budget in the first quarters of 2020 to the states compared to in 2019; however, only 8% was spent [34]. Even though Brazil was able to increase the number of ICU beds per 10,000 and the number of health care workers in 2020 compared to 2019, it did not handle the COVID-19 pandemic well, resulting in some of the greatest numbers of confirmed cases and deaths in the world [34,35]. Studies showed the second wave was more severe than the first one, with more hospitalizations and deaths [36,37]. Furthermore, patients affected during the second wave appeared to suffer more severe symptoms since more younger individuals who needed invasive ventilatory support died and were more hypoxemic, and more invasive ventilatory support was also required [36,37]. Several factors might have contributed to these differences. For instance, during the first wave, the most prevalent strains of SARS-CoV-2 were B.1.1.28 and B.1.1.3. In contrast, the second wave included the P.1 and P.2 strains (or Gamma strains) first described in the city of Manaus in Amazonas. It coincided with an increase in hospitalizations [37,38].
The Gamma strains might be more transmissible than non-Gamma strains and more deadly, especially to younger people, with a higher need for hospitalizations, ICU, and mechanical ventilation [39,40,41]. In the case of Manaus city, where the health care system collapsed in the second wave due to a lack of oxygen cylinders, there were increased numbers of hospitalizations and deaths and several allegations of corruption [33], showing that Brazil was not prepared for the second wave and that the first wave drained most of the resources available.
In the study period, we observed a correlation between case fatality rate and HDI. Brazilian states and the Federal District with a lower HDI score were likely to present a higher case fatality rate in individuals with and without comorbidities and among individuals who needed IMV, compared to individuals in states and the Federal District with a higher HDI score, which is similar to the findings in previous studies [20,21,42]. The present study follows a new Brazilian study that evaluated 203 cities in the São Paulo state, encompassing nearly 93% of the state’s population, in which HDI-GNI and HDI-LE were associated with COVID-19 mortality [43].
In addition, an Italian study encompassing 20 Italian regions in an ecological study and a study containing 189 countries reported a positive correlation between HDI and the COVID-19 death rate [22,44]. Previous Brazilian studies reported the HDI to be significantly associated with cumulative COVID-19 cases and faster dissemination of the virus; that is, the higher the HDI, the more cases were reported and in a shorter period [8,45]. Interestingly, one Brazilian study, which evaluated the cities from the Ceará state, also observed a positive correlation between HDI and the incidence of COVID-19, which is maybe related to worse sanitary conditions and the fact that Fortaleza, the capital of Ceará state, is a hub that attracts tourism, mainly from Europe, which could, ultimately, have enhanced the dissemination of SARS-CoV-2 [14,46].
Additionally, the associated results from the more effective health care systems, responsible for identifying early, asymptomatic, and subclinical cases, stated that individuals who live in high-HDI regions are more prone to have chronic diseases, which enhanced the COVID-19 mortality rate, dissemination rate, and the number of total cases [22,44]. In addition, areas with enhanced HDI were associated with more vaccines, ICU beds, ventilators, physicians, and nurses [47], thus, attenuating the impact of the COVID-19 case fatality rate. Unfortunately, regions with low HDI scores tend to have worse access to the health system and diagnostic tools for COVID-19, resulting in underreporting of COVID-19 cases and deaths, such as in Brazil, especially in regions where neglected peoples (e.g., indigenous peoples) live [44,48,49,50,51]. Countries with low HDI scores also tend to have worse surveillance systems, allowing an enhanced number of underreports and less detection of asymptomatic and mildly symptomatic patients [15], which was also the case in Brazil [49].
Since one of the components of the HDI is life expectancy, the higher the HDI-LE, the older the people who live in the region. Since the elderly are more likely to experience to severe cases of COVID-19, mainly due to the senescent immune system, we expected a positive correlation between death and HDI-LE similar to the literature [44,52,53,54]. Unfortunately, in Brazil, the mortality rate of COVID-19 individuals in the public health system is more significant than in the private health system [33,55]. Thus, individuals who live in low-HDI-GNI regions have a higher mortality since the low purchasing power makes it difficult to buy a private health plan. Finally, the Brazilian regions with lower HDI, mainly in the north and northeast, have the lowest average number of years of study, from 18 to 29 years, which is a risk factor for death from COVID-19 [56,57]. Previous studies observed low educational attainment as a risk factor for death from COVID-19 [56,58]. Maybe, these individuals with a higher level of education have a better understanding of the disease, making this group seek early medical care, which increases the diagnoses rates [44]. Additionally, a previous study demonstrated a significant association between cumulative death and HDI (an increase of 0.1 in the HDI was associated with a nearly 40-fold increased likelihood of death) [22]. However, we observed that an increase in HDI of 0.1 was associated with a lower case fatality rate, ranging from −7% to −13%. In the same way, a study analyzing BRICS (Brazil, Russia, India, China, and South Africa) observed a negative correlation between HDI and new COVID-19 cases [59], which is also in accordance with our study, and, perhaps, it shows a similar pattern between the socioeconomic characteristics from these countries.
On the other hand, in Europe, the HDI was not associated with COVID-19 whatsoever. A study conducted in Barcelona, Spain, and another one that evaluated the whole European Union did not find an association between this socioeconomic index and COVID-19 cases and deaths, demonstrating the socioeconomic differences between Brazil and European countries [60,61]. From a global perspective, however, the countries with higher HDI scores reported more COVID-19 deaths [15], perhaps due to their better surveillance system, which may provide fewer underreports, and, since one of the components of HDI is life expectancy, countries with more older people had a higher mortality rate, as aforementioned.
Furthermore, our study also follows previously published Brazilian studies [23,24,45,62]. These studies observed a higher infection and case fatality rate in some specific, low-HDI areas of Brazil, such as Pernambuco and Alagoas states. Additionally, one of the most critical risks factors for COVID-19 is HDI [23,24,45,62]. Since Brazil is a vast territory with culture and socioeconomic characteristics that vary from region to region, it is expected for the Brazilian regions to have different HDI scores. However, regions with the lowest HDI score were the most affected by COVID-19, perhaps due to inadequate access to quality health care [63], leading to insufficient access to complex treatments such as IMV. Furthermore, high-HDI regions seemed to adhere more to non-pharmaceutical interventions, such as social isolation, hand hygiene, and facial masks, decreasing infection and case fatality rate [42].
Although the HDI in Rio de Janeiro is higher than in Amazonas, a higher mortality rate was observed in Rio de Janeiro, which contrasts with our results. Nevertheless, the Rio de Janeiro state has the second highest absolute number of subnormal clusters, just behind São Paulo state, which has poor structure and is overcrowding which might contribute to the COVID-19 spread [10,11] and, ultimately, to COVID-19 deaths. Furthermore, a recent report observed higher levels of COVID-19 death underreporting in the Amazonas state when compared to Rio de Janeiro [64]. These factors might explain, at least in part, the enhanced mortality in Rio de Janeiro, even though it has a higher HDI score.
We also observed a positive association between the GINI index and mortality due to COVID-19, which is similar to the current literature [65,66,67]. Interestingly, studies also reported an association between the GINI index and the incidence of COVID-19 cases [68,69]. Perhaps this association is because lower-income individuals tended to remain in jobs with a higher risk of contracting COVID-19, such as restaurants and hotels, which require high levels of person-to-person contact. Additionally, these individuals were usually unable to work from home [66,70]. Finally, we observed a correlation between ICU beds and COVID-19 mortality, especially between the public ICU beds, following previous Brazilian studies [9,71]. In Brazil, nearly 80% of the population depends on the public health system; however, even before the COVID-19 pandemic, almost 75% of health institutions reported insufficient numbers of ICU beds [72]. This might have played a crucial role in the collapse of the public health system in Brazil. For instance, in 2021, most Brazilian states reported a critical situation of occupancy in the ICU, with nearly 100% of the beds occupied [73]; there were also reports of patients being treated in hallways while awaiting ICU vacancy [74].
The number of COVID-19 infections and case fatality rate were more associated with patterns of socioeconomic vulnerability than population age or prior comorbidities. Brazil reported the first COVID-19 cases in São Paulo and Rio de Janeiro. After, the number of deaths rapidly increased in states with socioeconomic vulnerabilities, especially in northern and northeastern states, due to vulnerability, poor governmental response, and enhanced structural inequalities, which resulted in worse outcomes in the COVID-19 pandemic [75]. Although Brazil faces a challenge with the COVID-19 pandemic, more target responses are necessary through epidemiological data analysis, alongside social-economical indexes, such as HDI, and governmental investments.
Limitations: We retrieved the data regarding the COVID-19 cases and death from an online dataset attributed to the Ministry of Health; thus, the data input might have errors made by the person who typed the information. Loss of data and a limited number of markers to be evaluated were limitations of the dataset. Brazil reported many underreported COVID-19 cases, which could bias our analysis. We only analyzed the Brazilian states, making it difficult to extrapolate our data to the municipalities. The GINI index might be underestimated due to a lack of proper data; some states have relatively low access to ICUs, and the total subnormal cluster results might also be lacking data since some places in Brazil are difficult to access, which could bias our analysis using these indexes, thus, highlighting the importance of HDI.

5. Conclusions

In Brazil, regions with a low HDI and components demonstrated a higher case fatality rate due to COVID-19, mainly in individuals who needed IMV, than those with higher HDI. Although other indexes studied, such as number of ICU beds and GINI index, were also associated with the COVID-19 case fatality rate, they were not as relevant as the HDI. Brazil is a vast territory comprising cultural, social, and economic diversity, which mirrors the diversity of the HDI scores. Brazil is a model for the study of HDI influence on COVID-19 pandemic measures such as the COVID-19 case fatality rate.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph19095306/s1, Supplementary Material S1: Results; Supplementary Material S2: Geospatial Graphs. References [9,25,26,29,31] are cited in the supplementary materials.

Author Contributions

F.E.V. and M.N.B., data collection. C.V.C.P., M.N.B., and F.A.L.M., data validation and statistical analysis. The authors wrote the manuscript, approved it, and agreed on its submission to the journal. All authors have read and agreed to the published version of the manuscript.

Funding

M.N.B., Fundação de Amparo à Pesquisa do Estado de São Paulo (São Paulo Research Foundation; #2021/05810-7). F.E.V., Fundação de Amparo à Pesquisa do Estado de São Paulo (São Paulo Research Foundation; #2021/08437-5).

Institutional Review Board Statement

Ethical review and approval were waived for this study because the data used in our study were made publicly available and do not contain patients’ data. The study was consent-free since it did not present risks to the research participants.

Informed Consent Statement

Patient consent was waived because the data used in our study were made publicly available and do not contain patients’ data. The study was consent-free since it did not present risks to the research participants.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sanyaolu, A.; Okorie, C.; Marinkovic, A.; Patidar, R.; Younis, K.; Desai, P.; Hosein, Z.; Padda, I.; Mangat, J.; Altaf, M. Comorbidity and Its Impact on Patients with COVID-19. SN Compr. Clin. Med. 2020, 2, 1069–1076. [Google Scholar] [CrossRef]
  2. Wang, B.; Li, R.; Lu, Z.; Huang, Y. Does Comorbidity Increase the Risk of Patients with COVID-19: Evidence from Meta-Analysis. Aging 2020, 12, 6049–6057. [Google Scholar] [CrossRef] [PubMed]
  3. Comber, L.; O Murchu, E.; Drummond, L.; Carty, P.G.; Walsh, K.A.; De Gascun, C.F.; Connolly, M.A.; Smith, S.M.; O’Neill, M.; Ryan, M.; et al. Airborne Transmission of SARS-CoV-2 via Aerosols. Rev. Med. Virol. 2021, 31, e2184. [Google Scholar] [CrossRef] [PubMed]
  4. Park, M.; Cook, A.R.; Lim, J.T.; Sun, Y.; Dickens, B.L. A Systematic Review of COVID-19 Epidemiology Based on Current Evidence. J Clin. Med. 2020, 9, 967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sun, J.; He, W.-T.; Wang, L.; Lai, A.; Ji, X.; Zhai, X.; Li, G.; Suchard, M.A.; Tian, J.; Zhou, J.; et al. COVID-19: Epidemiology, Evolution, and Cross-Disciplinary Perspectives. Trends Mol. Med. 2020, 26, 483–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sagar, A.D.; Najam, A. The Human Development Index: A Critical Review. Ecol. Econ. 1998, 25, 249–264. [Google Scholar] [CrossRef]
  7. Desai, M. Human Development: Concepts and Measurement. Eur. Econ. Rev. 1991, 35, 350–357. [Google Scholar] [CrossRef]
  8. Costa de Assis, S.J.; Lopes, J.M.; de Lima Filho, B.F.; José Bouzas Sanchis, G.; Sousa Rodrigues Guedes, T.; Limeira Cavalcanti, R.; Araujo, D.N.; José Sarmento da Nóbrega, A.; Barbosa Otoni Gonçalves Guedes, M.; Giuseppe Roncalli da Costa Oliveira, A. Dissemination of COVID-19 in Inland Cities of Northeastern Brazil. PLoS ONE 2021, 16, e0253171. [Google Scholar] [CrossRef]
  9. Palamim, C.V.C.; Marson, F.A.L. COVID-19—The Availability of ICU Beds in Brazil during the Onset of Pandemic. Ann. Glob. Health 2020, 86, 100. [Google Scholar] [CrossRef]
  10. Urban Poverty, Inclusive Cities and Housing|GPSC. Available online: https://www.thegpsc.org/knowledge-products/urban-poverty-and-housing (accessed on 21 April 2022).
  11. Hasan, S.M.; Das, S.; Hanifi, S.M.A.; Shafique, S.; Rasheed, S.; Reidpath, D.D. A Place-Based Analysis of COVID-19 Risk Factors in Bangladesh Urban Slums: A Secondary Analysis of World Bank Microdata. BMC Public Health 2021, 21, 502. [Google Scholar] [CrossRef]
  12. Barbosa, I.R.; Galvão, M.H.R.; de Souza, T.A.; Gomes, S.M.; de Almeida Medeiros, A.; de Lima, K.C. Incidence of and Mortality from COVID-19 in the Older Brazilian Population and Its Relationship with Contextual Indicators: An Ecological Study. Rev. Bras. Geriatr. Gerontol. 2020, 23. [Google Scholar] [CrossRef]
  13. Carneiro, F.F.; Franco Netto, G.; Corvalan, C.; de Freitas, C.M.; Sales, L.B.F. Saúde ambiental e desigualdades: Construindo indicadores para o desenvolvimento sustentável. Ciência Saúde Coletiva 2012, 17, 1419–1425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Maciel, J.A.C.; Castro-Silva, I.I.; de Farias, M.R. Initial Analysis of the Spatial Correlation between the Incidence of COVID-19 and Human Development in the Municipalities of the State of Ceará in Brazil. Rev. Bras. Epidemiol. 2020, 23. [Google Scholar] [CrossRef]
  15. Thazhathedath Hariharan, H.; Surendran, A.T.; Haridasan, R.K.; Venkitaraman, S.; Robert, D.; Narayanan, S.P.; Mammen, P.C.; Siddharth, S.R.; Kuriakose, S.L. Global COVID-19 Transmission and Mortality—Influence of Human Development, Climate, and Climate Variability on Early Phase of the Pandemic. Geohealth 2021, 5, e2020GH000378. [Google Scholar] [CrossRef]
  16. Bonanad, C.; García-Blas, S.; Tarazona-Santabalbina, F.; Sanchis, J.; Bertomeu-González, V.; Fácila, L.; Ariza, A.; Núñez, J.; Cordero, A. The Effect of Age on Mortality in Patients with COVID-19: A Meta-Analysis with 611,583 Subjects. J. Am. Med. Dir. Assoc. 2020, 21, 915–918. [Google Scholar] [CrossRef] [PubMed]
  17. Pförtner, T.-K.; Hower, K.I. Educational Inequalities in Risk Perception, Perceived Effectiveness, Trust and Preventive Behaviour in the Onset of the COVID-19 Pandemic in Germany. Public Health 2022, 206, 83–86. [Google Scholar] [CrossRef] [PubMed]
  18. Roser, M.; Human Development Index (HDI). Our World Data 2014. Available online: https://ourworldindata.org/human-development-index (accessed on 21 March 2022).
  19. Human Development Data Center|Human Development Reports. Available online: http://hdr.undp.org/en/data (accessed on 5 June 2021).
  20. Troumbis, A.Y. Testing the Socioeconomic Determinants of COVID-19 Pandemic Hypothesis with Aggregated Human Development Index. J Epidemiol Community Health 2020, 75, 414–415. [Google Scholar] [CrossRef]
  21. Khalatbari-Soltani, S.; Cumming, R.C.; Delpierre, C.; Kelly-Irving, M. Importance of Collecting Data on Socioeconomic Determinants from the Early Stage of the COVID-19 Outbreak Onwards. J. Epidemiol. Community Health 2020, 74, 620–623. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, K.; He, M.; Zhuang, Z.; He, D.; Li, H. Unexpected Positive Correlation between Human Development Index and Risk of Infections and Deaths of COVID-19 in Italy. One Health 2020, 10, 100174. [Google Scholar] [CrossRef]
  23. Baggio, J.A.O.; Machado, M.F.; do Carmo, R.F.; da Costa Armstrong, A.; dos Santos, A.D.; de Souza, C.D.F. COVID-19 in Brazil: Spatial Risk, Social Vulnerability, Human Development, Clinical Manifestations and Predictors of Mortality—A Retrospective Study with Data from 59 695 Individuals. Epidemiol. Infect. 2021, 149, e100. [Google Scholar] [CrossRef]
  24. de Santana Cabral Silva, A.P.; de Holanda, E.R.; de Abreu, P.D.; de Arruda Freitas, M.V. COVID-19 in Children in the State of Pernambuco: Spatial Analysis of Confirmed Severe Cases and the Human Development Index. Rev. Soc. Bras. Med. Trop. 2021, 54, e0782-2020. [Google Scholar] [CrossRef]
  25. Atlas Brasil. Available online: http://www.atlasbrasil.org.br/ranking (accessed on 5 June 2021).
  26. Bem Vindo-OPENDATASUS. Available online: https://opendatasus.saude.gov.br/ (accessed on 21 April 2022).
  27. Human Development Report 2019. 2019. Available online: http://www.atlasbrasil.org.br (accessed on 21 March 2022).
  28. Brasil|Cidades e Estados|IBGE. Available online: https://www.ibge.gov.br/cidades-e-estados (accessed on 5 June 2021).
  29. IBGE|Portal Do IBGE|IBGE. Available online: https://www.ibge.gov.br/pt/inicio.html (accessed on 3 October 2021).
  30. Bureau, U.C. Gini Index. Available online: https://www.census.gov/topics/income-poverty/income-inequality/about/metrics/gini-index.html (accessed on 19 April 2022).
  31. Available online: https://portal.cfm.org.br/ (accessed on 21 March 2022).
  32. Brazil: WHO Coronavirus Disease (COVID-19) Dashboard with Vaccination Data. Available online: https://covid19.who.int (accessed on 24 March 2022).
  33. Boschiero, M.N.; Palamim, C.V.C.; Ortega, M.M.; Mauch, R.M.; Marson, F.A.L. One Year of Coronavirus Disease 2019 (COVID-19) in Brazil: A Political and Social Overview. Ann. Glob. Health 2021, 87, 44. [Google Scholar] [CrossRef] [PubMed]
  34. Bigoni, A.; Malik, A.M.; Tasca, R.; Carrera, M.B.M.; Schiesari, L.M.C.; Gambardella, D.D.; Massuda, A. Brazil’s Health System Functionality amidst of the COVID-19 Pandemic: An Analysis of Resilience. Lancet Reg. Health-Am. 2022, 10, 100222. [Google Scholar] [CrossRef] [PubMed]
  35. WHO Coronavirus (COVID-19) Dashboard|WHO Coronavirus (COVID-19) Dashboard with Vaccination Data. Available online: https://covid19.who.int/ (accessed on 26 April 2021).
  36. Bastos, L.S.; Ranzani, O.T.; Souza, T.M.L.; Hamacher, S.; Bozza, F.A. COVID-19 Hospital Admissions: Brazil’s First and Second Waves Compared. Lancet Respir. Med. 2021, 9, e82–e83. [Google Scholar] [CrossRef]
  37. Zeiser, F.A.; Donida, B.; da Costa, C.A.; de Oliveira Ramos, G.; Scherer, J.N.; Barcellos, N.T.; Alegretti, A.P.; Ikeda, M.L.R.; Müller, A.P.W.C.; Bohn, H.C.; et al. First and Second COVID-19 Waves in Brazil: A Cross-Sectional Study of Patients’ Characteristics Related to Hospitalization and in-Hospital Mortality. Lancet Reg. Health-Am. 2022, 6, 100107. [Google Scholar] [CrossRef] [PubMed]
  38. Sabino, E.C.; Buss, L.F.; Carvalho, M.P.S.; Prete, C.A.; Crispim, M.A.E.; Fraiji, N.A.; Pereira, R.H.M.; Parag, K.V.; da Silva Peixoto, P.; Kraemer, M.U.G.; et al. Resurgence of COVID-19 in Manaus, Brazil, despite High Seroprevalence. Lancet 2021, 397, 452–455. [Google Scholar] [CrossRef]
  39. Banho, C.A.; Sacchetto, L.; Campos, G.R.F.; Bittar, C.; Possebon, F.S.; Ullmann, L.S.; de Carvalho Marques, B.; da Silva, G.C.D.; Moraes, M.M.; Parra, M.C.P.; et al. Impact of SARS-CoV-2 Gamma Lineage Introduction and COVID-19 Vaccination on the Epidemiological Landscape of a Brazilian City. Commun. Med. 2022, 2, 41. [Google Scholar] [CrossRef]
  40. Faria, N.R.; Mellan, T.A.; Whittaker, C.; Claro, I.M.; Candido, D.D.S.; Mishra, S.; Crispim, M.A.E.; Sales, F.C.S.; Hawryluk, I.; McCrone, J.T.; et al. Genomics and Epidemiology of the P.1 SARS-CoV-2 Lineage in Manaus, Brazil. Science 2021, 372, 815–821. [Google Scholar] [CrossRef] [PubMed]
  41. Funk, T.; Pharris, A.; Spiteri, G.; Bundle, N.; Melidou, A.; Carr, M.; Gonzalez, G.; Garcia-Leon, A.; Crispie, F.; O’Connor, L.; et al. Characteristics of SARS-CoV-2 Variants of Concern B.1.1.7, B.1.351 or P.1: Data from Seven EU/EEA Countries, Weeks 38/2020 to 10/2021. Eurosurveillance 2021, 26, 2100348. [Google Scholar] [CrossRef]
  42. Gozzi, N.; Tizzoni, M.; Chinazzi, M.; Ferres, L.; Vespignani, A.; Perra, N. Estimating the Effect of Social Inequalities on the Mitigation of COVID-19 across Communities in Santiago de Chile. Nat. Commun. 2021, 12, 2429. [Google Scholar] [CrossRef]
  43. de Souza, R.C.; Almeida, E.R.M.; Fortaleza, C.M.C.B.; Miot, H.A. Factors Associated with COVID-19 Mortality in Municipalities in the State of São Paulo (Brazil): An Ecological Study. Rev. Soc. Bras. Med. Trop. 2022, 55. [Google Scholar] [CrossRef] [PubMed]
  44. Shahbazi, F.; Khazaei, S. Socio-Economic Inequality in Global Incidence and Mortality Rates from Coronavirus Disease 2019: An Ecological Study. New Microbes New Infect. 2020, 38, 100762. [Google Scholar] [CrossRef] [PubMed]
  45. Martins, L.D.; da Silva, I.; Batista, W.V.; de Fátima Andrade, M.; de Freitas, E.D.; Martins, J.A. How Socio-Economic and Atmospheric Variables Impact COVID-19 and Influenza Outbreaks in Tropical and Subtropical Regions of Brazil. Environ. Res. 2020, 191, 110184. [Google Scholar] [CrossRef] [PubMed]
  46. Hub: Volume de Passageiros em Conexão Triplica em Fortaleza-Negócios. Available online: https://diariodonordeste.verdesmares.com.br/negocios/hub-volume-de-passageiros-em-conexao-triplica-em-fortaleza-1.2004425 (accessed on 26 March 2022).
  47. Galvan, D.; Effting, L.; Cremasco, H.; Adam Conte-Junior, C. Can Socioeconomic, Health, and Safety Data Explain the Spread of COVID-19 Outbreak on Brazilian Federative Units? Int. J. Environ. Res. Public Health 2020, 17, 8921. [Google Scholar] [CrossRef] [PubMed]
  48. Hallal, P.C.; Hartwig, F.P.; Horta, B.L.; Silveira, M.F.; Struchiner, C.J.; Vidaletti, L.P.; Neumann, N.A.; Pellanda, L.C.; Dellagostin, O.A.; Burattini, M.N.; et al. SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two Successive Nationwide Serological Household Surveys. Lancet Glob. Health 2020, 8, e1390–e1398. [Google Scholar] [CrossRef]
  49. Carvalho, T.A.; Boschiero, M.N.; Marson, F.A.L. COVID-19 in Brazil: 150,000 Deaths and the Brazilian Underreporting. Diagn. Microbiol. Infect. Dis. 2021, 99, 115258. [Google Scholar] [CrossRef]
  50. Mendes, M.F.; Pereira, L.R.; Lima, T.M.; Melani, V.F.; Palamim, C.V.C.; Boschiero, M.N.; Marson, F.A.L. COVID-19 Pandemic Evolution in the Brazilian Indigenous Population. J. Racial Ethn. Health Disparities 2021, 1–17. [Google Scholar] [CrossRef] [PubMed]
  51. Sansone, N.M.S.; Boschiero, M.N.; Ortega, M.M.; Ribeiro, I.A.; Peixoto, A.O.; Mendes, R.T.; Marson, F.A.L. Severe Acute Respiratory Syndrome by SARS-CoV-2 Infection or Other Etiologic Agents among Brazilian Indigenous Population: An Observational Study from the First Year of Coronavirus Disease (COVID)-19 Pandemic. Lancet Reg. Health-Am. 2022, 8, 100177. [Google Scholar] [CrossRef]
  52. Choi, K.W.; Chau, T.N.; Tsang, O.; Tso, E.; Chiu, M.C.; Tong, W.L.; Lee, P.O.; Ng, T.K.; Ng, W.F.; Lee, K.C.; et al. Outcomes and Prognostic Factors in 267 Patients with Severe Acute Respiratory Syndrome in Hong Kong. Ann. Intern. Med. 2003, 139, 715–723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hong, K.-H.; Choi, J.-P.; Hong, S.-H.; Lee, J.; Kwon, J.-S.; Kim, S.-M.; Park, S.Y.; Rhee, J.-Y.; Kim, B.-N.; Choi, H.J.; et al. Predictors of Mortality in Middle East Respiratory Syndrome (MERS). Thorax 2018, 73, 286–289. [Google Scholar] [CrossRef] [PubMed]
  54. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; et al. Clinical Course and Risk Factors for Mortality of Adult Inpatients with COVID-19 in Wuhan, China: A Retrospective Cohort Study. Lancet 2020, 395, 1054–1062. [Google Scholar] [CrossRef]
  55. Benchmarking-COVID-19-UTIs Brasileiras. Available online: http://www.utisbrasileiras.com.br (accessed on 9 June 2021).
  56. Concepción-Zavaleta, M.J.; Coronado-Arroyo, J.C.; Zavaleta-Gutiérrez, F.E.; Concepción-Urteaga, L.A. Does Level of Education Influence Mortality of SARS-CoV-2 in a Developing Country? Int. J. Epidemiol. 2021, 49, 2091–2093. [Google Scholar] [CrossRef] [PubMed]
  57. Anuário Brasileiro da Educação Básica 2020–Escolaridade. Available online: https://www.moderna.com.br/anuario-educacao-basica/2020/escolaridade.html (accessed on 6 June 2021).
  58. Feldman, J.M.; Bassett, M.T. Variation in COVID-19 Mortality in the US by Race and Ethnicity and Educational Attainment. JAMA Netw. Open 2021, 4, e2135967. [Google Scholar] [CrossRef]
  59. Zhu, J.; Yan, W.; zhu, L.; Liu, J. COVID-19 Pandemic in BRICS Countries and Its Association with Socio-Economic and Demographic Characteristics, Health Vulnerability, Resources, and Policy Response. Infect. Dis. Poverty 2021, 10, 97. [Google Scholar] [CrossRef] [PubMed]
  60. López-Gay, A.; Spijker, J.; Cole, H.V.S.; Marques, A.G.; Triguero-Mas, M.; Anguelovski, I.; Marí-Dell’Olmo, M.; Módenes, J.A.; Álamo-Junquera, D.; López-Gallego, F.; et al. Sociodemographic Determinants of Intraurban Variations in COVID-19 Incidence: The Case of Barcelona. J. Epidemiol. Community Health 2021, 76, jech-2020-216325. [Google Scholar] [CrossRef]
  61. Cifuentes-Faura, J. Factors Influencing the COVID-19 Mortality Rate in the European Union: Importance of Medical Professionals. Public Health 2021, 200, 1–3. [Google Scholar] [CrossRef]
  62. Alves, J.G.; Lima, T.P. COVID-19 Lethality in Non-Elderly Individuals in Cities with Different Human Development Index. Trop. Dr. 2021, 51, 124–125. [Google Scholar] [CrossRef]
  63. Campos, A.C.V.; Borges, C.M.; Vargas, A.M.D.; Leles, C.R.; e Ferreira, E.F. Social and Health Indicators as a Measure of Access to Primary Heathcare in Brazil. Ciência Saúde Coletiva 2011, 16, 4349–4355. [Google Scholar] [CrossRef] [Green Version]
  64. Paixão, B.; Baroni, L.; Pedroso, M.; Salles, R.; Escobar, L.; de Sousa, C.; de Freitas Saldanha, R.; Soares, J.; Coutinho, R.; Porto, F.; et al. Estimation of COVID-19 under-Reporting in the Brazilian States Through SARI. New Gener. Comput. 2021, 39, 623–645. [Google Scholar] [CrossRef]
  65. de Sousa Pinto, A.; Rodrigues, C.A.; Nascimento, C.L.; da Cruz, L.A.; dos Santos, E.G.; Nunes, P.C.; Costa, M.G.R.; da Costa Rocha, M.O. COVID-19 Epidemic Curve in Brazil: A Sum of Multiple Epidemics, Whose Inequality and Population Density in the States Are Correlated with Growth Rate and Daily Acceleration. An Ecological Study. Rev. Soc. Bras. Med. Trop. 2022, 55, e0118-2021. [Google Scholar] [CrossRef]
  66. Tan, A.X.; Hinman, J.A.; Abdel Magid, H.S.; Nelson, L.M.; Odden, M.C. Association Between Income Inequality and County-Level COVID-19 Cases and Deaths in the US. JAMA Netw. Open 2021, 4, e218799. [Google Scholar] [CrossRef] [PubMed]
  67. Oronce, C.I.A.; Scannell, C.A.; Kawachi, I.; Tsugawa, Y. Association between State-Level Income Inequality and COVID-19 Cases and Mortality in the USA. J. Gen. Intern. Med. 2020, 35, 2791–2793. [Google Scholar] [CrossRef] [PubMed]
  68. Ganasegeran, K.; Jamil, M.F.A.; Appannan, M.R.; Ch’ng, A.S.H.; Looi, I.; Peariasamy, K.M. Spatial Dynamics and Multiscale Regression Modelling of Population Level Indicators for COVID-19 Spread in Malaysia. Int. J. Environ. Res. Public Health 2022, 19, 2082. [Google Scholar] [CrossRef] [PubMed]
  69. da Silva, C.F.A.; Silva, M.C.; dos Santos, A.M.; Rudke, A.P.; do Bonfim, C.V.; Portis, G.T.; de Almeida Junior, P.M.; de Santana Coutinho, M.B. Spatial Analysis of Socio-Economic Factors and Their Relationship with the Cases of COVID-19 in Pernambuco, Brazil. Trop. Med. Int. Health 2022, 27, 397–407. [Google Scholar] [CrossRef]
  70. Baker, M.G. Nonrelocatable Occupations at Increased Risk during Pandemics: United States, 2018. Am. J. Public Health 2020, 110, 1126–1132. [Google Scholar] [CrossRef]
  71. de Souza, C.D.F.; de Almeida Santana, G.B.; Leal, T.C.; de Paiva, J.P.S.; da Silva, L.F.; Santos, L.G.; Santos, V.S.; do Carmo, R.F. Relationship between the Intensive Care Unit Beds and Mortality by COVID-19 in Brazil. Rev. Assoc. Med. Bras. 2021, 67, 645–649. [Google Scholar] [CrossRef]
  72. Rache, B.; Rocha, R.; Nunes, L.; Spinola, P.; Malik, A.M.; Massuda, A. Necessidades de Infraestrutura do SUS em Preparo à COVID-19: Leitos de UTI, Respiradores e Ocupação Hospitalar; Instituto de Estudos para Políticas de Saúde (IEPS): São Paulo, Brazil, 2020. [Google Scholar]
  73. COVID-19: Cenário é de Piora na Ocupação de Leitos de UTI. Available online: https://portal.fiocruz.br/noticia/covid-19-cenario-e-de-piora-na-ocupacao-de-leitos-de-uti (accessed on 22 April 2022).
  74. da Silva, S.J.R.; Pena, L. Collapse of the Public Health System and the Emergence of New Variants during the Second Wave of the COVID-19 Pandemic in Brazil. One Health 2021, 13, 100287. [Google Scholar] [CrossRef]
  75. Rocha, R.; Atun, R.; Massuda, A.; Rache, B.; Spinola, P.; Nunes, L.; Lago, M.; Castro, M.C. Effect of Socioeconomic Inequalities and Vulnerabilities on Health-System Preparedness and Response to COVID-19 in Brazil: A Comprehensive Analysis. Lancet Glob. Health 2021, 9, e782–e792. [Google Scholar] [CrossRef]
Figure 1. Description of inclusion and exclusion criteria for the study participants and the distribution of coronavirus disease COVID-19 individuals with ventilatory support and outcomes (death—case fatality rate and clinical recovery). SARI, severe acute respiratory infection; SARS, severe acute respiratory syndrome confirmed by a positive SARS-CoV-2 real-time polymerase chain reaction (RT-PCR) test. We obtained the data from OpenDataSUS [26].
Figure 1. Description of inclusion and exclusion criteria for the study participants and the distribution of coronavirus disease COVID-19 individuals with ventilatory support and outcomes (death—case fatality rate and clinical recovery). SARI, severe acute respiratory infection; SARS, severe acute respiratory syndrome confirmed by a positive SARS-CoV-2 real-time polymerase chain reaction (RT-PCR) test. We obtained the data from OpenDataSUS [26].
Ijerph 19 05306 g001
Figure 2. Distribution of the Human Development Index (HDI) (including the indicators for the educational level (HDI-E), life expectancy (HDI-LE), and gross national income per capita (HDI-GNI)) and the case fatality rate for coronavirus disease COVID-19 according to the Brazilian states and the Federal District. We presented the Pearson correlation matrix to compare HDI (HDI-LE, HDI-E, and HDI-GNI) with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities. We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. WC, without comorbidities; %, percentage. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Figure 2. Distribution of the Human Development Index (HDI) (including the indicators for the educational level (HDI-E), life expectancy (HDI-LE), and gross national income per capita (HDI-GNI)) and the case fatality rate for coronavirus disease COVID-19 according to the Brazilian states and the Federal District. We presented the Pearson correlation matrix to compare HDI (HDI-LE, HDI-E, and HDI-GNI) with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities. We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. WC, without comorbidities; %, percentage. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Ijerph 19 05306 g002
Figure 3. Distribution of the Human Development Index (HDI) (including the indicators for the educational level (HDI-E), life expectancy (HDI-LE), and gross national income per capita (HDI-GNI)) and the case fatality rate for coronavirus disease COVID-19 in individuals who needed invasive mechanical ventilation (IMV) according to the Brazilian states and the Federal District. We presented the Pearson correlation matrix to compare HDI (HDI-LE, HDI-E, and HDI-GNI) with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities. We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. WC, without comorbidities; %, percentage. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Figure 3. Distribution of the Human Development Index (HDI) (including the indicators for the educational level (HDI-E), life expectancy (HDI-LE), and gross national income per capita (HDI-GNI)) and the case fatality rate for coronavirus disease COVID-19 in individuals who needed invasive mechanical ventilation (IMV) according to the Brazilian states and the Federal District. We presented the Pearson correlation matrix to compare HDI (HDI-LE, HDI-E, and HDI-GNI) with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities. We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. WC, without comorbidities; %, percentage. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Ijerph 19 05306 g003
Figure 4. Univariate regression analysis between the case fatality rate due to coronavirus disease COVID-19 and the Human Development Index (HDI). (A) Overall case fatality rate. (B) The overall case fatality rate in COVID-19 individuals without comorbidities. (C) The case fatality rate in individuals who received invasive mechanical ventilation (IMV). (D) The case fatality rate in individuals who needed IMV and did not have comorbidities. The Y represents the case fatality rate as a dependent marker, and the Y describes the HDI as an independent marker. AC, Acre; AL, Alagoas; AP, Amapá; AM, Amazonas; BA, Bahia; CE, Ceará; ES, Espírito Santo; FD, Federal District; GO, Goiás; MA, Maranhão; MT, Mato Grosso; MS, Mato Grosso do Sul; MG, Minas Gerais; PA, Pará; PB, Paraíba; PR, Paraná; PE, Pernambuco; PI, Piauí; RJ, Rio de Janeiro; RN, Rio Grande do Norte; RS, Rio Grande do Sul; RO; Rondônia; RR, Roraima; SC, Santa Catarina; SP, São Paulo; SE, Sergipe; TO, Tocantins. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Figure 4. Univariate regression analysis between the case fatality rate due to coronavirus disease COVID-19 and the Human Development Index (HDI). (A) Overall case fatality rate. (B) The overall case fatality rate in COVID-19 individuals without comorbidities. (C) The case fatality rate in individuals who received invasive mechanical ventilation (IMV). (D) The case fatality rate in individuals who needed IMV and did not have comorbidities. The Y represents the case fatality rate as a dependent marker, and the Y describes the HDI as an independent marker. AC, Acre; AL, Alagoas; AP, Amapá; AM, Amazonas; BA, Bahia; CE, Ceará; ES, Espírito Santo; FD, Federal District; GO, Goiás; MA, Maranhão; MT, Mato Grosso; MS, Mato Grosso do Sul; MG, Minas Gerais; PA, Pará; PB, Paraíba; PR, Paraná; PE, Pernambuco; PI, Piauí; RJ, Rio de Janeiro; RN, Rio Grande do Norte; RS, Rio Grande do Sul; RO; Rondônia; RR, Roraima; SC, Santa Catarina; SP, São Paulo; SE, Sergipe; TO, Tocantins. We obtained the data from OpenDataSUS [26] and from the AtlasBR [25].
Ijerph 19 05306 g004
Figure 5. Pearson correlation matrix between GINI coefficient and case fatality rate according to the presence of comorbidities and the need for invasive mechanical ventilation (IMV). We presented the Pearson correlation matrix to compare the GINI coefficient (*, GINI of household income per capita, at average prices for the year; **, of the average real monthly income of people aged 14 and over actually received in the reference month, for all jobs, at average prices for the year) with the case fatality rate for overall coronavirus disease COVID-19 individuals and COVID-19 individuals without comorbidities (WC). We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. We obtained the data from OpenDataSUS [26] and from the Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) [29].
Figure 5. Pearson correlation matrix between GINI coefficient and case fatality rate according to the presence of comorbidities and the need for invasive mechanical ventilation (IMV). We presented the Pearson correlation matrix to compare the GINI coefficient (*, GINI of household income per capita, at average prices for the year; **, of the average real monthly income of people aged 14 and over actually received in the reference month, for all jobs, at average prices for the year) with the case fatality rate for overall coronavirus disease COVID-19 individuals and COVID-19 individuals without comorbidities (WC). We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. We obtained the data from OpenDataSUS [26] and from the Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) [29].
Ijerph 19 05306 g005
Figure 6. Univariate regression analysis between the case fatality rate due to coronavirus disease COVID-19 and the GINI coefficient and intensive care unit (ICU) beds per 10,000 inhabitants in the public health system. (A) The case fatality rate in individuals who needed IMV and did not have comorbidities. (B) The case fatality rate in individuals who received invasive mechanical ventilation (IMV). The Y represents the case fatality rate as a dependent marker, and the Y describes the GINI coefficient and ICU beds per 10,000 inhabitants at the public health system as an independent marker. AC, Acre; AL, Alagoas; AP, Amapá; AM, Amazonas; BA, Bahia; CE, Ceará; ES, Espírito Santo; FD, Federal District; GO, Goiás; MA, Maranhão; MT, Mato Grosso; MS, Mato Grosso do Sul; MG, Minas Gerais; PA, Pará; PB, Paraíba; PR, Paraná; PE, Pernambuco; PI, Piauí; RJ, Rio de Janeiro; RN, Rio Grande do Norte; RS, Rio Grande do Sul; RO; Rondônia; RR, Roraima; SC, Santa Catarina; SP, São Paulo; SE, Sergipe; TO, Tocantins. We obtained the data from OpenDataSUS [26], the Federal Council of Medicine website, Palamim and Marson (2020) [9,31], and from the Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) [29].
Figure 6. Univariate regression analysis between the case fatality rate due to coronavirus disease COVID-19 and the GINI coefficient and intensive care unit (ICU) beds per 10,000 inhabitants in the public health system. (A) The case fatality rate in individuals who needed IMV and did not have comorbidities. (B) The case fatality rate in individuals who received invasive mechanical ventilation (IMV). The Y represents the case fatality rate as a dependent marker, and the Y describes the GINI coefficient and ICU beds per 10,000 inhabitants at the public health system as an independent marker. AC, Acre; AL, Alagoas; AP, Amapá; AM, Amazonas; BA, Bahia; CE, Ceará; ES, Espírito Santo; FD, Federal District; GO, Goiás; MA, Maranhão; MT, Mato Grosso; MS, Mato Grosso do Sul; MG, Minas Gerais; PA, Pará; PB, Paraíba; PR, Paraná; PE, Pernambuco; PI, Piauí; RJ, Rio de Janeiro; RN, Rio Grande do Norte; RS, Rio Grande do Sul; RO; Rondônia; RR, Roraima; SC, Santa Catarina; SP, São Paulo; SE, Sergipe; TO, Tocantins. We obtained the data from OpenDataSUS [26], the Federal Council of Medicine website, Palamim and Marson (2020) [9,31], and from the Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) [29].
Ijerph 19 05306 g006
Figure 7. Pearson correlation matrix between intensive care unit (ICU) beds per 10,000 inhabitants (overall ICU beds per 10,000 inhabitants, overall ICU beds per 10,000 inhabitants in the public health system, and overall ICU beds per 10,000 inhabitants in the private health system) and case fatality rate according to the presence of comorbidities and the need for invasive mechanical ventilation (IMV). We presented the Pearson correlation matrix to compare the number of ICU beds per 10,000 inhabitants with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities (WC). We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. We obtained the data from OpenDataSUS [26], the Federal Council of Medicine website, and Palamim and Marson (2020) [9,31].
Figure 7. Pearson correlation matrix between intensive care unit (ICU) beds per 10,000 inhabitants (overall ICU beds per 10,000 inhabitants, overall ICU beds per 10,000 inhabitants in the public health system, and overall ICU beds per 10,000 inhabitants in the private health system) and case fatality rate according to the presence of comorbidities and the need for invasive mechanical ventilation (IMV). We presented the Pearson correlation matrix to compare the number of ICU beds per 10,000 inhabitants with the case fatality rate for overall COVID-19 individuals and COVID-19 individuals without comorbidities (WC). We considered the following categorization for the Pearson correlation test: very high positive/negative correlation, 0.9 to 1.0; high positive/negative correlation, 0.7 to 0.9; moderate positive/negative correlation, 0.5 to 0.7; low positive/negative correlation, 0.30 to 0.50; negligible correlation, 0.00 to 0.30. We presented an alpha error of 0.05 in all statistical analyses. We presented the case fatality rate as a percentage. We obtained the data from OpenDataSUS [26], the Federal Council of Medicine website, and Palamim and Marson (2020) [9,31].
Ijerph 19 05306 g007
Table 1. Description of case fatality rate and cases of coronavirus disease COVID-19 in Brazil during the first year of the pandemic regarding the presence of comorbidities and the need for invasive mechanical ventilation (IMV) *.
Table 1. Description of case fatality rate and cases of coronavirus disease COVID-19 in Brazil during the first year of the pandemic regarding the presence of comorbidities and the need for invasive mechanical ventilation (IMV) *.
States and the Federal DistrictWith and without ComorbiditiesIncrease in Cases of Deaths %Only without ComorbiditiesIncrease in Cases of Deaths %
All Individuals
-N (%)
IMV
-N (%)
All Individuals
-N (%)
IMV
-N (%)
Acre929/2349 (39.5)76/90 (84.4)44.9449/1273 (35.2)24/32 (75.0)39.8
Alagoas3496/8353 (41.8)1783/1958 (91.0)49.2793/2647 (30.0)394/436 (90.3)60.3
Amapá1001/2956 (33.8)628/769 (81.6)47.8290/1255 (23.1)181/244 (74.1)51.0
Amazonas12,042/27,845 (43.2)4142/4751 (87.1)44.04571/13,098 (34.9)1373/1616 (85.0)50.1
Bahia12,741/31,340 (40.6)6367/7315 (87.0)46.32385/8080 (29.5)1030/1237 (83.2)53.7
Ceará14,884/33,020 (45.0)5901/6991 (84.4)39.34005/10,907 (36.7)1334/1690 (78.9)42.2
Federal District6411/21,724 (29.5)3556/4356 (81.6)52.1989/5628 (17.7)516/702 (73.5)55.8
Espírito Santo4618/7668 (60.2)1900/2098 (90.5)30.31146/2191 (52.3)400/460 (87.0)34.7
Goiás11,551/29,458 (39.2)5341/6326 (84.4)34.73422/11,036 (31.0)1428/1814 (78.7)47.7
Maranhão4331/8402 (51.5)955/1107 (86.2)34.81403/3363 (41.7)255/310 (82.2)40.6
Mato Grosso3172/19,226 (16.4)992/1236 (80.2)63.8823/11,096 (7.4)231/322 (71.7)64.3
Mato Grosso do Sul4527/13,351 (34.0)2146/2590 (83.0)49.0619/4337 (14.3)261/423 (61.7)47.4
Minas Gerais25,833/73,165 (35.3)9061/11,535 (78.5)43.24598/20,325 (22.6)1461/2159 (67.7)45.1
Pará11,139/24,642 (45.2)3990/4357 (91.5)46.34057/11,082 (36.6)1293/1584 (81.6)45.0
Paraíba5822/13,271 (43.8)2950/3233 (91.2)47.31240/3717 (33.3)585/661 (88.5)55.2
Paraná15,619/48,978 (31.8)7801/9719 (80.2)48.43006/15,311 (19.6)1495/2085 (71.7)52.1
Pernambuco11,759/26,749 (43.9)1882/2337 (80.5)36.54155/11,191 (37.1)446/589 (76.1)39.0
Piauí3246/9021 (36.0)1934/2188 (88.3)52.3633/2774 (22.8)357/437 (81.7)58.9
Rio de Janeiro37,508/74,109 (50.6)10,567/12,224 (86.4)35.810,749/24,05 (44.6)2282/2770 (82.3)37.7
Rio Grande do Norte3828/8605 (44.4)1869/2115 (88.3)43.9643/1996 (32.2)294/356 (82.5)50.3
Rio Grande do Sul20,753/55,587 (37.3)10,813/13,349 (81.0)43.72273/13,458 (16.9)1152/1683 (68.4)51.5
Rondônia3523/7445 (47.3)2268/2566 (88.3)41.01475/3177 (46.4)933/1063 (87.7)41.3
Roraima1121/1957 (57.2)702/747 (94.0)36.8364/788 (46.2)226/248 (91.1)44.9
Santa Catarina11,068/33,210 (33.3)4770/6369 (74.9)41.62097/11,200 (18.7)869/1404 (61.9)43.2
São Paulo77,767/256,324 (30.3)30,297/40,358 (75.0)44.716,008/83,320 (19.2)5471/8393 (65.2)46.0
Sergipe3643/5806 (62.7)1835/2071 (88.6)25.91010/1847 (54.7)359/425 (84.5)29.8
Tocantins1836/3940 (46.6)612/735 (84.5)36.9520/1379 (37.7)127/169 (75.1)37.4
*, we described the deaths/(deaths and recovery cases) data. We calculated the case fatality rate increase by the difference between the % of deaths in individuals who needed IMV and the % of deaths among all individuals. We obtained the data from OpenDataSUS [26]. N, number of individuals; %, percentage.
Table 2. Description of GINI index *, intensive care unit (ICU) beds per 10,000 inhabitants **, and the number of occupied households in subnormal clusters *** in Brazil according to states and the Federal District.
Table 2. Description of GINI index *, intensive care unit (ICU) beds per 10,000 inhabitants **, and the number of occupied households in subnormal clusters *** in Brazil according to states and the Federal District.
States and Federal DistrictGINI a,bICU Beds Per 10,000 InhabitantsNumber of Occupied Households in Subnormal Clusters (%) d
TotalPublic Health SystemPrivate Health System (Beneficiaries) c
Acre0.5750.5280.900.713.5419,148 (8.53)
Alagoas0.5260.4801.450.865.1864,568 (6.68)
Amapá0.5600.4961.030.338.2736,835 (21.58)
Amazonas0.5720.5371.240.793.43393,955 (34.59)
Bahia0.5480.5371.320.646.49469,677 (10.62)
Ceará0.5530.5341.330.764.01243,848 (9.20)
Federal District0.5830.5533.390.898.7862,179 (6.65)
Espírito Santo0.5130.4682.721.195.60306,439 (26.10)
Goiás0.4740.4412.081.115.9235,801 (1.55)
Maranhão0.5280.5461.120.598.18144,625 (7.85)
Mato Grosso0.4570.4402.620.8910.6322,429 (1.99)
Mato Grosso do Sul0.4810.4701.780.944.716766 (0.74)
Minas Gerais0.5040.4862.061.303.14231,385 (3.43)
Pará0.5310.5271.180.576.42432,518 (19.68)
Paraíba0.5400.5381.510.945.4964,225 (5.07)
Paraná0.4850.4682.521.543.93135,188 (3.57)
Pernambuco0.5780.5511.961.096.31327,090 (10.55)
Piauí0.5460.5921.100.565.5750,382 (5.49)
Rio de Janeiro0.5240.4843.790.978.70717,326 (12.63)
Rio Grande do Norte0.5580.5491.710.945.2041,868 (3.97)
Rio Grande do Sul0.4860.4892.101.333.31133,021 (3.50)
Rondônia0.4780.4471.631.017.0823,236 (4.37)
Roraima0.5470.5320.920.576.293033 (2.12)
Santa Catarina0.4290.4171.581.032.6132,416 (1.46)
São Paulo0.5410.5332.631.193.801066,813 (7.09)
Sergipe0.5720.5601.481.013.4753,203 (7.37)
Tocantins0.4980.4751.430.868.379733 (2.14)
*, the GINI index was obtained from the Brazilian Institute of Geography and Statistics (IBGE) website (Instituto Brasileiro de Geografia e Estatística, in Portuguese) [29], and the calculation was performed in 2017; **, the number of occupied households in subnormal clusters was obtained from the IBGE website [29], and the calculation was performed in 2019; ***, the ICU beds distribution was obtained from the Federal Council of Medicine website and Palamim and Marson (2020) [9,31]. a, GINI index of household income per capita, at average prices for the year (first column); b, GINI index of the average real monthly income of people aged 14 and over actually received in the reference month, for all jobs, at average prices for the year (second column); c, the proportion represents only the number of ICU beds among individuals that have access for the private health system (beneficiaries). In the public health system, we considered all Brazilian individuals; d, we decided to present the absolute number of occupied households in subnormal clusters and the relative number (%) of occupied households in subnormal clusters, which represents the proportion of this type of occupation and the total number of households.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Palamim, C.V.C.; Boschiero, M.N.; Valencise, F.E.; Marson, F.A.L. Human Development Index Is Associated with COVID-19 Case Fatality Rate in Brazil: An Ecological Study. Int. J. Environ. Res. Public Health 2022, 19, 5306. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph19095306

AMA Style

Palamim CVC, Boschiero MN, Valencise FE, Marson FAL. Human Development Index Is Associated with COVID-19 Case Fatality Rate in Brazil: An Ecological Study. International Journal of Environmental Research and Public Health. 2022; 19(9):5306. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph19095306

Chicago/Turabian Style

Palamim, Camila Vantini Capasso, Matheus Negri Boschiero, Felipe Eduardo Valencise, and Fernando Augusto Lima Marson. 2022. "Human Development Index Is Associated with COVID-19 Case Fatality Rate in Brazil: An Ecological Study" International Journal of Environmental Research and Public Health 19, no. 9: 5306. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph19095306

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop