Using spatial analysis to demonstrate the heterogeneity of the cardiovascular drug-prescribing pattern in Taiwan
© Cheng et al; licensee BioMed Central Ltd. 2011
Received: 24 November 2010
Accepted: 24 May 2011
Published: 24 May 2011
Geographic Information Systems (GIS) combined with spatial analytical methods could be helpful in examining patterns of drug use. Little attention has been paid to geographic variation of cardiovascular prescription use in Taiwan. The main objective was to use local spatial association statistics to test whether or not the cardiovascular medication-prescribing pattern is homogenous across 352 townships in Taiwan.
The statistical methods used were the global measures of Moran's I and Local Indicators of Spatial Association (LISA). While Moran's I provides information on the overall spatial distribution of the data, LISA provides information on types of spatial association at the local level. LISA statistics can also be used to identify influential locations in spatial association analysis. The major classes of prescription cardiovascular drugs were taken from Taiwan's National Health Insurance Research Database (NHIRD), which has a coverage rate of over 97%. The dosage of each prescription was converted into defined daily doses to measure the consumption of each class of drugs. Data were analyzed with ArcGIS and GeoDa at the township level.
The LISA statistics showed an unusual use of cardiovascular medications in the southern townships with high local variation. Patterns of drug use also showed more low-low spatial clusters (cold spots) than high-high spatial clusters (hot spots), and those low-low associations were clustered in the rural areas.
The cardiovascular drug prescribing patterns were heterogeneous across Taiwan. In particular, a clear pattern of north-south disparity exists. Such spatial clustering helps prioritize the target areas that require better education concerning drug use.
Spatial analysis can identify inequities in disease distribution or in resources for healthcare utilization, resulting in financial incentives in payment policies [1–3]. Spatial analysis has taken giant steps forward in the recent decade, and Geographic Information Systems (GIS) software makes it possible to undertake a sophisticated visual approach to data analysis in medical issues [4, 5]. However, most studies using GIS have relied on their mapping capabilities rather than performance of statistical analyses [6–8]. Only when GIS are combined with spatial analytical methods can the result provide a helpful tool in the study of public health issues [9, 10].
Spatial analyses can be compared to societal norms or objectives. The traditional descriptive statistics fail to show the complexity of interactions between individual geographic units. In contrast, spatial analysis can reveal the co-variation of properties with a geo-space, since characteristics at proximal locations tend to be correlated, either positively or negatively, which leads to spatial autocorrelation . For example, a study shows that people are indirectly protected when more people in their neighborhood are vaccinated with cholera vaccine .
Another potential use of the spatial analysis is to examine patterns of drug use. Spatial analysis provides a method for delineating the complex nature of the associations between multiple levels of social structure and illicit drug use [13, 14]. These applications have important implications in policy-making, as the associations between geography and drug use could be useful for planning drug prevention programs.
Cardiovascular disease (CVD) remains the leading cause of death and disability, representing a major burden for health systems worldwide [15, 16], and impressive literature has documented geographic variation in the use of medication and procedures for treating cardiovascular disease in Western countries [17–21]. The cardiovascular medication-use patterns observed reflect the influence of a complex combination of demographic, social, economic, cultural, and environmental factors. Quality improvement programs based on geographic data may help to reduce the variation in quality of care .
A few previous studies have specifically examined the cardiovascular medication utilization in Taiwan [23, 24]; however, seldom is the spatial pattern of drug use addressed. Therefore, the objective of this study was to investigate prescribing patterns of cardiovascular medications based on the variations of geo-space in Taiwan. We hypothesized there would be significant geographic differences in prescribing cardiovascular medications among townships.
We conducted a cross-sectional study by using the non-sampled National Health Insurance Research Database (NHIRD) in Taiwan. This database is population-based and derived from the claims data of the National Health Insurance program, a mandatory-enrollment and single-payer system implemented in 1995 . Patients who had at least one outpatient prescription for CVD during the year 2004 were enrolled in this study. The dataset contained 22.13 million enrollees (about 97.5% of the total population in Taiwan) in 2004.
The data elements available for each patient in NHIRD contains minimal demographic characteristics, medical diagnoses (up to 5 for each admission, up to 3 for each ambulatory visit), procedures, expenditures, and all detailed prescriptions .
The prescriptions were coded according to the Anatomical Therapeutic Chemical coding system of the NHI Pharmaceutical subsidy, which was used as the interface for retrieving the pharmaceutical claims data. Five drug classes (and ATC numbers) were recorded, including agents acting on the rennin-angiotensin system (C09), called Group A hereafter; beta blocking agents (C07), called Group B hereafter; calcium channel blockers (C08), called Group C hereafter; diuretics (C03), called Group D hereafter; and antihypertensives (C02), called Group O hereafter. We presented the proportion of use of these cardiovascular drug classes among all patients based on geographic location. We used the "defined daily dose" (DDD) as a measure unit for drug consumption, which was developed by the World Health Organization (WHO) Collaborating Center for Drug Statistics Methodology (WHO, 2004). The DDD of a pharmaceutical substance represents the assumed average maintenance dose per day for a drug, when used for its main indication in adults. For example, the recommended maintenance dose of captopril, an angiotensin converting enzyme inhibitor, for hypertension is 50 mg to 75 mg per day, the DDD of this drug was define as 50 mg from WHO. The consumption of drugs was expressed as DDDs/1000 inhabitants per day, which may serve as an estimation of the proportion of the population receiving the drug treatment. The information on residents in each township in Taiwan was based on annual statistics of the Ministry of Interior .
We analyzed consumption of all classes of cardiovascular medications by DDD/1000 inhabitants per day across 352 townships. We presented the descriptive statistics of drug consumption by metropolitan areas, including the mean (± standard deviation). Coefficient of variation (COV, standard deviation divided by mean) is used to describe the variations of consumption in various areas. Bivariate analysis using the Student's t test was conducted to compare the south and north regions of Taiwan; a two-sided p value < 0.05 was considered statistically significant. Data were extracted from NHIRD, and descriptive analysis were performed with SAS version 9.1 (SAS Institute Inc., Cary, North Carolina)
where n indicates the total number of locations (352 townships), x i denotes the value of the variable of interest, X, at location i, x j denotes the observation at neighboring locations j, and is the sample average of X. w ij is the so-called spatial weights matrix (connectivity matrix), which defines spatial interaction across study regions. In general, w ij = 1 if location i and location j are neighboring, (share a common boundary); otherwise, w ij = 0. In the LISA analysis, if the test statistic is not significant at any sensible level, no spatial pattern is present in the areas; i.e., observations are spatially random. When it is significant, however, two possible patterns of local spatial clusters are likely to be exhibited: (1) When x i is higher than the average of the entire study area ( ) and so are its neighbors, a high-high (HH) association, a so-called hot spot, is indicated. (2) When both x i and its neighbors are lower than the average, the spatial tendency is low-low (LL), or a cold spot. The major limitation of Moran's I is that it only measures spatial clustering (spatial autocorrelation) at a global scale, but cannot detect important clusters at a local scale and determine the spatial patterns for the specific locations. To make spatial autocorrelation visually meaningful, the local Moran I is represented by cluster maps, in which the locations of significant spatial clusters are highlighted to identify the patterns of associations and p value < 0.01 was considered statistically significant.
Underlying spatial clusters can be geographic changes in risk factors, which introduce a degree of spatial autocorrelation in the outcome. In this study, spatial clusters of drug use reflect inefficiency in medical resources and utilization. Still, a variety of possible factors could lead to such inefficiency and they are location-specific. However, if there exist common factors creating diffusion or local clusters of medical inefficiency, it would be helpful and necessary to identify the locations of risk clusters and the extent of their clustering. Based on the sizes of risk areas and the significance of clusters, important risk factors shared by neighboring areas could be uncovered. Thus, from the perspective of policy making and research, priority should be given to these clustered areas.
We used the geographic information system software ArcGIS v9.2 to combine drug prescription data with Taiwan digital maps at the township level. The data was then processed by spatial analysis software GeoDa  to create the spatial weights matrix and to calculate the LISA value for each township.
Geographic map and geographic distribution of drugs consumption
Figure 1 illustrates the administrative division across 352 townships in Taiwan, including 5 major metropolitan areas, 2 secondary metropolitan areas, 206 rural areas, and 55 aboriginal townships.
The distribution of cardiovascular medications consumption by urbanicity, data presented as mean ± (SD)a
Taipei-Keelung metropolitan area (N = 44)
Taoyuan-Chungli metropolitan area (N = 12)
Hsinchu secondary metropolitan Area (N = 10)
Taichung-Changhwa metropolitan area (N = 18)
Chiayi secondary metropolitan area (N = 4)
Tainan metropolitan area (N = 14)
Kaohsiung metropolitan area (N = 28)
Rural townships area (N = 167)
Aboriginal townships area (N = 55)
Total (N = 352)
We did not observe any significant difference in cardiovascular medication consumption between northern and southern areas: Group A (21.2 ± 33.6, 17.8 ± 40.5; p = 0.67), Group B (11.8 ± 17.2, 10.1 ± 17.5; p = 0.51), Group C (20.8 ± 34.9, 19.0 ± 38.3; p = 0.86), Group D (6.3 ± 10.9, 5.6 ± 13.5; p = 0.83) and Group O (1.5 ± 3.1, 1.5 ± 3.4; p = 0.70), respectively.
LISA statistics for cardiovascular drugs
There were wide variations in utilization of all drug classes across the 352 townships in Taiwan. From the first glimpse of the LISA figures, under-utilization appeared to be the predominant spatial pattern for cardiovascular disorders across all drug classes: the total number of high-high associations (hot spots) amounts to 35, while low-low associations (cold spots) reaches to 85. Such spatial variations carried three statistical meanings: (1) variations were non-random; (2) variations were statistically significant; (3) variations exhibited effects of neighboring interactions.
North-south discrepancy distribution
Focusing on the west of Taiwan, where more than 97 percent of the population resides, patterns of drug use exhibited more high-high associations clustered in the southern (n = 30) than northern (n = 2) townships. Thus the southern region had higher-than-average utilization of cardiovascular drugs than other parts of Taiwan. However, at the same time, the number of low-low associations also presented in southern (n = 42) more than northern areas (n = 33).
The spatial analysis of five category cardiovascular medications by urbanicity, data presented as N (%)
Taipei-Keelung metropolitan area (N = 44)
Hsinchu secondary metropolitan Area (N = 10)
Tainan metropolitan area (N = 14)
Kaohsiung metropolitan area (N = 28)
Rural townships area (N = 167)
Aboriginal townships area (N = 55)
Total (N = 352)c
By drug class, both Group C and Group D had the highest amount of diverse clusters across the country, followed by Group B, Group A, and Group O. Hot spot areas concentrated mostly in the Kaohsiung and rural areas among all drug classes (Table 2). In contrast, cold spot areas showed up in rural and aboriginal townships areas, accounting for a particularly high portion of group A, followed by group B, group C, group O, and group D.
Our study showed a variability of cardiovascular prescription patterns with regard to geographic differences in Taiwan. The analysis of both HH and LL associations strongly indicated that there exists a discrepancy in drug use between the north and the south of Taiwan, and the southern areas exhibited spatial clustering of unusual utilization in cardiovascular drugs. There are two possible explanations for this finding. First, the disease prevalence may be different between the north and south in Taiwan as a result of urbanization levels. The prevalence of cardiovascular disease, such as hypertension and coronary heart disease, varies between the urban and rural areas. Some countries reported the prevalence was higher for people living in the urban than rural areas [31, 32], while others, including Taiwan, presented evidence to the contrary [33, 34]. Second, there were urban-rural differences in the appropriateness of disease treatment and the inclination of patients to receive medical care. Previous studies have documented regional variation in quality of care in cardiovascular disease [22, 35]. Patients hospitalized for acute myocardial infarction (AMI) in rural hospitals were less likely to receive recommended intervention than in urban hospitals. In addition, inappropriate or lower utilization of health care could be associated with lower socioeconomic status of patients in rural areas [36, 37].
Utilization of cardiovascular drugs appeared to be associated with different classes of pharmacological drugs among geographic regions. We found that underutilization of all classes of drugs was more likely to be observed than overutilization, and the former was scattered mainly in the southern, rural, and aboriginal townships. (Figure 2) It seemed to suggest a wide variation in prescribing behavior. Drug Group A, Group B, and Group C were generally the most widely used drugs for the treatment of cardiovascular diseases ; however, there were large variations across townships with regard to the use of drugs. A study conducted in France reported the socio-economic status at the area level influenced the utilization of specialty care . Differences between general practitioners and other specialties have been described in relation to the use of cardiovascular medicines [40–42]. For example, compared with cardiologists, family practitioners and general internists probably underutilize ACE inhibitors due to less knowledge about heart failure and poor adherence to guidelines . Among seventeen townships where ACE inhibitors were underutilized, we found that a clinic for primary care is the major medical resource in these areas . In addition, there are three aboriginal townships have the High-High hotspot. They tend to over-utilize Group B, C, D and O but not ACEI inhibitors. This might be attributable to the fact that the aboriginal residents only have access to general practitioners and medical treatment by specialists is lacking. Accordingly, the spatial patterns of drug use observed in this study suggested that there was a wide variation in physicians' prescribing behavior across Taiwan, to the extent that medical resources were not consistently utilized among 352 townships.
Utilization of cardiovascular drugs also associated with the disease patterns appeared to be different between townships. Some studies showed that rural residents are at higher risk for the burden of major diseases such as cardiovascular disorder and the severity of their complications [33, 44]. In contrast, other studies showed that higher prevalence rates of cardiovascular risk factors were found in urban than rural areas due to diet and stressful life styles [31, 45]. Drug choice should be based on clinical features of patients. For instance, rennin-angiotensin system agents have greatly improved clinical outcomes in patients with renal complications,  or diuretics are first-line agents for elderly hypertensive patients without other complications .
Spatial analytical methods recognize that spatial concordance between two proximal entities is closely connected by intrinsic causal relationship or hidden confounders. Although the geographical variations in drug use were widely documented, the concept of spatial dependence has received little attention in pharmacoepidemiology. It was our intention in this study to establish that "spatial perspectives in pharmacoepidemiology" would be an effective method in detecting the variation of drug utilization. That is, unusual drug prescriptions are not randomly located across locations but tend to be concentrated within certain areas.
LISAs are simply local derivations or disaggregation of global measures of spatial autocorrelation, and allow one to detail the information of local variations in spatial autocorrelation (e.g., hot spots and cold spots) for further study consideration. This method is useful to indicate spatial correlation of similar values at location i and its neighborhood. This means that location i and its neighborhood can both have values above the average value (hot spots), or both can have values below the average (cold spots). Such local variations in spatial autocorrelation are statistically significant. Departure from the classical assumption of statistical independence, the spatial clusters represent that data are not randomly distributed over the global region and can occur because there is a contagious factor to the underlying process that is commonly shared by the clustered region. Lack of the factor that triggers spatial dependence does not give rise to spatial clustering.
The primary strength of this study was that we incorporated GIS and tools for spatial analysis to explore the role of geography in drug utilization. In so doing, we were able to divert our examination from the traditional framework, which only considers individual decision-making, into the spatial realm. Epidemiologic investigations are commonly involved with data from contiguous or proximal geographic areas, where neighboring locations may share similar prescribing patterns due to, for example, inhabitants' lifestyle behaviors or environmental exposures. Although the trend in utilization of cardiovascular drugs in Taiwan was investigated, little was known regarding its geographical pattern across townships. The spatial analysis presented here may add essential information to pharmacoepidemiological studies and policy-making. On the one hand, spatial analysis allows us to examine spatial processes underlying the prescribing behavior of physicians and provides the first inquiry into whether spatial components should be accounted for in modeling to avoid erroneous inference . On the other hand, given the limited resources, the government can focus their policy efforts to target areas to improve drug utilization. For the reasons described above, our approach should be applicable worldwide. Developing countries, for which medical resources are increasingly disproportionately distributed, may see more distinct prescribing patterns than those for Taiwan or other developed countries.
Although spatial analysis with Moran's I is useful for this study, it does have limitations of which all should be aware. First, Moran's I (and other global indicators) is not powerful enough to differentiate well between a random pattern and a pattern without substantial spatial variations. Second, Moran's I is conditional, tied to how we define localities through a spatial weights matrix w ij . Different definitions of the matrix can lead to different values of spatial dependence. Finally, the measures of spatial dependence are subject to the modifiable areal unit problem. The problem arises due to the arbitrary choice and the modifiable nature of areal units in spatial analysis, in the sense that they can be aggregated to form units of different sizes of spatial arrangements. The degree of spatial association, as measured by Moran's I, changes with different levels of aggregation in areal units. While we are aware of the limitations of the spatial autocorrelation index, our results in this paper can be interpreted as a broad indication of the presence and magnitude of spatial dependence.
Other potential limitations of our study also should be noted. First, the utilization pattern of cardiovascular medications could be different as urbanization changes. Second, we only analyzed the classes of cardiovascular drugs based on the first three digits of ATC which have DDD assigned. Although we were able to depict the spatial patterns of drug use, whether the knowledge of new drugs was attributable to the geographic clustering requires a detailed classification of drug classes based on ATC. Third, we might underestimate the utilization of group A, B, D because complex medication (i.e., a compound made of two or more ingredients) without DDD assigned were not estimated in this study. The most complex medication are combinations with group D (i.e., group A plus group D or group B plus group D) to improve patient's compliance. Fourth, this study only analyzed one year of data, but the prescribing pattern could extend across years. Fifth, the cardiovascular medication pattern did not correspond to disease distribution. Most patients need two more classes of medication to control their disease, except simple essential hypertension. Finally, the fact that we have observed significant spatial autocorrelation in our data does not determine whether these are true spatial effects or are spurious in the sense that they can be attributed entirely to patterns in other variables, such as income or individual characteristics (such as age, gender, etc). If we control for all of these other factors and the spatial variable remains significant, then we have evidence that the patterns are consistent with a neighborhood effect and are not solely attributable to socioeconomic characteristics of these areas. Exploring theses sources of spatial patterns and the mechanisms that may be driving them is an obvious extension of the analysis presented here, and more detailed studies are necessary to better understand the factors that influence prescribing patterns.
The cardiovascular medication prescribing patterns were heterogeneous across Taiwan. In particular, a clear pattern of north-south disparity exists. The spatial variability of prescribing patterns documented in this study shows that priority-setting in research may heavily depend on the neighborhood association.
List of abbreviations
Acute myocardial infarction
Anatomical Therapeutic Chemical
Geographic Information Systems
Local Indicators of Spatial Association
National Health Insurance Research Database
World Health Organization
This research was funded by the Taiwan National Science Council (NSC 96-2320-B-006-028-MY3) and Multidisciplinary Center of Excellence for Clinical Trial and Research (DOH100-TD-B-111-002), Department of Health, Executive Yuan, Taiwan. The funding source had no role in the design, analysis, interpretation, or reporting of results or in the decision to submit the manuscript for publication.
- Allen NB, Holford TR, Bracken MB, Goldstein LB, Howard G, Wang Y, Lichtman JH: Geographic variation in one-year recurrent ischemic stroke rates for elderly Medicare beneficiaries in the USA. Neuroepidemiology. 2010, 34 (2): 123-129. 10.1159/000274804.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashton CM, Petersen NJ, Souchek J, Menke TJ, Yu HJ, Pietz K, Eigenbrodt ML, Barbour G, Kizer KW, Wray NP: Geographic variations in utilization rates in Veterans Affairs hospitals and clinics. N Engl J Med. 1999, 340 (1): 32-39. 10.1056/NEJM199901073400106.View ArticlePubMedGoogle Scholar
- Domanski M, Antman EM, McKinlay S, Varshavsky S, Platonov P, Assmann SF, Norman J: Geographic variability in patient characteristics, treatment and outcome in an International Trial of Magnesium in acute myocardial infarction. Control Clin Trials. 2004, 25 (6): 553-562. 10.1016/j.cct.2004.08.005.View ArticlePubMedGoogle Scholar
- Gatrell AC, Bailey TC: Interactive spatial data analysis in medical geography. Social Science & Medicine. 1996, 42 (6): 843-855. 10.1016/0277-9536(95)00183-2.View ArticleGoogle Scholar
- Parson RL, Frost JD: Interactive analysis of spatial subsurface data using GIS-based tool. Journal of Computing in Civil Engineering. 2000, 14 (4): 215-222. 10.1061/(ASCE)0887-3801(2000)14:4(215).View ArticleGoogle Scholar
- Gilberg K, Laouri M, Wade S, Isonaka S: Analysis of medication use patterns:apparent overuse of antibiotics and underuse of prescription drugs for asthma, depression, and CHF. J Manag Care Pharm. 2003, 9 (3): 232-237.PubMedGoogle Scholar
- Jia H, Ried LD, Wang X, Damush TM, Young LJ, Cameon RH, Williams LS: Geographic variation in poststroke depression among veterans with acute stroke. J Rehabil Res Dev. 2008, 45 (7): 1027-1035. 10.1682/JRRD.2007.06.0090.View ArticlePubMedGoogle Scholar
- Mathur AK, Ashby VB, Sands RL, Wolfe RA: Geographic variation in end-stage renal disease incidence and access to deceased donor kidney transplantation. Am J Transplant. 2010, 10 (4 Pt 2): 1069-1080.View ArticlePubMedGoogle Scholar
- Rezaeian M, Dunn G, St Leger S, Appleby L: Geographical epidemiology, spatial analysis and geographical information systems: a multidisciplinary glossary. J Epidemiol Community Health. 2007, 61 (2): 98-102. 10.1136/jech.2005.043117.View ArticlePubMedPubMed CentralGoogle Scholar
- Gesler W: The uses of spatial analysis in medical geography: a review. Soc Sci Med. 1986, 23 (10): 963-973. 10.1016/0277-9536(86)90253-4.View ArticlePubMedGoogle Scholar
- Anselin L: Local Indicators of Spatial Association - Lisa. Geographical Analysis. 1995, 27 (2): 93-115.View ArticleGoogle Scholar
- Emch M, Ali M, Root ED, Yunus M: Spatial and environmental connectivity analysis in a cholera vaccine trial. Soc Sci Med. 2009, 68 (4): 631-637. 10.1016/j.socscimed.2008.11.025.View ArticlePubMedPubMed CentralGoogle Scholar
- Paulozzi LJ, Ryan GW: Opioid analgesics and rates of fatal drug poisoning in the United States. Am J Prev Med. 2006, 31 (6): 506-511. 10.1016/j.amepre.2006.08.017.View ArticlePubMedGoogle Scholar
- Latkin C, Glass GE, Duncan T: Using geographic information systems to assess spatial patterns of drug use, selection bias and attrition among a sample of injection drug users. Drug Alcohol Depend. 1998, 50 (2): 167-175. 10.1016/S0376-8716(98)00025-8.View ArticlePubMedGoogle Scholar
- Wysowski DK, Governale LA, Swann J: Trends in outpatient prescription drug use and related costs in the US - 1998-2003. Pharmacoeconomics. 2006, 24 (3): 233-236. 10.2165/00019053-200624030-00003.View ArticlePubMedGoogle Scholar
- Jackevicius CA, Tu K, Filate WA, Brien SE, Tu JV, Canadian Cardiovascular Outcomes R: Trends in cardiovascular drug utilization and drug expenditures in Canada between 1996 and 2001. Canadian Journal of Cardiology. 2003, 19 (12): 1359-1366.PubMedGoogle Scholar
- Harris KM, Pastorius CA, Duval S, Harwood E, Henry TD, Carabello BA, Hirsch AT: Practice variation among cardiovascular physicians in management of patients with mitral regurgitation. Am J Cardiol. 2009, 103 (2): 255-261. 10.1016/j.amjcard.2008.09.065.View ArticlePubMedGoogle Scholar
- Kumar A, Fonarow GC, Eagle KA, Hirsch AT, Califf RM, Alberts MJ, Boden WE, Steg PG, Shao M, Bhatt DL, et al: Regional and practice variation in adherence to guideline recommendations for secondary and primary prevention among outpatients with atherothrombosis or risk factors in the United States: a report from the REACH Registry. Crit Pathw Cardiol. 2009, 8 (3): 104-111.View ArticlePubMedGoogle Scholar
- Laskey W, Spence N, Zhao X, Mayo R, Taylor R, Cannon CP, Hernandez AF, Peterson ED, Fonarow GC: Regional differences in quality of care and outcomes for the treatment of acute coronary syndromes: an analysis from the get with the guidelines coronary artery disease program. Crit Pathw Cardiol. 2010, 9 (1): 1-7.View ArticlePubMedGoogle Scholar
- Levin S, Welch VL, Bell RA, Casper ML: Geographic variation in cardiovascular disease risk factors among American Indians and comparisons with the corresponding state populations. Ethn Health. 2002, 7 (1): 57-67. 10.1080/13557850220146993.View ArticlePubMedGoogle Scholar
- O'Connor GT, Quinton HB, Traven ND, Ramunno LD, Dodds TA, Marciniak TA, Wennberg JE: Geographic variation in the treatment of acute myocardial infarction: the Cooperative Cardiovascular Project. JAMA. 1999, 281 (7): 627-633. 10.1001/jama.281.7.627.View ArticlePubMedGoogle Scholar
- Baldwin LM, MacLehose RF, Hart LG, Beaver SK, Every N, Chan L: Quality of care for acute myocardial infarction in rural and urban US hospitals. Journal of Rural Health. 2004, 20 (2): 99-108. 10.1111/j.1748-0361.2004.tb00015.x.View ArticlePubMedGoogle Scholar
- Lin TL, Yang YHK, Tang CH, Wu HM, Cheng HW: Impact of introduction of angiotensin II antagonist on the anti hypertensive utilization in Taiwan. Journal of Food and Drug Analysis. 2008, 16 (4): 21-29.Google Scholar
- Dai CL, Chu HF, Lin HD, Ng YY, Wu SC: Use of angiotensin converting enzyme inhibitors and angiotensin II receptor blockers in patients with type 2 diabetic chronic kidney disease. Taiwan Journal of Public Health. 2007, 26: 7-Google Scholar
- National Health Insurance Research Database. [http://w3.nhri.org.tw/nhird/en/Background.html]
- Hsiao FY, Yang CL, Huang YT, Huang WF: Using Taiwan's national health insurance research databases for pharmacoepidemiology research. Journal of Food and Drug Analysis. 2007, 15 (2): 99-108.Google Scholar
- Statistical Yearbook of Interior. [http://www.moi.gov.tw/english/index.aspx]
- Library information management. [http://gis.ascc.net/NDASupport/index-2.html]
- Anselin L, Florax RJGM: New directions in spatial econometrics. 1995, Berlin; New York: SpringerView ArticleGoogle Scholar
- Anselin L, Syabri I, Kho Y: GeoDa: An introduction to spatial data analysis. Geographical Analysis. 2006, 38 (1): 5-22. 10.1111/j.0016-7363.2005.00671.x.View ArticleGoogle Scholar
- Chadha SL, Gopinath N, Shekhawat S: Urban-rural differences in the prevalence of coronary heart disease and its risk factors in Delhi. Bull World Health Organ. 1997, 75 (1): 31-38.PubMedPubMed CentralGoogle Scholar
- Liu Z, Albanese E, Li S, Huang Y, Ferri CP, Yan F, Sousa R, Dang W, Prince M: Chronic disease prevalence and care among the elderly in urban and rural Beijing, China - a 10/66 Dementia Research Group cross-sectional survey. BMC Public Health. 2009, 9: 394-10.1186/1471-2458-9-394.View ArticlePubMedPubMed CentralGoogle Scholar
- Weeks WB, Wallace AE, Wang S, Lee A, Kazis LE: Rural-urban disparities in health-related quality of life within disease categories of veterans. Journal of Rural Health. 2006, 22 (3): 204-211. 10.1111/j.1748-0361.2006.00033.x.View ArticlePubMedGoogle Scholar
- Pan WH, Chang HY, Yeh WT, Hsiao SY, Hung YT: Prevalence, awareness, treatment and control of hypertension in Taiwan: results of Nutrition and Health Survey in Taiwan (NAHSIT) 1993-1996. Journal of Human Hypertension. 2001, 15 (11): 793-798. 10.1038/sj.jhh.1001268.View ArticlePubMedGoogle Scholar
- Sheikh K, Bullock C: Urban-rural differences in the quality of care for Medicare patients with acute myocardial infarction. Archives of Internal Medicine. 2001, 161 (5): 737-743. 10.1001/archinte.161.5.737.View ArticlePubMedGoogle Scholar
- Gornick ME, Eggers PW, Reilly TW, Mentnech RM, Fitterman LK, Kucken LE, Vladeck BC: Effects of race and income on mortality and use of services among Medicare beneficiaries. New England Journal of Medicine. 1996, 335 (11): 791-799. 10.1056/NEJM199609123351106.View ArticlePubMedGoogle Scholar
- Chan L, Hart LG, Goodman DC: Geographic access to health care for rural medicare beneficiaries. Journal of Rural Health. 2006, 22 (2): 140-146. 10.1111/j.1748-0361.2006.00022.x.View ArticlePubMedGoogle Scholar
- Chen CL, Chen LW, Yang WC: The influences of Taiwan's generic grouping price policy on drug prices and expenditures: Evidence from analysing the consumption of the three most-used classes of cardiovascular drugs. BMC Public Health. 2008, 8-Google Scholar
- Chaix B, Boelle PY, Guilbert P, Chauvin P: Area-level determinants of specialty care utilization in France: a multilevel analysis. Public Health. 2005, 119 (2): 97-104. 10.1016/j.puhe.2004.05.006.View ArticlePubMedGoogle Scholar
- Chin MH, Friedmann PD, Cassel CK, Lang RM: Differences in generalist and specialist physicians' knowledge and use of angiotensin-converting enzyme inhibitors for congestive heart failure. Journal of General Internal Medicine. 1997, 12 (9): 523-530. 10.1046/j.1525-1497.1997.07105.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Edep ME, Shah NB, Tateo IM, Massie BM: Differences between primary care physicians and cardiologists in management of congestive heart failure: Relation to practice guidelines. Journal of the American College of Cardiology. 1997, 30 (2): 518-526. 10.1016/S0735-1097(97)00176-9.View ArticlePubMedGoogle Scholar
- Baker DW, Hayes RP, Massie BM, Craig CA: Variations in family physicians' and cardiologists care for patients with heart failure. American Heart Journal. 1999, 138 (5): 826-834. 10.1016/S0002-8703(99)70006-2.View ArticlePubMedGoogle Scholar
- statistic for medical resource. [http://www.doh.gov.tw/CHT2006/DM/DM2_2.aspx?now_fod_list_no=9865&class_no=440&level_no=4]
- Jones A, Goza F: Rural, Suburban and Urban Differences in the Self-Diagnosis of Coronary Heart Disease in the United States. Journal of Biosocial Science. 2008, 40 (6): 895-909.View ArticlePubMedGoogle Scholar
- Chiu HC, Lee TK, Mau LW, Liu HW: Prevalence of cardiovascular risk factors of elderly persons in Taiwan. Kaohsiung J Med Sci. 2002, 18: 8-Google Scholar
- Ma TKW, Kam KKH, Yan BP, Lam YY: Renin-angiotensin-aldosterone system blockade for cardiovascular diseases: current status. British Journal of Pharmacology. 2010, 160 (6): 1273-1292.View ArticlePubMedPubMed CentralGoogle Scholar
- Chobanian AV, Bakris GL, Black HR, Cushman WC, Green LA, Izzo JL, Jones DW, Materson BJ, Oparil S, Wright JT, et al: Seventh Report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure. Hypertension. 2003, 42 (6): 1206-1252. 10.1161/01.HYP.0000107251.49515.c2.View ArticlePubMedGoogle Scholar
- Anselin L: Local indicators of spatial association - LISA. Geographical Analysis. 1995, 27: 22-Google Scholar
- Jerrett M, Burnett RT, Goldberg MS, Sears M, Krewski D, Catalan R, Kanaroglou P, Giovis C, Finkelstein N: Spatial analysis for environmental health research: Concepts, methods, and examples. Journal of Toxicology and Environmental Health-Part A. 2003, 66 (16-19): 1783-1810.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/11/380/prepub