For the annual yield there was one harvest per rotation, therefore the (lme) analysis considered
genotype, site and genotype x site only, and included block as a random effect. RLD was measured in two genotypes; the (lme) was used to uncover differences between genotype, site
and depth interval and test for interactions between these variables (genotype _ site _ depth), block and core number were both considered as random effects