我目前使用R中的lme
包构建了一个mlm。我试图从该模型中获取Cook距离的度量。我正在使用汽车包装中的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"))