Question: Filtering during RNA Seq analysis?
0
gravatar for aggregatibacter
22 months ago by
Bonn, Germany
aggregatibacter90 wrote:

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 COMMENTlink modified 22 months ago by i.sudbery2.3k • written 22 months ago by aggregatibacter90
4
gravatar for i.sudbery
22 months ago by
i.sudbery2.3k
Sheffield, UK
i.sudbery2.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 COMMENTlink written 22 months ago by i.sudbery2.3k
3
gravatar for Devon Ryan
22 months ago by
Devon Ryan81k
Freiburg, Germany
Devon Ryan81k 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 COMMENTlink written 22 months ago by Devon Ryan81k

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 REPLYlink written 22 months ago by aggregatibacter90
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: 978 users visited in the last hour