The association between genotypes and the Z2E(j) statistic was then evaluated with Bartlett's k-sample variance test [30], using the 'bartlett.test' function in R to find significant differences in the variance of Z2E(j) across each genotype using additive coding (i.e. AA, AB, BB). Resulting p-values were then adjusted for FDR using the 'p.adjust' function in R and the method of Benjamini and Hochberg [31].
