按升序或降序组织ranef.mer

Question

我试图弄清楚如何从一个只有随机截距和一个变量（性别）的简单lmer模型中组织ranef.mer随机效应列表。

fit.b <- lmer(Math ~ 1 + Sex + (1+Sex|SchoolID), data=pisa_com, REML=FALSE)

我使用qqmath绘制了随机效果，但我需要能够通过群集号（在本例中为学校）标记每个随机效果，或者组织ranef.mer输出。

Answer 1

昨晚解决了这个问题。 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")