我想用厨师在R的距离来检测离群值。
我正在使用r-package“betareg”(用于beta回归)。所以,很难使用普通的cooks.distance图。
但是使用r命令:cooks.distance(model)我得到一个带有每个观察点距离的矢量的答案。
mod1<- betareg(y~ a+b+c+d|a+b+c+d, data=d)
cooks.distance(mod1) #returns a vector
但我不知道如何将这个矢量转换成如下图所示(带有观察数字):
感谢您的帮助!
答案 0 :(得分:0)
cd <- cooks.distance(mod1)
plot(1:length(cd), cd, type = "h", ylim = range(cd)*c(1, 1.1))# stretch y-axis to make space for labels
threshold <- 0.25 # choose appropriate value
lab <- cd > thresh
myLabels <- 1:length(cd) # choose appropriate labels
text(which(lab), cd[lab], labels = myLabels[lab], pos = 3)# pos 3 puts labels above position