Question: deseq2 pre filter rna seq
0
gravatar for badribio
13 days ago by
badribio0
badribio0 wrote:

I see more and more papers these days that use a read cut off before deseq2 analysis, I wonder what is the rational and/or how do they arrive at how many reads is a right number..for example from one of the papers DESeq2 was applied to genes having more than 20 reads in half the number of samples ^ this study had n=8 4wt and 4ko. i could not comprehend how did they arrive at the number 20, even if they did what if the differences are genotype driven more over DESeq2 has independent filtering which takes care of low counts no?

Is there a way one can plot counts and determine what is a right threshold for pre filtering

Thanks

rnaseq deseq2 • 73 views
ADD COMMENTlink modified 13 days ago by zefrieira0 • written 13 days ago by badribio0

You'll likely get better answers if you employ some punctuation.

ADD REPLYlink written 13 days ago by swbarnes2220

excuse me, hope i did ok this time.

ADD REPLYlink written 13 days ago by badribio0
Answer: deseq2 pre filter rna seq
0
gravatar for zefrieira
13 days ago by
zefrieira0
zefrieira0 wrote:

DESeq2 filters genes based on their counts. Read the independent filtering section of the documentation.

ADD COMMENTlink written 13 days ago by zefrieira0

Correct. We do sometime perform minimal pre-filtering. Mostly this is to reduce unnecessary computation, whereas the independent filtering (or better yet, IHW) is for increasing power. There is no need to run DESeq2 for bulk RNA-seq on genes where only one or two samples have a single digit count and the rest have all 0's.

ADD REPLYlink written 13 days ago by Michael Love24k
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 16.09
Traffic: 335 users visited in the last hour