search for: covresult

...eterminant Estimator" # http://bioinfo.cipic.ucdavis.edu/publications/print_pub?pub_id=736&category=1 # Drop factors first factors <- names(x)[sapply(x,is.factor)] if(length(factors)>0) x <- x[-factors] # Get the robust location/scale estimates require(MASS) covResult <- cov.rob(x) # Calculate the mahalanobis distance for each datum distance <- mahalanobis(x,covResult$center,covResult$cov) n <- nrow(x) p <- ncol(x) h <- floor((n+p+1)/2) # Asymptotic chi-square method (page 11) # Often identifies too many points as outliers c...