Question: How to exercise RPKM cut-off
0
gravatar for Neu
2.4 years ago by
Neu10
India
Neu10 wrote:

Hi, I have 26 technical replicates. I have calculated RPKM values for each gene. Next, how to put the cut-off. Should I plot density plot of the log-trans RPKM value of all the genes and chose the peak as the cut off ? What is the standard procedure?

Thank you

sequencing rna-seq R • 1.0k views
ADD COMMENTlink modified 2.4 years ago by Macspider2.8k • written 2.4 years ago by Neu10
2

The standard procedure is not to use RPKM at all. There are much better ways to normalize read counts.

ADD REPLYlink written 2.4 years ago by Asaf5.6k
1
gravatar for Macspider
2.4 years ago by
Macspider2.8k
Vienna - BOKU
Macspider2.8k wrote:

So you have calculated RPKM for each gene for each sample. First, I would suggest you to use better normalization techniques, such as TPM ( http://www.rna-seqblog.com/rpkm-fpkm-and-tpm-clearly-explained/ ).

To set cutoffs, it always depends on your data. Making a density plot of the TPMs of each sample and comparing them will help you in understanding where the mean is, and therefore how much you lose with certain cutoffs. F.e. with FPKM a common practice was to select only FPKMs >= 1, with TPM as far as I know 1 works but requires a plot beforehand because it's a merely arbitrary number.

ADD COMMENTlink written 2.4 years ago by Macspider2.8k
1

Thank you so much. I plotted density of log2 values of TPM for each sample and selected 0.06 as the cut-off.

ADD REPLYlink written 2.4 years ago by Neu10
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: 712 users visited in the last hour