How to do DEGs based on CPM
1
0
Entering edit mode
li.zhaoz ▴ 10
@lizhaoz-22033
Last seen 3.6 years ago

I only have CPM value available. Just wondering if I could take the logCPM and put it input limma with trend = TRUE. That is what the voom paper descirbed if I get it right. Any advice? Thank you!

limma • 2.1k views
ADD COMMENT
5
Entering edit mode
@gordon-smyth
Last seen 11 hours ago
WEHI, Melbourne, Australia

See limma User's Guide, Section 15.4.

ADD COMMENT
2
Entering edit mode

Dear Gordon,

shouldn't the last chunk of code in Section 15.4 be (notice the trend=TRUE parm)?

Or, to give more weight to fold-changes in the gene ranking, one might use:

> fit <- lmFit(logCPM, design)
> fit <- treat(fit, lfc=log2(1.2), trend=TRUE)
> topTreat(fit, coef=ncol(design))

ADD REPLY
0
Entering edit mode

Yes, that is true. I will add the trend = TRUE.

ADD REPLY

Login before adding your answer.

Traffic: 853 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