[R-br] Anova, fatorial triplo, interacao dupla

Quarta Fevereiro 4 14:06:18 BRST 2015

Odirley,

Sua informação sobre as somas quadráticas "não baterem" em relação ao Excel
não parecem corretas, posto que os valores são iguais em ambas as tabelas
que você postou.

Em relação aos valores F, não é possível discutir porque você não postou a
linha de soma quadrática dos resíduos do Excel (e nem o nº de graus de
liberdade) .

2015-02-03 19:07 GMT-02:00 Odirley Campos <camposagro em yahoo.com.br>:

> # Boa noite! Novamente venho recorrer a lista para tirar duvidas sobre
> desdobramento na ANOVA.
> # Tenho resultados de materia seca de eucalipto de um experimento em
> arranjo fatorial triplo (solo; nc; np) com # UE ditribuidas em DBC.
> # Dados estão no: http://www.datafilehost.com/d/32efcd4b
> # Gostaria de corrigir os valores do Fcal no desdobramento das interações
> duplas solo:nc.
> # Na anova geral (avi) o GL e SQ do residuo são obtidos normalmente:
>
>
> avi<-aov(mspa~bloco + trat, data=a1)
> anova(avi)
>
>
> Response: mspa Df Sum Sq Mean Sq F value Pr(>F) bloco 2 2 0.993 0.1636
> 0.8493 trat 59 4963 84.119 13.8640 <2e-16 *** Residuals 118 716 6.067
>
>
> #refiz no excel e os resultados batem
>
>
>
> av1<-aov(mspa~bloco + solo*nc*np, data=a1) # com mais detalhes
> anova(av1)
>
>
> Response: mspa
>             Df  Sum Sq Mean Sq F value    Pr(>F)
> bloco        2    1.99    0.99  0.1636   0.84927
> solo         5 1987.74  397.55 65.5217 < 2.2e-16 ***
> nc           1    0.75    0.75  0.1232   0.72621
> np           4 2350.35  587.59 96.8431 < 2.2e-16 ***
> solo:nc      5   95.86   19.17  3.1599   0.01031 *
> solo:np     20  395.63   19.78  3.2603 3.261e-05 ***
> nc:np        4   10.54    2.64  0.4345   0.78348
> solo:nc:np  20  122.12    6.11  1.0064   0.46028
> Residuals  118  715.95    6.07
>
>
> # Após conferir no excel vi que: df, SS, MS, F e Pr estão corretos para av1
> # quando fui para o desdobramento conforme segue, obtive os resultados no R:
>
> av2 <- aov(mspa~bloco+solo/nc) #nc dentro solo
> anova(av2)
>
> names(coef(av2))
> length(names(coef(av2)))
>
>   # buscando as posiçoes por expressões regulares
>
> grep("ES", names(coef(av2))[9:14])
> grep("EU_0-10", names(coef(av2))[9:14])
> grep("EU_40-60", names(coef(av2))[9:14])
> grep("SE_0-10", names(coef(av2))[9:14])
> grep("SE_40-60", names(coef(av2))[9:14])
> grep("TM", names(coef(av2))[9:14])
>
>
> summary(av2, split=list("solo:nc"=list(
>   "ES"= c(1),
>   "EU_0-10"=c(2),
>   "EU_40-60"=c(3),
>   "SE_0-10"=c(4),
>   "SE_40-60"= c(5),
>   "TM"= c(6)
>   )))
>
>
> # saida da av2 abaixo:
>
>                      Df Sum Sq Mean Sq F value   Pr(>F)
> bloco                 2      2     1.0   0.046   0.9552
> solo                  5   1988   397.5  18.359 1.73e-14 ***
> solo:nc               6     97    16.1   0.744   0.6153
>   solo:nc: ES         1     11    10.6   0.488   0.4859
>   solo:nc: EU_0-10    1      0     0.5   0.022   0.8817
>   solo:nc: EU_40-60   1      2     2.1   0.096   0.7570
>   solo:nc: SE_0-10    1     81    81.0   3.741   0.0548 .
>   solo:nc: SE_40-60   1      0     0.0   0.000   0.9875
>   solo:nc: TM         1      2     2.5   0.114   0.7362
> Residuals           166   3595    21.7
>
>
> # vi que na parte que me interessa em destaque a SQ e QM estão corretas,
> mas o Fcal e P-valor estão errados pois o GL, SQ e QM do residuo não estão
> corretos.   Conforme conferi pelo excel:
> # saida do excel
> FV                            Sum Sq  Mean Sq         F value
> Pr(>F)
> solo:nc: ES                 10.6      10.6                1.740
> 0.1897
> solo:nc: EU_0-10          0.5       0.5                  0.079
> 0.2213
> solo:nc: EU_40-60        2.1       2.1                  0.343
> 0.4406
> solo:nc: SE_0-10          81.0     81.0                13.347
> 0.0004
> solo:nc: SE_40-60        0.0       0.0                  0.001
> 0.0236
> solo:nc: TM                 2.5        2.5                 0.406
> 0.4748
>
> # Existe alguma maneira de informar o df e Sum Sq corretos (que estão na
> avi ou av1) no procedimento da av2
> # de modo que o F e Pr sejam corrigidos na saida do comando?
> # Ou como eu poderia extrair a coluna do Sum Sq ou Mean Sq na forma de
> vector ou data.frame para fazer as correções manualmente no R?
>
>
> Odirley R. Campos
> Engenheiro Agrônomo UFV/MG
> Doutorando em Solos e Nutrição de Plantas UFV/MG
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150204/d1e7ae8b/attachment.html>