zdist:

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

Usage

zdist(m, pts = m, zloc = median, zscale = mad)

Arguments

m
pts
zloc
zscale

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 (m, pts = m, zloc = median, zscale = mad) 
{
    if (!is.matrix(m)) 
        stop("argument m should be a matrix")
    if (!is.matrix(pts)) 
        stop("argument pts should be a matrix")
    if (ncol(m) != ncol(pts)) 
        stop("Number of columns for m and pts are not equal")
    np <- ncol(m)
    val <- NA
    for (i in 1:nrow(pts)) {
        pval <- pts[i, ]
        START <- rep(1, np)/sqrt(np)
        temp <- nelderv2(m, np, FN = zdepth.sub, START = START, 
            zloc = zloc, zscale = zscale, pts = pval)
        temp <- temp/sqrt(sum(temp^2))
        y <- t(t(m) * temp)
        y <- apply(y, 1, sum)
        ppro <- sum(pval * temp)
        val[i] <- abs(ppro - zloc(y))/zscale(y)
    }
    val
  }

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