Thank you so much for the clarification.
• 0 views
•
link
written
by
camelest •
When is the correct timing to perform edgeR cpm calculation?
Hi, I have a question on edgeR cpm function. If we follow the user manual, the standard workflow should be
x <- read.delim("TableOfCounts.txt",row.names="Symbol")
group <- factor(c(1,1,2,2))
y <- DGEList(counts=x,group=group)
keep <- filterByExpr(y)
y <- y[keep,,keep.lib.sizes=FALSE]
y <- calcNormFactors(y) *1
design <- model.matrix(~group)
y <- estimateDisp(y,design) *2
When we want to get the count table of log2cpm by cpm() function, for example, for clustering or heatmap, as in
logcpm <- cpm(y, log=TRUE)
should it happen at point 1 just after calculating NormFactors or at point 2 after estimation of dispersion? Is it OK as long as it is after calcNormFactors()?
• 1,383 views
•
link
written
by
camelest •
1 answer
Yes, anytime after calcNormFactors(). Dispersion estimates don't effect the logCPM.
• 0 views
•
link
written
by
Gordon Smyth
Log in to answer this question.
More posts like this
-
How do I work out which condition the genes are down-regulated in?
written by Tota •
-
EdgeR analysis with CPM normalzed counts
written by silas008
-
Problems in differential expression analysis with edgeR
written by thancy •
-
Why 'filterByExpr' does not eliminate low read counts from the expression data?
written by svp
-
MDS Plot R
written by BioBaby •
-
How to perform multiple comparisons using edgeR?
written by silas008
-
edger differential expression analysis error
written by Janey •
-
R code for Differential Gene Expression
written by landrjos •
-
edgeR Factorial DE Design
written by hamindigoblue •
-
Using offsets in edgeR to correct GC content
written by Dawe
