[BioC] transcriptLocs2refLocs Bug when strand is negative

Cook, Malcolm MEC at stowers.org
Sat Sep 1 01:22:39 CEST 2012


Herve et al,

I believe that transcriptLocs2refLocs is implemented incorrectly for
negative strands.

The bug appears in the example of its use you provided in
https://stat.ethz.ch/pipermail/bioconductor/2011-March/038173.html

I think you will agree that the third result should be:
[[3]]
[1] 5899 6291

As a workaround, I have implemented a version in R using a new generic,
named `at`

https://gist.github.com/3560692


I would be happy to receive confirmation that this is a bug, and comments
from anyone on the approach to the workaround I have adopted, especially
as regards performance comparison with

Thanks,

Malcolm Cook
Stowers Institute for Medical Research



More information about the Bioconductor mailing list