我使用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。但是如果它不可能得到这样的输出,那么固定和随机的系数就会这样。
如果有人可以帮助我,我将不胜感激。
谢谢和问候,
答案 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
,则它们显然无法进入同一数据框。这就是使用数据框列表的原因,因此相同的结构可以概括为更复杂的模型。