Hello,
I am working with RNA-Seq data. I have the read counts for each, and before checking for differential expression, I want to normalize it using FPKM. In theory, I understand the difference between RPKM and FPKM, but could someone please tell me a formula to calculate the same?
Thank you.
•
written
by
vinayjrao
Although I would use DESeq2, Cuffdiff or other published and well established methods for finding DEGs. Here is a link where they explained the formula in details: RPKM, FPKM and TPM, clearly explained
Dear Rahul,
I have gone through the explanation provided on the link, but it still doesn't explain how you can relate read counts to fragment count. Is there any other way? As you said I could always use the tools, but I want to understand the method of calculating fpkm.
Rahul's point is if you want to do differential expression, let the software designed to do that handle the normalization. It does a better, more sophisticated method than using RPKM. You can probably get whatever software you use to give you the results of its normalization.
Dear Vinay,
once you've figured out the difference between read counts and fragments count, please share :) I need to use FPKM instead of RPKM since the experiment used was paired-end reads method.
Thank you.
Log in to answer this question.
An update (6th October 2018):
You should abandon RPKM / FPKM. They are not ideal where cross-sample differential expression analysis is your aim; indeed, they render samples incomparable via differential expression analysis:
Please read this: A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis
Also, by Harold Pimental: What the FPKM? A review of RNA-Seq expression units