在数据帧中提取固定效果和随机效果

时间:2016-03-21 16:31:36

标签: r mixed-models

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

E.g。

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

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

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

View(fixEffect)

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

  

视图错误:无法强制上课&#34;&#34; ranef.mer&#34;&#34;到data.frame

我实际上想要输出SAS,solutionF和solutionR。但是如果它不可能得到这样的输出,那么固定和随机的系数就会这样。

如果有人可以帮助我,我将不胜感激。

谢谢和问候,

1 个答案:

答案 0 :(得分: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,则它们显然无法进入同一数据框。这就是使用数据框列表的原因,因此相同的结构可以概括为更复杂的模型。