我正在r中绘制一个lme fit对象,并在图表上显示异常值id(studyID),但我想通过在plot对象中查找它们来自动访问这些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)
Run Code Online (Sandbox Code Playgroud)
我想要做的是访问plotObject的一些属性,该属性存储用于识别图表中由plot语句产生的异常值的id号.
谢谢.
给定输入向量(iv)
iv <- c(.10,.15,"hello","."," . ",". ")
Run Code Online (Sandbox Code Playgroud)
我正在使用:
out <- sub(regexp,NA,iv)
Run Code Online (Sandbox Code Playgroud)
我想要这样的输出向量:
.10,.15,"hello",NA,NA,NA
Run Code Online (Sandbox Code Playgroud)
但是,不知道如何形成正则表达式以获得我需要的东西.提前致谢.