-1

我正在使用 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     ###
4

1 回答 1

1

代替

y ~ 1 + Var1 + 1 | group,

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

如果您有多个变量,则将其写为:y ~ 1 + Var1 + Var2 + (1+Var1+Var2 | group)

于 2018-06-01T16:20:37.393 回答