box1way:

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

Usage

box1way(x, tr = 0.2, grp = c(1:length(x)))

Arguments

x
tr
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 (x, tr = 0.2, grp = c(1:length(x))) 
{
    J <- length(grp)
    print("The number of groups to be compared is")
    print(J)
    h <- vector("numeric", J)
    w <- vector("numeric", J)
    xbar <- vector("numeric", J)
    svec <- vector("numeric", J)
    for (j in 1:J) {
        xx <- !is.na(x[[j]])
        val <- x[[j]]
        x[[j]] <- val[xx]
        h[j] <- length(x[[grp[j]]]) - 2 * floor(tr * length(x[[grp[j]]]))
        svec[j] <- ((length(x[[grp[j]]]) - 1) * winvar(x[[grp[j]]], 
            tr))/(h[j] - 1)
        xbar[j] <- mean(x[[grp[j]]], tr)
    }
    xtil <- sum(h * xbar)/sum(h)
    fval <- h/sum(h)
    TEST <- sum(h * (xbar - xtil)^2)/sum((1 - fval) * svec)
    nu1 <- sum((1 - fval) * svec)
    nu1 <- nu1^2/((sum(svec * fval))^2 + sum(svec^2 * (1 - 2 * 
        fval)))
    nu2 <- (sum((1 - fval) * svec))^2/sum(svec^2 * (1 - fval)^2/(h - 
        1))
    sig <- 1 - pf(TEST, nu1, nu2)
    list(TEST = TEST, nu1 = nu1, nu2 = nu2, p.value = sig)
  }

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