Perhaps my understanding is wildly misguided, but I'm seriously confused about why people are still interested using FPKM values for cross sample gene expression analysis. My understanding is, that for some gene ${ i }$, its $FPKM$ is calculated as follows:

$$ FPKM_i \propto \frac{ \mathrm{fragments\ mapping \ to\ gene \ }i } { \sum_j \mathrm{ fragments\ mapping\ to \ gene\ } j } \bigg/ \mathrm{effective\_length}(i) $$

For comparing gene expression within a single sample, this doesn't seem too harmful to me. In fact, this seems like a good idea.

However, the use of $FPKM$ or even the famous $log_2(FPKM+1)$ in a similarity analysis seems to me plain wrong. True differences in RNA abundances of a couple of genes would significantly impact the FPKM values of the other genes leaving perhaps confused about the true differences between the samples.

Am I missing something here? Is it ever meaningful to use FPKM values for analyzing across samples?


Is it ever meaningful to use FPKM values for analyzing across samples?

One should never use FPKMs for anything important. They can occasionally be useful for plotting, but even in that case one needs to construct the FPKMs from properly normalized data and not use the original definition.

It's becoming ever rarer to see people using FPKMs these days, thankfully.

| improve this answer | |
  • $\begingroup$ Thanks, I also saw your response here biostars.org/p/317417 as the answer to my followup. $\endgroup$ – jaslibra Jul 3 '18 at 19:06
  • $\begingroup$ So what is a proper alternative? TMM? $\endgroup$ – LinuxBlanket Jul 10 '18 at 12:39
  • $\begingroup$ @LinuxBlanket TMM or another robust normalization method, yes. $\endgroup$ – Devon Ryan Jul 10 '18 at 12:40

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.