Diffbind: How to compare the DBs called with and without a blocking factor
Entering edit mode
zhangly811 • 0
Last seen 4.6 years ago

Hi, I'm doing a differential binding analysis between two experimental conditions across multiple (~5) cell lines. I'm interested in the highly conserved DBs across all cell lines as well as the DBs that are unique to each cell line. I have tried the blocking factor in diffbind, and got two lists of DBs from analysis with and without blocking factor. Is there built-in functions in Diffbind for me to know which DBs overlap between the two lists, and for the ones differ, from which cell line they are detected to be differentially bound? Thank you for your help.



diffbind differential binding analysis blocking • 540 views
Entering edit mode
Rory Stark ★ 4.4k
Last seen 5 days ago
CRUK, Cambridge, UK

You can use the dba.report() function in a special mode to get what you are interested in.

By setting bDB=TRUE, dba.report() will return a new DBA object containing different peaksets based on whether they are differentially bound and represent a gain or loss of binding. You can pass in a vector of contrasts to include differentially bound peaks from different comparisons, and a vector of analysis methods as well to compare sites from different analysis methods: eg. method=c(DBA_DESEQ2,DBA_DESEQ2_BLOCK).

The last code snippet on the help page shows an example -- I've highlighted the key line:

#Retrieve report-based DBA object
tamoxifen <- dba.contrast(tamoxifen, categories=DBA_CONDITION, block=tamoxifen$masks$MCF7)
tamoxifen <- dba.analyze(tamoxifen,bCorPlot=FALSE)
tamoxifen.DB <- dba.report(tamoxifen,method=c(DBA_DESEQ2,DBA_DESEQ2_BLOCK),
                          bDB=TRUE, bGain=TRUE, bLoss=TRUE, bAll=FALSE)
dba.plotVenn(tamoxifen.DB,1:4,label1="Single Factor GAIN",label2="Single Factor LOSS",
                             label3="Blocking Factor GAIN",label4="Blocking Factor LOSS")

This feature should enable you to get what you need!




Login before adding your answer.

Traffic: 403 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6