[R-br] Experimento fatorial. Desdobramento

tiago souza marçal tiagosouzamarcal em hotmail.com
Segunda Novembro 19 17:22:10 BRST 2012


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 em gmail.com
To: r-br em listas.c3sl.ufpr.br; camposagro em 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 em 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_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.



_______________________________________________
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/20121119/9a9cb5fe/attachment.html>


Mais detalhes sobre a lista de discussão R-br