我试图弄清楚如何从一个只有随机截距和一个变量（性别）的简单 lmer 模型中组织随机效应的 ranef.mer 列表。
fit.b <- lmer(Math ~ 1 + Sex + (1+Sex|SchoolID), data=pisa_com, REML=FALSE)
我已经使用 qqmath 绘制了随机效应，但我要么需要能够通过它们的簇号（在本例中为学校）来标记每个随机效应，要么组织 ranef.mer 输出。
问问题
734 次
1 回答
0
昨晚解决了这个问题。ranef.mer 可以强制转换为数据框。
我适合模型：
fit.b <- lmer(Math ~ 1 + Sex + (1+Sex|SchoolID), data=pisa_com, REML=FALSE)
然后通过包含识别变量将其强制转换为数据框
random.effects <- as.data.frame(ranef(fit.b)$SchoolID)
然后将其写入 .csv 以在 excel 中进行排序
write.csv(random.effects, file="~/folder/file.name.csv")
于 2013-11-07T15:12:58.513 回答