hierasymptkappa: some title

In WIAS-BERLIN/hierarchicalFDR: What the package does (short line)

Usage

hierasymptkappa(pValues, alpha = 0.05, kappa = 0.001)

Arguments

pValues
alpha
kappa

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 (pValues, alpha = 0.05, kappa = 0.001) 
{
    result <- list()
    if (is.list(pValues) == FALSE) {
        fpValues <- sort(pValues)
        m <- length(pValues)
        u <- as.integer(kappa * m) + 1
        conjpValues <- numeric()
        for (i in 1:(m - u + 1)) {
            locpvalue <- (fpValues[[(u - 1 + i)]] * (m - u + 
                1)/i)
            conjpValues <- c(conjpValues, locpvalue)
        }
        fpvalue <- min(conjpValues)
        if (fpvalue <= alpha * kappa) {
            hresult <- aorc(pValues, alpha, startIDX_SUD = u, 
                betaAdjustment = 0, silent = TRUE)
            result <- list(hresult)
        }
        else {
            result <- (list(fpValues != fpValues))
        }
    }
    else {
        k <- length(pValues)
        if ((1/k) <= kappa) {
            kappa <- (1/k)
        }
        for (i in 1:k) {
            if (length(pValues[[i]]) == 0) {
                lresult <- list(FALSE)
                result <- c(result, lresult)
                next
            }
            fpValues <- pValues[[i]]
            sortfpValues <- sort(fpValues)
            fm <- length(fpValues)
            fu <- as.integer(kappa * fm)
            if (fu < (kappa * fm)) {
                fu <- fu + 1
            }
            conjpValues <- numeric()
            for (i in 1:(fm - fu + 1)) {
                locpvalue <- (sortfpValues[[(fu - 1 + i)]] * 
                  (fm - fu + 1)/i)
                conjpValues <- c(conjpValues, locpvalue)
            }
            fpvalue <- min(conjpValues)
            if (fpvalue <= (alpha * kappa)) {
                lresult <- aorc(fpValues, alpha, startIDX_SUD = fu, 
                  betaAdjustment = 0, silent = TRUE)
                lresult <- list(lresult$rejected)
            }
            else {
                lresult <- list(fpValues != fpValues)
            }
            result <- c(result, lresult)
        }
    }
    return(result)
  }

WIAS-BERLIN/hierarchicalFDR documentation built on Dec. 30, 2019, 11:49 p.m.