Olá Tiago eu uso este programa para fazer o teste de médias.

# Teste Tukey
library(agricolae)

# Tratameto
df<- df.residual(anava1)
MSerror<-(deviance(anava1))/df
comparison <- HSD.test(QEA,Tratamento,df,MSerror, group=TRUE)

# Interação
df<- df.residual(mod)
MSerror<-(deviance(mod))/df
comparison <- HSD.test(QEA,Tratamento:Bloco,df,MSerror, group=TRUE

Daniele Portela de Oliveira

Mestranda Produção Animal - UEM

Cel: 44 9949 3525

67 9617 0965

De: Thiago De paula protásio <depaulaprotasio@yahoo.com.br>
Para: r-br <r-br@listas.c3sl.ufpr.br>
Enviadas: Sábado, 4 de Junho de 2011 11:42:18
Assunto: [R-br] Teste de Tukey! Experimento Desbalanceado!

Pessoal;

Venho novamente pedir a ajuda de vcs!

Tenho que analisar um experimento em blocos, mas os tratamentos possuem diferentes números de repetições. Estou utilizando a soma de quadrados do tipo III.

library(car)

anava1=aov(QEA~Tratamento+Bloco)

Anova(anava1,type="III")

Contudo, não estou conseguindo efetuar o teste de Tukey. Alguém pode me ajudar?

Muito obrigado!

Thiago de Paula Protásio
Acadêmico de Engenharia Florestal
Universidade Federal de Lavras
Laboratório de Energia da Biomassa Florestal
Laboratório de Biomateriais
(035) 9183-2246