我试图标签上使用卡方QQ图离群mvOutlier()
的功能MVN
包R
。
我设法通过标签识别异常值并获得它们的x坐标。我尝试使用 将前者放在图上text()
,但x和y坐标似乎被翻转了。
基于文档中的示例:
library(MVN)
data(iris)
versicolor <- iris[51:100, 1:3]
# Mahalanobis distance
result <- mvOutlier(versicolor, qqplot = TRUE, method = "quan")
labelsO<-rownames(result$outlier)[result$outlier[,2]==TRUE]
xcoord<-result$outlier[result$outlier[,2]==TRUE,1]
text(xcoord,label=labelsO)
Run Code Online (Sandbox Code Playgroud)
这会产生以下结果:
我也试过text(x = xcoord, y = xcoord,label = labelsO)
,当点靠近 y = x 线时很好,但当不满足正态性时可能会失败(并且点偏离这条线)。
有人可以建议如何访问卡方分位数或为什么函数的x坐标text()
似乎不服从输入参数。