3

我正在使用 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")

我得到以下情节: 在此处输入图像描述

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

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

预先感谢您的帮助。

4

1 回答 1

2

我认为这个index论点会起作用:

dotplot_diag(x = cooksd, index=levels(Subject), 
 data=Orthodont, cutoff = "internal", name = "cooks.distance") + 
 ylab("Cook's distance") + xlab("ID")
于 2014-05-30T18:22:50.560 回答