dotplot_diag for cooks distance(lmer model) - 主题值而不是索引

时间:2014-05-30 15:33:11

标签: r plot diagnostics lmer

我使用HLMdiag包中的工具来查看由lmer拟合的模型的诊断。我对绘制厨师距离的dotplot_diag函数有疑问。

以下是我正在做的事情的一个例子。

mod <- lmer(distance~age+Sex+(1|Subject), data=Orthodont)

cooksd <- cooks.distance(mod, group = "Subject")

dotplot_diag(x = cooksd, cutoff = "internal", name = "cooks.distance") + ylab("Cook's distance") + xlab("ID")

我得到以下情节: enter image description here

问题:有没有办法让有影响力的点标有Subject值而不是他们的索引?即它会显示类似M13或F10的内容(Subject变量中的值,见下文)。

head(unique(Orthodont$Subject))
[1] M01 M02 M03 M04 M05 M06

提前感谢您的帮助。

1 个答案:

答案 0 :(得分:2)

我认为index参数可行:

dotplot_diag(x = cooksd, index=levels(Subject), 
 data=Orthodont, cutoff = "internal", name = "cooks.distance") + 
 ylab("Cook's distance") + xlab("ID")