sedm:

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

Usage

sedm(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) 
{
    x <- as.matrix(x)
    if (is.na(y[1])) 
        y <- x[, 2]
    x <- x[, 1]
    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
    temp <- 0 - 2 * (v1 + v2 + v3 + v4)/(fx * fy) + q * (1 - 
        q)/fx^2 + q * (1 - q)/fy^2
    val <- sqrt(temp/n)
    val
  }

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