A total of 300 participants were included in this study, of which 202 (67.3%) were female and 98 (23.7%) were male. In terms of education level, 26 participants (8.7%) were under diploma, 69 (23.0%) had a diploma, 36 (12.0%) had an associate’s degree, 139 (46.3%) had a bachelor’s degree, 23 (7.7%) had a master’s degree, and 7 (2.3%) had a doctoral degree. Regarding marital status, 97 participants (32.3%) were single, 200 (66.7%) were married, and 3 (1.0%) were divorced. The age range of the participants was 18 - 65 years, with a mean (standard deviation) 36.16 ± 12.03.
To identify the underlying factors related to the RSS Scale, eigenvalues greater than one and the scree plot were used. A factor loading of 0.40 or higher was considered significant to confirm factor loading of an item. To assess the sampling adequacy, the Kaiser-Meyer-Olkin (KMO) measure of sampling adequacy was employed (
33). In assessing the EFA results, the KMO value was found to be 0.80, indicating an adequate sampling (> 0.70). Bartlett's test of sphericity was also significant (χ
2 = 1331.8, df = 91, P < 0.01). Examination of the scree plot suggested that factor loadings had a greater weight, and there was a sharp drop in eigenvalues after the fourth factor (
Figure 1).
Scree plot for factor analysis
Based on
Table 1, the items of the RSS Scale loaded onto four factors, overall accounting for 61.65% of the total variance.
| Factor | Eigenvalues | Explained Variance | Cumulative Variance | Mean ± SD | Skewness | Kurtosis |
|---|
| 1 | 4.4 | 31.47 | 31.47 | 2.43 ± 1.41 | 1.86 | 1.79 |
| 2 | 1.79 | 12.83 | 44.3 | 4.36 ± 1.52 | 1.55 | 1.67 |
| 3 | 1.35 | 9.67 | 53.97 | 1.01 ± 1.55 | 1.15 | 1.7 |
| 4 | 1.07 | 7.67 | 61.65 | 2.01 ± 1.97 | 0.73 | -0.29 |
Based on
Table 2, a varimax rotation of the items on the RSS Scale revealed a four-factor solution. The first factor consisted of items 1, 2, 3, 6, and 8, and was labeled “feeling religiously rejected”. The second factor, comprised of items 7, 9, 10, 11, and 12, was labeled “ambiguous religious/spiritual identity”. The third factor consisted of items 4 and 5, and was labeled “feeling threatened by the devil”. Finally, the fourth factor, consisted of items 13 and 14, was labeled “challenging the meaning of life”.
| Items | Factor Loading | Mean ± SD |
|---|
| First | Second | Third | Fourth |
|---|
| 1. I have felt that God has abandoned me. | 0.78 | -0.007 | 0.05 | 0.26 | 0.32 ± 0.76 |
| 2. I have felt that God is punishing me. | 0.73 | 0.04 | 0.20 | 0.27 | 0.62 ± 0.93 |
| 3. I have felt anger toward God. | 0.76 | 0.13 | 0.001 | 0.12 | 0.35 ± 0.74 |
| 4. I have felt harmed, mistreated, or insulted by religious/spiritual people. | 0.57 | 0.53 | -0.02 | -0.10 | 0.66 ± 1.09 |
| 5. I have felt anger toward organized religion. | 0.58 | 0.56 | -0.14 | -0.04 | 0.48 ± 1.02 |
| 6. I have had conflicts with others regarding religious issues. | 0.30 | 0.71 | -0.02 | -0.02 | 1.21 ± 1.07 |
| 7. I have struggled with myself regarding moral principles. | 0.06 | 0.70 | 0.12 | 0.17 | 0.94 ± 1.1 |
| 8. I have felt guilty for not adhering to my moral standards. | 0.02 | 0.53 | 0.22 | 0.23 | 1.05 ± 1.1 |
| 9. I have been distressed by doubting or questioning my faith or spirituality. | -0.15 | 0.46 | 0.17 | 0.28 | 0.63 ± 1.1 |
| 10. I have felt confused about my religious/spiritual beliefs. | 0.13 | 0.62 | 0.09 | 0.38 | 0.53 ± 0.88 |
| 11. I have worried that the devil or evil spirits caused my problems. | 0.01 | 0.09 | 0.89 | 0.04 | 0.59 ± 0.94 |
| 12. I have felt under attack by the devil or evil spirits. | 0.10 | 0.15 | 0.88 | 0.01 | 0.42 ± 0.74 |
| 13. I have questioned whether life really matters. | 0.31 | 0.15 | -0.03 | 0.74 | 1.29 ± 1.3 |
| 14. I have felt a lack of deeper meaning in my life. | 0.16 | 0.23 | 0.02 | 0.72 | 0.73 ± 1.004 |
The CFA was conducted to validate the findings of the EFA. The CFA deals with assessing the factor loading of each item, mutual correlations between factors, and fit indices to evaluate the overall model fit (
33). Maximum likelihood estimation (MLE) was employed to assess the model fit indices (
Figure 2).
Confirmatory factor analysis (CFA) of the four-factor model of the Religious/Spiritual Struggles Scale (initial model)
Based on the results, factor loadings of each item in its respective factor were significant at the 0.01 level (P < 0.001). Additionally, the mutual correlations or relationships between factors were significant at the 0.01 level (P < 0.001). Subsequently, the fit indices of the initial model were examined. The most prominent classification for organizing fit indices was proposed by Hair et al. (
34), who categorized fit indices into three groups: Absolute, relative, and parsimonious fit indices.
Based on
Table 3, the fit indices of the initial model, specifically the Goodness-of-Fit Index (GFI), root mean square error of approximation (RMSEA), Comparative Fit Index (CFI), Normed Fit Index (NFI), Incremental Fit Index (IFI), and Relative Fit Index (RFI), were unsatisfactory compared to the cut-off points and were below the acceptable range. To modify model fit, according to AMOS-24 suggestion, correlation coefficients (mutual pathway) were established between the errors of 1 and 2, 1 and 3, 2 and 3, and 9 and 10 covariances. The modified model results indicate that the NFI was marginally acceptable, with a difference of 0.02, and the RFI was also marginally acceptable, with a difference of 0.04. Other fit indices indicated improvement compared to the initial model and were within the acceptable range. The inflation of PCFI and PNFI is expected due to the complexity of the model (
Figure 3).
| Fit Indices | Initial Model | Modified Model | Cut-off Point (Norm) |
|---|
| Absolute | | | |
| χ2/df | 2.81 | 2.22 | 5 - 1 |
| GFI | 0.88 | 0.93 | > 90 |
| RMSEA | 0.09 | 0.06 | < 0.08 |
| CFI | 0.85 | 0.93 | > 0.90 |
| Relative | | | |
| NFI | 0.81 | 0.88 | > 0.90 |
| IFI | 0.85 | 0.93 | > 0.90 |
| RFI | 0.75 | 0.86 | > 0.90 |
| Parsimonious | | | |
| PCFI | 0.63 | 0.64 | < 0.50 |
| PNFI | 0.65 | 0.68 | < 0.50 |
Abbreviations: GFI, Goodness Fit Index; RMSEA, root mean square error of approximation; CFI, Comparative Fit Index; NFI, Normed Fit Index; IFI, Incremental Fit Index; RFI, Relative Fit Index; PCFI, Parsimonious Comparative Fit Index; PNFI, Parsimonious Normed Fit Index.
Confirmatory factor analysis (CFA) of the four-factor model of the Religious/Spiritual Struggles Scale (modified model)
Based on
Table 4, the correlation coefficients of the factors of feeling religiously rejected, ambiguous religious/spiritual identity, challenging the meaning of life, and the total score with spiritual well-being were -0.64, -0.27, -0.48, and -0.55, respectively, which were significant at the P < 0.01 level. The correlation coefficients of the factors of ambiguous religious/spiritual identity, challenging the meaning of life, and the total score with perceived stress were 0.14, 0.14, and 0.20, respectively, which were significant at the P < 0.05 level (P < 0.1). The Cronbach's alpha value, exceeding 0.70, indicated satisfactory internal consistency among the scale items. As shown in
Table 4, the Cronbach's alpha values for the factors and items ranged from 0.70 to 0.82. The CR values, exceeding 0.70, indicated satisfactory convergent validity (
35). In
Table 4, these values ranged from 0.80 to 0.91. Furthermore, the mean variance extracted values, exceeding 0.50, indicated high convergent validity of the scale, meaning that the factor in question is highly correlated with the total scale (
36). According to the results in the
Table 4, these values ranged from 0.51 to 0.83.
| Factors | Divergent Validity (Spiritual Well-being) | Convergent Validity (Perceived Stress) | Cronbach's Alpha | CR | AVE |
|---|
| Feeling religiously rejected | -0.64 a | 0.09 | 0.79 | 0.86 | 0.55 |
| Ambiguous religious/spiritual identity | -0.27 a | 0.14 b | 0.70 | 0.80 | 0.51 |
| Feeling threatened by the devil | -0.03 | 0.09 | 0.79 | 0.91 | 0.83 |
| Challenging the meaning of life | -0.48 a | 0.14 b | 0.70 | 0.84 | 0.73 |
| Total score | -0.55 a | 0.20 b | 0.82 | - | - |
Abbreviation: AVE, average variance extracted; CR, composite reliability.
a P < 0.05.
b P < 0.01.