The final model was built on the basis of the stepwise forward procedure (all two way interactions were tested), always according to a better Akaike’s Information Criterion (AIC) value and considering that the model with the lower AIC fit the data better [20].