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.
DesKinA<-summary(m1, split=list("A:K"=KinA))
DesKinA
DesAinK<-summary(m2, split=list("A:K"=AinK))
DesAinK
Outra forma seria pelo método dos efeitos aninhados
#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