9.3 Multivariate Ausreiss ererkennung

		library(chemometrics) data(glass) X=glass[,c("MgO","Cl")] data(glass.grp) library(robustbase) X_mcd=covMcd(X) pdf("outlier1.pdf",width=9,height=4.5) par(mfrow=c(1,2)) par(mar=c(4,4,2,2)) drawMahal(X,center=apply(X,2,mean),covariance=cov(X), quantile=0.975,pch=glass.grp,cex.lab=1.2) drawMahal(X,center=X_mcd$center,covariance=X_mcd$cov, quantile=0.975,pch=glass.grp,cex.lab=1.2) dev.off()
		library(chemometrics) data(glass) X=glass[,c("MgO","Cl")] data(glass.grp) pdf("outlier2.pdf",width=9,height=4.5) par(mar=c(4,4,2,2)) set.seed(100) res <- Moutlier(X,quantile=0.975,pch=glass.grp,cex.lab=1.2) dev.off()