查找固定效应不只是为了拦截

Question

我正在使用 lmer 函数来收集随机效应和固定效应。当运行我的函数时，我能够按组获得所有系数的随机效应。但是，当我搜索固定效果时，我只得到它的截距。如何获得所有系数的固定效应？

这是我的例子：

form_lmer <- "y  ~   1 + Var1 + 1 | group"
fm3<- lmer(as.formula(form_lmer), df, REML = FALSE, verbose = TRUE)
final<-ranef(fm3)
head(final)

$group
      (Intercept)                    Var1
500  0.0429171100           -0.0114512020
501 -0.0483822296            0.0793030797
502  0.0368572644           -0.0060213102
504  0.0390749512           -0.0075918832
505  0.0608323336           -0.0294764946
506  0.0030031361            0.0278525937

final2<-fixef(fm3)
final2
(Intercept) 
  0.9716914

我希望固定效果的最终输出是这样的：

new_output
(Intercept)    Var1
  0.9716914     ###

Answer 1

而不是

y ~ 1 + Var1 + 1 | group、

写

y ~ 1 + Var1 + (1+Var1 | group)。

如果你有多个变量，那么写成：y ~ 1 + Var1 + Var2 + (1+Var1+Var2 | group)

查找固定效应不只是为了拦截

Find Fixed Effects not just for intercept

r

lme4

random-effects