iband:

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

Usage

iband(x, alpha = 0.05, plotit = TRUE, sm = TRUE, SEED = TRUE, nboot = 500, grp = c(1:4), xlab = "First Difference", ylab = "Delta", crit = NA, print.all = FALSE, plot.op = FALSE)

Arguments

x
alpha
plotit
sm
SEED
nboot
grp
xlab
ylab
crit
print.all
plot.op

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, alpha = 0.05, plotit = TRUE, sm = TRUE, SEED = TRUE, 
    nboot = 500, grp = c(1:4), xlab = "First Difference", ylab = "Delta", 
    crit = NA, print.all = FALSE, plot.op = FALSE) 
{
    if (!is.list(x) && !is.matrix(x)) 
        stop("store data in list mode or a matrix")
    if (SEED) 
        set.seed(2)
    if (is.matrix(x)) 
        x <- listm(x)
    if (length(x) < 4) 
        stop("There must be at least 4 groups")
    for (j in 1:length(x)) x[[j]] = elimna(x[[j]])
    if (length(grp) != 4) 
        stop("The argument grp must have 4 values")
    x <- x[grp]
    n <- c(length(x[[1]]), length(x[[2]]), length(x[[3]]), length(x[[4]]))
    vals <- NA
    y <- list()
    if (is.na(crit)) {
        print("Approximating critical value. Please wait.")
        for (i in 1:nboot) {
            for (j in 1:4) y[[j]] <- rnorm(n[j])
            temp <- ks.test(outer(y[[1]], y[[2]], FUN = "-"), 
                outer(y[[3]], y[[4]], FUN = "-"))
            vals[i] <- temp[1]$statistic
        }
        vals <- sort(vals)
        ic <- (1 - alpha) * nboot
        crit <- vals[ic]
    }
    if (plot.op) {
        plotit <- F
        g2plot(v1, v2)
    }
    output <- sband(outer(x[[1]], x[[2]], FUN = "-"), outer(x[[3]], 
        x[[4]], FUN = "-"), plotit = plotit, crit = crit, flag = FALSE, 
        sm = sm, xlab = xlab, ylab = ylab)
    if (!print.all) {
        numsig <- output$numsig
        ks.test.stat <- ks.test(outer(x[[1]], x[[2]], FUN = "-"), 
            outer(x[[3]], x[[4]], FUN = "-"))$statistic
        output <- matrix(c(numsig, crit, ks.test.stat), ncol = 1)
        dimnames(output) <- list(c("number sig", "critical value", 
            "KS test statistics"), NULL)
    }
    output
  }

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