This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CPM plot for individual genes

Hi!

I am working with my DEGs list from edgeR with TMM normalization. I have seen in many publications that researches plot log2(FPKM) values for individual DEGs to show the statistical significance. Regarding this, I would like to know if it is posible to plot the log2(CPM) values of control and treatment samples for individual genes, as I have used edgeR with TMM and not RPKM. Is there a function in edgeR to obtain this values for control and treatment samples of each gene, instead of the total value of the final report? Or do I need to start with featureCounts data and transform it somehow? Can someone help me?

Thanks in advance!!

Iraia

rna-seq cpm edger tmm
• 3,093 views • link
updated by ATpoint ATpoint • written by iraia.munoa iraia.munoa

Well I have been reading the link you showed me, but I am a little bite confused and don't understand what I need to run in edgeR. In the coments James Ashmore provides a cpm.DGEList function, do I need to run it like that? Or running cpm(y, log=TRUE, lib.size=libsize, prior.count=2) after the calcNormFactor() selecting TMM method is enough? Thanks again!

• 0 views • link
updated by ATpoint ATpoint • written by iraia.munoa iraia.munoa

The latter is what I would do.

• 0 views • link
written by ATpoint ATpoint

Well thanks!!! I will try!

• 0 views • link
written by iraia.munoa iraia.munoa

0 answers

No answers yet.

Log in to answer this question.

More posts like this