PePr icon indicating copy to clipboard operation
PePr copied to clipboard

Remove chip1 / chip2 intersecting regions

Open rcavalcante opened this issue 8 years ago • 3 comments

I've observed that chip1 peaks and chip2 peaks can overlap a little bit at the edges. We've discussed how that comes from the windowing approach, but would it be possible to add a post-processing step that removes the regions that are in the intersection of chip1 and chip2 peaks from both of the results files?

rcavalcante avatar Jun 15 '16 15:06 rcavalcante

Hi Ray, could you send me an example of it so I know what it is like?

On Wednesday, June 15, 2016, Raymond Cavalcante [email protected] wrote:

I've observed that chip1 peaks and chip2 peaks can overlap a little bit at the edges. We've discussed how that comes from the windowing approach, but would it be possible to add a post-processing step that removes the regions that are in the intersection of chip1 and chip2 peaks from both of the results files?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/shawnzhangyx/PePr/issues/16, or mute the thread https://github.com/notifications/unsubscribe/AFOG6KdZmN4tqSfL9iqFqAXv88m3-fDFks5qMBzfgaJpZM4I2fxR .

Sent from Gmail Mobile

shawnzhangyx avatar Jun 15 '16 16:06 shawnzhangyx

Hey Yanxiao,

I attached a zip of chip1 and chip2 peaks, plus the result of bedtools intersect -wa -wb -a chip1 -b chip2.

peak_overlaps.zip

rcavalcante avatar Jun 15 '16 21:06 rcavalcante

I find the easiest way is to use bedtools subtract -a chip1 -b chip2 to trim the edges.

shawnzhangyx avatar Jun 27 '16 16:06 shawnzhangyx