plot outliers in X and Y data

find_outliers.r

#The input is two colums: x and corresponding y values
require(MASS) ## for mvrnorm()
set.seed(1)
mine <- read.table("xy.txt")
mine <- data.frame(mine)
names(mine) <- c("X","Y")
plot(mine)
res <- resid(mod <- lm(Y ~ X, data = mine))
res.qt <- quantile(res, probs = c(0.001,0.999))
want <- which(res >= res.qt[1] & res <= res.qt[2])
pdf("out.pdf")
plot(mine, type = "n")
#mine[-want,]
points(mine[-want,], col = "black", pch = 21, bg = "black", cex = 0.8)
points(mine[want,], col = "red", pch = 21, bg = "red", cex = 0.8)
abline(mod, col = "blue", lwd = 2)
dev.off()