Pessoal, poque o resultado dos testes HSD.test e lsmeans diferem neste exemplo?

O resultado fornecido pelo lsmeas está em acordo com o material de onde tirei esse exemplo e que havia sido rodado no SAS. Já o HSD.test difere tanto no valor das médias estimadas como na significancia.

install.packages(c("agricolae",lsmeans"))
library(agricolae)
library(lsmeans)

COVARIANCIA<-read.csv("https://www.dropbox.com/s/0uvwqnnvsdsvphc/Covariancia.csv?raw=1",head=TRUE)

modeloCOVAR<-lm(GANHO~PESOINIC+DIETA,data=COVARIANCIA)
anova(modeloCOVAR)
HSD.test(modeloCOVAR,"DIETA",console=TRUE)
lsmeans(modeloCOVAR,pairwise~DIETA,adjust="tukey")

--
=========================================
Fernando Souza
Zootecnista, DSc. Produção e Alimentação Animal
Celular: (31)99796-8781 (Vivo)
E-mail:nandodesouza@gmail.com
Lattes: http://lattes.cnpq.br/6519538815038307
Blog: https://producaoanimalcomr.wordpress.com/
==========================================