[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

______________________________________________
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.