*2.4. Data Analyses*

To compare the effect of fertilizer treatment and year on grain yield, grain quality, and soil chemical composition, we used analysis of variance (ANOVA) and multivariate analysis of variance MANOVA, followed by a post-hoc analysis performed by Tukey's HSD test (Statistica 13.3, TIBCO Software, Palo Alto, CA, USA). To analyse and visualize the relationships between fertilizer treatments, years, and soil chemical properties, we used principal component analysis (PCA), performed byAnalyse-it software (Analyse-it 5.30, Analyse-it Software, Ltd., Leeds, UK).
