Сумма остатков в квадрате для множественного вливания с использованием пакета мышей в R
С mice
пакет, как мы проверяем сумму остатков в квадрате объединенного анализа?
library(mice)
imp <- mice(nhanes, seed = 24500)
fit <- with(imp, lm(chl ~ age + bmi))
pool(fit)
summary(pool(fit))
fit
содержит анализ для каждого входного набора данных и pool(fit)
объединенные результаты. Есть ли команда для проверки суммы остатков в квадрате по стандарту? lm
объект, что-то вроде residuals(pool(fit))
?
1 ответ
Вы можете извлечь остатки для каждого входного набора данных с помощью функции res= sapply(fit[[4]], остатки)
и один из вариантов (хотя я не думаю, что он был опубликован где-либо) - это усреднить невязки по множественным вмененным наборам данных, применить (res,1, среднее)