[R-br] Diferença resultados agricolae::HSD.test e lsmeas

[Walmes Zeviani]

Isso acontece porque a lsmeans() usa médias ajustadas para o efeito de peso
inicial enquanto que HSD.test() usa as médias amostrais que não são
ajustadas. Com a função doBy::LSmeans() pode-se verificar que as médias
ajustadas são para um animal hipotético de peso igual a média dos animais
do experimento. Quando fizer análise de covariância, o mais correto é usar
médias ajustadas. As médias amostrais só serão iguais a média ajustada
quando os efeitos dos fatores a serem marginalizados são ortogonais. Esse é
o caso do experimento em DBC balanceado onde as médias ajustadas e as
médias amostrais são iguais pois o efeito de bloco é ortogonal ao de
tratamentos.

library(multcomp)
library(doBy)

lsm <- LSmeans(modeloCOVAR, effect = "DIETA")
lsm$K

# Média do peso inicial.
mean(COVARIANCIA$PESOINIC)

# Média do ganho por dieta.
aggregate(GANHO ~ DIETA, data = COVARIANCIA, FUN = mean)

À disposição.
Walmes.

​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170131/1991e3f0/attachment.html>