bdm2way:

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

Usage

bdm2way(J, K, x, grp = c(1:p), p = J * K)

Arguments

J
K
x
grp
p

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, grp = c(1:p), p = J * K) 
{
    if (is.matrix(x)) 
        x <- listm(x)
    xx <- list()
    for (j in 1:p) xx[[j]] <- x[[grp[j]]]
    Ja <- matrix(1, J, J)
    Ia <- diag(1, J)
    Pa <- Ia - Ja/J
    Jb <- matrix(1, K, K)
    Ib <- diag(1, K)
    Pb <- Ib - Jb/K
    cona <- kron(Pa, Jb/K)
    conb <- kron(Ja/J, Pb)
    conab <- kron(Pa, Pb)
    outA <- bdms1(xx, cona)
    releff = matrix(outA$q.hat, nrow = J, ncol = K, byrow = TRUE)
    outB <- bdms1(xx, conb)
    outAB <- bdms1(xx, conab)
    list(p.valueA = outA$p.value, p.valueB = outB$p.value, p.valueAB = outAB$p.value, 
        Relative.Effects = releff)
  }

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