r - 访问LME图中的异常值ID

我在r中绘制了一个lme fit对象，并在图形上显示了异常ID(studyID)，但我想通过在图对象中查找它们来自动访问这些ID。我不知道该怎么做。我正在做很多分析，因此能够自动执行此操作而不是实际查看每个图形中的异常ID号将有所帮助。

这是我正在做的简化示例:

fit <- lme(dv ~ studyID + Gender + Group * DOP, random=~1|studyID, cor=corSymm(), na.action="na.omit", method="ML", data=x$data)

require (car)

plotObject <- plot(fit, resid(., type = "p") ~ fitted(.) | Group*DOP, abline = 0, id=.05)

我想要做的是访问plotObject的某些属性，该属性存储用于标识由plot语句导致的图形中的异常值的id号。

谢谢你。

最佳答案

我不知道此信息是否实际存储在绘图对象中，但是计算您自己很容易。从?plot.lme:



所以我会说类似

library(nlme)
fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age
which(abs(residuals(fm1,type="normalized"))>qnorm(0.975))
## M09 M09 M13
##  34  35  49
plot(fm1,id=.05)  ## for comparison

似乎可以解决问题。