mcpPV:

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

Usage

mcpPV(pv, alpha = 0.05, opt = 1)

Arguments

pv
alpha
opt

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 (pv, alpha = 0.05, opt = 1) 
{
    K = length(pv)
    pick = NULL
    v = order(pv)
    ic = 0
    for (i in K:1) {
        flag = TRUE
        if (opt == 1) {
            i2 = i * 2
            if (i == K) 
                res = (0 - 2) * sum(log(pv))
            if (i < K) 
                res = (0 - 2) * sum(log(pv[-pick]))
            pvF = 1 - pchisq(res, i2)
        }
        if (opt == 2) {
            if (i == K) 
                res = sum(qnorm(pv/2)^2)
            if (i < K) 
                res = sum(qnorm(pv[-pick]/2)^2)
            pvF = 1 - pchisq(res, i)
        }
        if (opt == 3) {
            if (i == K) 
                res = max(pv)
            if (i < K) 
                res = max(pv[-pick])
            pvF = pbeta(res, i, 1)
        }
        if (pvF > alpha) 
            flag = TRUE
        if (pvF <= alpha/(K + 1 - i)) {
            ic = ic + 1
            pick = c(pick, v[ic])
            flag = FALSE
            if (pv[v[ic]] > alpha) 
                flag = TRUE
        }
        if (flag) 
            break
    }
    Decision = rep("Not Sig", length(pv))
    if (!is.null(pick)) 
        Decision[pick] = "Reject"
    nsig = sum(length(pick))
    list(p.values = pv, Decisions = as.matrix(Decision), num.sig = nsig)
  }

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