R 数据帧中固定效应和随机效应的提取

我正在使用lme4包运行混合模型。我想在单独的数据集中提取固定效应结果和随机效应结果，以便我们可以使用它进行进一步的分析。但不幸的是，我不能

例如

---

mixed_result使用
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)

比如说，如果你的斜率也因
主题
而不同，那么它们将位于与
主题
水平截距相同的
主题
数据帧中。然而，如果截获量也因其他变量
组
（与
受试者
的级别数不同）而不同，则它们显然不能进入同一数据帧。这就是为什么要使用数据帧列表，以便相同的结构可以概括为更复杂的模型。
请提供一个具有示例数据和一些代码的可复制示例<代码>帮助（“lmer”）

演示了可复制示例的外观。感谢罗兰的评论。我用小玩具例子更新了这个问题。

fixEffect

是一个（命名的）向量，

randfeffect

是一个data.frames列表<代码>视图只能处理数据帧。您可以打印。由于我不使用SAS，我不知道输出应该是什么样子。谢谢你的回答Gregor！你的回答完全满足了我的要求。但是，对于固定效应和随机效应模型，有可能得到系数的p值吗。。。但不容易。我曾经使用

arm:：se.ranef

以相同的格式输出标准错误。

broom

包应该可以简化这一点，例如，

broom:：tidy（混合结果，effects=“ran\u模式”）

，但我现在在这方面运气不好。这可能是因为我的工作计算机有一个旧版本的R。你可以看到它应该如何工作。感谢Gregor的所有帮助！就我而言，“整洁”起作用了。

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

View(fixEffect)

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)