[R-br] Teste média Subdividida

[Walmes Zeviani]

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 em 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 em 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 em ufpr.br                      twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes    linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110322/9f03c507/attachment.html>