[R-br] Enc: desdobramento de interação numa análise conjunta
helcio duarte pereira
hhelciopassos em yahoo.com.br
Quinta Janeiro 29 21:11:17 BRST 2015
Em Terça-feira, 27 de Janeiro de 2015 1:31, helcio duarte pereira <hhelciopassos em yahoo.com.br> escreveu:
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?
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150129/48e9d309/attachment.html>
Mais detalhes sobre a lista de discussão R-br