rmanovab:

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

Usage

rmanovab(x, tr = 0.2, alpha = 0.05, grp = 0, nboot = 599)

Arguments

x
tr
alpha
grp
nboot

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, tr = 0.2, alpha = 0.05, grp = 0, nboot = 599) 
{
    if (is.data.frame(x)) 
        x = as.matrix(x)
    if (!is.list(x) && !is.matrix(x)) 
        stop("Data must be stored in a matrix or in list mode.")
    if (is.list(x)) 
        mat = matl(x)
    if (is.matrix(x)) {
        if (sum(grp) == 0) 
            grp <- c(1:ncol(x))
        mat <- x[, grp]
    }
    mat = elimna(mat)
    J <- ncol(mat)
    connum <- (J^2 - J)/2
    bvec <- matrix(0, connum, nboot)
    set.seed(2)
    print("Taking bootstrap samples. Please wait.")
    data <- matrix(sample(nrow(mat), size = nrow(mat) * nboot, 
        replace = TRUE), nrow = nboot)
    xcen <- matrix(0, nrow(mat), ncol(mat))
    for (j in 1:J) xcen[, j] <- mat[, j] - mean(mat[, j], tr)
    bvec <- apply(data, 1, tsubrmanovab, xcen, tr)
    icrit <- round((1 - alpha) * nboot)
    bvec <- sort(bvec)
    crit <- bvec[icrit]
    test <- rmanova(mat, tr, grp)$test
    list(teststat = test, crit = crit)
  }

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