apanova:

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

Usage

apanova(data, grp = 0)

Arguments

data
grp

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 (data, grp = 0) 
{
    if (is.list(data)) {
        x <- matrix(0, length(data[[1]]), length(data))
        for (j in 1:length(data)) x[, j] <- data[[j]]
    }
    if (is.matrix(data)) 
        x <- data
    if (sum(grp == 0)) 
        grp <- c(1:ncol(x))
    x <- x[, grp]
    J <- ncol(x)
    n <- nrow(x)
    if (n <= 20) 
        print("With n<=20, suggest using bprm")
    rm <- matrix(rank(x), n, J)
    rv <- apply(rm, 2, mean)
    sm <- (n - 1) * winall(rm, tr = 0)$cov/(n - J + 1)
    jm1 <- J - 1
    cv <- diag(1, jm1, J)
    for (i in 2:J) {
        k <- i - 1
        cv[k, i] <- -1
    }
    cr <- cv %*% rv
    ftest <- n * t(cr) %*% solve(cv %*% sm %*% t(cv)) %*% cr/(J - 
        1)
    df1 <- J - 1
    df2 <- (J - 1) * (n - 1)
    siglevel <- 1 - pf(ftest, df1, df2)
    list(FTEST = ftest, df1 = df1, df2 = df2, p.value = siglevel)
  }

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