I wanted something to extract @ranges from a GRanges object along with its
@seqnames, @strand, and @seqinfo. Essentially, everything but the mcols.
Does this make sense? Is there a lighter-weight way to avoid any copying
in-flight?
setMethod("granges", "GRanges", function(x) {
GRanges(seqnames=seqnames(x),
ranges=ranges(x),
strand=strand(x),
seqinfo=seqinfo(x))
})
The fact that I'm constructing an entire new GRanges makes me a little
queasy... that said, it has turned out to be useful when I just want a
short list of locations, as for debugging plotting functions, profile
plots, or what have you.
Statistics is the grammar of science.
Karl Pearson <http://en.wikipedia.org/wiki/The_Grammar_of_Science>
[[alternative HTML version deleted]]
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel
