Convert DiffBind output to standard BED or NarrowPeak format
2
0
Entering edit mode
3.7 years ago
rbronste ▴ 390

Trying to figure out how to convert a differentially called list from DiffBind into one of these standard formats?

DiffBind ChIP-Seq BED NarrowPeak • 2.0k views
ADD COMMENT
2
Entering edit mode
3.7 years ago
Rory Stark ★ 1.1k

You can generate bed files from a GRanges object (as returned by dba.report) using export.bed() from the rtracklayer package.

> library(rtracklayer)
> dbsites <- dba.report(myDBA, contrast=1, th=0.05)
> export.bed(dbsites,"dbsites.bed")
ADD COMMENT
0
Entering edit mode
2.2 years ago

I've a slightly more complicated version of this issue. I want to export a list of peaks that do not overlap within 4 groups (e.g. myDBA$onlyA or myDBA$AandC).

From dba.analyze I have 6 contrasts. Which options do I use for dba.report to generate a contrast resembling:

myDBA.OL <- dba.overlap(myDBA, myDBA$masks$Consensus)

myDBA.OL$onlyA

or

myDBA.OL$AandC

so I can export these subsets as .bed?

ADD COMMENT

Login before adding your answer.

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