getting DiffBind full result list
1
0
Entering edit mode
@rafaelmunozviana-16038
Last seen 6.5 years ago

Dear Bioconducters

I am running ChIPseq analyses via DiffBind 2.2.12 using DESeq option.
As indicated in the vignette I run dba.analyze() and the extract the results with dba.report()

I obtain the GR of the peaks that are significative changed, but I wonder if it is possible to get the data (Fold and FDR) of those that show no appreciable change

 

Kind regards

deseq2 diffbind chipseq • 863 views
ADD COMMENT
1
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 4 weeks ago
Cambridge, UK

You can get stats for all the peaks by setting th=1 in the call to dba.report().

It's worth noting that it is not necessarily the case that peaks with higher FDR values are not changing, rather that there is insufficient evidence to be confident that they did change (potentially due to e.g. inadequate power given the variance in binding values for that peak).

 

ADD COMMENT

Login before adding your answer.

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