bin_permute_all: bin_permute_all

In LyonsZhang/ProgPermute: Progressive permutation for a dynamic representation of the robustness of microbiome discoveries

Usage

bin_permute_all(variable, testdata, top_pm, zoomn, sequen = NULL, alpha = 0.05, psigalpha = 0.05, method = "kruskal", goodtype = "goodpvname", initial_fun = bin_true_initial, test_fun = bin_multitest, exchange_fun = bin_exchange, paral = "no")

Arguments

variable
testdata
top_pm
zoomn
sequen
alpha
psigalpha
method
goodtype
initial_fun
test_fun
exchange_fun
paral

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 (variable, testdata, top_pm, zoomn, sequen = NULL, alpha = 0.05, 
    psigalpha = 0.05, method = "kruskal", goodtype = "goodpvname", 
    initial_fun = bin_true_initial, test_fun = bin_multitest, 
    exchange_fun = bin_exchange, paral = "no") 
{
    trueinitial <- initial_fun(variable = variable, testdata = testdata, 
        top_pm = top_pm, method = method, goodtype = goodtype)
    select <- trueinitial$select
    selname <- trueinitial$selname
    position0 <- trueinitial$position0
    position1 <- trueinitial$position1
    levs <- trueinitial$levs
    n <- trueinitial$n
    N <- trueinitial$N
    testlevel <- trueinitial$testlevel
    if (is.null(sequen)) {
        sequen <- seq(0, N - 1, 1)
    }
    require(doParallel)
    ncores <- parallel::detectCores()
    cl <- parallel::makeCluster(ncores)
    doParallel::registerDoParallel(cl)
    out <- {
        foreach(k = sequen) %dopar% exchange_fun(k, top_pm = top_pm, 
            n = n, zoomn = zoomn, variable = variable, select = select, 
            selname = selname, testdata = testdata, testlevel = testlevel, 
            levs = levs, position0 = position0, position1 = position1, 
            alpha = alpha, psigalpha = psigalpha, method = method, 
            test_fun = test_fun, paral = paral)
    }
    parallel::stopCluster(cl)
    out <- do.call(Map, c(f = rbind, out))
    return(list(alldist = out, selname = selname))
  }

LyonsZhang/ProgPermute documentation built on July 22, 2020, 4:09 p.m.