Descriptive analytics
The final 572 participants were enrolled in this study, with a mean age of 13.77 years (SD = 3.23) and a mean BMI of 25.34 (SD = 3.53). The sample was predominantly male, with 355 male participants (62.06%). The distribution of participants across educational levels was as follows: 235 (41.08%) from primary schools, 159 (27.80%) from junior high schools, and 178 (31.11%) from high schools. Descriptive statistics for the study variables are presented in Table 1.
Table 1
Socio-demographic characteristics of participants (N = 572)
Variable
|
Categories
|
M
|
SD
|
N
|
Percentage
|
Gender
|
|
|
|
|
|
|
Male
|
|
|
355
|
62.06%
|
|
Female
|
|
|
217
|
37.94%
|
Educational Levels
|
|
|
|
|
|
|
Primary
|
|
|
235
|
41.08%
|
|
Junior
|
|
|
159
|
27.80%
|
|
High
|
|
|
178
|
31.11%
|
Age
|
|
13.77
|
3.23
|
|
|
BMI
|
|
25.34
|
3.53
|
|
|
Weight stigma
|
|
12.45
|
4.67
|
|
|
Internalized stigma
|
|
38.93
|
11.50
|
|
|
Restrained eating
|
|
27.71
|
8.06
|
|
|
Emotional eating
|
|
33.03
|
13.26
|
|
|
External eating
|
|
33.47
|
7.79
|
|
|
Depression
|
|
9.21
|
5.51
|
|
|
School bullying
|
|
2.44
|
3.71
|
|
|
Social anxiety
|
|
8.72
|
5.05
|
|
|
The correlation matrix in Table 2 presents the bivariate correlations among the study variables. All three dysregulated eating were significantly and positively related to weight stigma and internalized stigma (ps < 0.05). It is noteworthy that emotional eating was significantly and positively correlated with depression (r = 0.21, p < 0.001), while restrained eating (r = 0.07, p > 0.05) and external eating (r = 0.03, p > 0.05) were not significantly correlated with depression. The associations between dysregulated eating and school bullying also varied. Emotional eating (r = 0.20, p < 0.001) and external eating (r = 0.12, p < 0.01) were significantly and positively correlated with school bullying. However, restrained eating was not significantly correlated with school bullying (r = 0.06, p > 0.05). Additionally, weight stigma was positively correlated with internalized stigma (r = 0.56, p < 0.001), depression (r = 0.29, p < 0.001) and social anxiety (r = 0.45, p < 0.001).
Table 2
Intercorrelations among variables used in this study (N = 572)
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
1.Weight stigma
|
1
|
|
|
|
|
|
|
|
2.Internalized stigma
|
0.56***
|
1
|
|
|
|
|
|
|
3.Restrained eating
|
0.14**
|
0.30***
|
1
|
|
|
|
|
|
4.Emotional eating
|
0.36***
|
0.26***
|
0.18***
|
1
|
|
|
|
|
5.External eating
|
0.31***
|
0.15***
|
-0.01
|
0.61***
|
1
|
|
|
|
6.Depression
|
0.29***
|
0.41***
|
0.07
|
0.21***
|
0.03
|
1
|
|
|
7.School bullying
|
0.53***
|
0.37***
|
0.06
|
0.20***
|
0.12**
|
0.35***
|
1
|
|
8.Social anxiety
|
0.45***
|
0.51***
|
0.11**
|
0.31***
|
0.33***
|
0.49***
|
0.32***
|
1
|
Gender
|
-0.05
|
0.001
|
-0.05
|
0.05
|
-0.03
|
-0.01
|
-0.03
|
-0.09
|
Age
|
-0.05
|
0.01
|
0.06
|
-0.03
|
-0.03
|
0.02
|
-0.01
|
-0.06
|
Educational Levels
|
-0.09*
|
0.01
|
0.04
|
-0.04
|
-0.05
|
0.03
|
-0.01
|
-0.07
|
BMI
|
0.22***
|
0.20***
|
0.05
|
0.04
|
-0.01
|
0.11*
|
0.13**
|
0.11**
|
Note. * p < 0.05, ** p < 0.01, *** p < 0.001. |
Gaussian graphical models (undirected networks)
Figure 1A presents the Gaussian graphical model (GGM). The estimated GGM has 19 non-zero edges out of a possible 28(67.86%), with a mean edge weight of 0.10 (edge weight matrix was showed in Supplementary Table 1). The strongest edge in the network is between emotional eating and external eating (Edge weight = 0.57). Figure 1B displays the centrality indicators for each node in the GGM. Among the centrality indicators, weight stigma WS emerges as the most important node, with the highest strength (Strength = 0.91) and expected influence (Expected Influence = 1.15). The node predictability, represented by the blue rings around each node in Fig. 1A, indicates the proportion of variance in each node that can be explained by its neighboring nodes in the network. The node predictability values range from 0.14 for restrained eating to 0.50 for weight stigma.
The GGM demonstrates very good stability (refer to Supplementary Figs. 1 and 2), particularly for edge weights (CS (cor = 0.7) = 0.75), strength (CS (cor = 0.7) = 0.67) and expected influence (CS (cor = 0.7) = 0.75). The stability of closeness is acceptable (CS (cor = 0.7) = 0.44), while the stability of betweenness is insufficient (CS (cor = 0.7) = 0.21). However, as closeness and betweenness are not the primary focus of this study, their lower stability does not significantly impact the interpretation of the results.
Bayesian networks (directed networks)
Figure 2 depicts the Bayesian networks and the beta coefficients are included in Table 3. Restrained eating and internalized stigma were conditionally independent of each other given the other variables(restrained eating ⊥ internalized stigma). This means that once controlling for the other variables, internalized stigma is the only parent variable upon restrained eating. In other words, the Bayesian network analysis indicates that internalized stigma is the sole direct predictor of restrained eating (β = 0.21).
In this Bayesian network, school bullying emerges as the highest parent node, indicating that it is a key predictor of several other variables in the network. School bullying directly influences both weight stigma (β = 0.67) and depression (β = 0.34). Weight stigma, in turn, directly influences internalized stigma (β = 1.35) and emotional eating (β = 1.03). Depression, which is influenced by both school bullying and internalized stigma, also directly affects social anxiety (β = 0.32). Emotional eating, which is directly influenced by weight stigma, appears to have a direct effect on external eating (β = 0.36). This relationship suggests that children who engage in emotional eating may also be more likely to eat in response to external cues.
Weight stigma plays a pivotal role in the network, acting as a crucial mediator between school bullying and a range of negative outcomes. It directly influences internalized stigma and emotional eating, setting off a cascade of effects on other variables in the network. Under weight stigma, internalized stigma emerges as a key parent node, directly influencing three important variables: depression (β = 0.15), restrained eating (β = 0.21), and social anxiety (β = 0.14).
Finally, social anxiety emerges as the lowest child node in the network, being directly influenced by external eating (β = 0.18), depression (β = 0.32), and internalized stigma (β = 0.14). This finding highlights the complex interplay of psychological distress, dysregulated eating, and internalized stigma in the development of social anxiety among obese or overweight children.
Table 3
Arc frequency, arc direction, and beta coefficients for the Bayesian networks.
From
|
To
|
Arc frequency
|
Arc direction
|
Beta coefficient
|
Emotional eating
|
External eating
|
1.00
|
0.75
|
0.36
|
External eating
|
Social anxiety
|
0.99
|
0.83
|
0.18
|
Depression
|
Social anxiety
|
1.00
|
0.71
|
0.32
|
School bullying
|
Depression
|
0.98
|
0.80
|
0.34
|
School bullying
|
Weight stigma
|
1.00
|
0.50
|
0.67
|
Weight stigma
|
Emotional eating
|
0.93
|
0.78
|
1.03
|
Weight stigma
|
Internalized stigma
|
1.00
|
0.62
|
1.35
|
Internalized stigma
|
Restrained eating
|
1.00
|
0.51
|
0.21
|
Internalized stigma
|
Depression
|
0.87
|
0.67
|
0.15
|
Internalized stigma
|
Social anxiety
|
0.99
|
0.89
|
0.14
|
Note: arc frequency refers to the proportion of times that an arc appeared in the bootstrap samples, regardless of direction. Arc direction refers to the proportion of times that the arc pointed in that given direction. |