On Tue, Mar 16, 2010 at 12:34 PM, Ivan Gregoretti <[email protected]>wrote:
> Hello everybody, > > Say A, B, C, D, .... are all RangedData instances. > > How do you come up with the list of ranges in A that have at least > some overlap with B, C, D, ...? > > I want to calculate the multiple intersection ignoring the extent of > the overlap. > > ?findOverlaps does not hint how to recover the ranges from A. > > Something like: ranges(A) %in% union(union(ranges(B), ranges(C)), ranges(D)) That syntax is a little verbose. I propose adding some operators to Ranges, RangesList and RangedData objects. "+" for union() "-" for setdiff() "!" for gaps() Then we could have: A %in% (B + C + D) What do people think? > 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 > [[alternative HTML version deleted]] _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
