[R] Split Plot and Tukey

[Silvano]

Hi,

I did the analysis of variance of a split-plot and the 
effect of treatment was significant.
I would like compare treatment means using Tukey.
I can't extract the mean square to apply HSD.test to use in 
agricolae package.

anava = aov(ganhos ~ Blocos + Trat*Supl + 
Error(Blocos/Trat))
names(anava)
summary(anava)

require(agricolae)
HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

Thanks

--------------------------------------
Silvano Cesar da Costa
Departamento de Estatística
Universidade Estadual de Londrina
Fone: 3371-4346