Teste média Subdividida

Bom dia, Estou com dificuldade para realizar o teste de média para um experimento em parcelas subdivididas, conforme segue abaixo: summary(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB)) Error: Bloco Df Sum Sq Mean Sq Bloco 5 0.0091 0.00182 Error: Bloco:Amb Df Sum Sq Mean Sq F value Pr(>F) Amb 2 1.20123 0.60062 34.396 3.293e-05 *** Residuals 10 0.17462 0.01746 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Prof 3 0.077294 0.0257648 8.8629 9.972e-05 *** Amb:Prof 6 0.045089 0.0075148 2.5850 0.03082 * Residuals 45 0.130817 0.0029070 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Até essas análises os resultados estão ok, conforme confronto de resultados da literatura. Porém quando insiro o comando de para análise de médias, os seguintes erros são apresentados: comando: TukeyHSD - TukeyHSD(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB)) Erro: Erro em UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied to an object of class "c('aovlist', 'listof')" comando: LTukey (pacote Laercio) Erro em UseMethod("LTukey") : no applicable method for 'LTukey' applied to an object of class "c('aovlist', 'listof')" Como devo proceder para fazer a análise de médias? Obrigado Atenciosamente, Eng. Agrícola Rafael Cesar Tieppo UNEMAT - Campus Tangará da Serra - MT Departamento de Agronomia Fone:+55-65 3311-4922 site:http://www2.unemat.br/rafaeltieppo "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

Rafael, A função TukeyHSD() não interpreta (ou aplica o teste) para objetos da classe "aovList", apenas para objetos da classe "lm" ou "aov". Ou seja, modelos em que você declara na fórmula o termo "Error()" não podem ser usados pela TukeyHSD(). Duas opções: * Usar a função agricolae::HSD.test(), fornecendo, se for o caso, o QM e GL. Confira na documentação da função. Lembre-se que o desdobramento de parcela dentro de subparcela deve-se usar QM e GL obtidos por aproximação de Satterthwaite. * Usar o pacote ExpDes, não oficial desenvolvido pelo Professor Eric Batista Ferreira da Unifal-MG. Na página ( https://sites.google.com/site/ericbferreira/unifal/downloads-1) dele tem o *.zip para download. Veja na documentação e irá encontrar uma função para Parcela Subdividida. A saída da função se parece muito com a saída das análises do SISVAR. A segunda opção é mais simples, embora a primeira não seja difícil. À disposição. Walmes. Em 22 de março de 2011 12:39, Rafael Tieppo <rafaeltieppo@yahoo.com.br>escreveu:
Bom dia,
Estou com dificuldade para realizar o teste de média para um experimento em parcelas subdivididas, conforme segue abaixo:
summary(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB))
Error: Bloco Df Sum Sq Mean Sq Bloco 5 0.0091 0.00182
Error: Bloco:Amb Df Sum Sq Mean Sq F value Pr(>F) Amb 2 1.20123 0.60062 34.396 3.293e-05 *** Residuals 10 0.17462 0.01746 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within Df Sum Sq Mean Sq F value Pr(>F) Prof 3 0.077294 0.0257648 8.8629 9.972e-05 *** Amb:Prof 6 0.045089 0.0075148 2.5850 0.03082 * Residuals 45 0.130817 0.0029070 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Até essas análises os resultados estão ok, conforme confronto de resultados da literatura.
Porém quando insiro o comando de para análise de médias, os seguintes erros são apresentados:
comando: TukeyHSD - TukeyHSD(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB)) Erro: Erro em UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied to an object of class "c('aovlist', 'listof')"
comando: LTukey (pacote Laercio) Erro em UseMethod("LTukey") : no applicable method for 'LTukey' applied to an object of class "c('aovlist', 'listof')"
Como devo proceder para fazer a análise de médias?
Obrigado
Atenciosamente,
Eng. Agrícola Rafael Cesar Tieppo UNEMAT - Campus Tangará da Serra - MT Departamento de Agronomia Fone:+55-65 3311-4922 site: http://www2.unemat.br/rafaeltieppo
"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- ========================================================================== 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 ==========================================================================

Rafael, Na link abaixo você encontra a implementação da primeira opção do Walmes. Script: http://www.leg.ufpr.br/~eder/CursoRFito/SubDividida/Subdividida.R dados: http://www.leg.ufpr.br/~eder/CursoRFito/SubDividida/ Até mais ----- ------------------------------ Éder David Borges da Silva Eng° Agrônomo ederdbs@gmail.com eder@leg.ufpr.br "Tudo é relativo" Albert Einstein -- View this message in context: http://r-br.2285057.n4.nabble.com/R-br-Teste-media-Subdividida-tp3396775p339... Sent from the R-br mailing list archive at Nabble.com.
participantes (3)
-
Eder David Borges da Silva
-
Rafael Tieppo
-
Walmes Zeviani