cov2med:

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

Usage

cov2med(x, y = NA, q = 0.5)

Arguments

x
y
q

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 = NA, q = 0.5) 
{
    if (is.na(y[1])) 
        val <- qse(x, q = q, op = 3)^2
    if (!is.na(y[1])) {
        if (sum((x - y)^2) == 0) 
            val <- qse(x, q = q, op = 3)^2
        if (sum((x - y)^2) > 0) {
            n <- length(x)
            m <- floor(q * n + 0.5)
            yord <- sort(y)
            flag <- (y <= yord[m])
            xord <- sort(x)
            xq <- xord[m]
            yord <- sort(y)
            yq <- yord[m]
            flag1 <- (x <= xq)
            flag2 <- (y <= yq)
            A <- mean(flag1 * flag2)
            flag1 <- (x <= xq)
            flag2 <- (y > yq)
            B <- mean(flag1 * flag2)
            flag1 <- (x > xq)
            flag2 <- (y <= yq)
            C1 <- mean(flag1 * flag2)
            flag1 <- (x > xq)
            flag2 <- (y > yq)
            D1 <- mean(flag1 * flag2)
            fx <- akerd(x, pts = xq, plotit = FALSE, pyhat = T)
            fy <- akerd(y, pts = yq, plotit = FALSE, pyhat = T)
            v1 <- (q - 1)^2 * A
            v2 <- (q - 1) * q * B
            v3 <- (q - 1) * q * C1
            v4 <- q * q * D1
            val <- ((v1 + v2 + v3 + v4)/(fx * fy))/n
        }
    }
    val
  }

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