FDR / q-values / tail-area p-values via fdrtool — morie_fdr

Thin extender over fdrtool::fdrtool for the Strimmer (2008) shrinkage estimator of local and tail-area false-discovery rates, q-values, and the underlying null-distribution scale parameter.

Usage

morie_fdr_qvalues(x, statistic = "normal", ...)

Arguments

x: Numeric vector of test statistics or p-values, as appropriate to statistic.
statistic: Character; the type of statistic in x. One of "normal" (default), "correlation", "pvalue", or "studentt".
...: Further arguments forwarded to fdrtool::fdrtool (e.g. plot, color.figure, verbose, cutoff.method, pct0).

Value

A list with $method = "fdrtool::fdrtool" and $raw (the fdrtool return list with pval, qval, lfdr, and param).

Examples

if (FALSE) { # \dontrun{
  if (requireNamespace("fdrtool", quietly = TRUE)) {
    set.seed(1)
    x <- c(stats::rnorm(900), stats::rnorm(100, mean = 3))
    morie_fdr_qvalues(x, statistic = "normal")
  }
} # }