从lmList对象中提取RMSE

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

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

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

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

2011-03-04 ASC

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

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

2011-03-04 20:48:19

它的工作原理，我很接近得到它正确的......我不知道，这就是所谓的R.感谢西格玛。 – ASC 2011-03-05 05:37:30

我也用'ldply（lm.model，function（x）summary（x）$ sigma）'代替... – ASC 2011-03-05 05:49:25

回答