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