BonEV (version 1.0)

BonEV-package: BonEV: An Improved Multiple Testing Procedure for Controlling False Discovery Rates

Description

BonEV is an improved multiple testing procedure for controlling false discovery rates which is developed based on the Bonferroni procedure with integrated estimates from the Benjamini-Hochberg procedure and the Storey's q-value procedure. It controls false discovery rates through controlling the expected number of false discoveries.

Details

Package:

BonEV

Type:

Package

Version:

1.0.0

Date:

2015-02-10

Depends:

R (>= 3.2.0), qvalue

License:

GPL (>= 2)

See Also

The Bon_EV function defined in this package. The qvalue package.

Examples

library(qvalue)
data(hedenfalk)
summary(hedenfalk)
pvalues <- hedenfalk$p
adjp <- Bon_EV(pvalues, 0.05)
summary(adjp)
results <- cbind(adjp$raw_P_value, adjp$BH_adjp, adjp$Storey_adjp, adjp$Bon_EV_adjp)
results

##Compare with Benjamini-Hochberg and Storey's q-value procedures
sum(adjp$raw_P_value <= 0.05)
sum(adjp$BH_adjp <= 0.05)
sum(adjp$Storey_adjp <= 0.05)
sum(adjp$Bon_EV_adjp <= 0.05)