Table 1 presents demographic and clinical characteristics of participants in this research.
| Variables | Percentage (%) |
|---|
| Age (y) | |
| 20 - 30 | 35.3 |
| 31 - 45 | 64.7 |
| Education | |
| High school diploma or lower | 56.9 |
| University degree | 43.1 |
| Occupation | |
| Housewife | 70.6 |
| Employee | 29.4 |
| Having children | |
| No children | 69.3 |
| Have children | 21.7 |
| Psychological treatment | |
| No | 90.8 |
| Yes | 9.2 |
| Cause of infertility | |
| Known | 77.0 |
| Unknown | 33.0 |
| Duration of treatment (y) | |
| 1 - 5 | 80.0 |
| More than 5 | 20.0 |
| Type of infertility | |
| Primary | 48.0 |
| Secondary | 52.0 |
| Type of treatment | |
| Drug therapy | 49.5 |
| IVF | 23.4 |
| Stimulated ovulation | 11.0 |
| IUI | 21.0 |
| Economic status (million Tomans/mo) | |
| ≤ 10 | 89.0 |
| > 10 | 11.0 |
Based on
Table 1, most participants were aged 31 - 45 years, and over half held a high school diploma or lower education. The majority were housewives and childless, and most had not received psychological treatment. In most cases, the cause of infertility was known, and the majority had undergone treatment for one to five years. The types of infertility were nearly evenly distributed between primary and secondary. Drug therapy was the most common treatment, followed by IVF and IUI. Most participants belonged to a lower economic status, earning less than 10 million Tomans per month.
Table 2 presents means and standard deviations as the descriptive indices of and standardized (1 - 100) scores of the research variables.
| General Factors and Subscales | Scores in Terms of 100 | Mean ± Standard Deviation |
|---|
| Relational resilience | | |
| Partner | 76.6 | 33.37 ± 6.72 |
| Spirituality | 73.8 | 27.18 ± 5.02 |
| Union | 61.0 | 46.4 ± 9.29 |
| Actor | 80.0 | 35 ± 5.76 |
| Relational resilience | 70.8 | 141.96 ± 23.1 |
| Infertility stress | | |
| Social concern | 47.5 | 33.39 ± 8.4 |
| Sexual concern | 48.8 | 26.73 ± 7.27 |
| Relationship concern | 49.5 | 34.29 ± 7.4 |
| Need for parenthood | 66.0 | 44.27 ± 10.19 |
| Rejection of a child-free lifestyle | 81.0 | 39.9 ± 6.97 |
| Infertility stress | 58.7 | 178.61 ± 26.06 |
| Dyadic coping | | |
| Self-coping | 86.3 | 49.36 ± 5.88 |
| Partner coping | 79.9 | 47.43 ± 7.43 |
| Dyadic coping | 85.2 | 96.79 ± 11.35 |
According to
Table 2, within the relational resilience subscales, 'Actor' (standardized score = 80.0) exhibited the highest level, while 'Union' (61.0) was the lowest. In the infertility stress subscales, 'Rejection of a child-free lifestyle' (81.0) was the highest, and 'Social Concern' (47.5) was the lowest. For dyadic coping, 'Self-Coping' (86.3) had the highest standardized score, with 'Partner Coping' (79.9) being slightly lower. Standardized scores facilitate direct comparison across subscales, with higher values reflecting higher levels of the respective construct. In the following,
Table 3 presents the correlation coefficients of the research variables.
| Variables | Infertility Stress | Need for Parenthood | Relationship Concern | Sexual Concern | Social Concern | Dyadic Coping | Other Coping | Self-coping | Relational Resilience | Spirituality | Union | Partner | Actor |
|---|
| Infertility stress | 1 | - | - | - | - | - | - | - | - | - | - | - | - |
| Need for parenthood | 0.75 a | 1 | - | - | - | - | - | - | - | - | - | - | - |
| Relationship concern | 0.70 a | 0.40 a | 1 | - | - | - | - | - | - | - | - | - | - |
| Sexual concern | 0.68 a | 0.29 a | 0.50 a | 1 | - | - | - | - | - | - | - | - | - |
| Social concern | 0.64 a | 0.24 a | 0.42 a | 0.39 a | 1 | - | - | - | - | - | - | - | - |
| Dyadic coping | -0.29 a | -0.19 a | -0.21 a | -0.24 a | -0.14 b | 1 | - | - | - | - | - | - | - |
| Partner coping | -0.33 a | -0.26 a | -0.23 a | -0.24 a | -0.17 b | 0.88 a | 1 | - | - | - | - | - | - |
| Self-coping | -0.14 b | -0.05 | -0.10 | -0.17 a | -0.05 | 0.81 a | 0.44 a | 1 | - | - | - | - | - |
| Relational resilience | -0.34 a | -0.14 b | -0.36 a | -0.37 a | -0.21 a | 0.52 a | 0.53 a | 0.33 a | 1 | - | - | - | - |
| Spirituality | -0.11 | 0.02 | -0.21 a | -0.17 a | -0.09 | 0.23 a | 0.22 a | 0.16 b | 0.69 a | 1 | - | - | - |
| Union | -0.39 a | -0.23 a | -0.37 a | -0.39 a | -0.20 a | 0.48 a | 0.50 a | 0.30 a | 0.93 a | 0.51 a | 1 | - | - |
| Partner | -0.31 a | -0.15 b | -0.34 a | -0.33 a | -0.19 a | 0.53 a | 0.58 a | 0.29 a | 0.90 a | 0.51 a | 0.80 a | 1 | - |
| Actor | -0.25 a | -0.05 | -0.26 a | -0.31 a | -0.21 a | 0.49 a | 0.46 a | 0.36 a | 0.85 a | 0.48 a | 0.72 a | 0.69 a | 1 |
The results showed that infertility stress (social anxiety, sexual anxiety, and relationship anxiety) had a negative and significant relationship with dyadic coping and its subscales, except for the relationship of need for parenthood, relationship anxiety, and social anxiety with self-coping, which was not significant. The correlation between infertility stress and relational resilience was negative and significant. The relationship between conflict stress and relational resilience subscales (cohesion, partner, and actor) was also negative and significant, although its relationship with spirituality was not significant. All subscales of infertility stress had a negative and significant relationship with relational resilience. Relationship anxiety and sexual anxiety had a negative and significant relationship with all subscales of relational resilience. Social anxiety was negatively and significantly correlated with alliance, partner, and actor, but not with spirituality. The need for parenthood had a negative and significant relationship with alliance and partner, but not with spirituality or actor. Dyadic coping was positively and significantly correlated with relational resilience and all its subscales. All subscales of dyadic coping were positively and significantly associated with relational resilience and its subscales.
Before testing model fit, structural equation assumptions were assessed. The Kolmogorov-Smirnov and Shapiro-Wilk tests for the dependent variable, infertility stress, indicated that the data distribution was normal. Model fitting showed a good fit to the data.
According to
Table 4, the initial mediation model demonstrated an acceptable fit based on both Absolute and Comparative Fit indices (CFI, GFI, and TLI). However, the RFI was slightly below the recommended threshold of 0.90, suggesting that some modifications could improve the overall model fit.
| Modification | NFI | RFI | IFI | TLI | CFI | GFI | Pratio | RMSEA | X2/df |
|---|
| Before | 0.88 | 0.85 | 0.92 | 0.90 | 0.92 | 0.91 | 0.74 | 0.08 | 2.54 |
| After | 0.92 | 0.90 | 0.96 | 0.95 | 0.96 | 0.94 | 0.71 | 0.07 | 1.95 |
Abbreviations: TLI, Tucker-Lewis Index; CFI, Comparative Fit Index; GFI, Goodness of Fit Index.
During model evaluation, it was found that the 'Rejection of a child-free lifestyle' subscale showed a weak correlation with the total infertility stress score. Given this poor correlation, the subscale was removed from the model for greater parsimony and theoretical consistency. The modified model was re-estimated to assess improvements in relationships among the remaining variables. In the revised model, the mediating role of relational resilience in the association between dyadic coping and infertility stress was tested among women undergoing infertility treatment (
Figure 1). The modified model demonstrated favorable and improved fit indices compared to the initial model.
The final model represents a complete mediation model, with the non-significant direct path from dyadic coping to infertility stress removed. This modification yielded a higher overall fit, confirming that the effect of dyadic coping on infertility stress is fully mediated by relational resilience. In other words, dyadic coping indirectly influences infertility stress by enhancing relational resilience, which then reduces infertility-related stress.
Figure 1 illustrates the study variables and their respective subscales. Dyadic coping (including self-coping and partner coping) was positively and significantly associated with relationship resilience (union, partner, spirituality, and actor subscales). The correlation coefficient between these two variables was 0.68, indicating that higher dyadic coping is associated with greater relational resilience among women undergoing infertility treatment.
Conversely, the relationship between dyadic coping and infertility stress (including social concern, sexual concern, relationship concern, and need for parenthood) was negative but not significant (R = -0.08). Due to this lack of significance, the direct path was removed from the final model, resulting in a simpler and better-fitting structure. There was a negative and significant relationship between relationship resilience and infertility stress (R = -0.51).
Table 5 shows the standardized and unstandardized regression weights of the model.
| Regression Weights | Standardize Estimate | P | Unstandardized Estimate |
|---|
| Dyadic coping → Relational resilience | 0.67 | 0.001 | 3.52 |
| Relational resilience → Infertility stress | -0.51 | 0.001 | -0.280 |
| Infertility stress → Social concern | 0.54 | 0.001 | 1 |
| Relational resilience → Union | 0.90 | 0.001 | 1 |
| Relational resilience → Actor | 0.79 | 0.001 | 0.54 |
| Relational resilience → Spirituality | 0.56 | 0.001 | 0.33 |
| Relational resilience → Partner | 0.89 | 0.001 | 0.71 |
| Dyadic coping → Partner coping | 0.87 | 0.001 | 3.99 |
| Dyadic coping → Self-coping | 0.51 | 0.001 | 1.87 |
| Infertility stress → Relationship concern | 0.76 | 0.001 | 1.23 |
| Infertility stress → Sexual concern | 0.68 | 0.001 | 1.08 |
| Infertility stress → Need for parenthood | 0.47 | 0.001 | 1.05 |
The results indicated that dyadic coping had a positive and significant effect on relational resilience (β = 0.67, P < 0.001), and relational resilience had a negative and significant effect on infertility stress (β = -0.51, P < 0.001). Relational resilience was positively and significantly associated with all its subscales (union, actor, spirituality, and partner), while dyadic coping was positively and significantly related to its subscales, self-coping and partner coping. Furthermore, infertility stress was positively and significantly associated with all its components, including social concern, relationship concern, sexual concern, and need for parenthood. These findings support the good fit of the model and confirm the significance of the hypothesized relationships among the main research variables.
According to
Table 6, the standardized estimate of the indirect effect of dyadic coping on infertility stress, with lower and upper limits of -0.44 and -0.15, respectively, was -0.30* (based on bootstrap sampling, n = 2000, confidence level = 0.05, P = 0.008). This indicates a significant indirect effect of positive dyadic coping on infertility stress via relational resilience.
| Variables | Standard Estimate | Lower Limit | Upper Limit | P-Value |
|---|
| Indirect effect of dyadic coping on infertility stress | -0.30 a | -0.44 | -0.15 | 0.008 |