$conf
[1] -0.004942252 0.378251413 $out
[1] -3.315391 2.685922 -3.055717 2.571203
> boxplot.stats(x)$out
[1] -3.315391 2.685922 -3.055717 2.571203 > boxplot(x) > y<-rnorm(100)
> df<-data.frame(x,y) > rm(x,y) > head(df)
x y 1 -3.31539150 0.7619774 2 -0.04765067 -0.6404403 3 0.69720806 0.7645655 4 0.35979073 0.3131930 5 0.18644193 0.1709528 6 0.27493834 -0.8441813 > attach(df)
> # find the index of outliers from x
> (a <-which(x %in% boxplot.stats(x)$out)) [1] 1 33 64 74
> # find the index of outliers from y
> (b <-which(y %in% boxplot.stats(y)$out)) [1] 24 25 49 64 74 > detach(df)
> # outliers in both x and y > (outlier.list<-intersect(a,b)) [1] 64 74 > plot(df)
> points(df[outlier.list,],col=\