Warning: fopen(/home/virtual/epih/journal/upload/ip_log/ip_log_2024-03.txt): failed to open stream: Permission denied in /home/virtual/lib/view_data.php on line 83 Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 84 The association between smoking or passive smoking and cardiovascular diseases using a Bayesian hierarchical model: based on the 2008-2013 Korea Community Health Survey

The association between smoking or passive smoking and cardiovascular diseases using a Bayesian hierarchical model: based on the 2008-2013 Korea Community Health Survey

Article information

Epidemiol Health. 2017;39.e2017026
Publication date (electronic) : 2017 June 22
doi : https://doi.org/10.4178/epih.e2017026
Graduate School of Public Health, Seoul National University, Seoul, Korea
Correspondence: Ho Kim Graduated School of Public Health, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea E-mail: hokim@snu.ac.kr
Received 2017 January 2; Accepted 2017 January 17.

Abstract

OBJECTIVES

Smoking and passive smoking have been extensively reported as risk factors of cardiovascular morbidity and mortality. Despite the biological mechanisms underlying the impact of hazardous chemical substances contained in tobacco in cardiovascular diseases (CVD), studies investigating the association between smoking and passive smoking with morbidity are at an inchoate stage in Korea. Therefore, this study aimed to estimate the risks of smoking and passive smoking on cardiovascular morbidity at the national and regional levels.

METHODS

This study calculated sex-standardized and age-standardized prevalence of CVD and smoking indices in 253 community health centers (si/gun/gu) in Korea using the 2008-2013 Korea Community Health Survey data. Furthermore, a Bayesian hierarchical model was used to estimate the association of smoking and passive smoking with the prevalence of CVD from the national and regional community health centers.

RESULTS

At the national level, smoking was significantly associated with stroke (relative risk [RR], 1.060) and hypertension (RR, 1.016) prevalence, whilst passive smoking at home and work were also significantly associated with prevalence of stroke (RR, 1.037/1.013), angina (RR, 1.016/1.006), and hypertension (RR, 1.010/1.004). Furthermore, the effects of smoking and passive smoking were greater in urban-industrial areas than in rural areas.

CONCLUSIONS

The findings of this study would provide grounds for national policies that limit smoking and passive smoking, as well as regionally serve as the basis for region-specific healthcare policies in populations with high CVD vulnerability.

INTRODUCTION

Smoking and passive smoking have extensively been documented as risk factors of morbidity and mortality by a number of diseases [1-4]. Among more than 4,000 types of chemical substances contained in a cigarette, polycyclic aromatic hydrocarbons and oxidizing gases cause cardiac toxicity [1], having adverse effects on various cardiovascular diseases (CVD), including coronary artery disease, ischemic stroke, non-traumatic subarachnoid hemorrhage [5,6].

The World Health Organization reported that CVD-related death (1.69 million) is the leading cause of smoking-related death. In a study on 360 thousand US men, smoking has been identified as a key risk factor for lung cancer, chronic heart disease, and stroke [3]. In addition, a Norwegian study revealed that smokers are 2.74 times more likely to have a stroke and 6.74 times more likely to die from it than are non-smokers [4]. Similar to the effects of smoking, passive smoking has also been reported to be the cause of respiratory and CVD, various cancers, and premature deaths in non-smoking adults [7]. Other studies have suggested that secondhand smoking in non-smokers increases their risk of developing a heart disease by 25 to 30% [8-10].

In South Korea (hereafter Korea), several studies have reported that smoking increases the risk of CVD-related death in adults [11-13]. As of 2012, the number of smoking-related deaths was 4,148 for stroke and 3,858 for ischemic heart disease, which are ranked the second and third highest causes of death by an illness respectively [14]. Furthermore, it also has been reported that 26.7% of all CVD-related deaths in men and women adults in Korea are caused by smoking [13]. Although adult smoking rate in Korea has consistently been on a decline since 1998 [15], past smoking history may have an impact on death [11]. Moreover, medical expenses and patient population related to hypertension, heart disease, and cerebrovascular disease are continuously on the rise [16], suggesting that smoking is a health hazard in terms of CVD and mortality requiring much attention.

As shown here, though there have been much research investigating the association between smoking and death by CVD in Korea, studies on smoking and morbidity is relatively lacking. Hence, this study aimed at quantitatively analyzing the association between smoking and the prevalence of various CVD (i.e., hypertension, stroke, myocardial infarction, angina), and to evaluate the relative risk (RR) of smoking on the prevalence of CVD using data from a large-scale cross-sectional study with more than 1.5 million participants. We also examined the effects of passive smoking on the prevalence of CVD.

MATERIALS AND METHODS

Data

