Obrigado Walmes, Obrigado pela indicação. Mas acabei optando por realizar o cálculo dos quadrados médios combinados e graus de liberdade combinados manualmente mesmo, era isso que eu estava na dúvida... e de fato deu resultado diferente do método que estava de exemplo no pacote agricolae. Acredito eu que está certo como fiz. Se a alguém interessar, ou se disponibilizar de conferir, mando os dados e comandos: # Análise de um experimento em faixas com interação significativa: dados<- structure(list(local = c(1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4), trat = c(0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60), bloco = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3), yield = c(150, 163, 164, 171, 188, 174, 182, 182, 220, 170, 170, 182, 147, 141, 133, 140, 157, 166, 144, 143, 173, 165, 152, 170, 146, 140, 149, 150, 172, 158, 142, 160, 204, 163, 150, 175)), .Names = c("local", "trat", "bloco", "yield" ), row.names = c(NA, 36L), class = "data.frame") model<-with(dados,strip.plot(bloco, trat, local, yield)) # Se interação não significativa: comparison1<-with(dados,HSD.test(yield,trat,model$gl.a,model$Ea));comparison1 comparison2<-with(dados,HSD.test(yield,local,model$gl.b,model$Eb)); comparison2 # Se tiver interação significativa: QMea<- model$Ea ;Dfea<-model$gl.a QMeb<- model$Eb; Dfeb<-model$gl.b QMec<- model$Ec ; Dfec<-model$gl.c a <-3; b <-4 #Erro medio a e c QMeac <-(QMea+(b-1)*QMec)/b ##Graus de liberdade a e c Glac<-((QMea+(b-1)*QMec)^2) / ( ((QMea^2)/Dfea) + ((((b-1)*QMec)^2)/Dfec) ) #Erro medio b e c QMebc<-(QMeb+(a-1)*QMec)/a ##Graus de liberdade a e c Glbc<-((QMeb+(a-1)*QMec)^2) / ( ((QMeb^2)/Dfeb) + ((((a-1)*QMec)^2)/Dfec) ) #--------------------------------------------------------------------------------- # Comparação de médias por Tukey: ## Tratamentos dentro de locais: L1<-subset(dados,dados$local=="1") t1<-HSD.test(L1$yield,L1$trat,Glac,QMeac) t1 L2<-subset(dados,dados$local=="2") t2<-HSD.test(L2$yield,L2$trat,Glac,QMeac) t2 L3<-subset(dados,dados$local=="3") t3<-HSD.test(L3$yield,L3$trat,Glac,QMeac) t3 L4<-subset(dados,dados$local=="4") t4<-HSD.test(L4$yield,L4$trat,Glac,QMeac) t4 ## Locais dentro de tratamentos: A1<-subset(dados,dados$trat=="0") ta1<-HSD.test(A1$yield,A1$local,Glbc,QMebc) ta1 A2<-subset(dados,dados$trat=="30") ta2<-HSD.test(A2$yield,A2$local,Glbc,QMebc) ta2 A3<-subset(dados,dados$trat=="60") ta3<-HSD.test(A3$yield,A3$local,Glbc,QMebc) ta3 # FIM ? Abraços! .................. Me. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030 ________________________________ De: R-br <r-br-bounces@listas.c3sl.ufpr.br> em nome de r-br-request@listas.c3sl.ufpr.br <r-br-request@listas.c3sl.ufpr.br> Enviado: domingo, 12 de março de 2017 12:00 Para: r-br@listas.c3sl.ufpr.br Assunto: Digest R-br, volume 75, assunto 8 Enviar submissões para a lista de discussão R-br para r-br@listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request@listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner@listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Desdobramento de um experimento em faixas (Walmes Zeviani) ---------------------------------------------------------------------- Message: 1 Date: Sat, 11 Mar 2017 17:44:03 -0300 From: Walmes Zeviani <walmeszeviani@gmail.com> To: Maurício Sangiogo <ms_sangiogo@hotmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br@listas.c3sl.ufpr.br> Subject: Re: [R-br] Desdobramento de um experimento em faixas Message-ID: <CAFU=EkYuWZnFRExk8eNhcLSq3u1=OG0Yg=ZQUNeGoxQQc8obzQ@mail.gmail.com> Content-Type: text/plain; charset="utf-8" Verifique se essa postagem do histórico da lista serve: http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html. Caso R-br - [R-br] Anova em faixas<http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html> r-br.2285057.n4.nabble.com [R-br] Anova em faixas. Srs, Preciso de uma ajuda. Preciso fazer uma ANOVA de um experimento com três fatores, porém o problema está no delineamento que foi em faixas. Tenho o quadro da análise... não funcione, envie no corpo da mensagem os resultados que você tem como sendo os certos e o resultado que você está conseguindo. Não esqueça de disponibilizar acesso aos dados usados no corpo da mensagem (use dput()). À disposição. Walmes. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170311/02ec2ff1/attachment-0001.html> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). ------------------------------ Fim da Digest R-br, volume 75, assunto 8 ****************************************