Pessoal, estou tentando fazer um desdobramento de interação em uma análise conjunta, que consiste em 2 experimentos em DBC com 6 tratamentos e 5 blocos cada, mas estou tendo problemas apenas em um dos desdobramentos.
Ao desdobrar tratamentos dentro de experimentos da seguinte forma:
desd1<-aov(formula = AACPD~EXPERIMENTO + EXPERIMENTO/BLOCO + EXPERIMENTO/TRATAMENTO, data = dados)
anova(desd1)
summary.aov(object = desd1, split = list("EXPERIMENTO:TRATAMENTO"=list("dag"=c(1,3,5,7,9), "vitorinha"=c(2,4,6,8,10))))
não ocorre nenhum problema.
Porém, ao tentar fazer o desdobramento de experimentos dentro de tratamentos da seguinte maneira:
desd2<-aov(formula = AACPD~TRATAMENTO + TRATAMENTO/EXPERIMENTO + EXPERIMENTO/BLOCO, data = dados)
anova(desd2)
primeiramente o quadro de Anava fica assim:
Response: AACPD
FV Df Sum Sq Mean Sq F value Pr(>F)
TRATAMENTO 5 170,5 34,1 0,4086 0,83999
EXPERIMENTO 1 5094,5 5094,5 61,0295 1,423e-09 ***
TRATAMENTO:EXPERIMENTO 5 1249,2 249,8 2,9930 0,02184 *
EXPERIMENTO:BLOCO 8 995,2 124,4 1,4903 0,19143
Residuals 40 3339,0 83,5
E o seguinte desdobramento
summary.aov(object = desd2,split = list("TRATAMENTO:EXPERIMENTO"= list("N:S.30:07"=1, "nit+sulf."=2, "nitrato"=3, "sem.N"=4, "sulfato"=5, "ureia"=6)))
resulta em:
FV Df Sum Sq Mean Sq F value Pr(>F)
TRATAMENTO 5 171 34 0,409 0,839988
EXPERIMENTO 1 5095 5095 61,030 1,42e-09 ***
TRATAMENTO:EXPERIMENTO 5 1249 250 2,993 0,021835 *
TRATAMENTO:EXPERIMENTO: N:S 30:07 1 5 5 0,054 0,816936
TRATAMENTO:EXPERIMENTO: nit+sulf 1 99 99 1,184 0,282958
TRATAMENTO:EXPERIMENTO: nitrato 1 1116 1116 13,364 0,000738 ***
TRATAMENTO:EXPERIMENTO: sem N 1 19 19 0,226 0,637284
TRATAMENTO:EXPERIMENTO: sulfato 1 11 11 0,137 0,713210
TRATAMENTO:EXPERIMENTO: ureia 1
EXPERIMENTO:BLOCO 8 995 124 1,490 0,191434
Residuals 40 3339 83
o que não está correto!
Alguém sabe onde errei nas funções acima?