Data analysis
Three participants who did not attend more than one intervention session or did not do the assigned homework, were excluded from the data analysis. Descriptive statistics included frequency, percentage, mean, and Standard Deviation (SD). At first, the Kolmogorov–Smirnov test confirmed the normal distribution of the dependent variables. One-way ANOVA, chi-square, and Kruskal–Wallis tests were used to determine the homogeneity of the study groups by comparing the participants’ characteristics and scores. Also, an independent t-test was used to determine the probable differences between the participants who took the pre-test at baseline. In the post-test, One-way ANOVA was utilized to assess the differences among the four study groups regarding the mean scores of the study variables. Then, Tukey’s post-hoc test was employed for pairwise comparisons. In the groups taking part in the pre-test, changes in the outcome variable from the pre-test to the post-test were evaluated using paired t-test. All tests were two-tailed, and the statistical level was considered 0.05. All data analyses were done using the SPSS 22 software.