https://github.com/cran/robCompositions
Tip revision: 2cc9b1ae6488319f4ae1d1655243b80e72fdcbe8 authored by Matthias Templ on 04 May 2011, 00:00:00 UTC
version 1.5.0
version 1.5.0
Tip revision: 2cc9b1a
plot.outCoDa.R
plot.outCoDa <- function(x, y, ...){
# plot(1:length(x$mahalDist), x$mahalDist, ylab="Mahalanobis distance", xlab="Index", type="n", ylim=c(0, max(x$mahalDist)) )
# points(1:length(x$mahalDist[x$outlierIndex]), x$mahalDist[x$outlierIndex], pch=3, col="red")
# points(1:length(x$mahalDist[!x$outlierIndex]), x$mahalDist[!x$outlierIndex])
# abline(h=x$limit, lty=2)
# #print(x$mahalDist)
if(x$method =="standard") yl <- "Mahalanobis distance" else yl <- "Robust Mahalanobis distance"
plot(1:length(x$mahalDist), x$mahalDist, ylab = yl,
xlab = "Index", ylim = c(0, max(x$mahalDist)),
pch = as.numeric(x$outlierIndex)*2+1)
abline(h=x$limit, lty=2)
}