Извлечение RMSE из объекта lmList

Question

У меня есть следующий код для извлечения коэффициентов из объекта lmList.

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

Этот код работает нормально.Однако, это эффективный способ извлечения коэффициентов?

Мой главный вопрос: как мне извлечь RMSE из того же объекта lm.model?

Answer 1

В этом случае ldply не требуется для получения коэффициентов (вы можете просто использовать coef (lm.model)). Чтобы получить остаточную стандартную ошибку:

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