cooks距离点图(lmer模型)-主题值而非索引
我正在使用HLMdiag包中的工具查看由cooks距离点图(lmer模型)-主题值而非索引,r,plot,diagnostics,lmer,R,Plot,Diagnostics,Lmer,我正在使用HLMdiag包中的工具查看由lmer安装的型号的诊断。我有一个关于绘制cooks距离的dotplot\u diag函数的问题 这是我正在做的一个例子 mod <- lmer(distance~age+Sex+(1|Subject), data=Orthodont) cooksd <- cooks.distance(mod, group = "Subject") dotplot_diag(x = cooksd, cutoff = "internal", name = "
lmer
安装的型号的诊断。我有一个关于绘制cooks距离的dotplot\u 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")
提前感谢您的帮助。我认为
索引
参数会起作用:
dotplot_diag(x = cooksd, index=levels(Subject),
data=Orthodont, cutoff = "internal", name = "cooks.distance") +
ylab("Cook's distance") + xlab("ID")
我已经试过了(对不起,我应该提到它),但它不起作用:
级别错误(主题):找不到对象“主题”
。cooks.distance
函数的输出似乎并不“存储”Subject
值,而只是它们的索引,这使得返回并重新安装没有几个特定主题的模型不太方便。我想一个解决方案是在cooksd
中包含Subject
信息,但我不知道怎么做。你还记得data=orthont
参数吗?它非常有效,谢谢!我觉得自己太傻了,没有多加注意。
dotplot_diag(x = cooksd, index=levels(Subject),
data=Orthodont, cutoff = "internal", name = "cooks.distance") +
ylab("Cook's distance") + xlab("ID")