bwmedbmcp:

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

Usage

bwmedbmcp(J, K, x, JK = J * K, grp = c(1:JK), con = 0, alpha = 0.05, dif = FALSE, pool = FALSE, bop = FALSE, nboot = 100, SEED = TRUE)

Arguments

J
K
x
JK
grp
con
alpha
dif
pool
bop
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, x, JK = J * K, grp = c(1:JK), con = 0, alpha = 0.05, 
    dif = FALSE, pool = FALSE, bop = FALSE, nboot = 100, SEED = TRUE) 
{
    if (is.matrix(x)) {
        y <- list()
        for (j in 1:ncol(x)) y[[j]] <- x[, j]
        x <- y
    }
    JK <- J * K
    data <- list()
    for (j in 1:length(x)) {
        data[[j]] <- x[[grp[j]]]
    }
    x <- data
    if (pool) {
        data <- list()
        m1 <- matrix(c(1:JK), J, K, byrow = T)
        for (k in 1:K) {
            for (j in 1:J) {
                flag <- m1[j, k]
                if (j == 1) 
                  temp <- x[[flag]]
                if (j > 1) {
                  temp <- c(temp, x[[flag]])
                }
            }
            data[[k]] <- temp
        }
        print("Group numbers refer to levels of Factor B")
        if (!dif) 
            temp <- lincdm(data, con = con, alpha = alpha, nboot = nboot, 
                mop = bop)
        if (dif) 
            temp <- qdmcpdif(data, con = con, alpha = alpha)
        return(temp)
    }
    if (!pool) {
        mat <- matrix(c(1:JK), ncol = K, byrow = T)
        for (j in 1:J) {
            data <- list()
            ic <- 0
            for (k in 1:K) {
                ic <- ic + 1
                data[[ic]] <- x[[mat[j, k]]]
            }
            print(paste("For level ", j, " of Factor A:"))
            if (!dif) 
                temp <- lincdm(data, con = con, alpha = alpha, 
                  nboot = nboot, mop = bop)
            if (dif) 
                temp <- qdmcpdif(data, con = con, alpha = alpha)
            print(temp$test)
            print(temp$psihat)
        }
    }
  }

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