Ok Mesmo,assim vou dar uma olhada. obrigado pela sugestão Alisson. Att. Tiago. Date: Mon, 19 Nov 2012 18:15:39 -0800 From: alissonluc@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Experimento fatorial. Desdobramento desculpa, enviei errado. Alisson Lucrécio da Costa From: tiago souza marçal <tiagosouzamarcal@hotmail.com> To: r-br@listas.c3sl.ufpr.br Sent: Monday, November 19, 2012 9:09 PM Subject: Re: [R-br] Experimento fatorial. Desdobramento Estes são os dados. Qualquer duvida eu estou a disposição. K A bloc rg ts nvi nv 0 37.5 1 14.55 136 22 56 30 37.5 1 21.51 159 2 62 60 37.5 1 24.62 156 0 66 120 37.5 1 21.88 171 2 68 180 37.5 1 28.11 190 0 82 0 50 1 17.16 140 20 63 30 50 1 24.44 193 6 86 60 50 1 28.5 200 6 94 120 50 1 32.94 208 7 86 180 50 1 28.76 237 10 97 0 62.5 1 14.13 132 23 58 30 62.5 1 19.71 178 10 80 60 62.5 1 31.05 231 5 99 120 62.5 1 38.35 271 0 110 180 62.5 1 36.97 262 1 106 0 37.5 2 15.72 120 22 48 30 37.5 2 19.72 146 0 64 60 37.5 2 24.29 159 1 62 120 37.5 2 25.39 152 0 60 180 37.5 2 28.45 159 0 62 0 50 2 15.44 124 12 52 30 50 2 28.12 239 9 103 60 50 2 31.49 248 4 103 120 50 2 29.95 190 1 87 180 50 2 30.06 240 4 101 0 62.5 2 15.42 140 14 64 30 62.5 2 23.68 167 7 77 60 62.5 2 33.27 237 3 96 120 62.5 2 36.63 240 1 94 180 62.5 2 41.04 248 0 105 0 37.5 3 12.77 106 4 44 30 37.5 3 20.45 157 3 63 60 37.5 3 24.35 204 1 84 120 37.5 3 27.15 169 0 66 180 37.5 3 24.08 147 1 58 0 50 3 16.01 133 22 56 30 50 3 24.57 199 7 81 60 50 3 28.85 211 3 92 120 50 3 26.23 223 7 91 180 50 3 26.36 261 5 100 0 62.5 3 16.08 128 21 50 30 62.5 3 20.99 164 9 71 60 62.5 3 30.45 216 3 94 120 62.5 3 36.83 215 2 89 180 62.5 3 34.6 228 2 89 0 37.5 4 14.26 125 11 60 30 37.5 4 23.71 185 7 75 60 37.5 4 22.76 159 1 67 120 37.5 4 22.46 176 2 76 180 37.5 4 22.97 154 1 63 0 50 4 15.54 133 27 51 30 50 4 18.48 149 4 60 60 50 4 25.54 170 1 67 120 50 4 30.29 215 2 90 180 50 4 30.64 186 4 80 0 62.5 4 17.26 147 18 62 30 62.5 4 20.51 154 10 67 60 62.5 4 24.47 179 2 75 120 62.5 4 34.52 229 0 90 180 62.5 4 37.62 216 1 87 0 37.5 5 10.3 92 5 36 30 37.5 5 16.28 133 5 57 60 37.5 5 23.61 139 2 53 120 37.5 5 29.66 166 0 69 180 37.5 5 23.34 128 1 50 0 50 5 14.41 131 26 57 30 50 5 17.24 150 1 61 60 50 5 29.08 218 4 92 120 50 5 29.52 236 1 96 180 50 5 27.98 261 4 110 0 62.5 5 17.56 150 17 69 30 62.5 5 20.05 156 13 63 60 62.5 5 29.9 215 3 95 120 62.5 5 33.21 128 0 55 180 62.5 5 35.5 208 2 90 Att. Tiago. From: sznelwar@uol.com.br To: tiagosouzamarcal@hotmail.com Subject: Fw: [R-br] Experimento fatorial. Desdobramento Date: Mon, 19 Nov 2012 20:07:34 -0200 Poderia enviar este arquivo UFPR1.txt, eu não estou conseguindo entrar na página do LEG. 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)))) DesKinA m2 <- 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.