Как получить коэффициенты и их доверительные интервалы в моделях со смешанными эффектами?
Вlm
а такжеglm
модели, я использую функцииcoef
а такжеconfint
для достижения цели:
m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous
coef(m)
confint(m)
Теперь я добавил случайный эффект к модели - использовал модели смешанных эффектов, используяlmer
функция из пакета lme4. Но тогда функцииcoef
а такжеconfint
не работай больше на меня!
> 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
Я пытался гуглить и использовать документы, но безрезультатно. Пожалуйста, укажите мне в правильном направлении.
РЕДАКТИРОВАТЬ: я также думал, подходит ли этот вопрос большеhttps://stats.stackexchange.com/ но я считаю, что это скорее технический, чем статистический, поэтому я пришел к выводу, что это подходит лучше всего (SO) ... что вы думаете?