我使用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
提前感谢您的帮助。
答案 0 :(得分:2)
我认为index
参数可行:
dotplot_diag(x = cooksd, index=levels(Subject),
data=Orthodont, cutoff = "internal", name = "cooks.distance") +
ylab("Cook's distance") + xlab("ID")