Thanks to Gael, I updated the scoring function of GO terms to use the grand total of genes annotated to a given GO term as the denominator of the average F score for each GO term.
This means that genes absent from the expression dataset but annotated to the GO term on ensembl will be considered to have a non-significant ANOVA result. Works very well on the AlvMac data