我正在使用 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
预先感谢您的帮助。
最佳答案
我认为 index
参数会起作用:
dotplot_diag(x = cooksd, index=levels(Subject),
data=Orthodont, cutoff = "internal", name = "cooks.distance") +
ylab("Cook's distance") + xlab("ID")
关于r - 厨师距离的 dotplot_diag(lmer 模型)- 主题值而不是索引,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/23958236/