“…Given the small sample sizes used in this study, a non-parametric ANOVA equivalent to a two-way general linear model, the Scheirer-Ray-Hare test (Scheirer et al 1976, Bonaglia et al 2017, was used to test the effects of treatment, species, season, and year, as well as their interactions on A max , g s , C i , T r , LSP, and LCP. To further detect differences in these parameters among treatments, a non-parametric test (Kruskal-Wallis test) was performed in SPSS 16.0 (Chicago, Illinois, USA).…”