johansp:

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

Usage

johansp(cmat, vmean, vsqse, h, J, K)

Arguments

cmat
vmean
vsqse
h
J
K

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 (cmat, vmean, vsqse, h, J, K) 
{
    p <- J * K
    yvec <- matrix(vmean, length(vmean), 1)
    test <- cmat %*% vsqse %*% t(cmat)
    invc <- solve(test)
    test <- t(yvec) %*% t(cmat) %*% invc %*% cmat %*% yvec
    temp <- 0
    klim <- 1 - K
    kup <- 0
    for (j in 1:J) {
        klim <- klim + K
        kup <- kup + K
        Q <- matrix(0, p, p)
        for (k in klim:kup) Q[k, k] <- 1
        mtem <- vsqse %*% t(cmat) %*% invc %*% cmat %*% Q
        temp[j] <- (sum(diag(mtem %*% mtem)) + (sum(diag(mtem)))^2)/(h[j] - 
            1)
    }
    A <- 0.5 * sum(temp)
    df1 <- nrow(cmat)
    df2 <- nrow(cmat) * (nrow(cmat) + 2)/(3 * A)
    cval <- nrow(cmat) + 2 * A - 6 * A/(nrow(cmat) + 2)
    test <- test/cval
    sig <- 1 - pf(test, df1, df2)
    list(teststat = test[1], p.value = sig)
  }

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