more work on factorial.anova.Rmd