The mean age of the male and female students was 21.12 ± 1.43 and 20.73 ± 1.84 years, respectively. The students were divided into two groups based on their gender, including 178 males (48%) and 194 females (52%). Causal models are analyzed based on the correlation matrix.
Table 1 shows the participants' descriptive indices (mean, SD, minimum and maximum score), as well as the results of Pearson’s correlation-coefficient to evaluate the correlation between the predictive, mediating, and criterion variables. Accordingly, all the variables demonstrated significant correlations (
Table 1). SEM was employed to simultaneously test the model of the hypothesized correlations. After the descriptive examination of the data, the statistical samples were divided into two groups based on their gender.
| Variables | Mean ± SD | Max. | Min. | 1 | 2 | 3 |
|---|
| Female | | | | | | |
| 1- Perceived social support | 47.11 ± 5.11 | 37 | 52 | 1 | | |
| 2- Self-compassion | 38.11 ± 6.67 | 58 | 111 | 0.36 a | 1 | |
| 3- Psychosomatic symptoms | 43.32 ± 6.34 | 21 | 58 | -0.27 a | -0.54 a | 1 |
| Male | | | | | | |
| 1- Perceived social support | 51.13 ± 6.21 | 42 | 55 | 1 | | |
| 2- Self-compassion | 46.14 ± 7.49 | 69 | 113 | 0.32 a | 1 | |
| 3- Psychosomatic symptoms | 38.61 ± 4.59 | 23 | 55 | -0.28 a | -0.56 a | 1 |
| Total | | | | | | |
| 1- Perceived social support | 49.55 ± 6.09 | 37 | 55 | 1 | | |
| 2- Self-compassion | 41.95 ± 7.42 | 58 | 113 | 0.34 a | 1 | |
| 3- Psychosomatic symptoms | 51.06 ± 5.75 | 21 | 58 | -0.27 a | -0.55 a | 1 |
Prior to data analysis, the assumptions of multivariate normality, linearity, multicollinearity, and independence of errors were tested and confirmed. Skewness and kurtosis tests were also performed to examine the normality of data distribution. As the absolute magnitudes of the skewness and kurtosis coefficients were respectively < 3 and < 10 for all the variables, the assumption of the normality of data distribution was confirmed in both groups, as well as the entire sample population. To assess multicollinearity, the tolerance statistic and variance inflation factor (VIF) were used. For all the variables, the VIF was < 10, and the tolerance statistic was > 0.1, which justified the multicollinearity assumption. The assumption of the independence of errors was also tested using the Durbin-Watson statistic, which was estimated at 2.05 for women and 1.79 for men; as this value was within the range of 1.5 - 2.5, the assumption held.
In this study, three models were examined, including the general model, model of the female students, and model of the male students.
Table 2 shows the fitness of the proposed and final models based on goodness-of-fit indices. Accordingly, the proposed model did not have a proper fit based on the root mean square error of approximation indicator (
Table 2). Therefore, the model was modified by correlating the errors of the latent variables, and the results indicated the proper fit of the final model.
Figure 1 depicts the standard path coefficients in the final model for the entire sample population.
| Fit Indicators | χ2 | df | (χ2/df) | IFI | TLI | CFI | NFI | RMSEA | PCLOSE |
|---|
| Initial model | 26.88 | 8 | 3.36 | 0.95 | 0.91 | 0.95 | 0.93 | 0.11 | 0.01 |
| Final model | 18.71 | 7 | 2.67 | 0.97 | 0.97 | 0.97 | 0.95 | 0.08 | 0.06 |
The standard path coefficients in the final model for the entire sample
Table 3 shows the structural model, paths, and standard coefficients of the final model for the entire sample population. Based on standard parameter coefficients and the corresponding significance level (
Table 3), all the paths of the final model were significant. Therefore, a significant, negative correlation was observed between perceived social support and psychosomatic symptoms (β = -0.14; P = 0.03). Moreover, perceived social support had a significant, positive correlation with self-compassion (β = 0.34; P = 0.001). The correlation between self-compassion and psychosomatic symptoms was also significant and negative (β = -0.52; P = 0.001).
| Path | The Final Model |
|---|
| Path Type | β | P |
|---|
| Perceived social support to psychosomatic symptoms | Direct | -0.14 | 0.03 |
| Self-compassion to psychosomatic symptoms | Direct | -0.52 | 0.001 |
| Perceived social support to self-compassion | Direct | 0.34 | 0.001 |
In the current research, percentile bootstrap and bias-corrected bootstrap tests were performed to determine the significance of the indirect correlations between the variables and examine the mediating role of self-compassion (
Table 4). To this end, the total effect of the predictive variable on the criterion variable was initially determined in the absence of the mediator. The obtained results indicated a significant correlation between perceived social support and psychosomatic symptoms without the mediator. In the indirect path and with the mediator variable, the correlation between perceived social support and psychosomatic symptoms was considered significant as well (
Table 4). Therefore, it could be concluded that self-compassion is partially affected by the impact of perceived social support on psychosomatic symptoms with a slightly mediating role in this regard.
| Path | Indirect Effect | Total Effect |
|---|
| β | P | β | P |
|---|
| Perceived social support to psychosomatic symptoms through the mediating role of self-compassion | -0.19 | 0.001 | -0.27 | 0.01 |
In the present study, multigroup analysis was performed to compare the two models between the male and female students. The invariance of the measurement models in the female and male groups was evaluated by the comparison of the unrestricted model (model A) and the weight-restricted model (model B).
Table 5 shows the unrestricted path coefficients of the male and female groups. Accordingly, all the paths were significant in both groups.
| Path | Female Students | Male Students |
|---|
| β | P | β | P |
|---|
| Perceived social support to psychosomatic symptoms | -0.14 | 0.04 | -0.15 | 0.02 |
| Self-compassion to psychosomatic symptoms | -0.54 | 0.001 | -0.50 | 0.001 |
| Perceived social support to self-compassion | 0.35 | 0.001 | 0.33 | 0.001 |
Figures 2 and
3 illustrate the unrestricted model (A) in both groups.
Table 6 shows the comparison of the goodness-of-fit indicators of the unrestricted (A) and restricted (B) adjusted models in the two groups. Based on these indicators, the moderating role of gender was not confirmed. In other words, no difference was observed in any of the paths between the male and female students, and gender did not moderate this model.
The standard path coefficients in the final model for the female students
The standard path coefficients in the final model for the male students
| Fit Indicators | χ2 | df | GFI | CFI | RMSEA | ∆χ2 | ∆df |
|---|
| Entire sample | 1.87 | 1 | 0.99 | 0.99 | 0.05 | 0.28 | 3 |
| Unrestricted model (A) | 41.25 | 16 | 0.97 | 0.96 | 0.07 |
| Restricted model (B) | 41.53 | 19 | 0.97 | 0.97 | 0.06 |