Bom dia pessoal, Mais uma vez venho recorrer a vocês em relação a uma dúvida... Estou trabalhando com um experimento em DBC em faixas, porém tenho encontrado poucos exemplos em relação a comparação de médias no caso de interação. Encontrei um exemplo em uma aula na internet e um exemplo com os comandos no pacote agricolae (http://127.0.0.1:25234/library/agricolae/html/strip.plot.html), porém mesmo usando os mesmos valores a comparação de médias não está batendo entre o exemplo da aula e quando rodo usando os comandos do exemplo do agricolae. Alguém ai poderia conferir se esses comandos que envio logo a baixo estão corretos? Principalmente na parte referente a comparação de médias. (Seria um experimento em blocos completos no esquema em faixas.) Desde já muito obrigado! Abraços # Yield library(agricolae) data(huasahuasi) YIELD<-huasahuasi$YIELD market <- YIELD$y1da + YIELD$y2da non_market <- YIELD$y3da yield <- market + non_market model<-with(YIELD,strip.plot(block, clon, trt, yield)) comparison<-with(YIELD,LSD.test(yield,clon,model$gl.a,model$Ea)) comparison<-with(YIELD,LSD.test(yield,trt,model$gl.b,model$Eb)) # simple effects A<-model$data a<-nlevels(A$clon) b<-nlevels(A$trt) r<-nlevels(A$block) Ea<-model$Ea; Eb<-model$Ec; Ec<-model$Ec; gla<-model$gl.a; glb<-model$gl.b; glc<-model$gl.c; B <-tapply.stat(A[,4],A[,2:3],mean) std<-tapply.stat(A[,4],A[,2:3],function(x) sd(x)/sqrt(length(x))) B<-data.frame(B[,1:2],yield=B[,3],std=std[,3]) cmab<-(b-1)*Ec + Ea cmba<-(a-1)*Ec + Eb # order.group # Tukey ta <- qtukey(0.95,a,gla) tb <- qtukey(0.95,b,glb) tc <- qtukey(0.95,a*b,glc) tab<- ((b-1)*Ec*tc + Ea*ta)/cmab tba<- ((a-1)*Ec*tc + Eb*tb)/cmba # Comparison of clones by treatment groups<-by(B,B[,2], function(x) order.group(x$clon,x$yield,N=b*r, cmab,tab, std.err=x$std,parameter=0.5,console=FALSE)) groups # Comparison of treatments by clon groups<-by(B,B[,1], function(x) order.group(x$trt,x$yield,N=a*r,cmba,tba, std.err=x$std,parameter=0.5,console=FALSE)) groups .................. Me. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030