我正在使用 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
预先感谢您的帮助。