Jak uzyskać współczynniki i ich przedziały ufności w modelach efektów mieszanych?
Wlm
iglm
modele, używam funkcjicoef
iconfint
aby osiągnąć cel:
m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous
coef(m)
confint(m)
Teraz dodałem efekt losowy do modelu - użyłem mieszanych modeli efektów używająclmer
funkcja z pakietu lme4. Ale potem działacoef
iconfint
nie pracuj już dla mnie!
> 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
Próbowałem google i używać dokumentów, ale bez rezultatu. Proszę wskazać mi właściwy kierunek.
EDYCJA: Myślałem również, czy to pytanie bardziej pasujehttps://stats.stackexchange.com/ ale uważam to za bardziej techniczne niż statystyczne, więc doszedłem do wniosku, że pasuje tutaj najlepiej (SO) ... co o tym sądzisz?