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 !
#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