rdepth:

In musto101/wilcox_R: What the package does (short line)

Usage

rdepth(d, x, y, sortx = T)

Arguments

d
x
y
sortx

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (d, x, y, sortx = T) 
{
    if (!is.vector(x) || !is.vector(y)) 
        stop("x and y should be vectors")
    n <- length(x)
    if (n < 2) 
        stop("you need at least two observations")
    xy <- cbind(x, y)
    b <- d[1]
    a <- d[2]
    if (sortx) 
        xy <- xy[order(xy[, 1], xy[, 2]), ]
    res <- xy[, 2] - a * xy[, 1] - b
    res[abs(res) < 1e-07] <- 0
    posres <- res >= 0
    negres <- res <= 0
    lplus <- cumsum(posres)
    rplus <- lplus[n] - lplus
    lmin <- cumsum(negres)
    rmin <- lmin[n] - lmin
    depth <- pmin(lplus + rmin, rplus + lmin)
    min(depth)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.