mutoss (version 0.1-13)

BY: Benjamini-Yekutieli (2001) step-up procedure

Description

The Benjamini-Yekutieli step-up procedure is applied to pValues. The procedure ensures FDR control for any dependency structure.

Usage

BY(pValues, alpha, silent=FALSE)

Value

A list containing:

adjPValues

A numeric vector containing the adjusted pValues

criticalValues

A numeric vector containing critical values used in the step-up-down test

rejected

A logical vector indicating which hypotheses are rejected

errorControl

A Mutoss S4 class of type errorControl, containing the type of error controlled by the function and the level alpha.

Author

WerftWiebke

Arguments

pValues

The used unadjusted pValues.

alpha

The level at which the FDR shall be controlled.

silent

If true any output on the console will be suppressed.

Details

The critical values of the Benjamini-Yekutieli (BY) procedure are calculated by replacing the alpha of the Benjamini-Hochberg procedure by alpha/sum(1/1:m)), i.e., c(i)=i*alpha/(m*(sum(1/1:m))) for i=1,...,m. For large number m of hypotheses the critical values of the BY procedure and the BH procedure differ by a factor log(m). Benjamini and Yekutieli (2001) showed that this step-up procedure controls the FDR at level alpha*m/m0 for any dependency structure among the test statistics.

References

Benjamini, Y. and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Annals of Statistics, 29(4):1165-1188.

Examples

alpha <- 0.05
p <-c(runif(10, min=0, max=0.01), runif(10, min=0.9, max=1))
result <- BY(p, alpha)
result <- BY(p, alpha, silent=TRUE)