Filtering during RNA Seq analysis?

Question: Filtering during RNA Seq analysis?

22 months ago by

Bonn, Germany

Hi,

I have roughly 100 samples (2 groups) of noisy human tissue data. Thinking about filtering to improve the multiple testing situation. Thus far, I filter by cpm values (calculated using edgeR after counting with featurecounts), and run voom/limma afterwards. I have also stumbled over htsfilter, and thought about variance filtering.

Now I heard that filtering, at least by variance or htsfilter, was per se incompatible with limma or other empirial Bayes methods.

What do you experts think/suggest?

Many thanks!

rna-seq • 885 views

ADD COMMENT • link •

modified 22 months ago by i.sudbery • 2.3k • written 22 months ago by aggregatibacter • 90

22 months ago by

i.sudbery • 2.3k

Sheffield, UK

i.sudbery • 2.3k wrote:

You might also try "Independent Hypothesis Weighting" as implemented in the IHW Bioconductor package and described it this paper: http://www.ncbi.nlm.nih.gov/pubmed/27240256

Instead of filtering genes to conserve power, it weights them. It is compatible at least with DESeq2 and is used in their vignette (i suspect it can be used in conjunction with other packages as well).

ADD COMMENT • link written 22 months ago by i.sudbery • 2.3k

22 months ago by

Devon Ryan ♦ 81k

Freiburg, Germany

Devon Ryan ♦ 81k wrote:

The package you're looking for is genefilter. As an aside, that's used automatically by DESeq2.

BTW, filtering by variance isn't compatible with any RNAseq package that I'm aware of (variance is convoluted with the DE test).

ADD COMMENT • link written 22 months ago by Devon Ryan ♦ 81k

Thanks guys for your suggestions. I will definitely also try DESeq2. I have the feeling that it is quite en vogue at the moment, is that correct? Thus far, I primarily use limma, because I was brought up with arrays, and usually need to correct for random and fixed factors.

ADD REPLY • link written 22 months ago by aggregatibacter • 90

Please log in to add an answer.

Similar posts • Search »