R statistical environment version 3

Data were collected in a dedicated database. A descriptive analysis was performed for sample description and frequencies and for IES-R and HADS subscales results (HADS-A, HADS-D, and HADS-General). A correlation analysis was implemented to detect the univariate associations between sociodemographic characteristics, clinical data, and the IES-R score, as well as the subscales of the HADS. The correlation between the continuous variables was calculated using Pearson's Correlation Coefficient (r). Eta squared was used for comparison between categorical variables and continuous variables. Student t-test was performed to assess the difference between the mean value in two groups. Statistical significance was achieved at a p < 0.05. All statistical analyses were performed using the using the R statistical environment, version 3.5.2 (The R Foundation for Statistical Computing; Vienna, Austria).

Descriptive statistical analyses of study sample were performed, and their results were expressed as mean ± standard error of the mean (SEM) and frequencies (%), depending on the nature of variables. Comparisons of socio-demographic and clinical parameters between two groups of interest were carried out with the t-test (for continuous variables) and the Pearson χ² test (for categorical variables). Comparisons between more than two groups were studied through one-way analysis of variance (ANOVA) followed by Bonferroni post-hoc test, where necessary. The correlation between continuous variables was also considered and calculated using Pearson’s Correlation Coefficient (r). P-values lower than 0.05 were regarded as significant. All the analyses were performed using the R statistical environment, version 3.5.2 (The R Foundation for Statistical Computing; Vienna, Austria), NCSS 2007 Statistical Software, version 0.7.10 (NCSS, LLC Company) and GraphPad Prism, version 7.04 (GraphPad Software; San Diego, USA).