从lmList对象中提取RMSE

Question

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

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

此代码正常运行。 HOwever，这是一种提取系数的有效方法吗？

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

Answer 1

在这种情况下，ldply不需要获得系数（你可以简单地使用coef（lm.model））。 要获得残留标准错误：

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