Использование effect.lme (пакет автомобилей) с функцией cooks.distance с моделью lme из nlme в R - PullRequest
0 голосов
/ 07 октября 2018

Я в настоящее время построил млм, используя пакет lme в R. Я пытаюсь измерить расстояние повара от этой модели.Я использую функцию influence из пакета cars, которую в документации сказано, что вы можете прочитать в функцию 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"))
...