Predicting South Korea adolescents vulnerable to depressive disorder using Bayesian nomogram: A community-based cross-sectional study.

BACKGROUND: Although South Korea has developed and carried out evidence-based interventions and prevention programs to prevent depressive disorder in adolescents, the number of adolescents with depressive disorder has increased every year for the past 10 years. AIM: To develop a nomogram based on a naïve Bayesian algorithm by using epidemiological data on adolescents in South Korea and present baseline data for screening depressive disorder in adolescents.
METHODS: Epidemiological data from 2438 subjects who completed a brief symptom inventory questionnaire were used to develop a model based on a Bayesian nomogram for predicting depressive disorder in adolescents.
RESULTS: Physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents were influential factors on depressive disorder in adolescents. Among them, physical symptoms were the most influential.
CONCLUSION: Active intervention by periodically checking the emotional state of adolescents and offering individual counseling and in-depth psychological examinations when necessary are required to mitigate depressive disorder in adolescents. ©The Author(s) 2022. Published by Baishideng Publishing Group Inc. All rights reserved.

Core Tip: The early detection and prevention of depressive disorder in adolescents is important because it not only adversely affects interpersonal relationships and academic achievement but also increases the probability of other related mental illnesses such as panic disorder. We developed a nomogram for screening depressive disorder using epidemiological data on 2438 adolescents. Physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents were influential factors on depressive disorder in adolescents.

INTRODUCTION

Depressive disorder causes an emotional response that can make a person feel a sense of failure, a sense of loss, and/or a sense of worthlessness as a result of a negative perception of him/herself[1]. It is defined as a persistent feeling of sadness or hopelessness to the extent of not being able to maintain daily activities for 2 wk in the past year[1]. It has been reported that South Koreans experience depressive disorder most frequently during adolescence compared to other stages of life[2]. A national survey of South Korean adolescents reported that one in four males and one in three females were diagnosed with depressive disorder[3]. In particular, it has been reported that depression during adolescence increases rapidly after middle school[3,4], suggesting that the mental health of adolescents is at risk during this period.

Adolescence involves the most physical, mental, and social changes that occur in a human lifetime[5]. Adolescents experiencing depressive disorder are highly likely to self-torture and/or express delinquent and aggressive behavior[5]. Moreover, if depressive disorder is not identified and managed early, it may progress to become a chronic illness with depression likely to recur during a person’s lifetime[6]. The early detection and prevention of depressive disorder in adolescence is an important social issue because it not only adversely affects interpersonal relationships and academic achievement but also increases the probability of developing other mental illnesses such as panic disorder[7].

Although South Korea has developed and carried out evidence-based intervention and prevention programs to mitigate depressive disorder in adolescents[2], the number of adolescents with depressive disorder has increased every year for the past 10 years[2]. Consequently, it is necessary to identify the influential factors causing depression and develop a predictive model with high accuracy that can identify groups highly vulnerable to depressive disorder as soon as possible.

Recently, the naïve Bayesian nomogram has been used as a method for predicting groups at high risk of developing diseases[8,9]. One of the advantages of this method is that it presents the risk probability according to multiple risk factors of a disease visually so that clinicians can easily understand the results[10]. In this study, a nomogram based on a naïve Bayesian algorithm using epidemiological data on adolescents in South Korea was developed and baseline data for screening depressive disorder in adolescents is presented.

MATERIALS AND METHODS

Data source

This is a secondary data analysis study using raw data from the 2019 Korean Children Youth Panel Study (KCYPS) survey from March to June 2019 provided by the National Youth Policy Institute. The study was approved by the Research Ethics Review Board of the National Youth Policy Institute (No. KCYPS-2018).

