News:New in DiffBind 3.2: dba.plotProfile() - Peak profiles and profile heatmaps
0
0
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 9 hours ago
CRUK, Cambridge, UK

For DiffBind users, there is a new function in Bioconductor 3.13, dba.plotProfile(), that makes it straightforward to calculate peak profiles and generate complex heatmap plots.

dba.plotProfile() flexibly packages different sets of sites and samples from a DiffBind analysis for the Bioconductor package profileplyr (by Tom Carroll and Doug Barrows), generating profileplyr objects that can be customized using that package, including exporting to deepTools.

A tutorial in the form of a markdown workbook is available, showing how to use the dba.plotProfile() function. For example, using the vignette tamoxifen dataset, the following plot can be generated with one simple call:

profiles <- dba.plotProfile(tamoxifen, doPlot=TRUE)

Default profile plot from `dba.plotProfile()` using sample analysis

profileplyr ATACSeq DiffBind ComplexHeatmap ChIPSeq • 304 views
ADD COMMENT
0
Entering edit mode

This looks great! Is there a way to adjust the +/- 1500 window size?

ADD REPLY
1
Entering edit mode

Yes, you have full control over the windows by passing through parameters to profileplyr::BamBigwig_to_chipProfile(), such as style, distanceAround, nOfWindows, and bin_size.

If you look at the tutorial linked in the original post, Section 7 explains how to control the windows.

ADD REPLY

Login before adding your answer.

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