Indices' data were analyzed through linear or generalized linear mixed models taking into account the potential correlations between multiple measurements taken on the same individual. In particular, for the analysis of EHI and PI, ordinal logistic regression models were used. For the PA% index, which was a proportion ranging from 0% to 100%, a linear model was used after a logit transformation of the dependent variable. Model selection, including checks for potential interaction effects, for the final multivariable models was based on likelihood ratio tests.