Default function for mapping expected counts to dispersion.

The function was estimated from circadian RNA-seq data from mouse liver (PRJNA297287), using local regression in DESeq2. In a negative binomial distribution, \(variance = mean + mean^2 * dispersion\).

Usage

defaultDispFunc(x)

Format

An object of class function of length 1.

Arguments

x

Numeric vector of mean counts.

Value

Numeric vector of dispersions.

See also

simphony()

Examples

means = 2^(6:10)
dispersions = defaultDispFunc(means)