t1way.effect:

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

Usage

t1way.effect(x, tr = 0.2, grp = NA, MAT = FALSE, lev.col = 1, var.col = 2)

Arguments

x
tr
grp
MAT
lev.col
var.col

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 (x, tr = 0.2, grp = NA, MAT = FALSE, lev.col = 1, var.col = 2) 
{
    library(MASS)
    if (MAT) {
        if (!is.matrix(x)) 
            stop("With MAT=T, data must be stored in a matrix")
        if (length(lev.col) != 1) 
            stop("Argument lev.col should have 1 value")
        temp = selby(x, lev.col, var.col)
        x = temp$x
        grp2 = rank(temp$grpn)
        x = x[grp2]
    }
    if (is.matrix(x)) 
        x <- listm(x)
    if (is.na(sum(grp[1]))) 
        grp <- c(1:length(x))
    if (!is.list(x)) 
        stop("Data are not stored in a matrix or in list mode.")
    J <- length(grp)
    h <- vector("numeric", J)
    w <- vector("numeric", J)
    xbar <- vector("numeric", J)
    pts = NULL
    for (j in 1:J) {
        xx <- !is.na(x[[j]])
        val <- x[[j]]
        pts = c(pts, val)
        x[[j]] <- val[xx]
        h[j] <- length(x[[grp[j]]]) - 2 * floor(tr * length(x[[grp[j]]]))
        w[j] <- h[j] * (h[j] - 1)/((length(x[[grp[j]]]) - 1) * 
            winvar(x[[grp[j]]], tr))
        xbar[j] <- mean(x[[grp[j]]], tr)
    }
    u <- sum(w)
    xtil <- sum(w * xbar)/u
    A <- sum(w * (xbar - xtil)^2)/(J - 1)
    B <- 2 * (J - 2) * sum((1 - w/u)^2/(h - 1))/(J^2 - 1)
    TEST <- A/(B + 1)
    nu1 <- J - 1
    nu2 <- 1/(3 * sum((1 - w/u)^2/(h - 1))/(J^2 - 1))
    sig <- 1 - pf(TEST, nu1, nu2)
    top = var(xbar)
    bot = winvarN(pts, tr = tr)
    e.pow = top/bot
    if (e.pow >= 1) {
        v1 = NULL
        v2 = NULL
        for (j in 1:J) {
            v1 = c(v1, rep(xbar[j], length(x[[j]])))
            v2 = c(v2, x[[j]])
        }
        e.pow = wincor(v1, v2, tr = tr)$cor^2
    }
    list(TEST = TEST, nu1 = nu1, nu2 = nu2, p.value = sig, Var.Explained = e.pow, 
        Effect.Size = sqrt(e.pow))
  }

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