Uma outra maneira de fazer isto e bem mais fácil é utilizar o pacote ExpDes, ele já faz o desdobramento automaticamente quando dá interação e bem fácil de utilizar.
Atte..
#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 ‘ ’ 1Odirley Campos
Engenheiro Agrônomo UFV/MG
_______________________________________________
R-br mailing list
R-br@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.