dotplot_diag 用于烹饪距离（lmer 模型） - 主题值而不是索引答案

【问题标题】：dotplot_diag for cooks distance (lmer model) - subject value instead of indexdotplot_diag 用于烹饪距离（lmer 模型） - 主题值而不是索引
【发布时间】：2014-05-30 15:33:11
【问题描述】：

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

提前感谢您的帮助。

【问题讨论】：

标签： r plot diagnostics lmer

【解决方案1】：

我认为index 参数会起作用：

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

【讨论】：

我试过了（对不起，我应该提到它）但它不起作用：Error in levels(Subject) : object 'Subject' not found。 cooks.distance 函数的输出似乎没有“存储”Subject 值，而只是它们的索引，这使得在没有一些特定主题的情况下返回并重新调整模型不太方便。我想一个解决方案是在cooksd 中包含Subject 信息，但我不知道该怎么做。
你还记得data=Orthodont 参数吗？
效果很好，谢谢！不仔细注意我觉得很傻。