Ischemic heart disease (IHD) is projected to continue to be a major cause of death in 2030 . Although the advances in cardiology, IHD has become an epidemic disease and one of the main causes of hospitalizations and morbimortality worldwide . According to the World Health Organization, more than 9 million deaths from IHD were estimated in 2016 . In 2015, most IHD deaths occurred in Brazil, which was the upper-middle-income country with the second-highest number of deaths .
It’s well known that socioeconomically disadvantaged groups are at high risk of cardiovascular diseases with poor outcome . In addition, the inequalities in healthcare, major disparities in individual income and education , low accessibility to healthcare centers and limited health coverage offered to the patients  may be responsible for IHD mortality rates in Brazil.
Easy access to healthcare centers lowers the mortality and morbidity generally . However, some studies have shown that the healthcare resources are unequally distributed across regions in Brazil . Thus, geostatistical tools such as spatial analysis to verify the geographically variance of IHD risks could increase the optimization of healthcare resource allocations .
Brazil is a perfect setting for an analysis of the interaction between healthcare access, socioeconomic inequality, and outcomes. Exploring and mapping spatially the disease’s risk and mortality is determinant for its comprehension and proper healthcare organization and planning . Using tools based on geographic information systems, it is possible to identify patterns associated with the sociodemographic profiles, coverage of tests, procedures and accessibility for each area, which collaborates in taking measures and modifying health policy plans [11, 12].
Given this context, this study aims to analyze the spatial distribution of mortality rates due to IHD and its association with sociodemographic indicators, exam coverage, and accessibility to health in the State of Paraná, Brazil, in order to verify whether high mortality due to IHD is related to the disparities in the municipalities.
This is an ecological study using spatial analysis tools based on secondary data from IHD mortality data in Paraná State, from 2013 to 2017. For the assessment of methodological quality, we followed the recommendations from the guideline Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) .
According to Brazilian Institute of Geography and Statistics data (2015), Paraná is located in the southern region of Brazil, occupying an area of 199,880 km2, with latitude coordinates 22º30’58’’ and 26º43’00’’, and longitude coordinates 48º05’37’’ and 54º37’08’’ . Paraná has a total of 10,439,601 inhabitants, with the majority living in urban area (85.3%), and being the sixth most populous state in Brazil (5.5% of the total population) [14, 15]. This state is bordered by the countries Argentina and Paraguay and the Brazilian states of São Paulo, Mato Grosso do Sul, and Santa Catarina. Paraná is divided by 10 geographical mesoregions, composed of 399 municipalities. The state has a GINI Index of 0.60, with values, diversified by municipality, ranging from 0.33 to 0.66 (Figure 1) .
In the Brazilian Mortality Information System (SIM), all deaths from IHD aged between 30 and 79 years old registered with the ICD-10 (I20.0 to I25.9) in the years 2013 to 2017 [5, 17] were collected. The population data by city, also aged between 30 and 79 years old, were obtained from Brasilian Institute of Geography and Statistics (IBGE) . Thus, we calculated the age adjusted specific IHD mortality rate per municipality, according to the American Heart Association Guidelines .
In order to verify the association between IHD mortality rate and exam coverage indicators, sociodemographic variables, and accessibility to health, this study evaluated the following variables: aging rate, income, illiteracy rate, proportion of family health program coverage, catheterization rate, scintigraphy rate, angioplasty rate, chemical reperfusion rate, echocardiography rate, revascularization rate, exercise test rate, and accessibility to cardiologists, chemical reperfusion centers, and mechanical reperfusion centers . These variables were selected to represent characteristics that may affect IHD mortality rates . The variables of race, ethnicity, and gender were analyzed in the study but excluded because they did not have significance in multivariate regression and increased multicollinearity.
The health coverage data were obtained from the Brazilian National Register of Health Facilities (CNES) , while the indicator variables to create the examination rates were obtained from the Brazilian Hospital Information System (SIH) by the Microdatasus package, using the R studio software .
Cartographic base data for all municipalities in Paraná State in 2017 was obtained from the Paraná Land, Cartography and Geology Institute – ITCG (www.pr.gov.br/itcg) .
All data are publicly and freely available from the Brazilian Health System Informatics Department (DATASUS). Table 1 lists all datasets accessed and utilized in this study.
|DATASUS – Mortality Information System (SIM)||
National Register of Health Facilities
Paraná Land, Cartography and Geology Institute
Brazilian Institute of Geography and Statistics
United Nations Development Program
(microdatasus) DATASUS – Hospital information system – SIH
Portal of the Secretariat of Primary Health Care
Two Step Floating Catchment Area (2SFCA) on ArcGIS® software (ESRI Company, USA, 2020) was used to create three accessibility indexes in order to verify the geographical accessibility (availability and proximity) , to cardiologists, chemical reperfusion centers and mechanical reperfusion centers. In the first step, the capacity of each center was calculated using the total number of cardiologists and service capacity available in referral centers chemical and mechanical, inside a buffer with 60 km radius surrounding each cardiologists’ office and the referral centers were created to verify the population potentially coverage . For the access time, the distance of 60 km was evaluated, considering that the patient must reach the reperfusion centers within the recommended time of 60 minutes. In the first step, the capacity of each center is determined, while the second step consists of summing up the capacities within the buffer taking into account the overlapped health service area. Thus, the accessibility is the available amount of cardiologists and referral centers chemical and mechanical per local population added up within 60 km of each health service [25, 26, 27]. The result is an accessibility index to cardiologists, chemical reperfusion referral centers and mechanical reperfusion centers for each municipality. A higher accessibility index indicates more access to these resources for a given municipality.
Spatial analysis was applied to determine the existence of significant Global Spatial Autocorrelation (Moran’s I) and Local Index of Spatial Autocorrelation (LISA) to verify the influence of spatiality on IHD mortality rate using GeoDa software version 1.12.0 .
Moran’s I only verifies the presence of spatiality globally, meaning that local patterns of spatial association may be hidden . In order to avoid this, LISA was applied to identify the significant formation of clusters locally , categorized as: high-high (HH) clusters, in which they are a set of municipalities with high IHD mortality rates surrounded by other municipalities with high IHD mortality rates. On the other hand, low-low clusters (LL) are groups of municipalities with low IHD mortality rates surrounded by municipalities with low IHD mortality rates. The global and local spatial autocorrelation coefficients were considered significant when p < 0.05 .
To identify which sociodemographic, health coverage, procedures and health access indicators had the greatest geospatial impact on the spatial distribution of IHD mortality rate, we performed a multivariate spatial regression analysis using Ordinary Least Squares (OLS) and Geographically Weighted Regression (GWR). The OLS regression does not take spatial dependence into the analysis . On the other hand, GWR considers spatial dependency by performing several separated regressions, identifying significant geographic clusters within the studied area . Additionally, GWR produces an estimate for the association between IHD mortality rate and its predictor variables, from analysis of the local spatial variability for each municipality.
The selection of variables for the final model was based on the criteria of low multicollinearity, with a result of 29.7 (Conditional number test – Geoda), preventively showing that multicollinearity does not influence the regression results . Thus, initially, each significant variable in the global model (OLS) were tested in the local model (GWR), allowing to verify the presence of spatial non-stationarity, and that the correlations of the regression variables vary locally in relation to space, as indicated by the chloroplastic maps, indicating the differential spatial impact of each variable on IHD mortality rates.
The chosen model between OLS and GWR was the one that provided the best fit based on the p < 0.05, highest adjusted R2, lowest Akaike Information Criterion (AIC) and Moran’s I residuals. GeoDa software was utilized to generate the OLS model; R Studio software was used to execute the GWR model and QGIS software, version 2.14.9, was used to generate choropleth maps.
This study did not require ethical approval and consent form, since we used data from secondary sources in publicly available, government and online databases (https://datasus.saude.gov.br/informacoes-de-saude-tabnet/). Nevertheless, this study followed the Brazilian National Guidelines for research with human beings, according to Resolution 510/2016 of the National Health Council.
During the analyzed period, 22,920 deaths due to IHD occurred mainly among men (64.5%), aging between 70 and 79 years (36.7%), white (73.2%), married (49.1%), and with up to seven years of schooling (57.6%). Most deaths occurred in the hospital (49.5%), followed by deaths at home (30.2%).
Regarding spatial patterns of death distribution from IHD in the 399 cities in Parana State, on average 44.0/100,000 (+/– 10.6; Standard Deviation) inhabitants died per year from IHD. The 399 cities presented IHD mortality rates between 22.8 and 76.3 per 100,000 inhabitants, these cities being located at high rates (50.6 to 76.3 per 100,000 inhabitants) mainly in Central-West, NorthWest, North Pioneer, and Southeast regions of the Paraná state (Figure 2A).
Global Moran Univariate Analysis indicated the existence of a positive spatial autocorrelation (Moran’s I = 0.633, p < 0.01), demonstrating that municipalities with high mortality rates by IHD tend to be surrounded by neighboring towns with similar high rates.
According to Figure 2B, the LISA analysis detected six high-high (HH) pattern clusters, indicating cities with high IHD mortality rates close to neighbors also showing high IHD mortality rates, covering 16.5% of the municipalities in the state. It was also observed in six low-high (LH) pattern clusters, indicating cities with low IHD mortality rates with neighbors with high IHD mortality rates. Only one high-low cluster (HL) was identified in the North Central of Paraná, indicating cities with high IHD mortality rates with neighbors with low IHD mortality rates. Considering the two largest HH clusters, one included municipalities located in the Central-West and North Central, and the other included municipalities located in the North Pioneer and East-Central of the state.
The spatial regression was performed and from the initial set of 14 variables, 8 remained significant. Table 2 presents the results of the multivariate spatial regression analysis to identify variables that were correlated with municipalities’ IHD mortality rates. According to global OLS results, the following variables were removed (p > 0.05): accessibility to chemical reperfusion centers, proportion of family health program coverage, scintigraphy rate, chemical reperfusion rate, echocardiography rate, and ergometric (exercise) test rate (Table 2). In addition, according to global OLS results, the variables that remained in the modeling process (p < 0.05) were: accessibility to cardiologists, accessibility to hospitals with chemical reperfusion, aging rate, income ratio, illiteracy, catheterization rate, angioplasty rate, and revascularization rate (Table 2).
|Variables||Global OLS coefficient||Local coefficientof GWR|
|Est||Std Error||T value||VIF||Pr(>|t|)||Min||IstQu||Median||3rd||Max|
|Access index of Cardiologists||–3.70||1.31||–2.83||4.17||< 0.01**||–2.59||–9.83||–4.73||–6.45||1.34|
|Access of chemical reperfusion centers||3.51||9.20||3.82||1.18||< 0.01**||–7.13||2.08||3.32||5.54||1.21|
|Access of mechanical reperfusion centers||–2.46||2.58||–0.09||3.63||0.92||–6.00||–5.86||7.20||1.33||3.37|
|Aging rate||1.02||0.28||3.64||1.42||< 0.01**||–4.31||7.07||1.03||1.34||2.22|
|Income ratio||–1.88||0.07||–2.64||1.17||< 0.01**||–1.05||–3.05||–1.89||–4.36||–6.16|
|Proportion of Family health program coverage||–5.12||0.03||–0.19||1.29||0.85||–9.63||–3.75||1.03||4.49||1.19|
|Catheterization rate||–2.14||0.53||4.01||1.16||< 0.01**||–6.78||–2.44||–1.72||–9.86||2.08|
|Chemical reperfusion rate||–1.51||5.85||–0.26||1.13||0.80||–3.89||–7.89||1.01||9.41||5.18|
|Moran’s I Residuals||0.05||–||–||–||–||0.00||–||–||–||–|
|Moran’s p Residuals||0.00||–||–||–||–||0.08||–||–||–||–|
As indicated in Table 2, the GWR multivariate model, considering the significant variables, showed a significant improvement, presenting better AIC, Moran’s I residuals and adjusted R2 when compared with the OLS model, indicating that these variables correlate differently with IHD mortality rates in the geographical level.
The GWR model had better performance in understanding the associations between the variables and IHD mortality rate. However, the model adjustment (Adjusted R2) varied spatially (Figure 3). The North Central, North Pioneer, East-Central, Metropolitan and parts of the West and South-West mesoregions presented better model adjustment. The Northwest, Center-West, South-Central, and Southeast mesoregions presented lower values of GWR model adjustment, below 0.6 (Figure 3).
The spatial impact of each significant variable included in the GWR model on IHD mortality rate are indicated in the Figures 4, 5, and 6.
The GWR analysis results presented the distinctive spatial impact and pattern of each variable, influencing IHD mortality rates. As indicated in Figures 4, 5, and 6, some variables presented a homogeneous influence over the municipalities, while others presented more heterogeneous spatial influence. The following variables presented negative correlation with IHD mortality rate: access to cardiologists (most on Central-West and part of Northwest, North-Central, West and Southwest mesoregions) (Figure 4 A3), access to chemical reperfusion centers (most on Central-West and few on South-Central mesoregions) (Figure 4 B3), angioplasty rate (most on Northwest, Central-West, West, and North-Central; few on South-Central, Southeast and Metropolitan mesoregions) (Figure 6 B3) and revascularization rate (most on Northwest and part of North Central and Central-West mesoregions) (Figure 6 C3). These data indicated that in these regions, higher access to cardiologists, chemical perfusion centers, and higher angioplasty rate and revascularization rate are associated with lower IHD mortality rate.
In addition, aging rate (Figure 5 A3) and illiteracy rate (Figure 5 C3) presented positive correlation, while income ratio presented negative correlation (Figure 5 B3), with IHD mortality rate, in all significant municipalities and mesoregions, indicating that municipalities’ higher aging rate and illiteracy rate and lower income ratio are associated with higher IHD mortality rates.
Concomitantly, in other mesoregions, differential spatial impact was observed concerning some of the accessibility and intervention variables and the following variables presented positive correlation with IHD mortality rate: access to cardiologists (East-Central mesoregion) (Figure 4 A3), access to chemical reperfusion centers (most on North-Central and Southeast, part of Northwest and Metropolitan, and few on West, North Pioneer and South-Central) (Figure 4 B3), catheterization rate (all significant municipalities and mesoregions) (Figure 6 A3), angioplasty rate (North Pioneer) (Figure 6 B3), and revascularization rate (most on South-Central and East-Central, part of West and Southwest and few on Southeast, North-Central, North Pioneer, and Metropolitan) (Figure 6 C3). These data indicated that in some specific municipalities and mesoregions facilitated access to healthcare and cardiological interventions did not present favorable outcomes; on the contrary, surprisingly, they were correlated with higher IHD mortality rates, suggesting that the quantity and quality of access and services in these locations are not enough and need to be improved to prevent IHD mortality. Additionally, other unidentified factors may be influencing these outcomes, as addressed in the Discussion section.
In this study, the spatial distribution of IHD mortality rate by municipalities was evaluated in order to verify its association with sociodemographic, exam coverage, and access to health indicators in Parana state, located in Brazil.
The GWR analysis presented the best fitted multivariate model and in the preliminary analysis showed that sex/gender and race/ethnicities were not significantly related with IHD mortality (p > 0.05) rate and were not included in the final model. This result is different from the study of Mozaffarian et al. (2015) who found that men and afro-descendants presented a higher mortality rate due to IHD in the United States .
The reason for such discrepancy is not known. However, the Parana state was colonized predominantly by Europeans, consisting of a white population. Thus, since the proportion of afro-descendants in Parana’s general population is very low, one possibility is that its impact on the municipalities IHD mortality rate are also very low; alternatively, it suggests that this population of men and afro-descendants may be subjected to different IHD risk factors according to Brazil’s region.
The results in our study indicated that accessibility to cardiologists, accessibility to chemical reperfusion centers, and angioplasty, catheterization, and revascularization rates presented differential spatial association with municipalities’ IHD mortality rates. The study showed that in some municipalities and mesoregions, higher access to cardiologists, chemical perfusion centers, and higher angioplasty rate and revascularization rate were associated with lower IHD mortality rates . These findings are supported by the literature, because referring these patients to cardiologists care is associated with more frequent use of efficacious therapies and may improve the outcome .
On the contrary, in some specific municipalities and mesoregions, facilitated access to the cardiologists and chemical perfusion centers and to higher angioplasty, catheterization and revascularization rates did not present favorable outcomes, and this was correlated with higher IHD mortality rates. Although unexpected, these findings may indicate that in these specific municipalities and mesoregions, others factors may be related with patients’ poor IHD outcomes, despite the higher healthcare accessibility and higher rates of cardiological interventions. These other factors may be:
Regarding the sociodemographic variables, aging rate and illiteracy rate presented positive correlation with IHD mortality rate, while income ratio presented negative correlation. Other studies also indicated the same correlation in high- and middle-income countries [39, 40]. Lower health literacy is an invisible barrier to healthcare delivery, which has been associated with limited knowledge of health conditions, medications, poorer overall health status, higher healthcare costs, and mortality . Additionally, the study by Paasche-Orlow et al. (2005) corroborates our study because they showed that lower health literacy is more prevalent among older adults and individuals with less education . In addition, a similar study also analyzed income and IHD mortality rates and the mortality were higher in low middle income regions because these regions have more difficulties than high-income regions to find resources for the treatment of IHD patients . Therefore, illiteracy, aging, and low income lead to poor outcome of IDH patients.
Geographic disparity in cardiovascular disease interventions is common worldwide . A study in the United States of America found that < 50% of patients were transferred to hospitals for cardiac catheterization , increasing patients’ poor outcomes. The present study also indicated that low accessibility to healthcare and low rates of interventions are associated with higher IHD mortality rates in some mesoregions of Parana state.
In the present study, the results indicated that in some specific municipalities and mesoregions of Parana state improvement of health accessibility and exam coverage were not associated with decrease of IHD mortality rates, so future studies are needed to identify the possible other factors involved in these specific settings and scenarios.
However, the data of the present study support that the expansion of the hospitals with high complex procedures must be considered; in addition, data suggest that public health policies aiming to reduce the time of the initial cardiological treatment of IHD patients may reduce the poor outcomes of these patients and decrease the Parana state IDH mortality rates, considering the spatial differences and specificities of each mesoregion, socioeconomic factors, and healthcare accessibility scenarios.
There are some limitations that need to be considered in the present study such as the use of secondary data, which may present under-notifications of IHD mortality cases in Parana state. However, the data quality of the Mortality Information System obtained from the website of the Brazilian Ministry of Health has increased its credibility and has been improved. In addition, since the number of cases of IHD mortality are high the final effect of the under-notifications is possibly low.
For the spatial distribution of IHD mortality rates, the association of risk indicators had an important influence on the clusters found, and it was possible to estimate the differential spatial impact of level of accessibility, assistance of health services, and socioeconomic factors on these rates in Parana state mesoregions. The study indicated the complexity of the factors associated with IDH mortality rates unveiled by the spatial analysis, showing the distinctive scenarios of intervention, indicating the vulnerable areas where public interventions may occur, which may facilitate the allocation of health resources and application of more suited prevention policies.
Those interested in the specific data and codes used in this study can enter the corresponding link: https://figshare.com/search?q=10.6084%2Fm9.figshare.13042649.
The present study was carried out with the support of the Coordination for the Improvement of Higher Education Personnel – Brazil (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior Brasil – CAPES) – Financing Code 001. We would like to thank the Group of Technologies and Geoprocessing in Health (Grupo de Tecnologias e Geoprocessamento em Saúde – GETS), for providing support with geospatial analysis.
The authors have no competing interests to declare.
ACD wrote the paper, collected the data, and made the figures; LLS and OKN helped with the writing and review of the manuscript; RBP, YPT, MT and MB cleared the data for statistical analysis; JFHCS and PHI performed the statistical analysis; RLF expertise in cardiology provided substantial information for the discussion; CAS and JRNV revised the statistical analysis; and LA designed the study and provided the final approval of the manuscript.
Mathers CD, Loncar D. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med. 2006; 3(11): e442. DOI: https://doi.org/10.1371/journal.pmed.0030442
Henriques A, Araújo C, Viana M, Laszczynska O, Pereira M, Bennett K, et al. Disability-adjusted life years lost due to ischemic heart disease in mainland Portugal, 2013. Rev Port Cardiol. 2017; 36(4): 273–81. DOI: https://doi.org/10.1016/j.repc.2016.08.011
World Health Organization. Global Health Estimates 2016: Deaths by Cause, Age, Sex, by Country and by Region, 2000–2016. Geneva, Switzerland: World Health Organization; 2018. http://www.who.int/healthinfo/global_burden_disease/estimates/en/. Accessed November 16, 2018.
Nowbar AN, Gitto M, Howard JP, Francis DP, Al-Lamee R. Mortality from ischemic heart disease. Circ Cardiovasc Qual Outcomes. 2019; 12(6): e005375. DOI: https://doi.org/10.1161/CIRCOUTCOMES.118.005375
Havranek EP, Mujahid MS, Barr DA, et al. American Heart Association Council on Quality of Care and Outcomes Research, Council on Epidemiology and Prevention, Council on Cardiovascular and Stroke Nursing, Council on Lifestyle and Cardiometabolic Health, and Stroke Council. Social determinants of risk and outcomes for cardiovascular disease: A scientific statement from the American Heart Association. Circulation. 2015; 132: 873–898. DOI: https://doi.org/10.1161/CIR.0000000000000228
Baena CP, Chowdhury R, Schio NA, et al. Ischaemic heart disease deaths in Brazil: Current trends, regional disparities and future projections. Heart. 2013; 99: 1359–1364. DOI: https://doi.org/10.1136/heartjnl-2013-303617
Ribeiro AL, Duncan BB, Brant LC, Lotufo PA, Mill JG, Barreto SM. Cardiovascular health in Brazil: Trends and perspectives. Circulation. 2016; 133(4): 422–33. DOI: https://doi.org/10.1161/CIRCULATIONAHA.114.008727
Okuyama K, Akai K, Kijima T, Abe T, Isomura M, Nabika T. Effect of geographic accessibility to primary care on treatment status of hypertension. PLoS One. 2019; 14(3): e0213098. DOI: https://doi.org/10.1371/journal.pone.0213098
Hertz JT, Fu T, Vissoci JR, Rocha TAH, Carvalho E, et al. The distribution of cardiac diagnostic testing for acute coronary syndrome in the Brazilian healthcare system: A national geospatial evaluation of health access. PLoS One. 2019; 14(1): e0210502. DOI: https://doi.org/10.1371/journal.pone.0210502
Du Q, Zhang M, Li Y, Luan H, Liang S, Ren F. Spatial patterns of ischemic heart disease in Shenzhen, China: A Bayesian multi-disease modelling approach to inform health planning policies. Int J Environ Res Public Health. 2016; 13: 436. DOI: https://doi.org/10.3390/ijerph13040436
Rocha-Brischiliari SC, Andrade L, Nihei OK, Brischiliari A, Hortelan MS, Carvalho MDB, et al. Spatial distribution of breast cancer mortality: Socioeconomic disparities and access to treatment in the state of Parana, Brazil. PLoS One. 2018; 13(10). DOI: https://doi.org/10.1371/journal.pone.0205253
Silva LL, Valverde TL, de Carvalho AD, de Andrade L, de Barros CMD, Pelloso SM, Pedroso RB, De Lima SRB. Geospatial analysis and factors associated with poor outcome for multidrug-resistant tuberculosis. The International Journal of Tuberculosis and Lung Disease. 2020; 24(8): 844–846(3). DOI: https://doi.org/10.5588/ijtld.19.0692
Von EE, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: Guidelines for reporting observational studies. Lancet. 2007; 370: 1453–57. https://www.equator-network.org/reporting-guidelines/strobe/. DOI: https://doi.org/10.1016/S0140-6736(07)61602-X
Instituto Brasileiro de Geografia e Estatística [homepage]. Mapas: bases e referenciais, bases cartográficas, malhas digitais, https://mapas.ibge.gov.br/bases-e-referenciais/bases-cartograficas/malhas-digitais; 2015 [accessed 18 January 2018].
Instituto de Terras, Cartografia e Geologia do Paraná. Dados geoespacias de referência, http://www.itcg.pr.gov.br/modules/faq/category.php?categoryid=8; 2015 [accessed 20 January 2020].
Programa das Nações Unidas para o Desenvolvimento. Relatório do Desenvolvimento Humano 2013, https://www.br.undp.org/content/brazil/pt/home/library/idh/relatorios-de-desenvolvimento-humano/relatorio-do-desenvolvimento-humano-200012.html; 2017 [accessed 16 January 2020].
Ministério da Saúde. Sistema de Informações de Mortalidade do SUS (SIM/SUS): O Ministério, http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sim/cnv/qiPR.def;  [accessed 22 January 2020].
Alrawashdeh A, Nehme Z, Williams B, et al. Emergency medical service delays in ST-elevation myocardial infarction: A meta-analysis. Heart. 2020; 106: 365–373. DOI: https://doi.org/10.1136/heartjnl-2019-315034
Ministério da Saúde. Sistema de Informações Hospitalares do SUS (SIH/SUS). [S.l.]: O Ministério, http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sih/cnv/qiPR.def;  [accessed 22 January 2020].
Atlas do Desenvolvimento Humano no Brasil. Base de Dados: Censos 2000 e 2010, http://www.atlasbrasil.org.br/2013/pt/home/; [s.d.] [accessed 20 January 2020].
CNES—Cadastro Nacional de Estabelecimentos de Saude: Departamento de Informatica do SUS; 2017. http://datasus.saude.gov.br/sistemas-e-aplicativos/cadastros-nacionais/cnes. [accessed 20 January 2020].
Saldanha RF, Bastos RR, Barcellos C. Microdatasus: Pacote para download e pré-processamento de microdados do Departamento de Informática do SUS (DATASUS). Cad Saúde Pública [Internet]. 2019; 35(9): e00032419. DOI: https://doi.org/10.1590/0102-311x00032419
Instituto de Terras, Cartografia e Geologia do Paraná. Dados geoespacias de referência, http://www.itcg.pr.gov.br/modules/faq/category.php?categoryid=8; 2015 [accessed 20 January 2020].
McGrail MR. Spatial accessibility of primary health care utilising the two-step floating catchment área method: An assessment of recente improvements. Int J Health Geographics. 2012; 11(50): 1–12. DOI: https://doi.org/10.1186/1476-072X-11-50
Vart P, Coresh J, Kwak L, Ballew SH, Heiss G, Matsushita K. Socioeconomic status and incidence of hospitalization with lower-extremity peripheral artery disease: Atherosclerosis risk in communities study. Journal of the American Heart Association. 2017; 6(8). DOI: https://doi.org/10.1161/JAHA.116.004995
Luo W, Wang F. Measures of spatial accessibility to health care in a GIS environment: synthesis and a case study in the Chicago region. Environ Plann B Plann Des. 2003; 30(6): 865–84. DOI: https://doi.org/10.1068/b29120
Bergamini M, Iora PH, Thiago Augusto HR, Tchuisseu YP, de Carvalho Dutra A, João Felipe Herman CS, et al. Mapping risk of ischemic heart disease using machine learning in a Brazilian state. PLoS One. 2020 12; 15(12). DOI: https://doi.org/10.1371/journal.pone.0243558
Dezman Z, Andrade L, Vissoci JR, El-Gabri D, Johnson A, Hirshon JM, et al. Hotspots and causes of motor vehicle crashes in Baltimore, Maryland: A geospatial analysis of five years of police crash and census data. Injury. 2016; 47(11): 2450–8. DOI: https://doi.org/10.1016/j.injury.2016.09.002
Fontes KB, Jacinto Alarcão AC, Nihei OK, Pelloso SM, Andrade L, Barros Carvalho MD. Regional disparities in the intimate partner sexual violence rate against women in Paraná State, Brazil, 2009–2014: An ecological study. BMJ Open. 2018; 8(2): e018437. DOI: https://doi.org/10.1136/bmjopen-2017-018437
Andrade L, Lynch C, Carvalho E, Rodrigues CG, Vissoci JRN, Passos GF, et al. System dynamics modeling in the evaluation of delays of care in ST-segment elevation myocardial infarction patients within a tiered health system. PLoS One. 2014; 9(7). DOI: https://doi.org/10.1371/journal.pone.0103577
Anselin L, Syabri I, Kho Y. Geoda: An introduction to Spatial Data Analysis. Geographical Analysis. 2006; 38: 05–22. DOI: https://doi.org/10.1111/j.0016-7363.2005.00671.x
Lu B, Harris P, Charlton M, Brunsdon C. The GWmodel R package: Further topics for exploring spatial heterogeneity using geographically weighted models. Geo-spatial Information Science. 2014; 17(2): 85–101. DOI: https://doi.org/10.1080/10095020.2014.917453
Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, et al. American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Heart disease and stroke statistics–2015 update: A report from the American Heart Association. Circulation. 2015; 131: e29–e322. DOI: https://doi.org/10.1161/CIR.0000000000000152
Auerbach AD, Hamel MB, Califf RM, et al. Patient characteristics associated with care by a cardiologist among adults hospitalized with severe congestive heart failure. SUPPORT Investigators. Study to understand prognoses and preferences for outcomes and risks of treatments. J Am Coll Cardiol. 2000; 36: 2119–25. DOI: https://doi.org/10.1016/S0735-1097(00)01005-6
Kalogeris T, Baines CP, Krenz M, Korthuis RJ. Cell biology of ischemia/reperfusion injury. Int Rev Cell Mol Biol. 2012; 298: 229–317. DOI: https://doi.org/10.1016/B978-0-12-394309-5.00006-7
Darling C, Saczynski JS, McManus DD, Lessard D, Spencer FA, Goldberg RJ. Delayed hospital presentation in acute decompensated heart failure: Clinical and patient reported factors. Heart Lung. 2013; 42(4): 281–286. DOI: https://doi.org/10.1016/j.hrtlng.2013.01.007
Gravely-Witte S, Jurgens CY, Tamim H, Grace SL. Length of delay in seeking medical care by patients with heart failure symptoms and the role of symptom-related factors: A narrative review. Eur J Heart Fail. 2010; 12(10): 1122–1129. DOI: https://doi.org/10.1093/eurjhf/hfq122
Goldberg RJ, Steg PG, Sadiq I, et al. Extent of, and factors associated with, delay to hospital presentation in patients with acute coronary disease (the GRACE registry). Am J Cardiol. 2002; 89(7): 791–796. DOI: https://doi.org/10.1016/S0002-9149(02)02186-0
Pednekar MS, Gupta R, Gupta PC. Illiteracy, low educational status, and cardiovascular mortality in India. BMC Public Health. 2011; 11: 567. Published 2011 Jul 15. DOI: https://doi.org/10.1186/1471-2458-11-567
Marmot M. Health in an unequal world. Lancet. 2006; 368: 2081–94. DOI: https://doi.org/10.1016/S0140-6736(06)69746-8
Magnani JW, Mujahid MS, Aronow HD, et al. Health literacy and cardiovascular disease: fundamental relevance to primary and secondary prevention: A scientific statement from the American Heart Association. Circulation. 2018; 138. DOI: https://doi.org/10.1161/CIR.0000000000000579
Paasche-Orlow MK, Parker RM, Gazmararian JA, Nielsen-Bohlman LT, Rudd RR. The prevalence of limited health literacy. J Gen Intern Med. 2005; 20(2): 175–184. DOI: https://doi.org/10.1111/j.1525-1497.2005.40245.x
Finegold JA, Asaria P, Francis DP. Mortality from ischaemic heart disease by country, region, and age: statistics from World Health Organization and United Nations. Int J Cardiol. 2013; 168(2): 934–945. DOI: https://doi.org/10.1016/j.ijcard.2012.10.046
Ferreira-Gonzalez I, Permanyer-Miralda G, Heras M, et al. Patterns of use and effectiveness of early invasive strategy in non-ST-segment elevation acute coronary syndromes: An assessment by propensity score. Am Heart J. 2008; 156: 946–53, 953.e2. DOI: https://doi.org/10.1016/j.ahj.2008.06.032
Roe MT, Chen AY, Delong ER, et al. Patterns of transfer for patients with non-ST-segment elevation acute coronary syndrome from community to tertiary care hospitals. Am Heart J. 2008; 156: 185–192. DOI: https://doi.org/10.1016/j.ahj.2008.01.033