hc4qtest:

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

Usage

hc4qtest(x, y, k, nboot = 500, SEED = TRUE)

Arguments

x
y
k
nboot
SEED

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, k, nboot = 500, SEED = TRUE) 
{
    if (SEED) 
        set.seed(2)
    x <- as.matrix(x)
    temp <- cbind(x, y)
    temp <- elimna(temp)
    pval <- ncol(temp) - 1
    x <- temp[, 1:pval]
    y <- temp[, pval + 1]
    x <- as.matrix(x)
    p <- ncol(x)
    pp <- p + 1
    temp <- lsfit(x, y)
    yhat <- mean(y)
    res <- y - yhat
    s <- lsfitNci4(x, y)$cov[-1, -1]
    s <- as.matrix(s)
    si <- s[k, k]
    b <- temp$coef[2:pp]
    qtest <- b[k]/sqrt(si)
    data <- matrix(runif(length(y) * nboot), nrow = nboot)
    data <- (data - 0.5) * sqrt(12)
    rvalb <- apply(data, 1, lsqtest4, yhat, res, x, k)
    sum <- sum(abs(rvalb) >= abs(qtest[1]))
    p.val <- sum/nboot
    list(p.value = p.val)
  }

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