Question: Applications of TPM normalisation
0
gravatar for rmf
15 months ago by
rmf950
rmf950 wrote:

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)))
}
ADD COMMENTlink modified 15 months ago • written 15 months ago by rmf950
0
gravatar for Bastien Hervé
15 months ago by
Bastien Hervé4.7k
Karolinska Institutet, Sweden
Bastien Hervé4.7k wrote:

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 COMMENTlink modified 15 months ago • written 15 months ago by Bastien Hervé4.7k
Please log in to add an answer.

Help
Access

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