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 <[email protected]> em nome de
[email protected] <[email protected]>
Enviado: domingo, 12 de março de 2017 12:00
Para: [email protected]
Assunto: Digest R-br, volume 75, assunto 8
Enviar submissões para a lista de discussão R-br para
[email protected]
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
[email protected]
Você poderá entrar em contato com a pessoa que gerencia a lista pelo
endereço
[email protected]
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 <[email protected]>
To: Maurício Sangiogo <[email protected]>, a lista Brasileira
oficial de discussão do programa R. <[email protected]>
Subject: Re: [R-br] Desdobramento de um experimento em faixas
Message-ID:
<CAFU=EkYuWZnFRExk8eNhcLSq3u1=OG0Yg=zqunegoxqqc8o...@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
[email protected]
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
****************************************
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo
m�nimo reproduz�vel.
