[R-br] Experimento fatorial. Desdobramento

[FHRB Toledo]

Odirley,

Os comando que o prof Walmes postou na wiki do LEG não são exatamente
transferidos aos seus dados... Observe:

Veja os resultados da sua anova, que você chama de "geral", temos 1 grau de
liberdade para o fator A e também um grau de liberdade para o fator K
dentro de A.

Em seguida você especifica 4 termos (coeficientes do modelo) para cada
fator K dentro de A!

Veja por exemplo qual a saída da função coef(m1). Nela será mostrado quais
coeficientes correspondem a quais termos do modelo (pelo nome do
coeficiente em cima). O que você deve por como os número no split da função
summary são justamente essas posições.

Com certeza essa não é a melhor explicação que você vai ver por aqui,
quanto ao seu thread, mas...

att,
FH

2012/11/16 Odirley Campos <camposagro em yahoo.com.br>

> #Pessoal! Boa noite !
> #Estou tentando sem sucesso fazer o desdobramento de interações com R.
> Cheguei a baixar alguns arquivos da net no endereço "
> http://www.leg.ufpr.br/doku.php/ridiculas#desdobramento_de_interacao_em_experimento_fatorial"
> . Porêm tive o mesmo problema que tive com meus dados.
> #Usei os mesmos comandos retirados da net para os dados tirados de lá. A
> anova "geral" saiu normalmente:
>
> >rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T) #após
> armezenar arquivo no meu PC#
> >rend
> >m1 <- aov(rg~bloc+A/K, data=rend)
> >summary(m1) #Anova geral saiu norma sem problemas#
>             Df Sum Sq Mean Sq F value  Pr(>F)
> bloc         1   51.3    51.3   3.585  0.0624 .
> A            1  495.6   495.6  34.633 1.2e-07 ***
> A:K          1 2213.7  2213.7 154.712 < 2e-16 ***
> Residuals   71 1015.9    14.3
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> >summary(m1,split=list("A:K"=list("A-37.5"=c(1,4,7,10),"A-50.0"=c(2,5,8,11),"A-62.5"=c(3,6,9,12))))
> #Mas o desdobramento não mostra Quadrado médio nem mesmo, F-value ou
> mesmo o P-value. O que eu posso ter feito de errado??
>               Df Sum Sq Mean Sq F value  Pr(>F)
> bloc           1   51.3    51.3   3.585  0.0624 .
> A              1  495.6   495.6  34.633 1.2e-07 ***
> A:K            1 2213.7  2213.7 154.712 < 2e-16 ***
>   A:K: A-37.5  4
>   A:K: A-50.0  4
>   A:K: A-62.5  4
> Residuals     71 1015.9    14.3
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> Odirley Campos
> Engenheiro Agrônomo 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/20121116/0da02675/attachment.html>