How to use Diffbind to quantify signal in a set of custom peaks

Question

[enanoide]

Hi, I would like to use Diffbind to quantify signal at a custom set of regions/peaks Therefore, I would like to leave the set of peaks "untouched", avoiding the generation of consensus peaks or merging, etc. I've found several posts on this issue, and tried many things, but Im still confused and the peaks seem to be always being merged to some degree whatever I try.

Could anyone point me to the right direction on how use Diffbind in this way? Any insight will be much aporeciated.

Best wishes

Answer 1

[LChart]

I assume your custom set of peaks have overlapping or book-ended peaks like:

    |---------------------------|
              peak1
                           |----------------------------------------------------|
                                                        peak2

and this is what you mean by "merging"? Have you tried setting DBA$config$mergeOverlap to 1e9?

[enanoide]

Thanks for the answer. I have tried to do dba twice, the second after setting that option to a number larger than the peak. It does not seem to work. Ill try with 1e9, thanks. I should do the dba twice, right?