med.effect:

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

Usage

med.effect(x, y, HD = TRUE, eq.var = FALSE, nboot = 100, loc.fun = median, varfun = pbvar)

Arguments

x
y
HD
eq.var
nboot
loc.fun
varfun

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, HD = TRUE, eq.var = FALSE, nboot = 100, loc.fun = median, 
    varfun = pbvar) 
{
    x = elimna(x)
    y = elimna(y)
    if (HD) {
        e1 = hd(x)
        e2 = hd(y)
    }
    if (!HD) {
        e1 = median(x)
        e2 = median(y)
    }
    if (eq.var) {
        s1sq = varfun(x)
        s2sq = varfun(y)
        spsq <- (n1 - 1) * s1sq + (n2 - 1) * s2sq
        sp <- sqrt(spsq/(n1 + n2 - 2))
        dval = (e1 - e2)/sp
    }
    if (!eq.var) {
        n1 = length(x)
        n2 = length(x)
        if (n1 == n2) 
            dval = var(c(e1, e2))/varfun(c(x, y))
        if (n1 != n2) {
            N = min(c(n1, n2))
            vals = 0
            for (i in 1:nboot) vals[i] = med.effect.sub(sample(x, 
                N), sample(y, N), HD = HD, varfun = varfun)
            dval = loc.fun(vals)
        }
    }
    dval = sqrt(dval)
    dval
  }

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