about function parameters
1
0
Entering edit mode
1006920318 • 0
@1006920318-23535
Last seen 4.0 years ago
write.table(res$meta.analysis$mu.hat, file="res$mu.hat.txt", sep="\t")
write.table(res$meta.analysis$pval, file="res$pval.txt", sep="\t")
write.table(res$meta.analysis$FDR, file="res$FDR.txt", sep="\t")
FDR <- res$meta.analysis$FDR
FDR_1E10 <- FDR[FDR[,1] < 1E-10,]
write.table(FDR_1E10, file="FDR_1E10.txt", sep="\t")

i don't know how to use function FDR,and FDR_1E10 is numeric empty

software error • 412 views
ADD COMMENT
0
Entering edit mode

This is not a Bioconductor question. I suggest you post this over at biostars.org and invest some effort to explain what you are actually doing. Code can be formatted with the 10101 button.

ADD REPLY
0
Entering edit mode
Kevin Blighe ★ 3.9k
@kevin
Last seen 24 days ago
Republic of Ireland

I think that you may just need to do

FDR_1E10 <- FDR[FDR < 1E-10]

write.table(data.frame(FDR_1E10 = FDR_1E10), file = 'FDR_1E10.txt', sep = '\t')

Also, as ATpoint mentions, for general bioinformatics questions, try http://www.biostars.org in the future.

Kevin

ADD COMMENT

Login before adding your answer.

Traffic: 330 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6