yuenv2:

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

Usage

yuenv2(x, y = NULL, tr = 0.2, alpha = 0.05, plotit = FALSE, plotfun = splot, op = TRUE, VL = TRUE, cor.op = FALSE, loc.fun = median, xlab = "Groups", ylab = "", PB = FALSE, nboot = 100, SEED = TRUE)

Arguments

x
y
tr
alpha
plotit
plotfun
op
VL
cor.op
loc.fun
xlab
ylab
PB
nboot
SEED

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, y = NULL, tr = 0.2, alpha = 0.05, plotit = FALSE, 
    plotfun = splot, op = TRUE, VL = TRUE, cor.op = FALSE, loc.fun = median, 
    xlab = "Groups", ylab = "", PB = FALSE, nboot = 100, SEED = TRUE) 
{
    if (tr == 0.5) 
        stop("Use medpb to compare medians.")
    if (tr > 0.5) 
        stop("Can't have tr>.5")
    if (is.null(y)) {
        if (is.matrix(x) || is.data.frame(x)) {
            y = x[, 2]
            x = x[, 1]
        }
        if (is.list(x)) {
            y = x[[2]]
            x = x[[1]]
        }
    }
    library(MASS)
    if (SEED) 
        set.seed(2)
    x <- x[!is.na(x)]
    y <- y[!is.na(y)]
    n1 = length(x)
    n2 = length(y)
    h1 <- length(x) - 2 * floor(tr * length(x))
    h2 <- length(y) - 2 * floor(tr * length(y))
    q1 <- (length(x) - 1) * winvar(x, tr)/(h1 * (h1 - 1))
    q2 <- (length(y) - 1) * winvar(y, tr)/(h2 * (h2 - 1))
    df <- (q1 + q2)^2/((q1^2/(h1 - 1)) + (q2^2/(h2 - 1)))
    crit <- qt(1 - alpha/2, df)
    m1 = mean(x, tr)
    m2 = mean(y, tr)
    mbar = (m1 + m2)/2
    dif = m1 - m2
    low <- dif - crit * sqrt(q1 + q2)
    up <- dif + crit * sqrt(q1 + q2)
    test <- abs(dif/sqrt(q1 + q2))
    yuen <- 2 * (1 - pt(test, df))
    xx = c(rep(1, length(x)), rep(2, length(y)))
    if (h1 == h2) {
        pts = c(x, y)
        top = var(c(m1, m2))
        if (!PB) {
            if (tr == 0) 
                cterm = 1
            if (tr > 0) 
                cterm = area(dnormvar, qnorm(tr), qnorm(1 - tr)) + 
                  2 * (qnorm(tr)^2) * tr
            bot = winvar(pts, tr = tr)/cterm
            e.pow = top/bot
            if (e.pow > 1) {
                x0 = c(rep(1, length(x)), rep(2, length(y)))
                y0 = c(x, y)
                e.pow = wincor(x0, y0, tr = tr)$cor^2
            }
        }
        if (PB) {
            bot = pbvar(pts)
            e.pow = top/bot
        }
    }
    if (n1 != n2) {
        N = min(c(n1, n2))
        vals = 0
        for (i in 1:nboot) vals[i] = yuen.effect(sample(x, N), 
            sample(y, N), tr = tr)$Var.Explained
        e.pow = loc.fun(vals)
    }
    if (plotit) {
        plot(xx, pts, xlab = xlab, ylab = ylab)
        if (op) 
            points(c(1, 2), c(m1, m2))
        if (VL) 
            lines(c(1, 2), c(m1, m2))
    }
    list(ci = c(low, up), n1 = n1, n2 = n2, p.value = yuen, dif = dif, 
        se = sqrt(q1 + q2), teststat = test, crit = crit, df = df, 
        Var.Explained = e.pow, Effect.Size = sqrt(e.pow))
  }

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