scorciMC:

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

Usage

scorciMC(x, y, nboot = 1000, alpha = 0.05, SEED = TRUE, plotit = TRUE, STAND = TRUE, corfun = pcor, cop = 3, ...)

Arguments

x
y
nboot
alpha
SEED
plotit
STAND
corfun
cop
...

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, nboot = 1000, alpha = 0.05, SEED = TRUE, plotit = TRUE, 
    STAND = TRUE, corfun = pcor, cop = 3, ...) 
{
    library(parallel)
    if (ncol(as.matrix(x)) != 1) 
        stop("x should be a single vector")
    m1 = cbind(x, y)
    m1 <- elimna(m1)
    nval = nrow(m1)
    x <- m1[, 1]
    y <- m1[, 2]
    est <- scor(x, y, plotit = plotit, STAND = STAND, corfun = corfun, 
        cop = cop, ...)$cor
    if (SEED) 
        set.seed(2)
    data <- matrix(sample(length(y), size = length(y) * nboot, 
        replace = TRUE), nrow = nboot)
    data = listm(t(data))
    bvec <- mclapply(data, scorsubMC, x, y, STAND = STAND, corfun = corfun, 
        cop = cop, ...)
    bvec = matl(bvec)
    bvec = as.vector(bvec)
    ihi <- floor((1 - alpha/2) * nboot + 0.5)
    ilow <- floor((alpha/2) * nboot + 0.5)
    bsort <- sort(bvec)
    corci <- 1
    corci[1] <- bsort[ilow]
    corci[2] <- bsort[ihi]
    phat <- sum(bvec < 0)/nboot
    sig <- 2 * min(phat, 1 - phat)
    list(cor.ci = corci, p.value = sig, cor.est = est)
  }

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