<div> </div>
<div><span style="color: #0000ff; font-size: medium;">Tentei rodar e não consegui, logo de saida não reconhece o comando strip.plot. Por acaso precisa de alguma bibliotéca?</span></div>
<div>
<div style="font-size: 12.8px;"> </div>
<br />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.<br /><br /><br />Acredito eu que está certo como fiz.<br /><br /><br /><br />Se a alguém interessar, ou se disponibilizar de conferir, mando os dados e comandos:<br /><br /><br /># Análise de um experimento em faixas com interação significativa:<br /><br />dados<- structure(list(local = c(1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,<br />1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1,<br />2, 3, 4), trat = c(0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60,<br />0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30,<br />30, 30, 60, 60, 60, 60), bloco = c(1, 1, 1, 1, 1, 1, 1, 1, 1,<br />1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3,<br />3, 3, 3, 3, 3, 3), yield = c(150, 163, 164, 171, 188, 174, 182,<br />182, 220, 170, 170, 182, 147, 141, 133, 140, 157, 166, 144, 143,<br />173, 165, 152, 170, 146, 140, 149, 150, 172, 158, 142, 160, 204,<br />163, 150, 175)), .Names = c("local", "trat", "bloco", "yield"<br />), row.names = c(NA, 36L), class = "data.frame")<br /><br />model<-with(dados,strip.plot(bloco, trat, local, yield))<br /><br /># Se interação não significativa:<br /><br />comparison1<-with(dados,HSD.test(yield,trat,model$gl.a,model$Ea));comparison1<br /><br />comparison2<-with(dados,HSD.test(yield,local,model$gl.b,model$Eb)); comparison2<br /><br /># Se tiver interação significativa:<br /><br />QMea<- model$Ea ;Dfea<-model$gl.a<br />QMeb<- model$Eb; Dfeb<-model$gl.b<br />QMec<- model$Ec ; Dfec<-model$gl.c<br /><br />a <-3; b <-4<br /><br />#Erro medio a e c<br />QMeac <-(QMea+(b-1)*QMec)/b<br /><br />##Graus de liberdade a e c<br />Glac<-((QMea+(b-1)*QMec)^2) / ( ((QMea^2)/Dfea) + ((((b-1)*QMec)^2)/Dfec) )<br /><br />#Erro medio b e c<br />QMebc<-(QMeb+(a-1)*QMec)/a<br /><br />##Graus de liberdade a e c<br />Glbc<-((QMeb+(a-1)*QMec)^2) / ( ((QMeb^2)/Dfeb) + ((((a-1)*QMec)^2)/Dfec) )<br /><br />#---------------------------------------------------------------------------------<br /># Comparação de médias por Tukey:<br /><br />## Tratamentos dentro de locais:<br /><br />L1<-subset(dados,dados$local=="1")<br />t1<-HSD.test(L1$yield,L1$trat,Glac,QMeac)<br />t1<br /><br />L2<-subset(dados,dados$local=="2")<br />t2<-HSD.test(L2$yield,L2$trat,Glac,QMeac)<br />t2<br /><br />L3<-subset(dados,dados$local=="3")<br />t3<-HSD.test(L3$yield,L3$trat,Glac,QMeac)<br />t3<br /><br />L4<-subset(dados,dados$local=="4")<br />t4<-HSD.test(L4$yield,L4$trat,Glac,QMeac)<br />t4<br /><br />## Locais dentro de tratamentos:<br /><br />A1<-subset(dados,dados$trat=="0")<br />ta1<-HSD.test(A1$yield,A1$local,Glbc,QMebc)<br />ta1<br /><br />A2<-subset(dados,dados$trat=="30")<br />ta2<-HSD.test(A2$yield,A2$local,Glbc,QMebc)<br />ta2<br /><br />A3<-subset(dados,dados$trat=="60")<br />ta3<-HSD.test(A3$yield,A3$local,Glbc,QMebc)<br />ta3<br /><br /># FIM ?<br /><br /><br />Abraços!<br /><br /><br /><br /><br />..................<br />Me. Eng. Agr. Maurício Sangiogo<br />TEL:(55) 96822030<br /><br /><br /><br />________________________________<br /><br /></div>