Methods for Function withinLaneNormalization in Package EDASeq
Within-lane normalization for GC-content (or other lane-specific) bias.
A numeric matrix representing the counts or a
A numeric vector representing the covariate to normalize for (if
Method used to normalized. See the details section and the reference below for details.
Should the normalized value be returned as an offset leaving the original counts unchanged?
The number of bins used to stratify the covariate for
If TRUE the normalization returns rounded values (pseudo-counts). Ignored if offset=TRUE.
This method implements four normalizations described in Risso et al. (2011).
loess normalization transforms the data by regressing the counts on
y and subtracting the loess fit from the counts to remove the dependence.
full normalizations are based on the stratification of the genes based on
y. Once the genes are stratified in
num.bins strata, the methods work as follows.
scales the data to have the same median in each bin.
the same but with the upper quartile.
forces the distribution of each stratum to be the same using a non linear full quantile normalization, in the spirit of the one used in microarrays.
signature(x = "matrix", y = "numeric")
It returns a matrix with the normalized counts if
offset=FALSEor with the offset if
signature(x = "SeqExpressionSet", y = "character")
It returns a
SeqExpressionSetwith the normalized counts in the
normalizedCountsslot and with the offset in the
D. Risso, K. Schwartz, G. Sherlock and S. Dudoit (2011). GC-Content Normalization for RNA-Seq Data. Manuscript in Preparation.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
library(yeastRNASeq) data(geneLevelData) data(yeastGC) sub <- intersect(rownames(geneLevelData), names(yeastGC)) mat <- as.matrix(geneLevelData[sub, ]) data <- newSeqExpressionSet(mat, phenoData=AnnotatedDataFrame( data.frame(conditions=factor(c("mut", "mut", "wt", "wt")), row.names=colnames(geneLevelData))), featureData=AnnotatedDataFrame(data.frame(gc=yeastGC[sub]))) norm <- withinLaneNormalization(data, "gc", which="full", offset=FALSE)
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.