[R] tukey.1

[Raphael Fraser]

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