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
•
written
by
Neu •
1 answer
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.
Log in to answer this question.
-
RPKM value interpretation
written by ymd56731 •
-
ChIP-Seq density plot between two groups of genes
written by hemantcnaik •
-
Calculate each value based on spearman correlation
written by Learner •
-
Meaningful DE using rpkm values in a RNA-Seq analysis
written by Seq225
-
Z-score for differentially expressed genes
written by Fra11 •
-
How to produce a unimodal distribution of genes expressed genes
written by mail2steff
-
Should I get the mean of 2 log2 transformed RPKM values?
written by BioinfGuru
-
Threshold for calling DNAse-seq peaks ON or OFF
written by Bioinformatist Newbie
-
Which method is good for getting logFC for pairwise comparsion of replicates?
written by mike
-
Rnaseq Differential Expression
written by Ashutosh Pandey
The standard procedure is not to use RPKM at all. There are much better ways to normalize read counts.