MBL.p.adjust {MHTdiscrete}

The adjusted p-values for Modified Benjamini-Liu (BL) step-down FDR controlling procedure.

Description

The function for calculating the adjusted p-values based on original available p-values and all attaianble p-values.

Usage

MBL.p.adjust(p, p.set, alpha, make.decision)

Arguments

p

numeric vector of p-values (possibly with NAs). Any other R is coerced by as.numeric. Same as in p.adjust.

p.set

a list of numeric vectors, where each vector is the vector of all attainable p-values containing the available p-value for the corresponding hypothesis.

alpha

significant level used to compare with adjusted p-values to make decisions, the default value is 0.05.

make.decision

logical; if TRUE, then the output include the decision rules compared adjusted p-values with significant level \alpha

Value

A numeric vector of the adjusted p-values (of the same length as p).

Note

The MBL procedure for discrete data controls FDR under the specific dependence assumption where the joint distribution of statistics from true nulls are independent of the joint distribution of statistics from false nulls.

Author(s)

Yalin Zhu

References

Benjamini, Y., and Liu, W. (1999). A step-down multiple hypotheses testing procedure that controls the false discovery rate under independence. Journal of Statistical Planning and Inference, 82: 163-170.

See Also

MBH.p.adjust, MBY.p.adjust

Examples

p <- c(pbinom(1,8,0.5),pbinom(1,5,0.75),pbinom(1,6,0.6))
p.set <-list(pbinom(0:8,8,0.5),pbinom(0:5,5,0.75),pbinom(0:6,6,0.6))
MBL.p.adjust(p,p.set)