1

lmList オブジェクトから係数を抽出する次のコードがあります。

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

このコードは正常に動作しています。しかし、これは係数を抽出する効率的な方法ですか?

私の主な質問は、同じ lm.model オブジェクトから RMSE を抽出するにはどうすればよいですか?

4

1 に答える 1

1

この場合、係数を取得するために ldply は必要ありません (単純に coef(lm.model) を使用できます)。残差標準誤差を取得するには:

ldply(slot(lm.model,".Data"),function(x) summary(x)$sigma)
于 2011-03-04T20:48:19.280 に答える