[Bioc-sig-seq] ChIPpeakAnno Error when Strand is *

[Dario Strbenac]

It has become standard for Bioconductor packages to use * to represent either strand. When I try to use annotatePeakInBatch with RangedData that has * in it, it fails. Could this be repaired ?

Minimal example :

library(ChIPpeakAnno)

peaksGR <- GRanges("chr1", IRanges(seq(100, 300, 100), width = 20))
genesGR <- GRanges("chr1", IRanges(200, 500), '+')

peaksRD <- as(peaksGR, "RangedData")
genesRD <- as(genesGR, "RangedData")

annotatePeakInBatch(peaksRD, AnnotationData = genesRD)

--------------------------------------
Dario Strbenac
Research Assistant
Cancer Epigenetics
Garvan Institute of Medical Research
Darlinghurst NSW 2010
Australia