The survey method for KCYPS is presented in Cho et al[11] (2018). Briefly, the KCYPS sampled 7th-grade students attending 162 middle schools across South Korea using a stratified multi-stage cluster sampling method. Schools were selected according to the probability proportional to the size sampling method for 27 clusters across 16 metropolitan cities, small and medium-sized cities, and rural areas. After checking the information on the number of 7th-grade classes and the number of students in each class at each school, samples were extracted by randomly selecting classes. The KCYPS collected data using a tablet-assisted personal interview method to compensate for the quality deterioration caused by existing questionnaire input errors or logical errors and to increase the accuracy and efficiency of the survey. In the present study, we analyzed 2438 subjects after excluding 152 cases with missing values in the depressive disorder screening part among 2590 people who completed the KCYPS questionnaire in 2019.

Measurements

Depression, the outcome variable, was defined by using ten items for measuring depression in the brief symptom inventory (BSI) (1984)[12], which was adapted for the South Korean population by standardizing the Symptom Checklist-90-Revision[13]. The BSI is a self-reporting test with each item being measured on a 4-point scale. Moreover, the total score ranges from 10 to 40 points. A higher score indicates more severe depression. Referring to Byeon et al[14] (2015), the threshold for depression in this study was 24 points, corresponding to 1 standard deviation (less than the 16th percentile). AORN J reported that Cronbach’s α (a measurement of reliability) for the BSI was 0.904 (0.882 in the present study)[15].

Explanatory variables included gender, environmental factors (number of siblings and mean conversation time with parents during weekdays), and personal factors (satisfaction with academic achievement, satisfaction with school life, mean sleeping hours during weekdays, social withdrawal, aggression, attention, and physical symptoms). The definitions of the explanatory variables are provided in Table 1.

Table 1

Measurements of explanatory variables

Classification	Variable	Characteristics
Sociodemographic factors	Gender	Male or female
	Number of siblings (including the subject)	1 person, 2 people, 3 people, or 4 people or more
Environmental factors	Mean conversation time with parents per day	< 30 min, ≥ 30 min and < 1 h, ≥ 1 h and < 2 h, ≥ 2 h and < 3 h, or ≥ 3 h
Personal factors	Satisfaction with academic achievement	Dissatisfied, not dissatisfied or satisfied, or satisfied
	Satisfaction with school life	Dissatisfied, not dissatisfied or satisfied, or satisfied
	Mean sleeping hours per day	< 5 h, 6 h, 7 h, 8 h, 9 h, or ≥ 10 h
	Social withdrawal	Continuous variable
	Aggression	Continuous variable
	Attention	Continuous variable
	Physical symptoms	Continuous variable

Social withdrawal was measured by using five items from the Behavior Problem Scale for Children and Adolescence (BPSCA) developed by Kim and Kim[16] (1998) after excluding items overlapping with other sub-domains. Each item was measured on a 4-point scale with the total score ranging from 5 to 20 points. A higher score indicates more severe social withdrawal. AORN J reported that Cronbach’s α for the tool was 0.850 (0.894 in the present study)[15].

Aggression was measured by using the Emotional or Behavioral Problems Scale (EPS) developed by Cho and Lim[17] (2003). Six items were used and each item was measured on a 4-point scale with the total score ranging from 6 to 24 points, and a higher score indicates a more aggressive condition. Cho and Lim[17] (2003) reported that the Cronbach’s α of the tool was 0.760 (0.809 in the present study).

Attention problems were measured by using 7 items in the EPS[17]. Each item was measured on a 4-point scale with the total score ranging from 7 to 28 points. A higher score indicates more severe attention problems. Kim and Song[18] (2014) reported Cronbach's α for the tool was 0.791 (0.813 in the present study).

Physical symptoms are when a person perceives that he or she is frequently ill or tired such as chest tightness or stomach discomfort without a pathological cause. These were measured by using eight items in the EPS[17]. Each item was measured on a 4-point scale with the total score ranging from 8 to 32 points. A higher score indicates more severe physical symptoms. Choi et al[19] (2017) reported Cronbach’s α for the tool was 0.87 (0.858 in the present study).

Developing the naïve Bayes nomogram for predicting adolescents vulnerable to depressive disorder

