Использование effect.lme (пакет автомобилей) с функцией cooks.distance с моделью lme из nlme в R
В настоящее время я построил MLM, используя lme
пакет в R. Я пытаюсь измерить расстояние повара от этой модели. Я использую influence
функция из пакета автомобилей, которую в документации сказано, что вы можете прочитать в cooks.distance
функция. Тем не менее, когда я положил результат этой функции в cooks.distance
Функция возвращает значения в диапазоне от 4,27 до 1609,96.
acc_2_d_l2b <- lme(Acc~load1+load2+PSWQd_Acc_C_p58, random=~1|Subject,
weights=varIdent(form=~1|trial), data=final_data_long_p58,
method="REML", na.action="na.exclude")
b <- influence(acc_2_d_l2b, group=c("Subject"), obs=TRUE)
a <- cooks.distance(b, group=c("Subject"))