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.