在R中使用来自带有nme的lme模型的带cooks.distance函数的Impact.lme(汽车包装)

时间:2018-10-06 21:44:21

标签: r nlme mlm

我目前使用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"))

0 个答案:

没有答案