# 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.1897solo:nc: EU_0-10 0.5 0.5 0.079 0.2213solo:nc: EU_40-60 2.1 2.1 0.343 0.4406solo:nc: SE_0-10 81.0 81.0 13.347 0.0004solo:nc: SE_40-60 0.0 0.0 0.001 0.0236solo: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/MGDoutorando em Solos e Nutrição de Plantas UFV/MG