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

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

Thiago, Tem um procedimento para isso nas Ridículas do LEG http://www.leg.ufpr.br/doku.php/ridiculas#experimento_com_dois_fatores_de_ef... Não se pode usar as médias amostrais como estimadoras das médias populacionais porque devido à não ortogonalidade elas são viesadas. Por isso o uso da agricolae::HSD.test() não é recomendado. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (3)
-
Daniele Portela de oliveira
-
Thiago De paula protásio
-
Walmes Zeviani