This study was a classical twin study. The study population was all twins living in Isfahan in the summer of 2021. In total, 174 participants 16 to 50 years old, twins enrolled in this sub-study. Subjects were obtained from Isfahan Twin Registry (ITR) (22). At the beginning of the online research questionnaire, a consent form was included that participants could only participate in the project after reading the consent form, if they agreed with the research conditions. Isfahan Twins Registry research has already been approved by the Research Ethics Committee of Isfahan University of Medical Sciences. The ethical code of this research is IR.mui.med.rec.1399.169.
The inclusion criteria were being at age range pairs being accessible, and living together for at least 10 years. Reluctance to participate in the research, absence of one of the twins, and not living together for more than 10 years were among the exclusion criteria of the research. In order to collect data first, a public call was announced by the ITR, then the responders who met the inclusion criteria were included in the research, and questionnaires were distributed among them. The research questionnaires included the consent form, the zygotic determination questionnaire, and the Hewitt and Flett Multidimensional Perfectionism Scale (HFMPS).
Instruments
Hewitt and Flett Multidimensional Perfectionism scale (HFMPS)
The multidimensional scale of Perfectionism was created by Hewitt and Flett in 1991. This scale consists of 30 items that measure the three components of self-oriented perfectionism, other-oriented perfectionism, and socially prescribed perfectionism. Each component consists of 10 statements that are answered by a 5-point Likert scale, ranging from "Strongly disagree" to "Strongly agree". In a validation study, in a sample of 480 Iranian students, Besharat (23) in 2005 reported adequate and satisfactory internal consistency, the Cronbach’s Alfa for which was 0.89 for self-oriented perfectionism, 0.83 for other-oriented perfectionism, and 0.78 for socially prescribed (23). In a more recent study the Cronbach’s Alfa of SOP, OOP, and SPP reported 0.86, 0.82, and 0.87 respectively (24). Also In this study, the Cronbach’s Alfa of SOP, OOP, and SPP were 0.82, 0.70, and 0.86.
Zygosity questionnaire
This questionnaire was developed by Song et al. (22) in 2010. They first collected the genetic information of 2060 twins (1030 pairs). Then the research questionnaire was completed by the twins. Finally, it was found that the accuracy of this questionnaire to determine the identical twins was 97.2% and to determine the fraternal twins was 95%. This questionnaire identifies the identical and fraternal twins by using three questions. The first question determines whether the parents of the twins confuse them with each other. The second question determines whether twins are mistaken for teachers and their friends. And the final question is whether strangers confuse twins (22).
Statistical analysis
After this step, data were entered into SPSS software. Descriptive analysis was performed in this software and then the data were transferred to R software.
The Kolmogorov-Smirnov test was then applied to the data to confirm the hypothesis of homogeneity of distribution. Two tests of intra-pair correlation coefficient and homogeneity of variance were performed. ACE structural equation model used to calculate the share of three components of heritability (additive and non-additive genetics), common environment, and specific environment. The calculation of these components is done by the correlation between MZ (monozygotic) and DZ (dizygotic) twins. The packages that have been used in R software were lme4, mets, psych, and nortest. In R software, first, the coefficient of intra-pair correlation was determined. Intra-pair correlation measures the correlation of pairs in each zygosity group. If the coefficient of intra-pair correlation in MZ twins is greater than that of DZ twins, it is estimated that the effect of genetics on the phenotype is greater than the environment. If the coefficient of intra-pair correlation is equal to or higher in the fraternal twins, then it is estimated that the influence of the environment on that phenotype is greater than the impact of genetics. In the next step, the F test was performed to compare the pairs of MZ and DZ twins. The purpose of this step was to ensure homogeneity of variance among the twins in each zygosity group. After this model, ACE was performed on the data. Also, after the ACE model AE, CE, E, ADE, and DE models were performed on them, respectively. After performing the ACE model, model fit was determined using Akaike's Information Criterion (AIC) and Bayesian Information Criterion (BIC). The AIC is used to select the best fit and the BIC is used to select the model that has the least amount of lost information. Better fitting models have smaller values. Also, after this test, the − 2L of models were compared with the reference models (ACE & ADE) using chi-square. It is recommended that models AE, CE, and E be compared with the ACE model and the DE model with the ADE model (15).