[BioC] qvalue ( ) ERROR: The estimated pi0 <= 0....

[Saurin Jani]

Hi,

I am trying to find qvalue and get FDR based on
pvalues. I am dong this below: 

> scores <- esApply(esetSub,1,function(x) {
           tmp <- t.test(x)
           c(tmp$statistic,tmp$p.value)
           });

> scores <- t(scores);
> colnames(scores) <- c("t.stat","p.value");
> temp <- scores[,"p.value"] < 0.05;
> esetSub2 <- esetSub[temp,];
> temp1 <- scores[temp,];
> pvalues <- temp1[,"p.value"];

> a <- qvalue(pvalues);
[1] "ERROR: The estimated pi0 <= 0. Check that you
have 
valid p-values or use another lambda method."


Thank you,
Saurin



		
__________________________________ 

Check out the new Yahoo! Front Page. 
www.yahoo.com