Learn R Programming

MAnorm2 (version 1.2.2)

inv.trigamma: Inversion of Trigamma Function

Description

inv.trigamma implements the Newton iteration for solving, given x, the equation for y: trigamma(y) = x. See appendix of the limma paper (see "References") for a theoretical deduction of the method.

Usage

inv.trigamma(x, eps = 1e-08)

Value

The solution, which is also a positive numeric scalar.

Arguments

x

A positive numeric scalar.

eps

The required precision of the solution.

References

Smyth, G.K., Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol, 2004. 3: p. Article3.

See Also

trigamma for the trigamma function.

Examples

Run this code
x <- trigamma(1:6)
vapply(x, inv.trigamma, numeric(1))

Run the code above in your browser using DataLab