Question: Limma filtering based on treatment group
0
gravatar for jamie.gearing
22 months ago by
jamie.gearing40 wrote:

Hello,
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. https://support.bioconductor.org/p/52762/) 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 • 347 views
ADD COMMENTlink modified 22 months ago by Steve Lianoglou12k • written 22 months ago by jamie.gearing40
Answer: Limma filtering based on treatment group
2
gravatar for Steve Lianoglou
22 months ago by
Denali
Steve Lianoglou12k wrote:

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

http://m.pnas.org/content/107/21/9546.long

 

 

ADD COMMENTlink written 22 months ago by Steve Lianoglou12k

Very good. Thanks Steve.

ADD REPLYlink written 22 months ago by jamie.gearing40
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: 213 users visited in the last hour