How to retrieve peak summit from DiffBind peaks
1
0
Entering edit mode
fbleao • 0
@7a40ce29
Last seen 4 months ago
United States

Hi,

I don't have experience with R, so I ran Diffbind using Galaxy. I didn't specify the Summit option, so all my peaks have a different width.

Now I want to get the summits of the differential peaks to run motif analysis. Is there an easy way to do this from the Diffbind output files?

Thank you

DiffBind • 248 views
ADD COMMENT
2
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 1 day ago
CRUK, Cambridge, UK

You can get the summits for your consensus set by running dba.count() again with summits=TRUE, and retrieving a GRanges object with the summit position indicated for each sample:

myDBA <- dba.count(myDBA, summits=TRUE, score=DBA_SCORE_SUMMITS_POS)
summits <- dba.peakset(myDBA, bRetrieve=TRUE)
ADD COMMENT

Login before adding your answer.

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