indt0:

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

Usage

indt0(x, y, nboot = 500, alpha = 0.05, flag = 1, SEED = TRUE)

Arguments

x
y
nboot
alpha
flag
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, nboot = 500, alpha = 0.05, flag = 1, 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)
    mflag <- matrix(NA, nrow = length(y), ncol = length(y))
    for (j in 1:length(y)) {
        for (k in 1:length(y)) {
            mflag[j, k] <- (sum(x[j, ] <= x[k, ]) == ncol(x))
        }
    }
    yhat <- 0
    res <- y - yhat
    print("Taking bootstrap sample, please wait.")
    data <- matrix(runif(length(y) * nboot), nrow = nboot)
    data <- (data - 0.5) * sqrt(12)
    rvalb <- apply(data, 1, indt0sub, yhat, res, mflag, x, tr)
    rvalb <- rvalb/sqrt(length(y))
    dstatb <- apply(abs(rvalb), 2, max)
    wstatb <- apply(rvalb^2, 2, mean)
    mstatb <- apply(abs(rvalb), 2, median)
    dstatb <- sort(dstatb)
    wstatb <- sort(wstatb)
    mstatb <- sort(mstatb)
    v <- c(rep(1, length(y)))
    rval <- indt0sub(v, yhat, res, mflag, x, tr)
    rval <- rval/sqrt(length(y))
    dstat <- NA
    wstat <- NA
    critd <- NA
    critw <- NA
    ib <- round(nboot * (1 - alpha))
    if (flag == 1 || flag == 3) {
        dstat <- max(abs(rval))
        critd <- dstatb[ib]
    }
    if (flag == 2 || flag == 3) {
        wstat <- mean(rval^2)
        critw <- wstatb[ib]
    }
    list(dstat = dstat, wstat = wstat, critd = critd, critw = critw)
  }

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