Experimento fatorial. Desdobramento
#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_expe..." . 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
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@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_expe..." . 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@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.
Vc pode enviar um CRM que te ajudo a fazer o fatorial, tenho um script ponto mas como usei muita funções anonima acho que vai ser dificil de vc entender. Alisson Lucrécio da Costa ________________________________ From: FHRB Toledo <fernandohtoledo@gmail.com> To: r-br@listas.c3sl.ufpr.br; Odirley Campos <camposagro@yahoo.com.br> Sent: Friday, November 16, 2012 9:52 PM Subject: Re: [R-br] Experimento fatorial. Desdobramento 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@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_expe..." . 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@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.
_______________________________________________ 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.
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.. Em 16 de novembro de 2012 21:41, Odirley Campos <camposagro@yahoo.com.br>escreveu:
#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_expe..." . 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@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.
Só completando, o pacote ExpDes é a maneira mais fácil, mas caso você queira fazer da forma que esta fazendo é só utilizar os seguintes comandos: rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T)attach(rend)A<-factor(A)K<-factor(K)bloc<-factor(bloc)m1 <- aov(rg~bloc+A/K)KinA<-sapply(paste("A",levels(A),sep=""),simplify=F,grep,x=names(coef(m1)[m1$assign==3]))DesKinA<-summary(m1, split=list("A:K"=KinA))DesKinA m2 <- aov(rg~bloc+K/A)AinK<-sapply(paste("K",levels(K),sep=""),simplify=F,grep,x=names(coef(m2)[m2$assign==3]))DesAinK<-summary(m2, split=list("A:K"=AinK))DesAinK Outra forma seria pelo método dos efeitos aninhados rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T)attach(rend)A<-factor(A)K<-factor(K)bloc<-factor(bloc)m1 <- aov(rg~bloc+A/K)effects(m1) #A partir da matriz de efeitos você seleciona os efeitos de interesse.DesKinA<-summary(m1, split=list("A:K"=list(A37.5=c(1,4,7,10),A50=c(2,5,8,11),A62.5=c(3,6,9,12))))DesKinAm2 <- aov(rg~bloc+K/A)effects(m2)DesAinK<-summary(m2, split=list("K:A"=list(K0=c(1,6),K30=c(2,7),K60=c(3,8),K120=c(4,9),K180=c(5,10))))DesAinK Espero que tenhamos sanado sua duvida. Att. Tiago. Date: Sat, 17 Nov 2012 09:21:23 -0200 From: jhguilhen@gmail.com To: r-br@listas.c3sl.ufpr.br; camposagro@yahoo.com.br Subject: Re: [R-br] Experimento fatorial. Desdobramento 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.. Em 16 de novembro de 2012 21:41, Odirley Campos <camposagro@yahoo.com.br> escreveu: #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_expe..." . 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@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. _______________________________________________ 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.
participantes (5)
-
Alisson Lucrecio -
FHRB Toledo -
Odirley Campos -
tiago souza marçal -
Ze Henrique