MBL.p.adjust: The adjusted p-values for Modified Benjamini-Liu (BL)...

In allenzhuaz/MHTdiscrete: Multiple Hypotheses Testing for Discrete Data

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 α

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)

allenzhuaz/MHTdiscrete documentation built on July 12, 2019, 2:27 p.m.