bwamcp:

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

Usage

bwamcp(J, K, x, tr = 0.2, JK = J * K, grp = c(1:JK), alpha = 0.05, op = T)

Arguments

J
K
x
tr
JK
grp
alpha
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 (J, K, x, tr = 0.2, JK = J * K, grp = c(1:JK), alpha = 0.05, 
    op = T) 
{
    if (is.matrix(x)) {
        y <- list()
        for (j in 1:ncol(x)) y[[j]] <- x[, j]
        x <- y
    }
    JK <- J * K
    if (!op) {
        data <- list()
        for (j in 1:length(x)) {
            data[[j]] <- x[[grp[j]]]
        }
        x <- data
        data <- list()
        jp <- 1 - K
        kv <- 0
        for (j in 1:J) {
            jp <- jp + K
            for (k in 1:K) {
                kv <- kv + 1
                if (k == 1) 
                  temp <- x[[jp]]
                if (k > 1) 
                  temp <- c(temp, x[[kv]])
            }
            data[[j]] <- temp
        }
        print("Group numbers refer to levels of Factor A")
        temp <- lincon(data, tr = tr, alpha = alpha)
    }
    if (op) {
        MJK <- K * (J^2 - J)/2
        JK <- J * K
        MJ <- (J^2 - J)/2
        cont <- matrix(0, nrow = J, ncol = MJ)
        ic <- 0
        for (j in 1:J) {
            for (jj in 1:J) {
                if (j < jj) {
                  ic <- ic + 1
                  cont[j, ic] <- 1
                  cont[jj, ic] <- 0 - 1
                }
            }
        }
        tempv <- matrix(0, nrow = K - 1, ncol = MJ)
        con1 <- rbind(cont[1, ], tempv)
        for (j in 2:J) {
            con2 <- rbind(cont[j, ], tempv)
            con1 <- rbind(con1, con2)
        }
        con <- con1
        if (K > 1) {
            for (k in 2:K) {
                con1 <- push(con1)
                con <- cbind(con, con1)
            }
        }
        print("Contrast Matrix Used:")
        print(con)
        temp <- lincon(x, con = con, tr = tr, alpha = alpha)
    }
    temp
  }

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