ESmainMCP:

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

Usage

ESmainMCP(J, K, x, tr = 0.2, nboot = 100, SEED = TRUE)

Arguments

J
K
x
tr
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, tr = 0.2, nboot = 100, SEED = TRUE) 
{
    if (is.matrix(x)) 
        x = listm(x)
    x = lapply(x, elimna)
    con = con2way(J, K)
    conA = con$conA
    FA = matrix(NA, nrow = ncol(conA), ncol = 3)
    ic = 0
    for (jj in 1:J) {
        for (jjj in 1:J) {
            if (jj < jjj) {
                ic = ic + 1
                FA[ic, 1] = jj
                FA[ic, 2] = jjj
            }
        }
    }
    for (j in 1:ncol(conA)) {
        flag1 = (conA[, j] == 1)
        flagm1 = (conA[, j] == -1)
        x1 = as.vector(matl(x[flag1]))
        x2 = as.vector(matl(x[flagm1]))
        FA[j, 3] = yuenv2(x1, x2, tr = tr, nboot = nboot, SEED = SEED)$Effect.Size
    }
    dimnames(FA) <- list(NULL, c("Level", "Level", "Effect.Size"))
    conB = con$conB
    FB = matrix(NA, nrow = ncol(conB), ncol = 3)
    ic = 0
    for (jj in 1:K) {
        for (jjj in 1:K) {
            if (jj < jjj) {
                ic = ic + 1
                FB[ic, 1] = jj
                FB[ic, 2] = jjj
            }
        }
    }
    for (j in 1:ncol(conB)) {
        for (jj in 1:J) {
            for (jjj in 1:J) {
                if (jj < jjj) {
                }
            }
        }
        flag1 = (conB[, j] == 1)
        flagm1 = (conB[, j] == -1)
        x1 = as.vector(matl(x[flag1]))
        x2 = as.vector(matl(x[flagm1]))
        FB[j, 3] = yuenv2(x1, x2, tr = tr, nboot = nboot, SEED = SEED)$Effect.Size
    }
    dimnames(FB) <- list(NULL, c("Level", "Level", "Effect.Size"))
    list(Factor.A = FA, Factor.B = FB)
  }

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