Last data update: 2014.03.03

R: Benjamini & Hochberg (1995) method for p-values correction
BHcorrectionR Documentation

Benjamini & Hochberg (1995) method for p-values correction

Description

This function evaluate the p-values correction according to the Benjamini & Hochberg (1995) method.

Usage

BHcorrection(pvs, alpha = 0.05)

Arguments

pvs

list of p-values

alpha

level of significance of the test

Value

list of corrected p-values

Author(s)

Stefano M. Pagnotta and Michele Ceccarelli

References

Benjamini, Y. and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. Roy. Statist. Soc. Ser. B 57 289-300.

Examples

###################
data <- matrix(rnorm(50000), 10000, 5)
pvs <- apply(data, 1, function(xx) t.test(xx, alternative = "less")$p.value)
qvs <- BHcorrection(pvs)

Results