t2way.no.p:

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

Usage

t2way.no.p(J, K, x, tr = 0.2, grp = c(1:p), alpha = 0.05, p = J * K)

Arguments

J
K
x
tr
grp
alpha
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, tr = 0.2, grp = c(1:p), alpha = 0.05, p = J * 
    K) 
{
    if (is.data.frame(x)) 
        x = as.matrix(x)
    if (is.matrix(x)) 
        x <- listm(x)
    if (!is.list(x)) 
        stop("Data are not stored in a matrix or in list mode")
    if (p != length(x)) {
        print("Warning: The number of groups in your data is not equal to JK")
    }
    for (j in 1:p) x[[j]] <- elimna(x[[j]])
    xbar <- 0
    h <- 0
    d <- 0
    R <- 0
    W <- 0
    d <- 0
    r <- 0
    w <- 0
    nuhat <- 0
    omegahat <- 0
    DROW <- 0
    DCOL <- 0
    xtil <- matrix(0, J, K)
    aval <- matrix(0, J, K)
    for (j in 1:p) {
        xbar[j] <- mean(x[[grp[j]]], tr)
        h[j] <- length(x[[grp[j]]]) - 2 * floor(tr * length(x[[grp[j]]]))
        d[j] <- (length(x[[grp[j]]]) - 1) * winvar(x[[grp[j]]], 
            tr)/(h[j] * (h[j] - 1))
    }
    d <- matrix(d, J, K, byrow = T)
    xbar <- matrix(xbar, J, K, byrow = T)
    h <- matrix(h, J, K, byrow = T)
    for (j in 1:J) {
        R[j] <- sum(xbar[j, ])
        nuhat[j] <- (sum(d[j, ]))^2/sum(d[j, ]^2/(h[j, ] - 1))
        r[j] <- 1/sum(d[j, ])
        DROW[j] <- sum(1/d[j, ])
    }
    for (k in 1:K) {
        W[k] <- sum(xbar[, k])
        omegahat[k] <- (sum(d[, k]))^2/sum(d[, k]^2/(h[, k] - 
            1))
        w[k] <- 1/sum(d[, k])
        DCOL[k] <- sum(1/d[, k])
    }
    D <- 1/d
    for (j in 1:J) {
        for (k in 1:K) {
            xtil[j, k] <- sum(D[, k] * xbar[, k]/DCOL[k]) + sum(D[j, 
                ] * xbar[j, ]/DROW[j]) - sum(D * xbar/sum(D))
            aval[j, k] <- (1 - D[j, k] * (1/sum(D[j, ]) + 1/sum(D[, 
                k]) - 1/sum(D)))^2/(h[j, k] - 3)
        }
    }
    Rhat <- sum(r * R)/sum(r)
    What <- sum(w * W)/sum(w)
    Ba <- sum((1 - r/sum(r))^2/nuhat)
    Bb <- sum((1 - w/sum(w))^2/omegahat)
    Va <- sum(r * (R - Rhat)^2)/((J - 1) * (1 + 2 * (J - 2) * 
        Ba/(J^2 - 1)))
    Vb <- sum(w * (W - What)^2)/((K - 1) * (1 + 2 * (K - 2) * 
        Bb/(K^2 - 1)))
    nu2 <- (J^2 - 1)/(3 * Ba)
    sig.A <- 1 - pf(Va, J - 1, nu2)
    nu2 <- (K^2 - 1)/(3 * Bb)
    sig.B <- 1 - pf(Vb, K - 1, nu2)
    Vab <- sum(D * (xbar - xtil)^2)
    dfinter <- (J - 1) * (K - 1)
    crit <- qchisq(1 - alpha, dfinter)
    hc <- (crit/(2 * dfinter)) * (1 + (3 * crit)/(dfinter + 2)) * 
        sum(aval)
    adcrit <- crit + hc
    list(Qa = Va, sig.A = sig.A, Qb = Vb, sig.B = sig.B, Qab = Vab, 
        critinter = adcrit)
  }

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