DqdifMC:

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

Usage

DqdifMC(x, y = NULL, q = 0.25, nboot = 1000, plotit = TRUE, xlab = "Group 1 - Group 2", SEED = TRUE, alpha = 0.05)

Arguments

x
y
q
nboot
plotit
xlab
SEED
alpha

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, q = 0.25, nboot = 1000, plotit = TRUE, 
    xlab = "Group 1 - Group 2", SEED = TRUE, alpha = 0.05) 
{
    library(parallel)
    if (SEED) 
        set.seed(2)
    if (q >= 0.5) 
        stop("q should be less than .5")
    if (!is.null(y)) {
        xy = elimna(cbind(x, y))
        dif = xy[, 1] - xy[, 2]
    }
    if (is.null(y)) 
        dif = elimna(x)
    x = as.matrix(x)
    n = length(dif)
    if (plotit) 
        akerd(dif, xlab = xlab)
    bvec = NA
    data <- matrix(sample(n, size = n * nboot, replace = TRUE), 
        nrow = nboot)
    data = listm(t(data))
    bvec <- mclapply(data, difQMC_sub, dif, q, mc.preschedule = TRUE)
    bvec = matl(bvec)
    est1 = hd(dif, q = q)
    est2 = hd(dif, q = 1 - q)
    pv = mean(bvec < 0) + 0.5 * mean(bvec == 0)
    p = 2 * min(c(pv, 1 - pv))
    low <- round((alpha/2) * nboot) + 1
    up <- nboot - low
    sbvec = sort(bvec)
    ci = sbvec[low]
    ci[2] = sbvec[up]
    list(est.q = est1, est.1.minus.q = est2, conf.interval = ci, 
        p.value = p)
  }

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