Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/estimateGLMTrendedDisp.R
Estimates the abundance-dispersion trend by Cox-Reid approximate profile likelihood.
1 2 3 4 5 |
y |
a matrix of counts or a |
design |
numeric design matrix, as for |
method |
method (low-level function) used to estimated the trended dispersions.
Possible values are |
offset |
numeric scalar, vector or matrix giving the linear model offsets, as for |
AveLogCPM |
numeric vector giving average log2 counts per million for each gene. |
weights |
optional numeric matrix giving observation weights |
... |
other arguments are passed to lower-level functions |
Estimates the dispersion parameter for each gene with a trend that depends on the overall level of expression for that gene. This is done for a DGE dataset for general experimental designs by using Cox-Reid approximate conditional inference for a negative binomial generalized linear model for each gene with the unadjusted counts and design matrix provided.
The function provides an object-orientated interface to lower-level functions.
When the input object is a DGEList
, estimateGLMTrendedDisp
produces a DGEList
object, which contains the estimates of the trended dispersion parameter for the negative binomial model according to the method applied.
When the input object is a numeric matrix, it returns a vector of trended dispersion estimates calculated by one of the lower-level functions dispBinTrend
, dispCoxReidPowerTrend
and dispCoxReidSplineTrend
.
Gordon Smyth, Davis McCarthy, Yunshun Chen
Cox, DR, and Reid, N (1987). Parameter orthogonality and approximate conditional inference. Journal of the Royal Statistical Society Series B 49, 1-39.
McCarthy, DJ, Chen, Y, Smyth, GK (2012). Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Research 40, 4288-4297. https://doi.org/10.1093/nar/gks042
dispBinTrend
, dispCoxReidPowerTrend
and dispCoxReidSplineTrend
for details on how the calculations are done.
1 2 3 4 5 6 7 |
Loading required package: limma
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.