Data from the 2008-2013 Korea Community Health Survey (KCHS) involving 253 community health centers were used for this study [17]. The KCHS was conducted by the Korea Centers for Disease Control and Prevention, and data collected from adults aged 19 years or older via interviews. The sample was extracted from an average of 900 adults per community (si/gun/gu) based on the type of housing within each dong/eup/myeon. The primary sample region was obtained using a probability proportional to size systematic sampling, after which the secondary sample families were selected [18]. A total of 1,567,930 people were included in the KCHS data used in this study (2008: 200,800; 2009: 220,258; 2010: 230,712; 2011: 229,229; 2012: 229,226; 2013: 228,781), and the participants were randomly extracted every year. Using this data, we calculated the prevalence of 4 CVD (hypertension in adults ≥ 30 years, stroke in adults ≥ 50 years, myocardial infarction in adults ≥ 40 years, angina in adults ≥ 40 years), smoking-related regional parameters (current smoking rate, passive smoking rate at home, and passive smoking rate at work), body mass index (BMI), and monthly drinking rate by community (si/gun/gu), after reflecting gender/age distributions based on the year of study. The prevalence of CVD was measured, based on a self-reported response regarding a physician’s diagnosis, whilst BMI and monthly drinking rates were also measured with self-reported data.

Statistical analysis

A generalized linear mixed model (GLMM) was used to analyze the association among regional and national smoking indices and prevalence of CVD that were repeatedly measured from 2008-2013. The prevalence of 4 CVD (hypertension in adults ≥ 30 years, stroke in adults ≥ 50 years, myocardial infarction in adults ≥ 40 years, and angina in adults ≥ 40 years) was set as the response variables. A Shapiro-Wilk test revealed that the prevalence of all 4 diseases did not follow a normal distribution, but instead showed a gamma distribution. Hence, a GLMM that assumes gamma distribution for the response variables was used. A log function was used as the link function. Mean BMI and monthly drinking rate per community (si/gun/gu) for 2008-2013, were set as the respective confounding variables and covariate, and a linear adjustment was made for annual effects to consider the effects of time. Smoking indices (current smoking rate, passive smoking rate at home, passive smoking rate at work) were used as the response variables, and we predicted associations with response variables with equal confounding variables and covariates.

