Applications of TPM normalisation
1
0
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.
#'
#' https://support.bioconductor.org/p/91218/
#'
tpm <- function(counts,len) {
  x <- counts/len
  return(t(t(x)*1e6/colSums(x)))
}
RNA-Seq deseq2 normalisation dge • 799 views
ADD COMMENT
0
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 : https://support.bioconductor.org/p/91218/

Your function seems ok

See also : https://youtu.be/TTUrtCY2k-w

ADD COMMENT

Login before adding your answer.

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