Association between number of institutions with coronary computed tomography angiography and regional mortality ratio of acute myocardial infarction: a nationwide ecological study using a spatial Bayesian model
International Journal of Health Geographics volume 17, Article number: 13 (2018)
Coronary computed tomography angiography (CTA) has demonstrated high diagnostic accuracy for detection of coronary artery stenosis, and healthcare providers can detect coronary artery disease in earlier stages before it develops into more serious clinical conditions such as acute myocardial infarction (AMI). We hypothesized that the mortality ratio of AMI in regions with a higher density of coronary CTA is lower than that in regions with a lower density of coronary CTA.
This ecological and cross-sectional study using secondary data targeted all secondary medical service areas (SMSAs) in Japan (n = 349). We obtained the numbers of cardiologists, institutions with coronary CTA, and institutions with a cardiac catheterization laboratory (CCL) as medical resources, socioeconomic factors, lifestyle factors, exercise habit factors, and AMI mortality data from a Japanese national database. We evaluated the association between the number of these medical resources and the standardized mortality ratio (SMR) of AMI in each SMSA using a hierarchical Bayesian model accounting for spatial autocorrelation (i.e., a conditional autoregressive model). We assumed a Poisson distribution for the observed number of AMI-related deaths and set the expected number of AMI-related deaths as the offset variable.
The number of institutions with coronary CTA was negatively and significantly associated with the SMR of AMI (relative risk [RR] 0.900; 95% credible interval [CI] 0.848–0.953), while the SMR in each SMSA was not significantly associated with the number of either cardiologists (RR 0.997; 95% CI 0.988–1.004) or institutions with a CCL (RR 1.026; 95% CI 0.963–1.096).
We observed a significant association between the number of institutions with coronary CTA and the SMR of AMI. Effective allocation of coronary CTA in each region is recommended, and it would be important to clarify the standing position of coronary CTA in regional networking for AMI treatment in the future.
Coronary artery disease (CAD) is the leading cause of death worldwide . Coronary artery disease caused 7.2 million deaths in 2008 and 7.4 million deaths in 2012, which exceeded 10% of the total global mortality [1, 2]. In Japan, CAD was the second-highest cause of death in 2015 , and the number of CAD-induced deaths in Japan was the third highest among high-income countries . Therefore, both primary and secondary prevention of CAD is of great importance.
During the past decade, the diagnostic accuracy of coronary computed tomography angiography (CTA), which is a noninvasive and anatomical imaging technique for detection of CAD, has been improved [4,5,6]. Using coronary CTA, health providers can detect preclinical CAD in advance and administer treatment before the development of a more critical condition such as acute myocardial infarction (AMI). For example, in the updated National Institute for Health and Care Excellence (NICE) guidelines, coronary CTA is the first-line investigation for patients presenting with new-onset chest pain due to suspected CAD . The updated NICE guideline CG95 poses a far wider and more uniform geographical delivery of coronary CTA as a major challenge . Considering that construction of regional networks is internationally recommended for treatment of AMI , evidence of how to allocate coronary CTA in each regional is important. Moreover, from the viewpoint of public health, appropriate allocation of coronary CTA seems of great importance to efficiently reduce the number of AMI-related deaths in each region. However, no study has evaluated the association between the geographical distribution of coronary CTA and AMI-related death by region. We hypothesized that the mortality ratio of AMI in regions with a higher density of coronary CTA would be lower than that in regions with a lower density of coronary CTA. In the present study, we evaluated the nationwide association between the geographical distribution of coronary CTA and the mortality ratio of AMI.
This was an observational, cross-sectional study using secondary data. This study was also a nationwide ecological study targeting the whole of Japan. Japan comprises 47 prefectures, and the Japanese government established subprefectural medical regions called secondary medical service areas (SMSAs) . An SMSA is defined as a medical unit that evaluates demand and supply of health resources for inpatient treatment, especially hospital care. We targeted all SMSAs in Japan according to a survey in 2011 (n = 349).
Geographical information, such as municipality boundary data, was obtained from the Municipality Map Maker for Web . We combined each municipality-level parameter to form the SMSA data because each SMSA consists of several municipalities. We used ArcGIS version 10.2.1 (ESRI Japan Inc., Tokyo, Japan).
The Japan Ministry of Health, Labour and Welfare (MHLW) conducts a biennial census survey of the working conditions of all physicians, dentists, and pharmacists. Details from this census survey were used to obtain the number of self-reported cardiologists in each SMSA. The Japan MHLW also conducts a detailed triennial survey of all medical institutions. The data from this survey were used to obtain the number of institutions with coronary CTA and the number with a cardiac catheterization laboratory (CCL) in each SMSA. Because of the survey questionnaire style, the institutions with coronary CTA included those with cardiac magnetic resonance imaging (MRI). However, we assumed that almost all institutions with cardiac MRI had also coronary CTA when compared with another data source from the Japanese Circulation Society . In addition, we regarded institutions with cardiac digital radiography as those with a CCL. We obtained permission from the MHLW to analyze these survey data. Both survey forms are available on the MHLW website .
Data regarding socioeconomics, lifestyle, and exercise habits were obtained from the national databases of the following surveys: Population Census, Survey of municipal tax situation, School Basic Survey, Comprehensive Survey of Living Conditions, and National Health and Nutrition Survey. All data are available from the national Japanese portal site, e-Stat .
The AMI mortality data were obtained from the official database of Vital Statistics, which covers all deaths in Japan in principle. In this study, we obtained the number of deaths classified by the International Classification of Diseases-10 codes I21 (AMI) and I22 (subsequent myocardial infarction).
We used the numbers of cardiologists, institutions with coronary CTA, and institutions with a CCL per 100,000 population. As socioeconomic covariates, we used available socioeconomic factors at an SMSA level from e-Stat that had been identified in previous studies [10, 15, 16]. These factors were the population density (people per hectare), average per capita income (in 1000 JYN), unemployment rate (%), proportion of university graduates per population aged > 15 years (%), and divorce rate (incidence per 1000 population). Because the population density distribution was extremely skewed (skewness, 3.12; kurtosis, 13.4), we took a logarithm of that variable. As covariates of lifestyle and exercise habits, we used data on smoking habits, drinking habits, and the estimated body mass index as lifestyle data and the average number of daily steps as the exercise habit datum. We defined a drinker as a person who drank at a frequency of “not less than 1 to 3 days per month,” and we defined a smoker as a person who “occasionally smoked.” Because a prefecture was the smallest unit serving as the target area for lifestyle and exercise habit data, these data at an SMSA level in the same prefecture took the same value.
Since the year in which these surveys are conducted differs, we obtained the data for the numbers of institutions with coronary CTA and with a CCL from the fiscal year 2011, the data for drinking habits from the fiscal year 2013, and all other data from the fiscal year 2010.
We calculated the ratios for AMI mortality based on standard mortality ratios (SMRs). The SMR is calculated by dividing the observed number by the expected number of AMI-related deaths. We calculated the ratios separately for men and women and for eight age subgroups: < 20, 21–30, 31–40, 41–50, 51–60, 61–70, 71–80, and > 80 years. We calculated the national AMI mortality rates in each subgroup, then multiplied the national AMI mortality rates by the population of each subgroup in each SMSA to obtain the expected number of AMI-related deaths.
Typically, when we handle spatial data, there is likely to be a correlation relating to the location of each region; i.e., a spatial autocorrelation. Moran’s I statistic was selected to quantify the degree of spatial autocorrelation. To model the spatial autocorrelation, we used a conditional autoregressive (CAR) model, which is a hierarchical Bayesian model . In the CAR model, random effects are set as error terms in general regression models by the class of CAR prior distribution, which is a type of Markov random field model. More concretely, random effects are represented by a set of full conditional distributions using a binary neighborhood matrix. If two regions are considered to be neighbors, their random effects are correlated. Additional details on the CAR model have been previously published . Among several kinds of CAR models, we chose the CAR Leroux model, which was reported in a simulation-based study to be consistently suitable for various spatial correlation scenarios . A Poisson distribution was assumed for the observed number of AMI-related deaths, and we set the expected number of AMI-related deaths as the offset variable in the CAR Leroux model. We used Markov chain Monte Carlo (MCMC) simulations with 110,000 iterations; the first 10,000 iterations were discarded as burn-in. We used Geweke’s diagnostic to check MCMC convergence . We estimated the relative risk (RR) and 95% Bayesian credible interval (CI) of each variable. According to a previous study , we considered that an association is not significant if the 95% CI of RR includes 1.
We used the widely applicable information criterion (WAIC) as a measure of the goodness-of-fit of the Bayesian statistical model; a model with a lower WAIC was considered the better-fit model . We compared the CAR Leroux model with a normal Poisson regression model without consideration for spatial autocorrelation using the WAIC.
We evaluated multicollinearity of covariates using the variance inflation factor (VIF) . Although the proportions of university graduates and the average per capita income were important factors, we removed them because they had the greatest VIF. All other variables had a VIF of < 2.5 and were entered into the CAR Leroux model.
The mean and standard deviation of the SMR of AMI were 1.039 and 0.313, respectively, indicating a wide range across SMSAs. Figure 1 shows the distribution of the SMR of AMI across SMSAs in Japan in 2011. Moran’s I statistic for the SMR was 0.306 (p < 0.001), which indicates strong spatial autocorrelation of the SMR in Japan. Figure 2 shows the number of institutions with coronary CTA per 100,000 population. Moran’s I statistic for the institutions with coronary CTA was 0.053 (p = 0.03), which indicates moderate and statistically significant spatial autocorrelation.
Table 1 shows descriptive data for the 349 SMSAs included in this study. The third quantiles of the numbers of all three medical resources (i.e., cardiologists, institutions with coronary CTA, and institutions with a CCL) were more than twice the first quantiles of those variables. This result indicates that these medical resources varied across SMSAs.
Results of CAR Leroux model
Table 2 shows the detailed results of the CAR Leroux model. The absolute value of Geweke’s diagnostic for each parameter was < 1.96, which means that all parameters used in the MCMC simulations converged .
Table 2 shows the significantly negative association between the number of institutions with coronary CTA and the SMR of AMI (RR = 0.900, 95% CI = 0.848–0.953). Across Japan, an increase of 1300 institutions with coronary CTA was associated with a 10% decrease in the SMR of AMI. In contrast, the SMR in SMSAs was not significantly associated with the number of either cardiologists or institutions with a CCL. A smoking habit was positively associated with the SMR of AMI, while a drinking habit was negatively associated. All other socioeconomic factors were not significantly associated with the SMR.
The WAIC was 5683.486 for the normal Poisson regression model without consideration for spatial autocorrelation and 2791.118 for the CAR Leroux model. Because the model with the lower WAIC was considered the better fit, this result indicates that the CAR Leroux model was much better fit than the normal Poisson regression model.
The main finding of this nationwide ecological study is the significant geographical association between the number of institutions with coronary CTA and the SMR of AMI. To the best of our knowledge, the current study is the first to clarify this association. This association can be said to be an important international finding in that it directly focuses on both coronary CTA allocation and regional clinical outcomes, considering that allocation of coronary CTA is receiving more attention in countries beyond Japan as in the updated NICE guidelines [7, 8]. This finding was obtained because of spatial statistics; application of the method to the topic is also a novelty of the current study.
Two major randomized controlled trials, PROMISE  and SCOT-HEART , showed several effects of coronary CTA on clinical outcomes. The former reported a reduction in death or nonfatal myocardial infarction in the coronary CTA group over the initial 12-month follow-up period (hazard ratio, 0.66; p = 0.049) compared with the functional testing group, although this difference was not recognized at the end of follow-up. The latter reported that the use of coronary CTA led to a 38% reduction in coronary heart disease-related death and nonfatal myocardial infarction over a median follow-up of 1.7 years compared with standard care alone. Although further investigations will be required to elucidate the association between coronary CTA and clinical outcomes , the current study demonstrated an impact of coronary CTA on clinical outcomes, particularly for AMI mortality, from the viewpoint of an ecological association.
The results of the current study suggest an important role of coronary CTA in specialized regionalized network systems. Compared with CCL, coronary CTA is easy to establish in rural medical areas and has been expected to be a gatekeeper for stratifying CAD patients . Considering this role of coronary CTA, how to allocate coronary CTA in each region is becoming increasingly more important. The construction of specific regional networks is reportedly essential for timely AMI treatment, particularly percutaneous coronary intervention [28,29,30,31]. There is a trend of treating AMI on regional networks worldwide . In contrast, the standing position of coronary CTA, a relatively new technique, remains unclear in regional networks. The findings of this study suggest that effective allocation of coronary CTA in regionalized networks would be important to improve regional outcomes for AMI. For example, it has been proposed that smaller facilities could be used as triage points , and coronary CTA equipment in these facilities might enhance triage performance.
Prediction of the future geographical distribution of coronary CTA is also important to effectively allocate coronary CTA. The spatial competition model would be helpful for this prediction . According to this model, even if the distribution of service resources is concentrated in cities with large populations for maximal profit, relocation of these resources toward smaller cities can occur because of the increasingly keener competition for profit in large cities. This model could be applied to medical devices such as CT, MRI, and positron emission tomography . Because coronary CTA is relatively new medical technique, its maldistribution could be improved in the future. For example, the spread of coronary CTA is still insufficient in the United Kingdom, and the maldistribution might be improved with an increase in coronary CTA . The frequency of using CT is higher in Japan than in other countries , and the results in Japan might be helpful for future allocation of coronary CTA in countries where coronary CTA is not sufficient.
While the number of institutions with coronary CTA was associated with the SMR, the number of institutions with a CCL was not significantly associated with the SMR. Considering that the number of coronary angiography (CAG) varies among medical institutions, we also evaluated the association between the number of CAG in each SMSA and the SMR of AMI. We found no significant association between them [see Additional file 1: Table A1]. This result could stem from the recent trend in which percutaneous coronary intervention is performed not only for surviving patients but also to improve patients’ quality of life, particularly patients with stable CAD as shown in the COURAGE trial .
The number of cardiologists was not significantly associated with the SMR of AMI. This result is consistent with a previous study in Japan . In contrast, however, a previous study in the US showed that the mortality risk in patients hospitalized for AMI in regions with a low density of cardiologists was higher than in regions with a high density of cardiologists . The discordance between our results and the study in the US may have two explanations. First, while the study in the US used patient-level data, the current study was based on SMSA-level data. Because of the difference in sample sizes, the previous study had higher statistical power to detect associations than did our study. Second, the previous study did not consider the influence of medical devices, particularly coronary CTA. As an additional analysis, we constructed a Bayesian model that excluded the number of institutions with coronary CTA and institutions with a CCL [see Additional file 1: Table A2]. This model showed a greater tendency toward a statistically significant negative association between the number of cardiologists and the SMR of AMI (RR 0.993; 95% CI 0.986–1.001) than the results shown in Table 2 (RR 0.997; 95% CI 0.988–1.004). If the previous study had considered the influence of coronary CTA, the association might be not significant.
One of the biggest advantages of the present study is its consideration of the spatial autocorrelation. Ignoring this parameter introduces a risk of incorrect statistical model estimation, particularly a risk of false-positive significance . As shown in Table 2, the posterior distribution of the spatial autocorrelation exhibited a moderate spatial autocorrelation (ρ 0.444; 95% CI 0.217–0.761). Given that the WAIC of the CAR Leroux model was approximately half that of the normal Poisson regression model, we should explicitly use spatial autocorrelation in statistical models that handle spatial data. This advantage of our approach is expandable in that it is applicable not only to the distribution of coronary CTA but also to the distribution of any medical resources.
This study has several limitations. First, it was an ecological analysis and did not account for patient-level factors that might contribute to AMI-related death, such as comorbidities, the severity of AMI in each patient, and other personal factors, resulting in possible ecological fallacy. Future researchers should also collect and analyze microdata regarding AMI-related deaths. Second, this study had a risk of information bias derived from the use of secondary data sources. For example, based on the survey questionnaire style, the number of institutions with coronary CTA included the number of institutions with cardiac MRI. In addition, we obtained the data regarding cardiologists by counting the number of cardiologists whose main medical department was “cardiology.” However, because of the questionnaire format, there might have been several physicians who responded that they were “internists” despite the fact that they were engaged in cardiology. Third, the results of this study might not be applicable to other countries. Notably, the frequency of using CT is higher in Japan than in other countries . In addition, because using CT is popular in Japan, much fewer institutions tend to perform functional stress testing compared to coronary CTA , while functional stress testing is a major non-invasive testing especially in the US. Therefore, we could not obtain detailed data of the number of functional stress testing institutions. Further data from other countries should be used to verify the generalizability of the results in this study.
Despite these limitations, this study provides the first evidence of an association between the geographical distribution of coronary CTA and regional AMI-related deaths. To the best of our knowledge, no other studies have assessed this association, and our findings will therefore help the government to determine how to efficiently allocate coronary CTA. We have enhanced the study’s accuracy by using spatial statistical models to consider the spatial autocorrelation.
The number of institutions with coronary CTA was significantly negatively associated with the SMR of AMI, while the SMR in each SMSA was not significantly associated with the number of either cardiologists or institutions with a CCL. Proactive allocation of coronary CTA is recommended, and it would be important to clarify the standing position of coronary CTA in regional networking for AMI treatment in the future.
acute myocardial infarction
coronary artery disease
cardiac catheterization laboratory
computed tomography angiography
Markov chain Monte Carlo
Ministry of Health, Labour and Welfare
magnetic resonance imaging
standardized mortality ratio
secondary medical service area
variance inflation factor
widely applicable information criterion
Finegold JA, Asaria P, Francis DP. Mortality from ischaemic heart disease by country, region, and age: statistics from World Health Organisation and United Nations. Int J Cardiol. 2013;168:934–45.
World Health Organization: Cardiovascular Diseases (CVDs): Fact Sheet No. 317. WHO. http://www.who.int/mediacentre/factsheets/fs317/en/. Accessed 13 Mar 2017.
Ministry of Health, Labour and Welfare: Overview of Vital Statistics in 2015. http://www.mhlw.go.jp/toukei/saikin/hw/jinkou/kakutei15/. Accessed 13 Mar 2017 (in Japanese).
Mowatt G, Cummins E, Waugh N, Walker S, Cook J, Jia X et al. Systematic review of the clinical effectiveness and cost-effectiveness of 64-slice or higher computed tomography angiography as an alternative to invasive coronary angiography in the investigation of coronary artery disease. Health Technol Assess. 2008;12:iii–iv, ix–143.
Budoff MJ, Dowe D, Jollis JG, Gitter M, Sutherland J, Halamert E, et al. Diagnostic performance of 64-multidetector row coronary computed tomographic angiography for evaluation of coronary artery stenosis in individuals without known coronary artery disease. Results from the prospective multicenter ACCURACY. J Am Coll Cardiol. 2008;52:1724–32.
Schroeder S, Achenbach S, Bengel F, Burgstahler C, Cademartiri F, De Feyter P, et al. Cardiac computed tomography: indications, applications, limitations, and training requirements: report of a writing group deployed by the Working Group Nuclear Cardiology and Cardiac CT of the European Society of Cardiology and the European Council of Nuclear Cardiology. Eur Heart J. 2008;29:531–56.
Moss AJ, Williams MC, Newby DE, Nicol ED. The updated NICE guidelines: cardiac CT as the first-line test for coronary artery disease. Curr Cardiovasc Imaging Rep. 2017;10:15.
Nicol E, Padley S, Roditi G, Roobottom C; on behalf of the British Society of Cardiovascular Imaging/British Society of Cardiovascular Computed Tomography. The challenge of national CT coronary angiography (CTCA) provision in response to NICE CG95 update 2016. http://www.bsci.org.uk/images/Documents/2016_BSCI_BSCCT_response_to_NICE_CG95.pdf. Accessed 24 Apr 2018.
Huber K, Gersh BJ, Goldstein P, Granger CB, Armstrong PW. The organization, function, and outcomes of ST-elevation myocardial infarction networks worldwide: current state, unmet needs and future directions. Eur Heart J. 2014;35:1526–32.
Park S, Lee J, Ikai H, Otsubo T, Imanaka Y. Decentralization and centralization of healthcare resources: investigating the associations of hospital competition and number of cardiologists per hospital with mortality and resource utilization in Japan. Health Policy. 2013;113:100–9.
Kirimura T, Nakaya T, Yano K. Building a spatio-temporal GIS database about boundaries of municipalities: municipality Map Maker for Web. Theory Appl GIS. 2011;19:83–92 (in Japanese).
The Japanese Circulation Society: Annual report of survey of clinical practice for cardiovascular diseases in 2011. http://www.j-circ.or.jp/jittai_chosa/jittai_chosa2011web.pdf. Accessed 22 Apr 2017 (in Japanese).
Ministry of Health, Labour and Welfare: Survey Questionnaires. http://www.mhlw.go.jp/toukei/chousahyo/index.html. Accessed 4 Apr 2017 (in Japanese).
General government statistics portal GL01010101. http://www.e-stat.go.jp/SG1/estat/eStatTopPortal.do. Accessed 2 Mar 2017 (in Japanese).
Kulkarni VT, Ross JS, Wang Y, Nallamothu BK, Spertus JA, Normand SLT, et al. Regional density of cardiologists and rates of mortality for acute myocardial infarction and heart failure. Circ Cardiovasc Qual Outcomes. 2013;6:352–9.
Dupre ME, George LK, Liu G, Peterson ED. Association between divorce and risks for acute myocardial infarction. Circ Cardiovasc Qual Outcomes. 2015;8:244–51.
MacNab YC. Hierarchical Bayesian modeling of spatially correlated health service outcome and utilization rates. Biometrics. 2003;59:305–16.
Lee D. A comparison of conditional autoregressive models used in Bayesian disease mapping. Spat Spatiotemporal Epidemiol. 2011;2:79–89.
Geweke J. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments. In: Bayesian statistics 4. Oxford: Oxford University Press; 1992. p. 169–93.
Helbich M, Leitner M, Kapusta ND. Lithium in drinking water and suicide mortality: interplay with lithium prescriptions. Br J Psychiatry. 2015;207:64–71.
Watanabe S. Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res. 2010;11:3571–94.
Allison PD. Multiple regression: a primer (research methods and statistics). 1st ed. Newbury Park: Pine Forge Press; 1999.
R Core Team. R: a language and environment for statistical computing. http://www.r-project.org/. Accessed 9 Apr 2017.
Douglas PS, Hoffmann U, Patel MR, Mark DB, Al-Khalidi HR, Cavanaugh B, et al. Outcomes of anatomical versus functional testing for coronary artery disease. N Engl J Med. 2015;372:1291–300.
The SCOT-HEART investigators. CT coronary angiography in patients with suspected angina due to coronary heart disease (SCOT-HEART): an open-label, parallel-group, multicentre trial. Lancet. 2015;385:2383–91.
Doris M, Newby DE. Coronary CT angiography as a diagnostic and prognostic tool: perspectives from the SCOT-HEART trial. Curr Cardiol Rep. 2016;18:18.
Maclachlan H, Thomas R, Langtree J, Hare C. Is there a role for a local inpatient CT coronary angiography service in selected patients with acute coronary syndrome? A cohort analysis of inpatient tertiary centre referrals for invasive coronary angiography. Open Heart. 2016;3:e000389.
Concannon TW, Kent DM, Normand SL, Newhouse JP, Griffith JL, Cohen J, et al. Comparative effectiveness of ST-segment-elevation myocardial infarction regionalization strategies. Circ Cardiovasc Qual Outcomes. 2010;3:506–13.
Manari A, Ortolani P, Guastaroba P, Casella G, Vignali L, Varani E, et al. Clinical impact of an inter-hospital transfer strategy in patients with ST-elevation myocardial infarction undergoing primary angioplasty: the Emilia-Romagna ST-segment elevation acute myocardial infarction network. Eur Heart J. 2008;29:1834–42.
Moyer P, Ornato JP, Brady WJ, Davis LL, Ghaemmaghami CA, Gibler WB, et al. Development of systems of care for ST-elevation myocardial infarction patients: the emergency medical services and emergency department perspective. Circulation. 2007;116:e43–8.
Solla DJF, De Mattos Paiva Filho I, Delisle JE, Braga AA, De Moura JB, De Moraes X, et al. Integrated regional networks for ST-segment-elevation myocardial infarction care in developing countries the experience of Salvador, Bahia, Brazil. Circ Cardiovasc Qual Outcomes. 2013;6:9–17.
Califf RM, Faxon DP. Need for centers to care for patients with acute coronary syndromes. Circulation. 2003;107:1467–70.
Newhouse JP. Geographic access to physician services. Annu Rev Public Health. 1990;11:207–30.
Matsumoto M, Koike S, Kashima S, Awai K. Geographic distribution of CT, MRI and PET devices in Japan: a longitudinal analysis based on national census data. PLoS ONE. 2015;10:e0126036.
OECD Health Statistics 2017. 2017. http://www.oecd.org/els/health-systems/health-data.htm. Accessed 26 Nov 2017.
Boden WE, O’Rourke RA, Teo KK, Hartigan PM, Maron DJ, Kostuk WJ, et al. Optimal medical therapy with or without PCI for stable coronary disease. N Engl J Med. 2007;356:1503–16.
Carl G, Kühn I. Analyzing spatial autocorrelation in species distributions using Gaussian and logit models. Ecol Model. 2007;207:159–70.
HK was responsible for the study concept and design. HK and SK were responsible for acquisition of data. All authors were responsible for analysis and interpretation of data. HK drafted the manuscript and SK, RS, and KO undertook for critical revision of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The data that support the findings of this study are available from the Japan Ministry of Health, Labour and Welfare but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Japan Ministry of Health, Labour and Welfare.
Consent for publication
Ethics and approval and consent to participate
The Ethics Committee of the Graduate School of Medicine and Faculty of Medicine at The University of Tokyo assessed and gave permission for this study (Assessment Number 10493).
This study was supported by a Health Labour Sciences Research Grant from the Ministry of Health, Labour and Welfare, Tokyo, Japan (H26 Research on Statistics and Information 001) and JSPS KAKENHI Grant Number JP16K00432.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kawaguchi, H., Koike, S., Sakurai, R. et al. Association between number of institutions with coronary computed tomography angiography and regional mortality ratio of acute myocardial infarction: a nationwide ecological study using a spatial Bayesian model. Int J Health Geogr 17, 13 (2018). https://doi.org/10.1186/s12942-018-0133-0