Question: Calculating RPKM for DESeq
gravatar for Whoknows
5.7 years ago by
Whoknows800 wrote:

Hi friends 

 DESeq works via count table and create DE genes by count fold change, but it's better to know about gene RPKM values.

Is there any possible way to use / show RPKM for DESeq's output? 

And how can i normalized those RPKM values?


rna-seq deseq tophat htseq • 6.8k views
ADD COMMENTlink modified 5.7 years ago by Devon Ryan95k • written 5.7 years ago by Whoknows800

Why is it better to know about RPKM values?  

ADD REPLYlink written 5.7 years ago by Sean Davis26k

 I think  fold-change where generated based on the count table could be a good measurement, but sometimes read count is not practical (for example those  genes with 0 count in one condition) and you have to use another measurement like RPKM and FPKM , Using RPKM make study more robust and confident based on the A.Mortazavi 2008 paper that he said  you can detect with 95% confidence if a transcript has RPKM >11.By the way, using RPKM or FPKM let you to decide for choosing only those genes  with an exact threshold based on your coverage , and other sequencing criterias. 

ADD REPLYlink modified 5.7 years ago • written 5.7 years ago by Whoknows800
gravatar for Devon Ryan
5.7 years ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k wrote:

Do not mix DESeq(2) and RPKMs. RPKMs are inferior to count data for a variety of reasons. Rethink why you want these values to begin with.

If you want to convert count data to RPKMs, take the normalized counts produced by DESeq2 (or edgeR, or any similar tool), divide them by the gene length in KB (typically one would just sum the exons) and divide by 1 million.

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Devon Ryan95k


I want to use rpkm() function in edgeR to do this but i am facing problems with gene length do you know what am I might be doing wrong?

ADD REPLYlink modified 5 months ago by RamRS27k • written 4.5 years ago by kanika.15180
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1793 users visited in the last hour