see
two_way_anova_factorial_anova