Como obter coeficientes e seus intervalos de confiança em modelos de efeitos mistos?
Emlm
eglm
modelos, eu uso funçõescoef
econfint
para alcançar o objetivo:
m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous
coef(m)
confint(m)
Agora eu adicionei efeito aleatório ao modelo - usei modelos de efeitos mistos usandolmer
função do pacote lme4. Mas então, funçõescoef
econfint
não trabalhe mais para mim!
> mix1 = lmer(resp ~ 0 + var1 + var1:var2 + (1|var3))
# var1, var3 categorical, var2 continuous
> coef(mix1)
Error in coef(mix1) : unable to align random and fixed effects
> confint(mix1)
Error: $ operator not defined for this S4 class
Eu tentei google e usar documentos, mas sem resultado. Por favor me aponte na direção correta.
EDIT: Eu também estava pensando se esta questão se encaixa mais parahttps://stats.stackexchange.com/ mas eu considero isso mais técnico do que estatístico, então eu concluí que se encaixa melhor aqui (SO) ... o que você acha?