¿Cómo obtener los coeficientes y sus intervalos de confianza en modelos de efectos mixtos?
Enlm
yglm
modelos, uso funcionescoef
yconfint
para lograr el objetivo:
m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous
coef(m)
confint(m)
Ahora agregué efectos aleatorios al modelo - usé modelos de efectos mixtos usandolmer
Función del paquete lme4. Pero entonces, las funciones.coef
yconfint
¡No trabajes más para mí!
> 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
Intenté buscar en Google y usar documentos pero sin resultado. Por favor, apúntame en la dirección correcta.
EDITAR: También estaba pensando si esta pregunta se ajusta más ahttps://stats.stackexchange.com/ pero lo considero más técnico que estadístico, así que llegué a la conclusión de que encaja mejor aquí (SO) ... ¿qué piensas?