I have been trying to do tukey's test to no avail. This is what I have and the error produced:
pen.df = data.frame(blend, treatment, y) source("tukey.1.r") tukey.1(aov.pen, pen.df) Error in tukey.1(aov.pen, pen.df) : the model must be two-way This is a two-way design. Therefore I am confused. Can anyone help? Raphael > pen.df blend treatment y 1 1 A 89 2 2 A 84 3 3 A 81 4 4 A 87 5 5 A 79 6 1 B 88 7 2 B 77 8 3 B 87 9 4 B 92 10 5 B 81 11 1 C 97 12 2 C 92 13 3 C 87 14 4 C 89 15 5 C 80 16 1 D 94 17 2 D 79 18 3 D 85 19 4 D 84 20 5 D 88 ______________________________________________ R-help@r-project.org mailing list 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.