Hi everybody,
I have two sets of ChIP-seq peaks: A and B.
How do I select the peaks in A that are not overlapped by any peak in B?
Now in more detail:
I do
library(rtracklayer)
A <- import('condition_A.bed', 'bed')
B <- import('condition_B.bed', 'bed')
they are both RangedData instances.
I need to compute A - B.
How do you guys do it?
Thank you,
Ivan
Ivan Gregoretti, PhD
National Institute of Diabetes and Digestive and Kidney Diseases
National Institutes of Health
5 Memorial Dr, Building 5, Room 205.
Bethesda, MD 20892. USA.
Phone: 1-301-496-1592
Fax: 1-301-496-9878
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
