mcp3med:

Usage Arguments Examples

Usage

1
mcp3med(J, K, L, x, tr = 0.2, con = 0, alpha = 0.05, grp = NA, op = FALSE)

Arguments

J
K
L
x
tr
con
alpha
grp
op

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
##---- 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, con = 0, alpha = 0.05, grp = NA, 
    op = FALSE) 
{
    if (is.data.frame(x)) 
        x = as.matrix(x)
    JKL <- J * K * L
    if (is.matrix(x)) 
        x <- listm(x)
    if (!is.na(grp[1])) {
        yy <- x
        x <- list()
        for (j in 1:length(grp)) x[[j]] <- yy[[grp[j]]]
    }
    if (!is.list(x)) 
        stop("Data must be stored in list mode or a matrix.")
    for (j in 1:JKL) {
        xx <- x[[j]]
        x[[j]] <- xx[!is.na(xx)]
    }
    if (JKL != length(x)) 
        warning("The number of groups does not match the number of contrast coefficients.")
    for (j in 1:JKL) {
        temp <- x[[j]]
        temp <- temp[!is.na(temp)]
        x[[j]] <- temp
    }
    temp <- con3way(J, K, L)
    conA <- temp$conA
    conB <- temp$conB
    conC <- temp$conC
    conAB <- temp$conAB
    conAC <- temp$conAC
    conBC <- temp$conBC
    conABC <- temp$conABC
    if (!op) {
        Factor.A <- msmed(x, con = conA, alpha = alpha)
        Factor.B <- msmed(x, con = conB, alpha = alpha)
        Factor.C <- msmed(x, con = conC, alpha = alpha)
        Factor.AB <- msmed(x, con = conAB, alpha = alpha)
        Factor.AC <- msmed(x, con = conAC, alpha = alpha)
        Factor.BC <- msmed(x, con = conBC, alpha = alpha)
        Factor.ABC <- msmed(x, con = conABC, alpha = alpha)
    }
    All.Tests <- NA
    if (op) {
        Factor.A <- NA
        Factor.B <- NA
        Factor.C <- NA
        Factor.AB <- NA
        Factor.AC <- NA
        Factor.BC <- NA
        Factor.ABC <- NA
        con <- cbind(conA, conB, conB, conAB, conAC, conBC, conABC)
        All.Tests <- msmed(x, con = con, alpha = alpha)
    }
    list(Factor.A = Factor.A, Factor.B = Factor.B, Factor.C = Factor.C, 
        Factor.AB = Factor.AB, Factor.AC = Factor.AC, Factor.BC = Factor.BC, 
        Factor.ABC = Factor.ABC, All.Tests = All.Tests, conA = conA, 
        conB = conB, conC = conC, conAB = conAB, conAC = conAC, 
        conBC = conBC, conABC = conABC)
  }

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