CPM plot for individual genes
0
0
Entering edit mode
5.6 years ago
iraia.munoa ▴ 130

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 • 2.3k views
ADD COMMENT
1
Entering edit mode

You can get the TMM-corrected CPM values from your analysis, see output TMM normalized counts with edgeR. No need to start over.

ADD REPLY
0
Entering edit mode

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!

ADD REPLY
1
Entering edit mode

The latter is what I would do.

ADD REPLY
0
Entering edit mode

Well thanks!!! I will try!

ADD REPLY

Login before adding your answer.

Traffic: 1497 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6