[BioC] GRanges setdiff With Features on Both Strands

[Steve Lianoglou]

Hi,

2012/1/23 Hervé Pagès <hpages at fhcrc.org>:
[snip]

>  union(`strand<-`(x, "*"), `strand<-`(y, "*"))
>  intersect(`strand<-`(x, "*"), `strand<-`(y, "*"))
>  setdiff(`strand<-`(x, "*"), `strand<-`(y, "*"))
>
> is just a cosmetic question but I agree that the output of the latter
> will look prettier (and might be slightly less confusing). So I can do
> this. Does anybody object?

Not sure who you're inviting to the tally, but I'm +1 on doing this
latter way when ignore.strand=TRUE, ie. setting strand to "*" then
doing the set algebra as is suggested in the quoted piece above ...

-steve

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact