This function (re)implements the empirical bayes shrinkage estimate of Smyth (2004), which is also implemented in the Limma package. This function is normally called from the function variance_adjust
, and is not normally intended for stand-alone use.
sigmashrink(s2, d)
"Standard" estimates of sigma^2
"Standard" degrees of freedom of the residuals
A list containing
Estimates of sigma^2 using the empirical bayes shrinkage method of Smyth (2004)
Estimate of degrees of freedom using the empirical bayes shrinkage method of Smyth (2004)
Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Smyth, 2004.
Using control genes to correct for unwanted variation in microarray data. Gagnon-Bartsch and Speed, 2012. Available at: http://biostatistics.oxfordjournals.org/content/13/3/539.full.
Removing Unwanted Variation from High Dimensional Data with Negative Controls. Gagnon-Bartsch, Jacob, and Speed, 2013. Available at: http://statistics.berkeley.edu/tech-reports/820.