R 找到固定效果,而不仅仅是拦截效果

R 找到固定效果,而不仅仅是拦截效果,r,lme4,random-effects,R,Lme4,Random Effects,我使用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

我使用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
而不是

y~1+Var1+1 |组

y~1+Var1+(1+Var1 |组)


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

通常的语法是使用括号表示随机效果:
(1 |组)
表示随机截取,或
(Var1 |组)
表示相关随机斜率/截取。因此,也许你想要
y~1+Var1+(Var1 | group)
当我这样做时,它确实会给我固定效果输出,但它会删除随机效果部分中的Var1。请随意编写up@AlexW ; 你能花点钱吗?OP使用随机坡度模型并不罕见;感谢您的回复:对不起,我似乎遗漏了一些东西:这个模型在重复测量/eg生长曲线模型中很常见。(而且,
?lmer
帮助中有几个例子:
lmer(反应~天+(天|主题),睡眠研究)
。这些似乎与op的做法一致。
new_output
(Intercept)    Var1
  0.9716914     ###