[Bioc-devel] ordering GenomicRanges without regard for strand

Hervé Pagès hpages at fhcrc.org
Tue Dec 17 08:26:39 CET 2013


OK. Let's go for unstrand() too. Will add this in the next few days.

H.


On 12/12/2013 05:07 PM, Kasper Daniel Hansen wrote:
> we should have the argument for sort, but I also find unstrand() to be a
> great idea.
>
>
> On Thu, Dec 12, 2013 at 7:19 PM, Hervé Pagès <hpages at fhcrc.org
> <mailto:hpages at fhcrc.org>> wrote:
>
>     Hi Michael,
>
>     Feels like Christmas is not far away ;-)
>
>     Right now only the "match" method in the GenomicRanges-comparison.R
>     family supports the 'ignore.strand' arg and I think it would be a good
>     thing to have other "comparison" methods support it. It's too bad
>     for order() though (but maybe there is a way to hack the definition of
>     the generic in BiocGeneric?). Should be easy for sort(). Will do.
>
>     Cheers,
>     H.
>
>
>
>     On 12/12/2013 01:44 PM, Michael Lawrence wrote:
>
>         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 at r-project.org <mailto:Bioc-devel at r-project.org>
>         mailing list
>         https://stat.ethz.ch/mailman/__listinfo/bioc-devel
>         <https://stat.ethz.ch/mailman/listinfo/bioc-devel>
>
>
>     --
>     Hervé Pagès
>
>     Program in Computational Biology
>     Division of Public Health Sciences
>     Fred Hutchinson Cancer Research Center
>     1100 Fairview Ave. N, M1-B514
>     P.O. Box 19024
>     Seattle, WA 98109-1024
>
>     E-mail: hpages at fhcrc.org <mailto:hpages at fhcrc.org>
>     Phone: (206) 667-5791 <tel:%28206%29%20667-5791>
>     Fax: (206) 667-1319 <tel:%28206%29%20667-1319>
>
>
>     _________________________________________________
>     Bioc-devel at r-project.org <mailto:Bioc-devel at r-project.org> mailing list
>     https://stat.ethz.ch/mailman/__listinfo/bioc-devel
>     <https://stat.ethz.ch/mailman/listinfo/bioc-devel>
>
>

-- 
Hervé Pagès

Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M1-B514
P.O. Box 19024
Seattle, WA 98109-1024

E-mail: hpages at fhcrc.org
Phone:  (206) 667-5791
Fax:    (206) 667-1319



More information about the Bioc-devel mailing list