If you don't mind, I think I'll wait till after you're finished testing before merging.
Please keep me up to date, Michael On Fri, Aug 2, 2013 at 7:43 AM, Hector Corrada Bravo <hcorr...@gmail.com>wrote: > Thanks again Kasper, > All of these are now fixed: > * empty subject seqlevel, > * out-of-order seqinfo's, > * show,GIntervalTree-method > > Michael, please merge from > https://github.com/hcorrada/IRanges [version 1.19.20] > https://github.com/hcorrada/GenomicRanges [version 1.13.36] > > I'll test the SummarizedExperiment stuff later today. > > Cheers, > Hector > > > > On Thu, Aug 1, 2013 at 5:21 PM, Kasper Daniel Hansen < > kasperdanielhan...@gmail.com> wrote: > >> If you're fixing tonight here is a possible related. I get an error when >> I do the findOverlaps and I have two GRanges with the same seqlevels but in >> different order (so, I guess, technically not the same seqlevels). If this >> is not supported, I suggest an upfront check. >> >> gr1 = GRanges(seqnames = c("chr1", "chr2"), ranges = IRanges(1:2, width = >> 1)) >> gr2 = gr1 >> seqlevels(gr2, force = TRUE) = c("chr2", "chr1") >> findOverlaps(GIntervalTree(gr1), GIntervalTree(gr2)) >> >> I also got a seqfault on a similar example (but with real data). I am >> trying to see if I can make a small reproducible example (which is how I >> got here), hypothesizing that it has something to do with the seqinfo slot >> in the two GRanges. >> >> Kasper >> >> >> >> >> >> >> On Thu, Aug 1, 2013 at 4:48 PM, Hector Corrada Bravo >> <hcorr...@gmail.com>wrote: >> >>> Thanks. I'll fix tonight. >>> >>> >>> On Thu, Aug 1, 2013 at 4:28 PM, Kasper Daniel Hansen < >>> kasperdanielhan...@gmail.com> wrote: >>> >>>> Bug: >>>> >>>> > library(GenomicRanges) ## 1.13.35 >>>> >>>> > gr = GRanges(seqnames = c("chr1", "chr2"), ranges = IRanges(1:2, >>>> width = 1)) >>>> >>>> > gr.tree <- GIntervalTree(gr) >>>> >>>> >>>> > findOverlaps(gr.tree, gr.tree) >>>> >>>> >>>> Hits of length 2 >>>> >>>> >>>> queryLength: 2 >>>> >>>> >>>> subjectLength: 2 >>>> >>>> >>>> queryHits subjectHits >>>> >>>> >>>> <integer> <integer> >>>> >>>> >>>> 1 1 1 >>>> >>>> >>>> 2 2 2 >>>> >>>> >>>> !> findOverlaps(gr.tree, gr.tree[1]) >>>> >>>> >>>> >>>> >>>> >>>> *** caught segfault *** >>>> >>>> >>>> address (nil), cause 'memory not mapped' >>>> >>>> >>>> >>>> >>>> >>>> Traceback: >>>> >>>> >>>> 1: .Call(.NAME, ..., PACKAGE = PACKAGE) >>>> >>>> >>>> 2: .Call2(fun, object@ptr, ..., PACKAGE = "IRanges") >>>> >>>> >>>> 3: .IntervalForestCall(subject, fun, query, partitionIndices, >>>> elementLengths(partitioning), query_ord) >>>> >>>> 4: .local(query, subject, maxgap, minoverlap, type, select, ...) >>>> >>>> >>>> 5: findOverlaps(qlist, subject@ranges, maxgap = maxgap, minoverlap = >>>> minoverlap, type = type, select = "all") >>>> >>>> 6: findOverlaps(qlist, subject@ranges, maxgap = maxgap, minoverlap = >>>> minoverlap, type = type, select = "all") >>>> >>>> 7: .local(query, subject, maxgap, minoverlap, type, select, ...) >>>> >>>> >>>> 8: findOverlaps(gr.tree, gr.tree[1]) >>>> >>>> >>>> 9: findOverlaps(gr.tree, gr.tree[1]) >>>> >>>> >>>> >>>> >>>> >>>> Possible actions: >>>> >>>> >>>> 1: abort (with core dump, if enabled) >>>> >>>> >>>> 2: normal R exit >>>> >>>> >>>> 3: exit R without saving workspace >>>> >>>> >>>> 4: exit R saving workspace >>>> >>>> >>>> Selection: >>>> >>>> >>>> On Thu, Aug 1, 2013 at 3:45 PM, Kasper Daniel Hansen < >>>> kasperdanielhan...@gmail.com> wrote: >>>> >>>>> Ah! Super-nice! I'll test as well, when I have time. >>>>> >>>>> Kasper >>>>> >>>>> >>>>> On Thu, Aug 1, 2013 at 3:34 PM, Hector Corrada Bravo < >>>>> hcorr...@gmail.com> wrote: >>>>> >>>>>> This is not tested, but... >>>>>> >>>>>> GIntervalTree extends GenomicRanges so, in principle, can be used in >>>>>> the rowData slot of SummarizedExperiment. I tried to make GIntervalTree >>>>>> support as much of the GenomicRanges interface as I could so this may >>>>>> work >>>>>> already for overlap queries where the SummarizedExperiment object is the >>>>>> subject. I'll test this soon and let you know. >>>>>> >>>>>> >>>>>> On Thu, Aug 1, 2013 at 3:12 PM, Kasper Daniel Hansen < >>>>>> kasperdanielhan...@gmail.com> wrote: >>>>>> >>>>>>> Now that we have GIntervalTree, I am really interested in a >>>>>>> SummarizedExperiment that stores its index in a slot, so it can be >>>>>>> indexed >>>>>>> once. I am sure other people are greedily eyeing this as well. >>>>>>> Either >>>>>>> changing the class or extending it. >>>>>>> >>>>>>> Unfortunately, I will not have time to work on this for some weeks, >>>>>>> but I >>>>>>> wanted to put the idea out there. >>>>>>> >>>>>>> Best, >>>>>>> Kasper >>>>>>> >>>>>>> [[alternative HTML version deleted]] >>>>>>> >>>>>>> _______________________________________________ >>>>>>> Bioc-devel@r-project.org mailing list >>>>>>> https://stat.ethz.ch/mailman/listinfo/bioc-devel >>>>>>> >>>>>>> >>>>>> >>>>> >>>> >>> >> > [[alternative HTML version deleted]] _______________________________________________ Bioc-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/bioc-devel