Applications of TPM normalisation
Entering edit mode
23 months ago
rmf ★ 1.1k

I use raw counts for DGE in DESeq2, edgeR etc. And I use VST for PCA and clustering.

Now, I am curious about TPM. What are the applications of a TPM normalised count table? What can it be used for? What are the drawbacks of TPM normalised counts?

enter image description here enter image description here enter image description here

I hope my TPM calculation is correct. Here is what I use.

#' @title Compute TPM from a read count matrix
#' @param counts A numeric data.frame of read counts with samples (columns) and genes (rows).
#' @param len A vector of gene cds length equal to number of rows of dfr.
tpm <- function(counts,len) {
  x <- counts/len
RNA-Seq deseq2 normalisation dge • 799 views
Entering edit mode
23 months ago

You can use TPM to compare gene count in the same sample, but it is not the best way to go to compare genes across sample, prefer RLE or TMM method for example.

This is one way to convert to TPM :

Your function seems ok

See also :


Login before adding your answer.

Traffic: 1872 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6