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

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

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

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

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


r plyr lme4
person ASC    schedule 04.03.2011    source источник

Ответы (1)


arrow_upward
1
arrow_downward

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

ldply(slot(lm.model,".Data"),function(x) summary(x)$sigma)
person Wojciech Sobala    schedule 04.03.2011
comment
Это работает, я был близок к правильному... Я не знал, что в R это называется сигмой. Спасибо. - person ASC; 05.03.2011