Furthermore, a Bayesian hierarchical model was used with a GLMM to measure the effects of smoking indices of each community (si/gun/gu) and adjust for the correlations according to distance. The associations between community (si/gun/gu)-specific smoking indices and CVD were examined with random intercept and random slope, and the prior distributions of random intercept and random slope were assumed to follow a normal distribution with mean zero and an assumption that the inverse of variance follows a gamma distribution (parameter= 1, 5*10-5). Spatial correlation was adjusted for, with random intercept using the Besag-York-Mollie method [19], and all random intercepts and random slopes were assumed to be independent. Based on existing studies, the prior distribution of fixed effects was assumed a multivariate normal, which assumes variance to be a Gaussian Markov random field [20], and the posterior credibility interval (CI) was estimated based on the prior distribution. All analyses were performed using the R and R version 3.2.2 (https://cran.r-project.org/bin/windows/base/old/3.2.2/) for integrated nested Laplace approximation (INLA) [21]. More detailed formula and prior conditions are reported in Appendix 1.

RESULTS

Table 1 shows the descriptive statistics for the variables used in this study. Hypertension over 30 years was the most prevalent (18.2%) throughout the 6 years of study (2008-2013), while angina over 40 years was the least prevalent (1.9%). The distributions of prevalence varied widely across the country. The average smoking rate as an interest variable during 6 years was 25.2% (14.9 to 33.4%), with an average passive smoking rate at home of 11.8% (2.3 to 37.3%) and an average passive smoking rate at work of 29.2% (6.0 to 76.7%).

Descriptive statistics for regional prevalence of CVD and confounders

Table 2 shows the top and bottom 5 community health centers in terms of the prevalence of CVD and smoking indices (6-year average). Hypertension over 30 years was generally prevalent in the Gangwon province (21.72 to 22.35%) and less prevalent in community health centers in the Gyeongsang province and Jeolla province (14.03 to 14.72%). Stroke over 50 years was generally prevalent in community health centers in the North Gyeongsang province (4.68 to 5.33%) and less prevalent in the South Jeolla province (1.93 to 2.12%). Myocardial infarction over 40 years tended to be prevalent in community health centers in Seoul and the Gyeonggi province (2.12 to 2.25%) and less prevalent in community health centers in North Jeolla province (0.78 to 0.83%). Angina over 40 years tended to be prevalent in community health centers in Busan (2.88 to 3.25%) and less prevalent in the North Jeolla province (1.03 to 1.08%). In terms of smoking parameters, smoking rate did not vary by region. Community health centers in Jeju Island showed higher passive smoking rates at home (19.88 to 21.80%) while community health centers in Ulsan showed higher passive smoking rates at work (40.93 to 41.70%).

Higher and lower 5 community health centers (si/gun/gu) of cardiovascular disease prevalence

The national average RR of smoking indices on the prevalence of CVD is shown in Table 3, whilst the regional average (RR) for top and bottom 5 regions is shown in Table 4. All RR shown in Tables 3 and 4 represent RR when smoking parameters increase by 5%. As shown in Table 3, smoking rate was most highly associated with stroke (RR, 1.060; 95% CI, 1.022 to 1.100), and it was also significantly associated with hypertension (RR, 1.016; 95% CI, 1.004 to 1.029). Passive smoking at home was most strongly associated with a rise in the prevalence of stroke (RR, 1.037; 95% CI, 1.023 to 1.051), and its RR for hypertension, stroke, and angina (but not myocardial infarction) was statistically significant (RR, 1.010 to 1.037). Passive smoking at work had the greatest effect on stroke (RR, 1.013; 95% CI, 1.007 to 1.019) and also had a significant effect on hypertension (RR, 1.004; 95% CI, 1.002 to 1.006).

Relative risk of smoking/passive smoking (at home)/passive smoking (at work) by cardiovascular diseases (per 5% smoking or passive smoking rate)

Higher and lower 5 community health centers (si/gun/gu) of relative risk (RR) to smoking/passive smoking (at home)/passive smoking (at work)

Table 4 shows the top and bottom three community health centers in terms of the RRs of smoking indices, whilst Figures 1-3 are national maps showing the effects of community-specific smoking rate and passive smoking rate at home and work on the prevalence of CVD. RR of smoking rate for hypertension was high in community health centers in Gyeonggi and Gangwon provinces (RR, 1.042 to 1.046) and low in community health centers in North Jeolla and Gyeonsang provinces (RR, 0.985 to 0.996). Furthermore, smoking was strongly associated with stroke in community health centers in North Gyeongsang province (RR, 1.095 to 1.099) but weakly associated with stroke in community health centers in Busan and Ulsan (RR, 1.033 to 1.034). The RR of smoking for myocardial infarction was high in community health centers in Daejeon (RR, 1.029 to 1.029), but there were no regional differences in terms of low RR (RR, 0.973 to 0.983). The effects of smoking on angina was high in the Gyeonggi region (RR, 1.040 to 1.035) and low in the southern region of the Korean peninsula (North Gyeonsang, North Jeolla, and Jeju Island) (RR, 0.975 to 0.989). The distributions of RR of passive smoking at home and work were similar to those of smoking, and the community rankings were similar as well. Exposure to passive smoking at home (RR, 1.003 to 1.037) had greater impact on all 4 CVD than did exposure to passive smoking at work (RR, 0.998 to 1.013). In addition, the RR of all smoking indices tended to be higher in the Seoul metropolitan area (Figures 1-3).

Figure 1.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and smoking.

Figure 2.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and passive smoking at home.

Figure 3.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and passive smoking at work.

DISCUSSION

This study analyzed the associations between regional smoking indices (smoking rate and passive smoking rate at home and work) and 4 CVD (hypertension in adults ≥ 30 years, stroke in adults ≥ 50 years, myocardial infarction in adults ≥ 40 years, angina in adults ≥ 40 years) with the KCHS data using a Bayesian hierarchical model. The results showed that smoking rate was significantly associated with stroke (RR, 1.060) and hypertension (RR, 1.016), and that passive smoking was significantly associated with stroke (RR, 1.013 to 1.037), angina (RR, 1.006 to 1.016), and hypertension (RR, 1.004 to 1.010) at the national level. Furthermore, exposure to passive smoking at home (RR, 1.003 to 1.037) had a greater impact on CVD than did exposure to passive smoking at work (RR, 0.998 to 1.013).

In addition to the national-level analysis, this study also analyzed the effects of smoking regionally in 253 regions, based on the assumption that specific regional features (geographical location, environment, income, population matrix) may have varying effects on the prevalence of CVD and may alter the effects of smoking on CVD. Our findings showed that RR of smoking varied across regions, suggesting that the effects of smoking on CVD varied across regions. The effects of smoking and passive smoking on all CVD were high in urban-industrial regions.

Nicotine, the primary ingredient of tobacco, is a sympathetic stimulant that acts on the central and peripheral nervous systems to stimulate the release of catecholamine and other neurotransmitters. It leads to cardiovascular consequences, such as elevating the heart rate, blood pressure, and cardiac output. Furthermore, it leads to a rise in low-density lipoprotein and decrease in high-density lipoprotein by mobilizing free fatty acids, thereby intensifying vasoconstriction and accelerating the progression of blood epithelial cell injury and atherosclerosis [1]. Carbon monoxide that is inhaled through smoking and passive smoking binds with hemoglobin to induce hypoxia, which increases the number of red blood cells and subsequently, blood viscosity thereby, having a direct impact on thrombosis and atherosclerosis [22]. Through these actions, smoking not only causes coronary artery diseases but also induces structural damage to arterial walls, and has also been associated with ischemic stroke caused by atherosclerosis and nontraumatic subarachnoid hemorrhage caused by formation and rupture of aneurysm [23].

This study presents data that support these biological causes. Although smoking indices were not statistically significantly associated with myocardial infarction and angina, the RR of smoking was above 1 for all CVDs at the national level (RR, 1.004 to 1.010). These results are in line with those reported by US and Norwegian studies on the relationship between smoking and heart diseases [3,4]. According to a foreign 20-year retrospective study, the greater the degree of exposure to passive smoking (cotinine concentration) was, the greater the risk for coronary artery disease and stroke [24]. Another study involving experimental and control groups revealed that the risk for myocardial infarction increases by 68% among passive smokers [25]. Furthermore, a meta-analysis of studies in 6 countries found that the RR for CVDs was between 1.2 and 1.3 [26]. In one study, passive smoking was found to be more of a health-hazard than direct smoking because smoke from the tip of a cigarette is more harmful than filtered smoke [27]. Our study also found similar results, and our results showed that the RR of passive smoking at home on angina (RR, 1.016) was higher than that of direct smoking (RR, 1.007).

Nevertheless, this study has a few limitations. First, this study merely analyzed associations and not the causal relationships between smoking indices and CVD. It is difficult to consider the temporal relations of smoking indices and CVD with KCHS data because this was a cross-sectional study [18]. In addition, this study did not consider risk factors at the individual level, so the results should be interpreted as a regional association and not as causal relations. Second, we were able to predict regions that were vulnerable to smoking and passive smoking but could not identify the biological and epidemiological causes of their vulnerability. More specifically, this study was able to verify that people with CVD in urban regions (including industrial regions) are relatively more vulnerable to smoking and passive smoking, but we could not shed light on the fundamental hypothesis and cause after adjusting for gender and age. We expect follow-up studies to evaluate the cause and to explain these results via diverse investigations and advanced analytical methods. Finally, this study did not consider social indicators. CVD may be influenced by social factors, such as income level and environment, but this study could not take these into considerations due to a lack of appropriate data. However, we partially compensated for this by examining potential regional features through performing a longitudinal analysis using a random effect model [28,29].

Despite a few limitations, this study offers several benefits. First, the large sample size of more than 1.5 million subjects increased the accuracy of our analysis in identifying the associations between CVD and smoking indices. Particularly, most of the 4 CVD studied were strongly associated with the smoking indices, with stroke having the greatest association with smoking and passive smoking (RR, 1.013 to 1.060). Second, by using an advanced statistical method, we were able to predict RRs not only at the national but also regional levels. The INLA method [20] predicts post-distributions as does the Markov Chain Monte Carlo method that is widely used in Bayesian statistical analysis, but the INLA is quicker and is easier for expanding distributions [30]. This study is valuable in that it simultaneously analyzed national and regional level risks by performing a Bayesian hierarchical analysis, considering regional correlations with the INLA.

Furthermore, the findings of this study can be used as the grounds for public healthcare policies that limit smoking and passive smoking nationwide to promote people’s health. Particularly, they will contribute to devising region-specific healthcare policies for regions populated by groups of people who may be vulnerable to CVD. In other words, there is a need to implement more active measures to limit smoking and to respond to the problems of indoor passive smoking in regions with high populations of people vulnerable to CVD, such as the elderly and financially unstable groups.

Acknowledgements

This work was supported by the Korea Ministry of Environment as Climate Change Correspondence Program (grant no. 20140013 10007).

Notes

The authors have no conflicts of interest to declare for this study.

SUPPLEMENTARY MATERIAL

Supplementary material (Korean version) is available at http://www.e-epih.org/.

References

1. Bullen C. Impact of tobacco smoking and smoking cessation on cardiovascular risk and disease. Expert Rev Cardiovasc Ther 2008;6:883–895.
2. Centers for Disease Control and Prevention (CDC). Smoking-attributable mortality, years of potential life lost, and productivity losses--United States, 2000-2004. MMWR Morb Mortal Wkly Rep 2008;57:1226–1228.
3. Kuller LH, Ockene JK, Meilahn E, Wentworth DN, Svendsen KH, Neaton JD. Cigarette smoking and mortality. MRFIT Research Group. Prev Med 1991;20:638–654.
4. Håheim LL, Holme I, Hjermann I, Leren P. Risk factors of stroke incidence and mortality. A 12-year follow-up of the Oslo Study. Stroke 1993;24:1484–1489.
5. Kannel WB, D’Agostino RB, Belanger AJ. Fibrinogen, cigarette smoking, and risk of cardiovascular disease: insights from the Framingham Study. Am Heart J 1987;113:1006–1010.
6. Ockene IS, Miller NH. Cigarette smoking, cardiovascular disease, and stroke: a statement for healthcare professionals from the American Heart Association. American Heart Association Task Force on Risk Reduction. Circulation 1997;96:3243–3247.
7. International Agency for Research on Cancer (IARC). IARC monographs on the evaluation of carcinogenic risks to humans. 2004. [cited 2017 Jun 20]. Available from: http://monographs.iarc.fr/ENG/Monographs/vol106/index.php.
8. Glantz SA, Parmley WW. Passive smoking and heart disease. Epidemiology, physiology, and biochemistry. Circulation 1991;83:1–12.
9. Law MR, Morris JK, Wald NJ. Environmental tobacco smoke exposure and ischaemic heart disease: an evaluation of the evidence. BMJ 1997;315:973–980.
10. Thun M, Henley J, Apicella L. Epidemiologic studies of fatal and nonfatal cardiovascular disease and ETS exposure from spousal smoking. Environ Health Perspect 1999;107 Suppl 6:841–846.
11. Lee EH, Park SK, Ko KP, Cho IS, Chang SH, Shin HR, et al. Cigarette smoking and mortality in the Korean Multi-center Cancer Cohort (KMCC) study. J Prev Med Public Health 2010;43:151–158. (Korean).
12. Meng KH. Smoking-attributable mortality among Korean adults. Korean J Epidemiol 1988;10:138–145. (Korean).
13. Jee SH, Lee JK, Kim IS. Smoking-attributable mortality among Korean adults: 1981-2003. Korean J Epidemiol 2006;28:92–99. (Korean).
14. Jung KJ, Yun YD, Baek SJ, Jee SH, Kim IS. Smoking-attributable mortality among Korean adults, 2012. J Korea Soc Health Inform Stat 2013;38:36–48. (Korean).
15. Korean Statistical Information Service. Current smoking rate 2008-2016. [cited 2017 Jul 20]. Available from: http://kosis.kr/statHtml/statHtml.do?orgId=117&tblId=DT_H_SM&conn_path=I2&language=en.
16. Kim J, Son M. National health insurance statistical yearbook Seoul: Health Insurance Review and Assessment Service and National Health Insurance Service; 2014. p. 62–67. (Korean).
17. Kim YT, Choi BY, Lee KO, Kim H, Chun JH, Kim SY, et al. Overview of Korean Community Health Survey. J Korean Med Assoc 2012;55:74–83. (Korean).
18. Lee WH, Choo JY, Son JY, Kim H. Association between long-term exposure to air pollutants and prevalence of cardiovascular disease in 108 South Korean communities in 2008-2010: a cross-sectional study. Sci Total Environ 2016;565:271–278.
19. Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math 1991;43:1–20.
20. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Series B Stat Methodol 2009;71:319–392.
21. Martino S, Rue H. Implementing approximate Bayesian inference using integrated nested Laplace approximation: a manual for the INLA program. 2008. [cited 2017 Jun 20]. Available from: http://www.bias-project.org.uk/GMRFCourse/inla-program.pdf.
22. Ambrose JA, Barua RS. The pathophysiology of cigarette smoking and cardiovascular disease: an update. J Am Coll Cardiol 2004;43:1731–1737.
23. Kurth T, Kase CS, Berger K, Gaziano JM, Cook NR, Buring JE. Smoking and risk of hemorrhagic stroke in women. Stroke 2003;34:2792–2795.
24. Whincup PH, Gilg JA, Emberson JR, Jarvis MJ, Feyerabend C, Bryant A, et al. Passive smoking and risk of coronary heart disease and stroke: prospective study with cotinine measurement. BMJ 2004;329:200–205.
25. Wells AJ. Heart disease from passive smoking in the workplace. J Am Coll Cardiol 1998;31:1–9.
26. He J, Vupputuri S, Allen K, Prerost MR, Hughes J, Whelton PK. Passive smoking and the risk of coronary heart disease--a metaanalysis of epidemiologic studies. N Engl J Med 1999;340:920–926.
27. Dobson R. Smoke from cigarette tip is more toxic than main inhaled smoke. BMJ 2005;331:1425.
28. Laird NM, Ware JH. Random-effects models for longitudinal data. Biometrics 1982;38:963–974.
29. Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika 1986;73:13–22.
30. Blangiardo M, Cameletti M. Spatial and Spatio-temporal Bayesian models with R-INLA New York: John Wiley & Sons; 2015.

Appendices

Appendix 1. Statistical analysis formulation

Gamma distribution

f(yi)=1Γ(v)(vyμ)vexp(-vyμ),y0,μ>0,v>0Yi~Γ(μ,v)log(v)~logGamma(1, 0.01)

where Yi is a observed values with i= 1, 2, … N. We assumed Gamma distribution as a distribution of Yi with parameter μ and ν. And we set the prior distribution of log(ν) as logGamma distribution with parameter 1, and 0.01 each.

Integrated nested Laplace approximation (INLA)

p(yβ, ψ) = i=1np(yiβi, ψ)p(βψ)=(2π)-n2Q(ψ)12exp(-12β'Q(ψ)β)

The INLA algorithm, proposed by Rue al. (2009) [1] and Martino & Rue [2] is a deterministic algorithm for Bayesian inference rather than simulation based, MCMC (Markov Chain Monte Carlo). β is a set of regression 16 coefficient parameter and Q(ψ) is a sparse precision matrix known as Gaussian Markov random field (GMRF).

Regression model

log(μij)~(β0+b0j+b2j)+(β1+b1j)Xij+β-0,1Zb0j~Normal(0,σb02),b1j~Normal(0,σb12)1/σv02~logGamma(1,0.00005),1/σv12~logGamma(1,0.00005)

β0 is an intercept of regression model and X is an interesting variable (smoking, indirect smoking indices, each). Bayesian hierarchical model is used with hyper parameters;ψ=(σv02, σv12 etc).

Besag-York-Mollie (BYM) model

b2j=uj+vjuju-j~Normal(1Ni kajkuk,sj2),vj=Normal(0,σv2)

Considering N areas, each characterized by a set of neighbors Nj. And ajk is 1 if areas j and k are neighbors and 0 otherwise.

sj2=σu2/Nj is the variance for the same area. vj is the area-specific effect modeled as exchangeable.

References

1. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Series B Stat Methodol 2009;71:319–392. 2. Martino S, Rue H. Implementing approximate Bayesian inference using integrated nested Laplace approximation: a manual for the INLA program. 2008. [cited 2017 Jun 20]. Available from: http://www.bias-project.org.uk/GMRFCourse/inla-program.pdf.

Article information Continued

Figure 1.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and smoking.

Figure 2.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and passive smoking at home.

Figure 3.

Geographical distribution of relative risks (RR). The associations between 4 cardiovascular diseases (A) hypertension over 30 years, (B) stoke over 50 years, (C) myocardial infarction over 40 years, and (D) angina over 40 years and passive smoking at work.

Table 1.

Descriptive statistics for regional prevalence of CVD and confounders

Diseases Period Mean (SD) Min 25% Median 75% Max
CVD (yr) Hypertension (> 30) 2008-2013 18.2 (2.4) 10.0 16.5 18.3 19.9 25.6
Stroke (> 50) 3.3 (1.1) 0.4 2.6 3.2 4.0 7.7
Myocardial infarction (> 40) 1.4 (0.6) 0.1 1.0 1.4 1.8 4.6
Angina (>40) 1.9 (0.7) 0.2 1.4 1.8 2.3 4.4
Confounders (%) Smoking rate 2008-2013 25.2 (2.9) 14.9 23.3 25.2 27.2 33.4
Passive smoking (at home) rate 2009-2011, 2013 11.8 (4.1) 2.3 9.1 11.3 14.1 37.3
Passive smoking (at work) rate 2009-2011, 2013 29.2 (9.3) 6.0 23.8 28.6 34.9 76.7
Drinking rate per month 2008-2013 56.5 (5.4) 32.9 53.2 57.1 60.4 68.8
Body mass index 2008-2013 23.0 (0.7) 22.2 22.8 23.0 23.1 34.5

CVD, cardiovascular disease; SD, standard deviation; Min, minimum; Max, maximum.

Table 2.

Higher and lower 5 community health centers (si/gun/gu) of cardiovascular disease prevalence

Si/do Si/gun/gu Mean Standard deviation
Prevalence (%)
 Hypertension (>30 yr) Higher 5 communities Gangwon Sokcho 22.35 1.74
Incheon Dong 22.10 2.21
Gangwon Goseong 21.98 0.63
Gangwon Samcheok 21.88 1.05
Gangwon Yeongwol 21.72 2.13
Lower 5 communities Gyeongnam Hapcheon 14.03 1.50
Jeonbuk Imsil 14.05 2.85
Gyeongbuk Seongju 14.48 1.78
Jeonnam Naju 14.50 1.13
Gyeongnam Sancheong 14.72 1.15
 Stroke (> 50 yr) Higher 5 communities Gyeongbuk Gumi-Gumi 5.33 1.71
Chungnam Asan 4.88 1.36
Gangwon Yeongwol 4.75 1.79
Gyeongbuk Gimcheon 4.72 0.81
Gyeongbuk Andong 4.68 0.41
Lower 5 community health centers Busan Suyeong 1.87 0.97
Jeonnam Jangheung 1.93 0.75
Jeonnam Yeongam 1.97 0.76
Jeonnam Wando 2.00 0.58
Jeonnam Gangjin 2.12 0.38
 Myocardial infarction (> 40 yr) Higher 5 community health centers Jeonnam Damyang 2.38 1.24
Jeju Jeju-Jeju 2.28 0.50
Gyunggi Goyang-Ilsanseo 2.25 1.18
Seoul Nowon 2.20 0.84
Seoul Dongjak 2.12 0.44
Lower 5 community health centers Gyeongnam Hamyang 0.68 0.27
Daegu Dalseong 0.72 0.40
Gyeongbuk Sungju 0.72 0.26
Jeonbuk Sunchang 0.78 0.12
Jeonbuk Jangsu 0.83 0.32
 Angina (> 40 yr) Higher 5 community health centers Busan Yeonje 3.25 0.49
Gwangju Nam 2.97 0.69
Chungbuk Cheonan 2.92 0.80
Busan Dong 2.90 0.18
Busan Gijang 2.88 0.96
Lower 5 community health centers Gyeongbuk Sungju 0.87 0.12
Jeonbuk Jinan 1.03 0.36
Ulsan Dong 1.07 0.34
Jeonbuk Wanju 1.08 0.58
Chungbuk Jecheon 1.20 0.70
Smoking index (%)
 Smoking rate Higher 5 community health centers Gangwon Taebaek 31.62 1.61
Gyeonggi Dongducheon 29.72 1.77
Chungbuk Eumseong 29.67 2.80
Busan Jung 29.37 2.19
Gyeonggi Bucheon -Ojeong 29.35 1.19
Lower 5 community health centers Gyeonggi Gwacheon 17.37 1.68
Gyeonggi Seongnam-Bundang 17.50 2.01
Seoul Seocho 19.57 1.24
Gyeonggi Yongin-Suji 20.00 1.70
Chungnam Gyeryong 20.77 1.58
 Passive smoking (at home) rate Higher 5 community health centers Incheon Ongjin 24.45 8.65
Jeju Dongbu 21.80 3.56
Jeju Seogwipo (west) 19.90 6.24
Jeju Jeju-Seo(west) 19.88 3.74
Jeonbuk Gimje 19.10 8.13
Lower 5 community health centers Busan Suyeong 6.28 2.47
Gyeongnam Yangsan 6.28 0.75
Jeonnam Gurye 6.73 3.40
Gyeonggi Gwacheon 6.78 2.44
Gangwon Wonju 6.78 2.42
 Passive smoking (at work) rate Higher 5 community health centers Chungnam Dangjin 42.83 24.56
Ulsan Nam 41.70 10.22
Ulsan Jung 41.40 9.40
Ulsan Ulju 40.93 12.45
Gyeongnam Changwon-Masan 40.83 5.25
Lower 5 community health centers Gyeongbuk Uiseong 15.70 7.03
Jeonnam Gurye 15.83 8.81
Jeonnam Goheung 15.85 13.30
Jeonbuk Jangsu 17.58 13.55
Jeonnam Hampyeong 18.55 4.71

Table 3.

Relative risk of smoking/passive smoking (at home)/passive smoking (at work) by cardiovascular diseases (per 5% smoking or passive smoking rate)

Smoking Exposure to passive smoking
Home Work
Hypertension (> 30 yr) 1.016 (1.004, 1.029) 1.010 (1.006, 1.014) 1.004 (1.002, 1.006)
Stroke (> 50 yr) 1.060 (1.022, 1.100) 1.037 (1.023, 1.051) 1.013 (1.007, 1.019)
Myocardial infarction (> 40 yr) 1.004 (0.958, 1.051) 1.003 (0.986, 1.021) 0.998 (0.991, 1.006)
Angina (> 40 yr) 1.007 (0.966, 1.049) 1.016 (1.001, 1.032) 1.006 (1.000, 1.013)

Values are presented relative risk (95% credibility interval).

Table 4.

Higher and lower 5 community health centers (si/gun/gu) of relative risk (RR) to smoking/passive smoking (at home)/passive smoking (at work)

Smoking
Exposure passive smoking
Home
Work
Si/do Si/gun/gu RR Si/do Si/gu RR Si/do Si/gu RR
Hypertension (>30yr) Higher 5 community health centers Gyeonggi Gapyeong 1.046 Gyeonggi Gapyeong 1.036 Gangwon Samcheok 1.030
Gangwon Samcheok 1.045 Incheon Dong 1.036 Gyeonggi Gapyeong 1.030
Incheon Dong 1.045 Gangwon Samcheok 1.035 Incheon Dong 1.030
Gangwon Sokcho 1.042 Gangwon Sokcho 1.033 Chungcheongnam Yeongi 1.027
Gyeonggi Paju 1.042 Chungnam Yeongi 1.032 Gangwon Sokcho 1.027
Lower 5 community health centers Gyeongnam Hapcheon 0.985 Gyeongnam Hapcheon 0.977 Gyeongnam Hapcheon 0.972
Busan Gangseou 0.991 Busan Gangseou 0.984 Busan Gangseou 0.979
Jeonbuk Imsil 0.994 Jeonbuk Imsil 0.987 Jeonbuk Imsil 0.982
Jeonbuk Gochang 0.996 Jeonbuk Gochang 0.988 Jeonbuk Gochang 0.984
Gyeongbuk Seongju 0.996 Jeju Seogwipo (west) 0.989 Gyeongbuk Seongju 0.985
Stroke (>50yr) Higher 5 community health centers Gyeongbuk Gumi-Gumi 1.099 Gyeongbuk Gumi-Gumi 1.068 Gyeongbuk Gumi-gumi 1.040
Daejeon daedeok 1.099 Gyeonggi Gapyeong 1.066 Gyeonggi Gapyeong 1.040
Gyeonggi Gapyeong 1.098 Daejeon Daedeok 1.066 Daejeon daedeok 1.039
Gyeongbuk Cheongdo 1.095 Daegu Dalseong 1.064 Gyeongbuk Cheongdo 1.038
Jeonnam Asan 1.095 Jeonnam Asan 1.063 Jeonnam Asan 1.037
Higher 5 community health centers Gyeongnam Hadong 1.029 Gyeongnam Hadong 1.005 Gyeongnam Hadong 0.986
Chungnam Wando 1.032 Incheon Ongjin 1.010 Ulsan Buk 0.990
Busan Suyeong 1.033 Chungnam Wando 1.011 Busan Suyeong 0.991
Ulsan Dong 1.034 Ulsan Buk 1.012 Chungnam Wando 0.991
Ulsan Buk 1.034 Busan Suyeong 1.013 Ulsan Dong 0.992
Myocardial infarction (> 40 yr) Higher 5 community health centers Jeonnam Yeosu 1.040 Jeonnam Yeosu 1.039 Jeonnam Yeosu 1.035
Daejeon Dong 1.029 Daejeon Dong 1.029 Daejeon Dong 1.024
Daejeon Daedeok 1.029 Daejeon Daedeok 1.028 Daejeon Daedeok 1.024
Chungnam Gongju 1.027 Chungnam Gongju 1.026 Chungnam Gongju 1.022
Jeju Jeju-Jeju 1.026 Jeju Jeju-Jeju 1.025 Jeonnam Gokseong 1.021
Lower 5 community health centers Gyeongbuk Seongju 0.973 Gyeongbuk Seongju 0.972 Gyeongbuk Seongju 0.967
Daegu Dalseong 0.974 Daegu Dalseong 0.973 Daegu Dalseong 0.969
Gyeongnam Hamyang 0.977 Gyeongnam Hamyang 0.977 Gyeongnam Hamyang 0.972
Gyeonggi Sungnam-Sujung 0.980 Gyeonggi Sungnam-Sujeong 0.980 Gyeonggi Sungnam-Sujeong 0.975
Chungnam Yesan 0.983 Busan Gangseo 0.983 Chungnam Yesan 0.978
Angina (> 40 yr) Higher 5 community health centers Chungnam Cheonan 1.047 Chungnam Cheonan 1.053 Chungnam Cheonan 1.042
Jeonbuk Gimje 1.041 Jeonbuk Gimje 1.047 Jeonbuk Gimje 1.038
Gyeonggi Gapyeong 1.040 Gyeonggi Gapyeong 1.047 Gyeonggi Gapyeong 1.037
Busan Gijang 1.038 Busan Gijang 1.045 Busan Gijang 1.034
Gyeonggi Ansung 1.035 Gyeonggi Ansung 1.042 Gyeonggi Ansung 1.032
Lower 5 community health centers Gyeongbuk Seongju 0.975 Gyeongbuk Seongju 0.985 Gyeongbuk Seongju 0.976
Jeonbuk Wanju 0.985 Jeju Jeju (west) 0.994 Jeonbuk Wanju 0.985
Gyeonggi Sungnam-Sujung 0.986 Jeonbuk Wanju 0.995 Gyeonggi Sungnam-Sujung 0.986
Jeju Jeju (west) 0.986 Gyeonggi Sungnam-Sujung 0.996 Jeju Jeju (west) 0.986
Ulsan Dong 0.989 Incheon Ganghwa 0.999 Incheon Ganghwa 0.989