factorial_anova improved (still work to do)