EdgeR cpm values after batch effects

[danielfortin86]

Is there a way to obtain the CPM values from EdgeR after batch correction? The "cpm" function only returns values after library size normalization. When I look at the result of glmFit there are "abundance" "fitted.values" and "counts". The description in the documentation (?glmFit) is not entirely clear to me. Are the "fitted.values" cpm values after batch correction or should I be doing something else?

Thanks!

[Gordon Smyth]

The answer is probably yes, but this is quite a technical question that can probably only be answered by the package authors, so it would best be sent to the main edgeR help list. See "how to get help" in the edgeR User's Guide.

[gfmgfm]

CPM after batch effect correction in edgeR

Hello,

I wonder if you figured out how to print CPM after batch effect correction. I am trying to do exactly the same.

Many thanks.

[Gordon Smyth]

You could try:

library(edgeR)
logCPM <- cpm(y, log=TRUE, prior.count=3)
logCPM <- removeBatchEffect(y, batch=batch)

Otherwise, follow the advice given in my earlier post.

[gfmgfm]

Thank you. This is very helpful!