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



________________________________