mvoutlier (version 2.1.1)

chisq.plot: Chi-Square Plot

Description

The function chisq.plot plots the ordered robust mahalanobis distances of the data against the quantiles of the Chi-squared distribution. By user interaction this plotting is iterated each time leaving out the observation with the greatest distance.

Usage

chisq.plot(x, quan=1/2, ask=TRUE, ...)

Arguments

x

matrix or data.frame containing the data

quan

amount of observations which are used for mcd estimations. has to be between 0.5 and 1, default ist 0.5

ask

logical. specifies whether user interacton is allowed or not. default is TRUE

...

additional graphical parameters

Value

outliers

indices of the outliers that are removed by left-click on the plotting device.

Details

The function chisq.plot plots the ordered robust mahalanobis distances of the data against the quantiles of the Chi-squared distribution. If the data is normal distributed these values should approximately correspond to each other, so outliers can be detected visually. By user interaction this procedure is repeated, each time leaving out the observation with the greatest distance (the number of the observation is printed on the console). This method can be seen as an iterative deletion of outliers until a straight line appears.

References

R.G. Garrett (1989). The chi-square plot: a tools for multivariate outlier recognition. Journal of Geochemical Exploration, 32 (1/3), 319-341.

Examples

# NOT RUN {
data(humus)
res <-chisq.plot(log(humus[,c("Co","Cu","Ni")]))
res$outliers # these are the potential outliers
# }