qcorp1:

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

Usage

qcorp1(x, y, q = 0.5, xout = FALSE, outfun = outpro, plotit = FALSE, ...)

Arguments

x
y
q
xout
outfun
plotit
...

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 (x, y, q = 0.5, xout = FALSE, outfun = outpro, plotit = FALSE, 
    ...) 
{
    X = cbind(x, y)
    X = elimna(X)
    x <- as.matrix(x)
    p = ncol(x)
    x = X[, 1:p]
    p1 = p + 1
    y = X[, p1]
    if (xout) {
        x <- as.matrix(x)
        flag <- outfun(x, plotit = plotit, ...)$keep
        x <- x[flag, ]
        y <- y[flag]
        x <- as.matrix(x)
        X = cbind(x, y)
    }
    est = qreg(x, y, q = q)$coef
    top = qreg.sub(X, est, qval = q)
    null = hd(y, q)
    v = c(null, rep(0, p))
    bot = qreg.sub(X, v, qval = q)
    ce = sqrt(1 - top/bot)
    if (p == 1) 
        ce = sign(est[2]) * ce
    list(cor = ce)
  }

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