Storey q-value procedure (adaptive FDR)

Estimates the proportion of true null hypotheses (pi0) and tightens the BH thresholds by that factor. Delegates to qvalue::qvalue (Bioconductor) when installed; otherwise falls back to an inline Storey-style cutoff so the wrapper keeps working on CRAN-only installs.

Usage

storey_q(p_values, alpha = 0.05, lambda_param = 0.5, labels = NULL)

Arguments

p_values

Numeric vector of raw p-values.

alpha

Significance level.

lambda_param

Tuning parameter in (0, 1) for the pi0 estimator.

labels

Optional character vector of test labels.