Soap ANOVA with par(ps=12)