bbwtrimbt:

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

Usage

bbwtrimbt(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)) 
        data = as.matrix(x)
    if (is.matrix(x)) {
        y <- list()
        for (j in 1:ncol(x)) y[[j]] <- x[, j]
        x <- y
    }
    p <- J * K * L
    if (p > length(x)) 
        stop("JKL is less than the Number of groups")
    JK = J * K
    v <- matrix(0, p, p)
    data <- list()
    xx = list()
    for (j in 1:length(x)) {
        data[[j]] <- x[[grp[j]]]
        xx[[j]] = x[[grp[j]]]
        data[[j]] = data[[j]] - mean(data[[j]], tr = tr)
    }
    test.stat = bbwtrim(J, K, L, xx, tr = tr)
    x <- data
    if (SEED) 
        set.seed(2)
    testA = NA
    testB = NA
    testC = NA
    testAB = NA
    testAC = NA
    testBC = NA
    testABC = NA
    bsam = list()
    bdat = list()
    aboot = NA
    bboot = NA
    cboot = NA
    abboot = NA
    acboot = NA
    bcboot = NA
    abcboot = NA
    nvec = NA
    for (j in 1:JK) {
        nvec[j] = length(x[[j]])
        for (ib in 1:nboot) {
            ilow <- 1 - L
            iup = 0
            for (j in 1:JK) {
                ilow <- ilow + L
                iup = iup + L
                nv = length(x[[ilow]])
                bdat[[j]] = sample(nv, size = nv, replace = T)
                for (k in ilow:iup) {
                  bsam[[k]] = x[[k]][bdat[[j]]]
                }
            }
            temp = bbwtrim(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.