question about DiffBind
0
0
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 5 weeks ago
Cambridge, UK
Hi Rong- You are correct that the peak merging logic results in wider peak areas that may incorporate multiple binding sites and/or make it more difficult to determine the precise binding site. This is true of the peak calling step as well. In our experience, the scientific questions we are actually asking do not require the precise binding location in the differential analysis step -- we are first looking for places where there is a significant difference in binding levels. It is only after that, and only in a minority of cases, that we look to narrow the region down to a specific binding site of interest (eg, to find enriched motifs). Tehre are a numebr of tools available to find more narrow, precise binding sites in a set of candidate regions. The most recent release of DiffBind does, however, include some features to better a) focus on the likely binding site and b) work with more uniform, narrower peak widths. It accomplishes this using the "summits" parameter in dba.count. This parameter allows you to identify the "summit" for each peak interval (the point of highest pileup), and optionally re-center the peaks around this point (eg, taking 200bp upstream and downstream of the summit for each merged peak; summits=200). If there were several "true" binding sites in a merged peak, you are likely to end up with only one of them, but the results are more consistent overall. Cheers- Rory On 12/05/2014 17:40, "Rong Cong" <rong.cong at="" inserm.fr=""> wrote: >Hi Rory, >Thank you very much for your patient and efficient reply to my naive >question. The short version answer made me quite clear of the idea. > >Indeed, overlapping and merging peaks may form a quite wider peak to be >compared among different peaksets, but for biologists, the wider peak >often has no meaning to us(some biologically functional elements are >quite short),therefore is it possible for DiffBind to adjust some >parameters to narrow a little the overlapped and merged peak for >comparison? (Again, I am not a bioinformatician, sorry for this naive >question.) >Looking forward to your reply. >All the best, >Rong >
DiffBind DiffBind • 1.2k views
ADD COMMENT

Login before adding your answer.

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