It is very common to want to order a GenomicRanges by *position* but not
strand. Would be nice to have an ignore.strand argument for ordering
GenomicRanges. I realize that this is not possible by the order generic.
But it would be possible for sort().
I also realize that someone could do this:
gr2 <- gr
strand(gr2) <- "*"
gr[order(gr2)]
But that is long-winded. Perhaps there could be an unstrand()?
gr[order(unstrand(gr))]
But this is still better:
sort(gr, ignore.strand=TRUE)
Sound like a good thing?
Michael
[[alternative HTML version deleted]]
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel
