Nothing
outlier <- function(x, ...) UseMethod("outlier")
outlier.RRF <- function(x, ...) {
if (!inherits(x, "RRF")) stop("x is not a RRF object")
if (x$type == "regression") stop("no outlier measure for regression")
if (is.null(x$proximity)) stop("no proximity measures available")
outlier.default(x$proximity, x$y)
}
outlier.default <- function(x, cls=NULL, ...) {
if (nrow(x) != ncol(x)) stop ("x must be a square matrix")
n <- nrow(x)
if (is.null(cls)) cls <- rep(1, n)
cls <- factor(cls)
lvl <- levels(cls)
cls.n <- table(cls)[lvl]
id <- if (is.null(rownames(x))) 1:n else rownames(x)
outlier <- structure(rep(NA, n), names=id)
for (i in lvl) {
iclass <- cls == i
out <- rowSums(x[iclass, iclass]^2)
out <- n / ifelse(out == 0, 1, out)
out <- (out - median(out)) / mad(out)
outlier[iclass] <- out
}
outlier
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.