[Bioc-devel] How best to remap S4Vectors::Hits indices?
p@rik@heet@n@nd@ @ending from uconn@edu
Wed May 23 01:57:37 CEST 2018
I'm working on a package that does some trivial GRanges position
classifications; primarily to standardize nomenclature according to the
literature in workflows.
The API for S4Vectors::Hits() generally doesn't seem amenable to modify
Hits objects, except for the remapHits() feature (which I see underneath
the covers really generates a new Hits object).
I was hoping someone could take a quick look at a short function I'm using
to subset and reindex Hits in the da_tss() function:
Yes, to illustrate the problem I'm having, I've directly used the @-style
S4 access which is, of course, a terrible thing to do because it defeats
the purpose of S4 object validation, which is why I'm e-mailing the list
for an alternative. I feel like casting to something like a data.frame,
changing the indices, and changing back to Hits would be wasteful and
improperly using the Bioconductor framework?
Here are the corresponding tests that run the da_tss() function:
What it comes down to is this:
I want to compare a subset of GRanges for hits, but revert to the original
GRanges indices when returning the results.
Thanks for any advice!
[[alternative HTML version deleted]]
More information about the Bioc-devel