boa noite,
estou precisando verificar o efeito linear e quadrático para um modelo usando a lsmeans. Imagina neste exemplo DIETA sendo quantitativo 10, 20 e 30. Como eu poderia programar tal análise usando a lsmeans?
Grato
Em terça-feira, 31 de janeiro de 2017 13:17:52 BRST, Walmes Zeviani via R-br <r-br@listas.c3sl.ufpr.br> escreveu:


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.


_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br

https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.