Two-stage FDR control for basin-specific differential abundance scans

Two-stage FDR control for basin-specific differential abundance scans

Usage

two.stage.basin.fdr(
  pval.mat,
  q = 0.1,
  basin.global.method = c("simes", "fisher")
)

Arguments

pval.mat

Numeric matrix of p-values with rows = taxa and cols = basins.

q

Target FDR level (e.g., 0.10).

basin.global.method

Method to aggregate basin-level p-values. Currently "simes" or "fisher".

Value

List with selected.basins, basin.pvals, within.q, and within.adj.p (list).