Dear All, Suppose I do :-
head(warpbreaks) model1<- aov(breaks ~ wool*tension,data = warpbreaks) summary(model1) There is significant interaction. So I need to test for simple effects of wool at each level of tension and vice versa. I can do a subset and then do a one way anova for each level of tension and wool. My query is that is there any command programmed which can do this automatically? Best Regards, Ashim [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.