How to obtain FDR across all regions investigated in DiffBind
2
0
Entering edit mode
29 days ago

Hi guys, I'm interested to make a Volcano plot and using different color i want to show which genes are UP and which are DOWN. To do this I need to use the logFc and FDR relative to ALL investigated regions and not just of DEGs. While logFC was found in in this way:

dbObj_region <- dba.analyze(dbObj_region, method=DBA_EDGER)
b <- dbObj_region$contrasts
all_site <- as.data.frame(b[[1]][5]$edgeR$LRT$table)

I can't find FDR value of all sites. Can you help me??

Thanks a lot Giuseppe

DiffBind • 189 views
ADD COMMENT
1
Entering edit mode
29 days ago
Rory Stark ★ 1.0k

The "correct" way to retrieve results in DiffBind is via the dba.report() function. Setting th=1 returns all results with FDR <= 1.0, which is all of them.

ADD COMMENT
0
Entering edit mode
29 days ago
seidel 7.7k

Could it be as simple as:

p.adjust(all_site$PValue, method="fdr")

? Your code isn't quite clear (to me), but if edgeR is used in the process, and I'm guessing that last line returns a typical edgeR table containing all the things investigated, then a typical way to adjust p-values (i.e. examine FDR) is to use the p.adjust function).

ADD COMMENT
0
Entering edit mode

Good solution... Thanks :)

Giuseppe

ADD REPLY

Login before adding your answer.

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