R help - Feb 2010 - Split Plot and Tukey

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

Hi:
this should work, (but you don't specify what is "ganhos" but
I'm guessing
is the right object...)

anava = aov(ganhos ~ Blocos + Trat*Supl +
Error(Blocos/Trat))
df<-summary(anava)[[2]][[1]][2,1]
MSerror <- summary(anava)[[2]][[1]][2,3]
HSD.test(ganhos, Trat, df , MSerror, alpha = 0.05)

Gabriela

______________________________
Lic. Mar?a Gabriela Cendoya
Mag?ster en Biometr?a
Profesor Adjunto
C?tedra de Estad?stica y Dise?o
Facultad de Ciencias Agrarias
Universidad Nacional de Mar del Plata
______________________________



----- Original Message ----- 
From: "Silvano" <silvano at uel.br>
To: <r-help at r-project.org>
Sent: Wednesday, February 17, 2010 9:15 AM
Subject: [R] Split Plot and Tukey


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

___________________________________________________________________________

Aviso:
====
El contenido del presente e-mail y sus posibles adjuntos pertenecen al INTA y
pueden contener informaci?n confidencial. Si usted no es el destinatario
original de este mensaje y por este medio pudo acceder a dicha informaci?n, por
favor solicitamos contactar al remitente y eliminar el mensaje de inmediato. Se
encuentra prohibida la divulgaci?n, copia, distribuci?n o cualquier otro uso de
la informaci?n contenida en el presente e-mail por parte de personas distintas
al destinatario.


This e-mail contents and its possible attachments belong to INTA and may contain
confidential information. If this message was not originally addressed to you,
but you have accessed to such information by this means, please contact the
sender and eliminate this message immediately. Circulation, copy, distribution,
or any other use of the information contained in this e-mail is not allowed on
part of those different from the addressee.


Antes de imprimir este mensaje, aseg?rese de que sea necesario. Proteger el
medio ambiente est? tambi?n en su mano.

Please look at the maiz example, the last example in ?MMC
in the HH package, to see how to display pairwiuse contrasts in
split-plot experiments.
Look at the entire example, all the way to the end of the help file.

Rich