logBH: BH correction on log-p-values

In MarioniLab/scran: Methods for Single-Cell RNA-Seq Data Analysis

BH correction on log-p-values

Description

Perform a Benjamini-Hochberg correction on log-transformed p-values to get log-adjusted p-values, without the loss of precision from undoing and redoing the log-transformations.

Usage

.logBH(log.p.val)

Arguments

log.p.val

Numeric vector of log-transformed p-values.

Value

A numeric vector of the same length as log.p.val containing log-transformed BH-corrected p-values.

Author(s)

Aaron Lun

Examples

log.p.values <- log(runif(1000))
obs <- .logBH(log.p.values)
head(obs)

ref <- log(p.adjust(exp(log.p.values), method="BH"))
head(ref)

MarioniLab/scran documentation built on Sept. 7, 2024, 6:25 a.m.