ancGLOB_pv:

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

Usage

ancGLOB_pv(n1, n2, est = tmean, fr1 = 0.8, fr2 = 0.8, nboot = 500, SEED = TRUE, iter = 1000, nmin = 12, MC = TRUE, alpha = 0.05, PRM = FALSE, pts = NULL, ...)

Arguments

n1
n2
est
fr1
fr2
nboot
SEED
iter
nmin
MC
alpha
PRM
pts
...

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 (n1, n2, est = tmean, fr1 = 0.8, fr2 = 0.8, nboot = 500, 
    SEED = TRUE, iter = 1000, nmin = 12, MC = TRUE, alpha = 0.05, 
    PRM = FALSE, pts = NULL, ...) 
{
    if (SEED) 
        set.seed(45)
    bvec = list()
    np1 = min(c(n1, n2)) + 1
    nmax = max(c(n1, n2))
    for (i in 1:iter) {
        bvec[[i]] = rmul(nmax, p = 4)
        if (n1 != n2) 
            bvec[[i]][np1:nmax, 1:2] = NA
    }
    if (MC) {
        library(parallel)
        prm = mclapply(bvec, ancGLOB_sub2, fr1 = fr1, fr2 = fr2, 
            est = est, SEED = SEED, nboot = nboot, pts = pts, 
            ...)
    }
    if (!MC) 
        prm = lapply(bvec, ancGLOB_sub2, fr1 = fr1, fr2 = fr2, 
            est = est, SEED = SEED, nboot = nboot, pts = pts, 
            ...)
    prm = elimna(as.vector(matl(prm)))
    ef.iter = length(prm)
    p.crit = hd(prm, alpha)
    prm = sort(elimna(prm))
    if (!PRM) 
        prm = NULL
    list(p.crit = p.crit, prm = prm, ef.iter = ef.iter)
  }

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