edgeR
1
0
Entering edit mode
@da307361
Last seen 2.3 years ago
Turkey

Hi, in edgeR how can I write the FDR values to DE file?

et <- exactTest(dgelist)
fdr <- p.adjust(et$table$PValue, method="BH")

I know exactTest doesnt produce FDR values but with fdr <- p.adjust(et$table$PValue, method="BH") i think i am calculating fdr values but how can I include them to DE file?

rnaseq edgeR • 817 views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode
@gordon-smyth
Last seen 5 minutes ago
WEHI, Melbourne, Australia

The documentation shows that you should use:

topTags(et)

The output not only includes FDR values but also ranks the genes in order of DE significance.

ADD COMMENT

Login before adding your answer.

Traffic: 868 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