提取 Dataframe 中的固定效应和随机效应

Question

我将 lme4 包用于运行混合模型。我想在单独的数据集中提取固定效应结果和随机效应结果，以便我们可以将其用于进一步分析。可惜我做不到。

例如

mixed_result<- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)

我尝试使用以下方法提取固定效应和随机效应：

fixEffect<-fixef(mixed_result)
randEffect<-ranef(mixed_result)

View(fixEffect)

我分别尝试了 fixef 和 ranef 的固定效果和随机效果，并尝试使用它的结果创建数据集。但它给了我以下错误：

Error in View : cannot coerce class ""ranef.mer"" to a data.frame

我实际上想要我们在 SAS 、 solutionF 和 solutionR 中获得的输出。但是如果不可能得到这样的输出，固定和随机的系数就可以了。

如果有人能帮助我，我将不胜感激。

感谢和问候，

Answer 1

使用str查看对象的结构。

str(fixEffect)
# named vector, can probably be coerced to data.frame

View(as.data.frame(fixEffect))
# works just fine


str(randEffect)
# list of data frames (well, list of one data frame in this case)

View(randEffect$Subject)

如果你有，比方说，斜率也变化 Subject，它们将进入与 Subject 级别截距相同的 Subject 数据帧。但是，如果截距也因其他变量 group 而异，级别数与 Subject 不同，则它们显然不能进入同一数据帧。这就是使用数据框列表的原因，以便相同的结构可以推广到更复杂的模型。

提取 Dataframe 中的固定效应和随机效应

Extract Fixed Effect and Random Effect in Dataframe

r

mixed-models