We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Risk factors for interstitial lung disease: a 9-year Nationwide population-based study
BMC Pulmonary Medicinevolume 18, Article number: 96 (2018)
Understanding the risk factors that are associated with the development of interstitial lung disease might have an important role in understanding the pathogenetic mechanism of interstitial lung disease as well as prevention. We aimed to determine independent risk factors of interstitial lung disease development.
This was a retrospective cohort study with nationwide population-based 9-year longitudinal data. We selected subjects who were aged > 40 years at cohort entry and with a self-reported history of cigarette smoking. Cases were selected based on International Classification of Diseases codes. A cohort of 312,519 subjects were followed until December 2013. We used Cox regression analysis to calculate the hazard ratios (HRs) for interstitial lung disease development.
Interstitial lung disease developed in 1972 of the 312,519 subjects during the 9-year period. Smoking (HR: 1.2; 95% confidence interval [CI]: 1.1–1.4), hepatitis C (HR: 1.6; 95% CI: 1.1–2.3), history of tuberculosis (HR: 1.5; 95% CI: 1.1–1.9), history of pneumonia (HR: 1.6; 95% CI: 1.3–2.0), and chronic obstructive pulmonary disease (HR: 1.8; 95% CI: 1.6–2.1), men (HR: 1.9; 95% CI: 1.7–2.1) were significantly associated with the development of interstitial lung disease. The risk of interstitial lung disease development increases with age, and the risk was 6.9 times higher (95% CI: 5.9–8.0) in those aged over 70 than in their forties.
Smoking, hepatitis C, history of tuberculosis, history of pneumonia, chronic obstructive pulmonary disease, male sex, and older age were significantly associated with interstitial lung disease development.
If bilateral reticular or reticulonodular opacities are found on chest radiography, interstitial lung disease (ILD) is suspected. ILD is not an uncommon disease . In 2012, prevalence estimates of ILD with fibrosis ranged from 42.7–63 per 100.000 population in USA, and 1.25–23.4 per 100.000 population in Europe .
Patients with ILD are often asymptomatic until the lesion progresses significantly. Interstitial pulmonary abnormality (ILA), which is considered to be an early lesion of interstitial lung disease, has a higher mortality rate than patients without ILA  as well as ILD . Furthermore, in patients with ILD, a study using health insurance claim data found that lung cancer incidence is higher than chronic obstructive pulmonary disease (COPD) . ILD could have significant impact on health.
Three different population-based studies determined that cigarette smoking is a risk factor for lung parenchymal as well as interstitial abnormalities in addition to airway abnormalities [6,7,8]. However, previous studies examined risk factors from a limited sample size, and most of them used case-control methodology and did not consider many of the potential confounding variables as risk factors for ILD.
Determining the risk factors that are associated with the development of ILD might have an important role in understanding the pathogenetic mechanism of ILD, early diagnosis, adequate treatment and prevention.
The purpose of this study was to identify the independent risk factors for the development of ILD in a 9-year follow-up longitudinal population-based study.
The National Health Insurance Service (NHI) covers more than 99% of all Korean residents and includes all health claim data, including diagnostic codes, procedures, prescription drugs, patient personal information, and hospital information. There is one health insurance system with a unique resident registration number for each citizen; therefore, duplication of subjects can be avoided. This study used data from the National Health Insurance Service-National Sample Cohort(NHIS-NSC) 2002–2013 , which was released by the KNHIS in 2015. It includes all medical claims filed from January 2002 to December 2013 for 1,099,094 nationally representative randomly selected subjects, accounting for approximately 2.2% of the entire population in the KNHIS in 2002. The data were produced by the KNHIS using a systematic sampling method to generate a representative sample of all 46,605,433 Korean residents in 2002. The cohort population underwent biennial medical evaluations through the NHI Corporation between January 1, 2002, and December 31, 2013.
This was a nationwide population-based 9-year longitudinal study. We included subjects with a self-reported smoking history, who were ≥ 40 years old, and with co-morbidities diagnosed before the index date. Health examination data confirmed self-reported cigarette smoking history. The health examination data every 2 years was linked with the cohort data.
We did not include patients with ILD from 2002 to 2004 to exclude preexisting cases of ILD during the two medical evaluations and 1 year of follow-up. We also excluded patients with ILD who did not visit the clinic 30 days after the index date. The final sample included 312,519 subjects. Then, we identified patients with newly diagnosed ILD between January 2005 and December 2013. Among the 312,519 subjects, the control group was selected by excluding those who had ILD between 2005 and 2013.
Definition of interstitial lung disease
Cases of ILD were selected based on International Classification of Disease-10 (ICD-10) code J84 for other interstitial lung diseases, excluding drug-induced interstitial lung disorders, interstitial emphysema, and lung diseases caused by external agents. Connective tissue disease-associated ILD, hypersensitivity pneumonitis, and sarcoidosis were excluded.
Comorbidities diagnosed before the index date that could be associated with an increased risk of lung fibrosis were identified using ICD-10, including COPD, hepatitis C, gastroesophageal reflux disorder (GERD), and diabetes [10,11,12,13].
Baseline characteristics (including age, sex, and comorbidities) for cases and controls are summarized using descriptive statistics such as proportion. A chi-squared test was used to compare frequencies of risk factors between ILD and the control group. Cox proportional hazards regression models were used to evaluate the risk factors for ILD and analyze the associations between ILD and different variables and comorbidities. The final multivariate models included age, sex, smoking status (former or current smoker vs. never smoker), household income, and comorbidities such as hepatitis C, herpes, tuberculosis, pneumonia, GERD, COPD, diabetes, and hepatitis B. Risk factor models for ILD were selected according to sex and smoking as sensitivity analysis. Model selection method was forward stepwise procedure using likelihood ratio test with p-value < 0.05 as entry criterion, and p-value ≥0.10 as removal criterion.
A P value < 0.05 considered to be statistically significant. All statistical analyses were performed using SAS V.9.2 (SAS Institute, Cary, North Carolina, USA).
Incidence and baseline characteristics
The final sample included 312,519 subjects, of which 1972 developed ILD during the 9-year study period (Fig. 1). ILD incidence was 70.1 cases per 100,000 person-year.
All subjects were tracked by December 31, 2013, Follow up duration was median 65.6 months (interquartile range: 35.5, 89.0) in the ILD group and median 107.5 months (interquartile range; 100.2, 109.5) in the non-ILD group.
The ILD group had a higher percentage of men than the control. Compared with the control, subjects with ILD were older, and more likely to be smokers. The ILD group was more likely to have comorbidities such as respiratory diseases (tuberculosis and pneumonia), diabetes, chronic renal failure, malignancy, GERD, hepatitis C, and COPD than the control (Table 1).
Risk factors for developing ILD
Based on a multivariate Cox regression analysis of all variables (Table 2), smoking was significantly associated with the development of ILD (HR: 1.2; 95% CI: 1.1–1.4).
Co-morbidities such as hepatitis C (HR: 1.6; 95% CI: 1.1–2.3), history of tuberculosis (HR: 1.5; 95% CI: 1.1–1.9), history of pneumonia (HR: 1.6; 95% CI: 1.3–2.0), and COPD (HR: 1.8; 95% CI: 1.6–2.1) were significantly associated with the development of ILD. Men were almost twice as likely to associated with the development of ILD as women.
The risk of ILD development was 1.9 times in the 50s, 4.1 times in the 60s, and 6.9 times in the 70s, which increased sharply with age. In multivariate analysis, age was found to be most associated to ILD development showing dose response. All the risk factors included in the model showed relatively narrow confidence intervals.
There were no variables that changed the direction of the hazard ratio in the univariate and multivariate analyses except hepatitis B, which was statistically not significant. The hazard ratios of COPD and history of tuberculosis were reduced by almost half after multivariate analysis. History of pneumonia, hepatitis C and diabetes showed also reduced hazard ratios in multivariate analysis. Smoking was associated with the occurrence of ILD, but its magnitude was relatively small compared to other variables.
Multivariate analysis stratified by smoking showed same direction and similar magnitude of hazard ratios in all variables. However, hazard ratios stratified by sex were somewhat different in several variables. As a result of stratified analysis on the basis of gender, the hazard ratios of COPD according to man and woman were similar 1.8 and 1.9, respectively. In the case of GERD, multivariate analysis of male subjects showed that GERD was not a significant risk factor, but the hazard ratio for women was 1.3, which was consistent with the 1.3 observed in the univariate analysis (Table 3). GERD, Hepatitis C are significant risk factors in females, in contrast malignancy and diabetes are significant risk factors in males.
In the present study, the development of ILD was associated with older age, male sex, cigarette smoking, hepatitis C, history of tuberculosis, history of pneumonia, and COPD.
Smokers were at greater risk of developing ILD than non-smokers (HR: 1.2). Similarly, a previous study showed that smoking might contribute to the development of ILD by fibrosis (Odds Ratio: 1.6) . The findings of three different studies also support a strong association between ILA and exposure to tobacco smoke and smoking status [6,7,8]. In this study, the risk of developing ILD was 1.4 in smokers compared to non-smokers in the univariate analysis, but decreased 1.2-fold in multivariate analysis. The lowest risk among the variables associated with ILD occurrence is low, indicating that smoking is less involved in ILD development.
COPD was associated with the development of ILD in the present study. A 2012 review of the pathogenesis of IPF and COPD showed similarities between the basic pathogenic mechanisms involved in the development of either emphysema or fibrosis . Coexisting pulmonary fibrosis and emphysema is now a distinct entity , and studies show that the pathologic changes associated with these coexisting entities are mostly found in smokers [16,17,18,19]. In this study, the risk of developing ILD in COPD was 3.7 in univariate analysis but decreased 1.8-fold in multivariate analysis. This is presumably due to the control of cigarette smoking variables.
Hepatitis C was another risk factor associated with the development of ILD in the present study. There have been conflicting results regarding the prevalence of anti-hepatitis C virus (HCV) antibody in patients with IPF [10, 20, 21]. However, a 2008 study with HCV-infected patients and HBV-infected controls showed that ILD with fibrosis developed at a significantly greater rate in the HCV group than in the HBV group . In 2002, Idilman et al. reported an increased bronchoalveolar lavage neutrophil count in individuals with hepatitis C, suggesting an inflammatory reaction in the lungs leading to fibrotic changes . The findings of these studies suggest that systemic factors stimulating fibrosis, such as HCV infection, may affect the development of lung fibrosis.
Local factors such as history of pneumonia or tuberculosis may be associated with lung fibrosis. Nonresolving pneumonia may result in organizing pneumonia commonly in bacterial infections . Pneumonia due to mycoplasma and Legionnaires’ disease has been mostly implicated with development of pulmonary fibrosis [25, 26]. However, the occurrence of organization in cases of pneumonia is more common than expected. In 1952, Auerbach et al. studied the material from 307 necropsies and found organization in 38 cases . In 1989, Shachor et al. found that the incidence of tuberculosis in subjects with ILD was 4.5 times higher than that of the general population . Dheda et al. studied lung remodeling and fibrosis associated with lung injury from tuberculosis infection. Lung remodeling can result in extensive fibrosis and may be interstitial . Therefore, ILD may develop due to tuberculosis or other lung infections.
GERD is a well-known risk factor for IPF [11, 30,31,32]. However, our data did not show a significant association between GERD and ILD. In patients with esophagitis, the associated odds ratio of pulmonary fibrosis was 1.3–1.6 [30, 33], but there was no significant association between reflux esophagitis and pulmonary fibrosis in this study (Table 2). However, GERD were observed only in women as a significant risk factor for the development of ILD (Table 3). In the present study, ICD-10 code defined a wider range of ILD. Therefore, in this study, it is estimated that mild cases of pulmonary fibrosis are more involved than previous studies, and that there was no significant association between reflux esophagitis and ILD in men.
Diabetes is also prevalent with IPF [12, 33]. Although diabetes was an important risk factor for developing ILD in the univariate analysis, the significance reduced to borderline in the multivariate analysis (HR 1.1, P = 0.06).
Our study shows that men and older individuals have a higher risk of developing ILD. Two other studies conducted with subjects with interstitial abnormalities found that they were significantly older [7, 8]. Several studies have shown that the incidence of ILD with fibrosis is higher in men and increases with advanced age [34,35,36,37]. In our study, the HR for developing ILD in those aged ≥70 years was almost 7 times higher than for those aged 40–49 years. We suggest that ILD would be a result of the aging process.
We found reduction of risk estimates from the univariate to the multivariate analysis of some variables such as COPD and GERD. We did sensitivity analysis based on gender. Both males and females have similar risk estimates for ILD development in relation to COPD. However, GERD is significant risk factor in females but not in males.
ILD is known as non-homogeneous diseases. Although assessing risk factors for developing specific forms of ILD would be ideal, it might be difficult performing such kinds of study because surgical lung biopsy was rarely performed. In addition, assessing risk factors for the ILD development, the control group should be a general population. Therefore, we assess risk factors for developing ILD by population-based cohort database, though it does not provide information on specific forms of ILD.
The incidence of ILD was 70.1 cases per 100,000 per year, which is higher than previously reported . However, the present incidence was calculated for individuals > 40 years old.
Limitations of the present study
The most important limitation of the present study is that the diagnoses of ILD and other comorbidities were defined based on ICD codes, which may be inaccurate compared to the diagnoses obtained from a medical chart. Underreporting of asymptomatic ILD or misclassification was also possible.
The validity of the medical insurance claims data for ILD has not been determined in Korea. This database consists of random samples of national insurance claim data without identification numbers. Therefore, it was impossible to validate individual cases through a chart review.
We may have underestimated ILD incidence because of inaccurate ILD data. However, previous studies of ILD incidence using data from the Health Insurance Review and Assessment Service of Korea  reported similar results to those of previous studies [35, 36, 39]. The incidence rate of ILD was reportedly 48.5 per 100,000 person-years in Korea based on all claims data from 2008 to 2012 . In the present study, the incidence rate was 70.1 per 100,000 patients from 2005 to 2013, which appears reasonable because we excluded subjects aged < 40 years.
The present study did not take into consideration of relevant risk factor such as occupational and environmental exposure. However, annual income could be a weak proxy of the occupational exposure .
The present study may be affected by selection bias because the controls were identified based on medical claims. Thus, the controls were more likely to have comorbidities than controls selected from the general population. Although we excluded patients with a diagnosis of ILD between January 2002 and December 2004 to exclude pre-existing ILD diagnosed before the first year of our study (2005), patients with ILD could be miscounted as new cases if the patient did not require medical care between 2002 and 2004, which may have caused inaccuracies. In addition, other risk factors, such as pulmonary function and high-resolution computed tomography findings for ILD, could not be evaluated because of the NHIS-NSC 2002–2013 primarily included medical claims.
The population is aging, and the incidence of ILD increases with advancing age. Aged men with the history of previous lung infection or COPD could be a clinical marker of suspect ILD. The demonstration of the association of these clinical attributes with ILD may encourage the search for these factors in health examination and stimulate the translational research on ILD development.
Chronic obstructive pulmonary disease
Gastroesophageal reflux disorder
International classification of disease
Interstitial lung disease
Idiopathic pulmonary fibrosis
National health insurance
Coultas DB, Zumwalt RE, Black WC, Sobonya RE. The epidemiology of interstitial lung diseases. Am J Respir Crit Care Med. 1994;150(4):967–72.
Nalysnyk L, Cid-Ruzafa J, Rotella P, Esser D. Incidence and prevalence of idiopathic pulmonary fibrosis: review of the literature. Eur Respir Rev. 2012;21(126):355–61.
Putman RK, Hatabu H, Araki T, Gudmundsson G, Gao W, Nishino M, Okajima Y, Dupuis J, Latourelle JC, Cho MH, et al. Association between interstitial lung abnormalities and all-cause mortality. Jama. 2016;315(7):672–81.
Choi WI, Park SH, Dauti S, Park BJ, Lee CW. Interstitial lung disease and risk of mortality: 11-year nationwide population-based study. Int J Tuberc Lung Dis. 2018;22(1):100–5.
Choi WI, Park SH, Park BJ, Lee CW. Interstitial lung disease and lung Cancer development: a 5-year Nationwide population-based study. Cancer Res Treat. 2017;50(2):374–81.
Lederer DJ, Enright PL, Kawut SM, Hoffman EA, Hunninghake G, van Beek EJ, Austin JH, Jiang R, Lovasi GS, Barr RG. Cigarette smoking is associated with subclinical parenchymal lung disease: the multi-ethnic study of atherosclerosis (MESA)-lung study. Am J Respir Crit Care Med. 2009;180(5):407–14.
Washko GR, Hunninghake GM, Fernandez IE, Nishino M, Okajima Y, Yamashiro T, Ross JC, Estepar RS, Lynch DA, Brehm JM, et al. Lung volumes and emphysema in smokers with interstitial lung abnormalities. N Engl J Med. 2011;364(10):897–906.
Hunninghake GM, Hatabu H, Okajima Y, Gao W, Dupuis J, Latourelle JC, Nishino M, Araki T, Zazueta OE, Kurugol S, et al. MUC5B promoter polymorphism and interstitial lung abnormalities. N Engl J Med. 2013;368(23):2192–200.
Lee J, Lee JS, Park SH, Shin SA, Kim K. Cohort profile: the National Health Insurance Service-National Sample Cohort (NHIS-NSC), South Korea. Int J Epidemiol. 2017;46(2):e15.
Ueda T, Ohta K, Suzuki N, Yamaguchi M, Hirai K, Horiuchi T, Watanabe J, Miyamoto T, Ito K. Idiopathic pulmonary fibrosis and high prevalence of serum antibodies to hepatitis C virus. Am Rev Respir Dis. 1992;146(1):266–8.
Tobin RW, Pope CE 2nd, Pellegrini CA, Emond MJ, Sillery J, Raghu G. Increased prevalence of gastroesophageal reflux in patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 1998;158(6):1804–8.
Enomoto T, Usuki J, Azuma A, Nakagawa T, Kudoh S. Diabetes mellitus may increase risk for idiopathic pulmonary fibrosis. Chest. 2003;123(6):2007–11.
Chilosi M, Poletti V, Rossi A. The pathogenesis of COPD and IPF: distinct horns of the same devil? Respir Res. 2012;13:3.
Baumgartner KB, Samet JM, Stidley CA, Colby TV, Waldron JA. Cigarette smoking: a risk factor for idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 1997;155(1):242–8.
Travis WD, Costabel U, Hansell DM, King TE Jr, Lynch DA, Nicholson AG, Ryerson CJ, Ryu JH, Selman M, Wells AU, et al. An official American Thoracic Society/European Respiratory Society statement: update of the international multidisciplinary classification of the idiopathic interstitial pneumonias. Am J Respir Crit Care Med. 2013;188(6):733–48.
Auerbach O, Garfinkel L, Hammond EC. Relation of smoking and age to findings in lung parenchyma: a microscopic study. Chest. 1974;65(1):29–35.
Katzenstein AL, Mukhopadhyay S, Zanardi C, Dexter E. Clinically occult interstitial fibrosis in smokers: classification and significance of a surprisingly common finding in lobectomy specimens. Hum Pathol. 2010;41(3):316–25.
Cottin V, Nunes H, Brillet PY, Delaval P, Devouassoux G, Tillie-Leblond I, Israel-Biet D, Court-Fortune I, Valeyre D, Cordier JF, et al. Combined pulmonary fibrosis and emphysema: a distinct underrecognised entity. Eur Respir J. 2005;26(4):586–93.
Inomata M, Ikushima S, Awano N, Kondoh K, Satake K, Masuo M, Kusunoki Y, Moriya A, Kamiya H, Ando T, et al. An autopsy study of combined pulmonary fibrosis and emphysema: correlations among clinical, radiological, and pathological features. BMC Pulm Med. 2014;14:104.
Irving WL, Day S, Johnston ID. Idiopathic pulmonary fibrosis and hepatitis C virus infection. Am Rev Respir Dis. 1993;148(6 Pt 1):1683–4.
Meliconi R, Andreone P, Fasano L, Galli S, Pacilli A, Miniero R, Fabbri M, Solforosi L, Bernardi M. Incidence of hepatitis C virus infection in Italian patients with idiopathic pulmonary fibrosis. Thorax. 1996;51(3):315–7.
Arase Y, Suzuki F, Suzuki Y, Akuta N, Kobayashi M, Kawamura Y, Yatsuji H, Sezaki H, Hosaka T, Hirakawa M, et al. Hepatitis C virus enhances incidence of idiopathic pulmonary fibrosis. World J Gastroenterol. 2008;14(38):5880–6.
Idilman R, Cetinkaya H, Savas I, Aslan N, Sak SD, Bastemir M, Sarioglu M, Soykan I, Bozdayi M, Colantoni A, et al. Bronchoalveolar lavage fluid analysis in individuals with chronic hepatitis C. J Med Virol. 2002;66(1):34–9.
Cordier JF, Cottin V, Lazor R, Thivolet-Bejui F. Many faces of bronchiolitis and organizing pneumonia. Semin Respir Crit Care Med. 2016;37(3):421–40.
Kaufman JM, Cuvelier CA, Van der Straeten M. Mycoplasma pneumonia with fulminant evolution into diffuse interstitial fibrosis. Thorax. 1980;35(2):140–4.
Chastre J, Raghu G, Soler P, Brun P, Basset F, Gibert C. Pulmonary fibrosis following pneumonia due to acute Legionnaires' disease. Clinical, ultrastructural, and immunofluorescent study. Chest. 1987;91(1):57–62.
Auerbach SH, Mims OM, Goodpasture EW. Pulmonary fibrosis secondary to pneumonia. Am J Pathol. 1952;28(1):69–87.
Shachor Y, Schindler D, Siegal A, Lieberman D, Mikulski Y, Bruderman I. Increased incidence of pulmonary tuberculosis in chronic interstitial lung disease. Thorax. 1989;44(2):151–3.
Dheda K, Booth H, Huggett JF, Johnson MA, Zumla A, Rook GA. Lung remodeling in pulmonary tuberculosis. J Infect Dis. 2005;192(7):1201–9.
el-Serag HB, Sonnenberg A. Comorbid occurrence of laryngeal or pulmonary disease with esophagitis in United States military veterans. Gastroenterology. 1997;113(3):755–60.
Lee JS, Collard HR, Anstrom KJ, Martinez FJ, Noth I, Roberts RS, Yow E, Raghu G, Investigators IP. Anti-acid treatment and disease progression in idiopathic pulmonary fibrosis: an analysis of data from three randomised controlled trials. Lancet Respir Med. 2013;1(5):369–76.
Lee JS, Ryu JH, Elicker BM, Lydell CP, Jones KD, Wolters PJ, King TE Jr, Collard HR. Gastroesophageal reflux therapy is associated with longer survival in patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2011;184(12):1390–4.
Gribbin J, Hubbard R, Smith C. Role of diabetes mellitus and gastro-oesophageal reflux in the aetiology of idiopathic pulmonary fibrosis. Respir Med. 2009;103(6):927–31.
Navaratnam V, Fleming KM, West J, Smith CJ, Jenkins RG, Fogarty A, Hubbard RB. The rising incidence of idiopathic pulmonary fibrosis in the UK. Thorax. 2011;66(6):462–7.
Esposito DB, Lanes S, Donneyong M, Holick CN, Lasky JA, Lederer D, Nathan SD, O'Quinn S, Parker J, Tran TN. Idiopathic pulmonary fibrosis in United States automated claims. Incidence, prevalence, and algorithm validation. Am J Respir Crit Care Med. 2015;192(10):1200–7.
Raghu G, Chen SY, Yeh WS, Maroni B, Li Q, Lee YC, Collard HR. Idiopathic pulmonary fibrosis in US Medicare beneficiaries aged 65 years and older: incidence, prevalence, and survival, 2001-11. Lancet Respir Med. 2014;2(7):566–72.
Raghu G, Chen SY, Hou Q, Yeh WS, Collard HR. Incidence and prevalence of idiopathic pulmonary fibrosis in US adults 18-64 years old. Eur Respir J. 2016;48(1):179–86.
Gjonbrataj J, Choi WI, Bahn YE, Rho BH, Lee JJ, Lee CW. Incidence of idiopathic pulmonary fibrosis in Korea based on the 2011 ATS/ERS/JRS/ALAT statement. Int J Tuberc Lung Dis. 2015;19(6):742–6.
Hutchinson J, Fogarty A, Hubbard R, McKeever T. Global incidence and mortality of idiopathic pulmonary fibrosis: a systematic review. Eur Respir J. 2015;46(3):795–806.
Evans GW, Kantrowitz E. Socioeconomic status and health: the potential role of environmental risk exposure. Annu Rev Public Health. 2002;23:303–31.
We would like to thanks to Ms. Byeong Ju Park for the preparing data sets.
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MSIP) (No. 2014R1A5A2010008).
Availability of data and materials
The datasets supporting the conclusion of this article is not available. The dataset comes from whole national population. National Health Insurance (NHI) service strictly controlled the database. NHI completely restricted the data sets used in this article open to the public.
Ethics approval and consent to participate
The current study was approved by the institutional review board at Dongsan Hospital, Keimyung University School of Medicine. The need for written informed consent was waived. This study was conducted in compliance with the Declaration of Helsinki.
The authors state that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.