Running the following code on the current GenomicRanges devel version (1.31.3):
library(GenomicRanges)
example(GenomicRanges, echo=FALSE)
olap <- findOverlaps(gr, gr)
class(sort(olap))
... gives a Hits object, while previously in the release (1.30.0) it gave a SortedByQueryHits object. Is this change intentional? I am curious because my InteractionSet package defines its own findOverlaps method, and I would like to know what the correct return class should be.
Running the following code on the current GenomicRanges devel version (1.31.3):
... gives a
Hits
object, while previously in the release (1.30.0) it gave aSortedByQueryHits
object. Is this change intentional? I am curious because my InteractionSet package defines its ownfindOverlaps
method, and I would like to know what the correct return class should be.