Biostar Beta. Not for public use.
how to retrieve binding affinity matrix in diffbind
0
Entering edit mode
19 months ago
hamaor • 0

hey everyone, I'm performing ChIP-Seq analysis using diffbind, after doing:

dba.count(peaksets, summits=250)

I would like to get the binding affinity Matrix into a text file.

help will be highly appriciated!!

ADD COMMENTlink
0
Entering edit mode
20 months ago
iamjli • 10

For completeness, this question has been answered here.

ADD COMMENTlink
0
Entering edit mode
16 months ago
Rory Stark • 510
University of Cambridge, Cancer Researc…
dba.peakset(myDBA, bRetrieve=TRUE, writeFile="BindingAffinityMatrix.csv")

The values will be the score used in dba.count(), the default is DBA_SCORE_TMM_MINUS_FULL (TMM normalised read counts after subtracting control reads, based on the number of reads in each bam file). To get, for example, raw read counts, use DBA_SCORE_READS. See help page for dba.count(). You can change the score without recounting:

myDBA <- dba.count(myDBA, peaks=NULL, score=DBA_SCORE_xxx)
ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1