A nomogram is used to visually present complex functions or calculations[20,21]. In particular, it is used as a method to visually present the diagnosis, recurrence, and survival prediction of a disease[20,21]. It is expressed graphically (Figure 1) in which a line is assigned to each attribute used as an input item and the possible value of the attribute is displayed on the line[22]. The score corresponding to the position of the attribute value becomes the individual score of the point displayed at the top.

Figure 1

An example of a nomogram[Citation: Byeon H. Developing a nomogram for predicting the depression of senior citizens living alone while focusing on perceived social support. World J Psychiatry 2021; 11: 1314-1327. Copyright ©The Authors 2021. Published by Baishideng Publishing Group Inc.

We used a naïve Bayes classifier as the algorithm to develop the nomogram. A naïve Bayes classifier model determines the probability for a specific class by applying the Bayesian theorem under the assumption that the attributes of data and events are independent of each other. When the attributes are assumed to be independent, the posterior probability indicating the probability that an object (belongs to class C can be expressed as follows:

Where c is the target class of the nomogram. However, if it is a class other than c, and P(|X) represents the probability that object X does not belong to class c[9], then the odds ratio (OR) for these two probabilities can be calculated as:

We calculated the general accuracy, precision, recall, F-1 score, the area under the curve (AUC), and calibration plot using leave-one-out cross-validation (LOOCV) of the developed Bayesian algorithm-based nomogram to validate its predictive performance.

Precision is defined as the proportion of classifications that are true actually being true:

Recall is defined as the ratio of the number of model predictions that are true over the number that are actually true:

Accuracy is an evaluation index that can most intuitively indicate the performance of a model:

However, since using accuracy alone to overcome bias due to data imbalance is limited, it is necessary to present the F-1 score as an additional predictive performance indicator to overcome bias.

The F-1-score is the harmonic mean of Precision and Recall; i.e.

AUC is an indicator used to evaluate the performance of a binary classifier. The maximum value is 1, and a value close to 1 means that the performance of the model is good (i.e., the recall is larger than the fall-out). All analyses were performed by using Python version 3.10.0 (https://www.python. org/downloads; accessed on November 28, 2021).

RESULTS

General characteristics of the subjects

Of the 2438 subjects, the majority were male (54.1%) in a family with two siblings including the subject (61%); 32.0% had 30 minutes or more but less than 1 h of mean conversation time with their parents during weekdays; 37.2% were neither satisfied nor dissatisfied with last semester's school performance, 44.2% were satisfied with last semester's school life, and 40.4% slept for an average of 8 h per day during weekdays (Table 2). Their aggression, attention problems, physical symptoms, and social withdrawal are presented in Figure 2.

Table 2

General characteristics of subjects (mean ± SD)

Characteristic	n	%
Depressive disorder
No	1999	82.0
Yes	439	18.0
Gender
Male	1318	54.1
Female	1120	45.9
Number of siblings (including the subject)
1 person	358	14.7
2 people	1487	61.0
3 people	515	21.1
4 people	78	3.2
Mean sleeping hours per day
< 5 h	63	2.6
6 h	236	9.7
7 h	600	24.6
8 h	986	40.4
9 h	454	18.6
≥ 10 h	99	4.1
Mean conversation time with parents per day
< 30 min	456	18.7
≥ 30 min and < 1 h	781	32.0
≥ 1 h and < 2 h	644	26.4
≥ 2 h and < 3 h	351	14.4
≥ 3 h	206	8.4
Satisfaction with academic achievement
Dissatisfied	577	23.9
Not dissatisfied or satisfied	906	37.6
Satisfied	928	38.5
Satisfaction with school life
Dissatisfied	144	5.9
Not dissatisfied or satisfied	616	25.4
Satisfied	1666	68.7
Attention	15.2 ± 3.9
Aggression	11.4 ± 3.5
Social withdrawal	10.6 ± 3.5
Physical symptoms	14.9 ± 4.8

Figure 2

Test results. A: Subject's attention test; B: Subject's aggression test; C: Subject's social withdrawal test; D: Subject's physical symptoms test.

General sample characteristics and prevalence of depressive disorder

The results of Chi-square tests showed significant differences (P < 0.05) between the groups with and without depressive disorder in gender, mean sleeping hours per day, mean conversation time with parents per day, satisfaction with academic achievement, satisfaction with school life, attention, aggression, social withdrawal, and physical symptoms (Table 3).

Table 3

Characteristics by prevalence of depressive disorder, n (%) (mean ± SD)

Characteristic	Depressive disorder		P value
	No (n = 1999)	Yes (n = 439)
Gender			< 0.001
Male	1119 (84.9)	199 (15.1)
Female	880 (78.6)	240 (21.4)
Number of siblings (including the subject)			0.671
1 person	301 (84.1)	57 (15.9)
2 people	1217 (81.8)	270 (18.2)
3 people	419 (81.4)	96 (18.6)
4 people	62 (79.5)	16 (20.5)
Mean sleeping hours per day			<0.001
< 5 h	44 (69.8)	19 (30.2)
6 h	191 (80.9)	45 (19.1)
7 h	512 (85.3)	88(14.7)
8 h	841 (85.3)	145 (14.7)
9 h	350 (77.1)	104 (22.9)
≥ 10 h	61 (61.6)	38 (38.4)
Mean conversation time with parents per day			< 0.001
< 30 min	240 (74.6)	116 (25.4)
≥ 30 min and < 1 h	645 (82.6)	136 (17.4)
≥ 1 h and < 2 h	539 (83.7)	105 (16.3)
≥ 2 h and < 3 h	293 (83.5)	58 (16.5)
≥ 3 h	182 (88.3)	24 (11.7)
Satisfaction with academic achievement			< 0.001
Dissatisfied	434 (75.2)	143 (24.8)
Not dissatisfied or satisfied	735 (81.1)	171 (18.9)
Satisfied	812 (87.5)	116 (12.5)
Satisfaction with school life			< 0.001
Dissatisfied	63 (43.8)	81 (56.3)
Not dissatisfied or satisfied	470 (76.3)	146 (23.7)
Satisfied	1457 (87.5)	209 (12.5)
Attention	14.6 ± 3.8	17.7 ± 3.1	< 0.001
Aggression	10.6 ± 3.1	15.1 ± 3.1	< 0.001
Social withdrawal	10.0 ± 3.3	13.6 ± 2.9	< 0.001
Physical symptoms	13.7 ± 4.2	20.2 ± 3.7	< 0.001

Correlations between variables

Correlation analysis results between the major variables used in this study are presented in Figure 3. Depressive disorder was significantly and positively correlated with attention, aggression, social withdrawal, and physical symptoms (P < 0.05).

Figure 3

Correlation between variables: Scatter matrix. depN: Depressive disorder; brot_n: Number of siblings; sleep_N: Mean sleeping hours per day; Pa_talk_N: Mean conversation time with parents per day; record_N: Satisfaction with academic achievement; sch_life_N: Satisfaction with school life; atten: Attention; attxk: Aggression; etop: Social withdrawal; phy: Physical symptoms.

Predicting the group of adolescents vulnerable to depressive disorder by using the Bayesian nomogram

Figure 4 shows the Bayesian nomogram for predicting the adolescent group vulnerable to depressive disorder. We developed a nomogram comprising seven variables with high importance: Physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents were the major influential factors associated with depression in adolescents. Physical symptoms comprised the most influential factor for predicting depression in this high-risk group. We predicted the depression risk of South Korean adolescents by using the developed nomogram (Figure 4). The high-risk group comprised those who received 15.5 points for physical symptoms (EPS test), 11.5 points for aggression (EPS test), 10.5 points for social withdrawal (BPSCA test), and 17.5 points for attention (EPS test) and were dissatisfied with their school life, slept 10 h or more per day on average, and talked with parents less than 30 min (84% of developing depression).

Figure 4

A model for predicting adolescent groups vulnerable to depressive disorder by using Bayesian nomograms. phy: Physical symptoms; attxk: Aggression; etop: Social withdrawal; atten: Attention; sch_life_N: Satisfaction with school life (1 = dissatisfied, 2 = not dissatisfied or satisfied, 3 = satisfied); Pa_talk_N: Mean conversation time with parents per day (1: < 30 min, 2: ≥ 30 min and < 1 h, 3: ≥ 1 h and < 2 h, 4: ≥ 2 h and < 3 h; 5: ≥ 3 h).

The predictive performance of the developed nomogram for predicting the adolescent group highly vulnerable to depressive disorder was validated by using AUC, F-1 score, accuracy, and a calibration plot. The results of the LOOCV evaluation show that the model had an AUC of 0.90 (Figure 5), F-1 score of 0.86, general accuracy of 0.85, precision of 0.88, and recall of 0.86. Adolescents with and without depressive disorder were compared by using a calibration plot (Figure 6) and Chi-square tests based on the predicted and observed probability, between which there was no significant difference (P = 0.683).

Figure 5

Receiver operating characteristic analysis result of the developed model.

Figure 6

Calibration plot comparing predicted to actual probability of depressive disorder.

DISCUSSION

This study was conducted to present baseline data for preventing depressive disorder in adolescents by identifying multiple influential risk factors. The results reveal that physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents were significant predictors. Among them, physical symptoms had the greatest influence on the depressive disorder of adolescents. The outcomes of numerous previous studies on variables associated with depression in adolescents identify peer relationships, the home environment, and the school environment as significant risk factors[23-26], which supports the findings of the present study.

From the perspective of the socioecological model, family, peer group, and school are three major domains directly affecting the mental health of adolescents[27,28]. Since risk and protective factors are generated in these three domains[27,28], the viewpoint of the socioecological model is useful for explaining the depressive disorder of adolescents as the outcome of multiple risk factors. Nevertheless, there were limitations in previous studies about explaining the relationship between multiple risk factors and depressive disorder[2,25]. First, researchers mainly used regression models as the method of exploring the risk factors associated with depressive disorder in adolescents. Although calculating ORs by using regression analysis is useful for identifying individual risk factors, its ability to identify complex multiple risk factors is limited. Second, only sociodemographic and environmental factors as risk factors associated with depressive disorder were identified in previous studies[23,25]. Indeed, comprehensive analysis of risk factors associated with depressive disorder in adolescents by using individual factors such as difficulty with attention, social withdrawal, as well as environmental factors, has still not been sufficiently conducted. Third, as normality and independence between the variables are assumed in regression analysis, it is difficult to draw accurate conclusions because the data on many diseases such as depressive disorder are unbalanced, thereby violating the normality assumption.

In summary, there are limitations when utilizing depression prediction models for adolescents based on regression analysis in the primary medical care environment because it is difficult to identify the complex relationships between multiple risk factors solely relying on ORs based on regression models. Thus, we identified an adolescent group highly vulnerable to depressive disorders by using multiple risk factors based on a Bayesian nomogram to overcome these limitations. Our results predicted that adolescents who received 15.5 points in physical symptoms, 11.5 points in aggression, 10.5 points in social withdrawal, and 17.5 points in attention and who were dissatisfied with their current school life, slept for 10 h or more per day on average, and talked with their parents less than 30 min have a depression risk of 84%. Therefore, communities and schools must continually monitor the high-risk group for the early identification and prevention of depressive disorder in adolescents with these multiple risk factors.

Another important finding of the present study is that physical symptoms in adolescents comprised the most influential risk factor in predicting depressive disorder. Ryu and Hong[29] (2019) also explored factors affecting depressive disorder in 1881 middle school students and confirmed that physical symptoms of adolescents comprised the main risk factor influencing depressive disorder. Choi et al[19] (2017) also revealed that physical symptoms and depressive disorder had a positive correlation in fourth graders.

Physical symptoms in adolescents, which are related to mental activities and the psychological state, are generally overlooked as an early symptom of depressive disorder because they cannot be found by internal or neurological examination, or even when a physical abnormality is found, the symptoms are insufficient for disease diagnosis. However, although depressive disorder of adolescents is similar to the adult psychopathology, unlike in adults, clinical characteristics are often accompanied by physical symptoms (e.g., fatigue, insomnia, muscle pain, and headache) and aggression[30]. In particular, Jung et al[31] (2004) reported that depressed people excessively focus on physical symptoms or amplify their bodily sensations. Therefore, frequent complaints by adolescents of physical symptoms without a known medical cause are likely to be early signs of depressive disorder, even when the physical symptoms seem superficially less severe. Consequently, the community and school must pay attention to them and actively intervene by periodically checking the emotional state of adolescents, as well as providing individual counseling and in-depth psychological testing.

The strength of the present study is that it identified the group at high risk of developing depressive disorder based on multiple risk factors by using epidemiological data on South Korean adolescents and provided evidence for the early screening and management of depression. However, it does have some limitations, with the first being that there could be more potential variables for depressive disorder in addition to the explanatory variables used in this study because we analyzed secondary data. Second, the results cannot be generalized for all high school students because we identified a high-risk group for depressive disorder in seventh graders only. Third, the variables used (including depressive disorder) were measured based on a self-report questionnaire. Thus, future studies are needed to identify groups at high risk of depressive disorder by integrating qualitative research methods such as Delphi analysis and in-depth interviews in addition to self-report questionnaires. Fourth, since the results were based on a cross-sectional approach, it is difficult to determine causal relationships. Hence, additional prospective cohort studies should be conducted to prove causality between the depressive disorder high-risk group and depressive disorder found in the present study.

CONCLUSION

We showed that physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents are influential factors associated with depressive disorder in adolescents. Among them, physical symptoms comprise the most influential factor in the prediction of depressive disorder. Therefore, periodically checking on the emotional state of adolescents is required, along with providing individual counseling and conducting in-depth psychological examinations when necessary. Moreover, longitudinal studies based on clinical depressive disorder data targeting depressive disorder in the high-risk group confirmed in this study should be conducted.

ARTICLE HIGHLIGHTS

Research background

Although South Korea has developed and carried out evidence-based intervention and prevention programs to mitigate depressive disorder in adolescents, the number of adolescents with depressive disorder has increased every year for the past 10 years. Consequently, it is necessary to identify the influential factors causing depression and develop a predictive model with high accuracy that can identify groups highly vulnerable to depressive disorder as soon as possible.

Research motivation

Recently, the naïve Bayesian nomogram has been used as a method of predicting groups at high risk of developing diseases. One of the advantages of this method is that it presents the risk probability according to multiple risk factors of a disease visually so that clinicians can easily understand the results.

Research objectives

In this study, a nomogram based on a naïve Bayesian algorithm using epidemiological data on adolescents in South Korea was developed and baseline data for screening depressive disorder in adolescents was presented.

Research methods

We used a naïve Bayes classifier as the algorithm to develop the nomogram. Also, we calculated the general accuracy, precision, recall, F-1 score, the area under the curve, and calibration plot using leave-one-out cross-validation of the developed Bayesian algorithm-based nomogram to validate its predictive performance.

Research results

We showed that physical symptoms, aggression, social withdrawal, attention, satisfaction with school life, mean sleeping hours, and conversation time with parents were influential factors associated with depressive disorder in adolescents.

Research conclusions

Periodically checking on the emotional state of adolescents is required, along with providing individual counseling and conducting in-depth psychological examinations when necessary.

Research perspectives

Longitudinal studies based on clinical depressive disorder data targeting depressive disorder in the high-risk group confirmed in this study should be conducted.