r - 从 lmList 对象中提取 RMSE

Question

我有以下代码从 lmList 对象中提取系数。

library(lme4)
library(plyr)
lm.model <- lmList(Y ~ X | eventID, df)
param <- ldply(lm.model, coef)

此代码工作正常。但是，这是提取系数的有效方法吗？

我的主要问题是，如何从同一个 lm.model 对象中提取 RMSE？

Answer 1

在这种情况下，不需要 ldply 来获取系数（您可以简单地使用 coef(lm.model)）。要获得残差标准误差：

ldply(slot(lm.model,".Data"),function(x) summary(x)$sigma)