Filtering during RNA Seq analysis?

Question: Filtering during RNA Seq analysis?

4.4 years 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 • 1.8k views

ADD COMMENT • link •

modified 4.4 years ago by i.sudbery ♦ 10k • written 4.4 years ago by aggregatibacter • 140

4.4 years ago by

i.sudbery ♦ 10k

Sheffield, UK

i.sudbery ♦ 10k 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 4.4 years ago by i.sudbery ♦ 10k

4.4 years ago by

Devon Ryan ♦ 98k

Freiburg, Germany

Devon Ryan ♦ 98k 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 4.4 years ago by Devon Ryan ♦ 98k

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 4.4 years ago by aggregatibacter • 140

Please log in to add an answer.

Similar posts • Search »