[BioC] LIMMA decideTests result zero from contrast matrix

Ekta Jain Ekta_Jain at jubilantbiosys.com
Wed May 9 10:31:47 CEST 2012


Dear Alex,
Thank you very much. It appears much clear now. I agree, I wouldn't want to apply the filter to unadjusted P values since it just shows that the results obtained are very insignificant. I did not understand in detail how the decideTests() works and was looking for some info regarding the same in case someone else had a similar issue. The venn Diagrams are not important if the data is bad, it is not going to lead anywhere. It was very intriguing when summary(decideTests) gave zeros for all contrasts and now I know exactly why.

Thank you kindly,
Ekta
-----Original Message-----
From: bioconductor-bounces at r-project.org [mailto:bioconductor-bounces at r-project.org] On Behalf Of Alex Gutteridge
Sent: 09 May 2012 13:50
To: bioconductor at r-project.org
Subject: Re: [BioC] LIMMA decideTests result zero from contrast matrix

On 09.05.2012 04:21, Ekta Jain wrote:
> Dear Jim,
> I did change things around when you pointed out for the first time.
> All I have been doing is
>> numGenes <- rownames(eset)
>> topTable(fit2, coef=1, adjust="BH", sort.by="B", number=numGenes)
> And
>> results<-decideTests(fit2, method ="global", lfc =0)
> As you mention in your email that "by default decideTests() uses a BH
> adjusted p-value to filter genes" so am i not applying the same
> adjustment for both the toptable() and decideTests() here?
>
> This is what I am clueless about since I still get zero genes. I
> cannot seem to figure out how to not let decideTests use a BH adjust
> for p-value. For the sake of detail, the code worked fine for all my
> cell lines treated with CPI since the p-values were not as bad as the
> ones for treatment with DMSO.
>
> Thank you,
> Ekta

I think it is the difference in the default P value cutoffs for 
topTable and decideTests that is confusing the issue. From the 
decideTests docs:

Usage:

      
decideTests(object,method="separate",adjust.method="BH",p.value=0.05,lfc=0)

Look at the default p.value cutoff - 0.05. Earlier in the thread you 
said your adjusted P values were ~0.9, hence *nothing* will come through 
the filter. The default for topTable is 1:

Usage:

      topTable(fit, coef=NULL, number=10, genelist=fit$genes, 
adjust.method="BH",
               sort.by="B", resort.by=NULL, p.value=1, lfc=0, 
confint=FALSE)

If you want to apply the filter to unadjusted P values the docs say 
this:

adjust.method: character string specifying p-value adjustment method.
           Possible values are ‘"none"’, ‘"BH"’, ‘"fdr"’ (equivalent to
           ‘"BH"’), ‘"BY"’ and ‘"holm"’. See ‘p.adjust’ for details.

Though I'm not sure why you would want to do this.

-- 
Alex Gutteridge

_______________________________________________
Bioconductor mailing list
Bioconductor at r-project.org
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
The information contained in this electronic message and in any attachments to this message is confidential, legally privileged and intended only for use by the person or entity to which this electronic message is addressed. If you are not the intended recipient, and have received this message in error, please notify the sender and system manager by return email and delete the message and its attachments and also you are hereby notified that any distribution, copying, review, retransmission, dissemination or other use of this electronic transmission or the information contained in it is strictly prohibited. Please note that any views or opinions presented in this email are solely those of the author and may not represent those of the Company or bind the Company. Any commitments made over e-mail are not financially binding on the company unless accompanied or followed by a valid purchase order. This message has been scanned for viruses and dangerous content by Mail Scanner, and is believed to be clean. The Company accepts no liability for any damage caused by any virus transmitted by this email.
www.jubl.com


More information about the Bioconductor mailing list