Statistical Analysis
Sample size was calculated considering the results obtained during our previous observationalThe results were summarized with counts and percentages for categorical variables and mean and SD for normally distributed quantitative variables. For variables with asymmetric distribution, we also reported the median and IQR. The data were analyzed with generalized linear models considering the analysis of variance for quantitative variables, binomial distribution for binary data, generalized logit for multinomial nonordinal data, and cumulative logit for multinomial ordinal data (eAppendix 2 in 