2.1 Participants
The 2018 Korea National Survey on Suicide (KNSS) [37], conducted by the Ministry of Health and Welfare of Korea, served as the basis for this study. This survey was conducted based on a complex sample design to ensure representativeness. Structured face-to-face interviews were conducted with 1,500 adults aged 19 to 75 between November 21 and December 17, 2018. The interviewers were trained by experts and supported during the survey. Ten households were randomly selected for each district, and an interviewer visited each household and conducted the survey. Following the interview, 30% of the participants underwent phone follow-up for quality control.
The Seoul National University Hospital Institutional Review Board approved and monitored this study (IRB No. 1405-019-577). Interviewers obtained informed consent from participants before commencing the interviews. This study adhered to the Strengthening the Reporting of Observational Studies in Epidemiology reporting guidelines [38] and was conducted in compliance with the Declaration of Helsinki.
2.2 Measurements
Demographic characteristics of participants, including age, gender, education, employment, and religion, were collected. Participants were also asked about any previous experience with suicidal ideation.
The Questionnaire on Attitudes Toward Suicide (ATTS) was employed to assess the normative attitudes of participants toward suicide. Developed and validated by Renberg and Jacobsson [22], the ATTS is suitable for large-scale studies and widely used in various countries [39]. The Korean version, created through the 2013 KNSS [40], was utilized in the 2018 KNSS. Each item of the ATTS was scored on a Likert scale of 1-5 points, where 1 indicated “disagreement” and 5 indicated “agreement.”
2.2.1 Factor structure of ATTS
Kim et al. previously reported a factor structure of ATTS based on the 2013 and 2018 KNSS [41]. In this study, we reevaluated the previous factor structure, as multiple methods for determining the number of factors indicated that a three-factor structure best fits the 2018 KNSS data alone (Supplementary Fig. 1). Unreliable and semantically redundant factors were removed or combined to create a more parsimonious structure while maintaining the basic framework of the original factor structure. The revised factor structure encompassed 3 factors and 15 items (Table 1). Through confirmatory factor analysis (Supplementary Fig. 2), we assessed the goodness of fit and acquired factor score. The revised factor structure had a good absolute fit and parsimony correction, and the comparative fit was within an acceptable or satisfactory range [42]. Acquisition of factor score used Barlett’s method.
Table 1. Revised factor structure for the Questionnaire on Attitudes Toward Suicide. Items were scored on a Likert scale ranging from 1 to 5, with 1 indicating disagreement and 5 indicating agreement. Reliability coefficients were calculated using Cronbach's alpha and omega total. English expressions for each item are shortened. SE standard error, α Cronbach’s alpha, ω Omega total.
Factor
|
Number
|
Items
|
Mean
|
SE
|
Reliability
|
Permissiveness (PA)
|
5
|
Suicide acceptable for incurable disease
|
2.6
|
0.041
|
α = 0.70 ω = 0.78
|
16
|
Situations where suicide is the only solution
|
2.6
|
0.044
|
18
|
Suicide as relief
|
2.7
|
0.044
|
20
|
Consider suicide if incurable disease—self
|
2.9
|
0.046
|
34
|
Right to commit suicide
|
2.5
|
0.042
|
36
|
Get help for suicide if incurable disease—self
|
3.1
|
0.040
|
Unjustified behavior (UB)
|
2
|
Suicide never justifiable
|
3.9
|
0.037
|
α = 0.55 ω = 0.61
|
3
|
Suicide is among the worst for relatives
|
4.3
|
0.032
|
19
|
Youth suicides are particularly puzzling
|
3.6
|
0.056
|
27
|
Express suicide wish without meaning it—self
|
3.5
|
0.039
|
Readiness to help/ Preventability (RP)
|
1
|
Always able to help
|
3.3
|
0.052
|
α = 0.50 ω = 0.57
|
9
|
Duty to restrain suicidal act
|
3.9
|
0.036
|
30
|
Ready to help a suicidal person—self
|
3.2
|
0.038
|
33
|
Suicide talkers are not always completers
|
3.7
|
0.036
|
37
|
Suicide preventable
|
3.9
|
0.032
|
The factors names, “Permissiveness (PA),” “Unjustified behavior (UB)” and “Readiness to help/Preventability (RP),” were derived from the original factor structure. PA semantically represents permissive and acceptable attitudes toward suicide, considering it a right or option in certain circumstances, such as incurable diseases. UB represents opposition to justifying suicide and regards it as morally reprehensible behavior that should be restrained. Intriguingly, this factor also included items that perceive suicide as puzzling or deceptive. RP represents proactive attitudes toward suicide prevention, characterized by a readiness to help and taking an active role in preventing suicide. Compared to PA, participants in the 2018 KNSS generally demonstrated greater agreement with items in UB and RP, and their standard errors were narrower for most items. The reliability of UB and RP, evaluated using Cronbach's alpha, is low; however, it surpasses the requisite value of 0.5 for preliminary or exploratory research, as proposed by Nunnally [43].
2.2.2 Suicide mortality
Microdata on suicide mortality was obtained through the MicroData Integrated System (MDIS) provided by Statistics Korea [44]. The MDIS provides individual information, including the deceased’s gender, age, address, and cause of death, as registered in the Causes of Death Statistics. Researchers remotely accessed and processed the data stored on the MDIS server, and Statistics Korea examined and approved the anonymization of the processed data when exported. This study included individuals who died between January 1 and December 31, 2018, aged 19 to 75 years, and whose causes of death corresponded to the X60-X84 (intentional self-harm) codes of the Korean Standard Classification of Diseases 7th edition [45].
2.3 Statistical analysis
The country was divided into 30 regional units for ecological analysis. Survey strata comprising 80 or more participants were reorganized into clusters of 40 to 70 participants based on geographical proximity, forming regional units (Supplementary Table 1). Basic characteristics, ATTS, and microdata on suicide mortality were aggregated by these regional units, utilizing survey weights based on a complex sample design. Comparisons were conducted between regions with high and low suicide rates to demonstrate differences in regional characteristics according to suicide rates. If the Shapiro-Wilk test rejected the normal distribution, Wilcoxon's test was employed for comparison; otherwise, Welch's t-test was utilized.
A negative binomial regression model was employed to investigate the effect of ATTS and potential covariates on suicide mortality. This model can be generalized to overdispersed count data in place of Poisson models and is widely used [46]. A univariate model (Modeluniv) and three multivariate models (Model 1 to Model 3) were constructed. In Model 1, adjustments were made for the proportions of females and adults over 60 years. Using the proportion of older adults as a covariate rather than the mean age allowed for consideration of South Korea's unique characteristics, where elderly suicide is notably prevalent, and its etiology appears to be differentiated [47,48]. In Models 2 and 3, one additional covariate from Model 1 was further adjusted by adding lifetime experience of suicidal ideation (Model 2) and college education (Model 3). The dependent variable was the incidence of suicide in each region, with the regional population set as an offset variable. The result of each model was presented as an incidence rate ratio (IRR), and model performance was assessed by Akaike Information Criteria (AIC) and Bayesian Information Criteria (BIC). Additionally, due to differences in the mechanisms resulting in suicide based on gender and age [49-52], stratified analyses were performed. Stratification was based on gender (male, female) and age group (-39, 40–59, 60- years).
R Statistical Software v4.2.2 [53] was used for all statistical analyses. Furthermore, package lavaan v0.6-13 [54] was utilized for confirmatory factor analysis, package survey v4.1-1 [55] for analysis based on the complex sample design, and package MASS v7.3-58.2 [56] for the negative binomial regression model. All statistical tests were two-tailed, and a P-value of 0.05 was considered the threshold for significance. Asterisks were used to indicate the significance level for enhanced readability: * for P < 0.05, ** for P < 0.01, and *** for P < 0.001.