wwwtrimbt:

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

Usage

wwwtrimbt(J, K, L, x, tr = 0.2, JKL = J * K * L, con = 0, alpha = 0.05, grp = c(1:JKL), nboot = 599, SEED = TRUE, ...)

Arguments

J
K
L
x
tr
JKL
con
alpha
grp
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 (J, K, L, x, tr = 0.2, JKL = J * K * L, con = 0, alpha = 0.05, 
    grp = c(1:JKL), nboot = 599, SEED = TRUE, ...) 
{
    if (is.data.frame(x)) 
        x = as.matrix(x)
    if (is.matrix(x)) {
        y <- list()
        for (j in 1:ncol(x)) y[[j]] <- x[, j]
        x = y
    }
    ncon = ncol(con)
    p <- J * K * L
    JKL = p
    if (p > length(x)) 
        stop("JKL is less than the Number of groups")
    JK = J * K
    KL = K * L
    data <- list()
    xx = list()
    for (j in 1:length(x)) {
        xx[[j]] = x[[grp[j]]]
        data[[j]] = xx[[j]] - mean(xx[[j]], tr = tr)
    }
    if (SEED) 
        set.seed(2)
    bsam = list()
    bdat = list()
    aboot = NA
    bboot = NA
    cboot = NA
    abboot = NA
    acboot = NA
    bcboot = NA
    abcboot = NA
    test.stat = wwwtrim(J, K, L, xx, tr = tr)
    nv = length(x[[1]])
    for (ib in 1:nboot) {
        bdat[[j]] = sample(nv, size = nv, replace = T)
        for (k in 1:JKL) bsam[[k]] = data[[k]][bdat[[j]]]
        temp = wwwtrim(J, K, L, bsam, tr = tr)
        aboot[ib] = temp$Qa
        bboot[ib] = temp$Qb
        cboot[ib] = temp$Qc
        acboot[ib] = temp$Qac
        bcboot[ib] = temp$Qbc
        abboot[ib] = temp$Qab
        abcboot[ib] = temp$Qabc
    }
    pbA = NA
    pbB = NA
    pbC = NA
    pbAB = NA
    pbAC = NA
    pbBC = NA
    pbABC = NA
    pbA = mean(test.stat$Qa[1, 1] < aboot)
    pbB = mean(test.stat$Qb[1, 1] < bboot)
    pbC = mean(test.stat$Qc[1, 1] < cboot)
    pbAB = mean(test.stat$Qab[1, 1] < abboot)
    pbAC = mean(test.stat$Qac[1, 1] < acboot)
    pbBC = mean(test.stat$Qbc[1, 1] < bcboot)
    pbABC = mean(test.stat$Qabc[1, 1] < abcboot)
    list(p.value.A = pbA, p.value.B = pbB, p.value.C = pbC, p.value.AB = pbAB, 
        p.value.AC = pbAC, p.value.BC = pbBC, p.value.ABC = pbABC)
  }

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