Question: Limma filtering based on treatment group
gravatar for jamie.gearing
2.1 years ago by
jamie.gearing40 wrote:

I have a question about filtering lowly expressed probes (or genes) when using limma.
For the example in section 17.4 of the limma user's guide, probes are called expressed if they exceed a cut-off in more than a given number of samples (equal to the size of the smallest treatment group):

isexpr <- rowSums(y$E > cutoff) >= 4​

I was wondering whether it would make sense to use the treatment group information to look for probes that are expressed in every sample in any treatment group (or perhaps a proportion of samples for any treatment). Something like this perhaps:

proportion <- 1.0
isexpr2 <- apply(y$E > cutoff, 1, function(z){
                           any(sapply(levels(Treatment), function(treat){
                             sum(z[Treatment == treat]) >= sum(Treatment == treat)*proportion

For this example, the normal approach yields 32754 expressed probes, whereas this yields 30840 probes.
I have seen other answers on this subject (e.g. warning against filtering based on variance because it would affect the limma algorithms, but I am not sure whether this is quite the same thing.
Would this be a technically valid? Even if it is, it may just be unnecessarily complicated.

limma filtering • 379 views
ADD COMMENTlink modified 2.1 years ago by Steve Lianoglou12k • written 2.1 years ago by jamie.gearing40
Answer: Limma filtering based on treatment group
gravatar for Steve Lianoglou
2.1 years ago by
Steve Lianoglou12k wrote:

No. Don't do that. Your filtering criterion needs to be independent of the test statistic:



ADD COMMENTlink written 2.1 years ago by Steve Lianoglou12k

Very good. Thanks Steve.

ADD REPLYlink written 2.1 years ago by jamie.gearing40
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 278 users visited in the last hour