Learn R Programming

mutoss (version 0.1-12)

oracleBH: Bejamini-Hochberg (2000) oracle linear step-up Procedure...

Description

Bejamini-Hochberg (2000) oracle linear step-up Procedure

Usage

oracleBH(pValues, alpha, pi0, silent=FALSE)

Value

A list containing:

adjPValues

A numeric vector containing the adjusted pValues

rejected

A logical vector indicating which hypotheses are rejected

criticalValues

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

errorControl

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

Author

HackNiklas

Arguments

pValues

pValues to be used

alpha

the level at which the FWER should be controlled

pi0

miraculousy known number of true null hypotheses

silent

Logical, if true any output on the console will be suppressed.

Details

Knowledge of the number of true null hypotheses (m0) can be very useful to improve upon the performance of the FDR controlling procedure. For the oracle linear step-up procedure we assume that m0 were given to us by an 'oracle', the linear step-up procedure with q0 = q*m/m0 would control the FDR at precisely the desired level q in the independent and continuous case, and would then be more powerful in rejecting hypotheses for which the alternative holds.

Examples

Run this code
pval <- c(runif(50), runif(50, 0, 0.01))
result <- oracleBH(pValues=pval,alpha=0.05,pi0=0.85)

Run the code above in your browser using DataLab