[Bioc-devel] flank function with start argument a named vector
Leonard Goldstein
goldstein.leonard at gene.com
Fri Mar 28 17:15:55 CET 2014
Dear Bioc developers,
I ran into problems when using the flank function with the start
argument set to a named vector.
This returns an IRanges with named starts and ends, which causes
problems downstream. Please see example below.
Thanks for your help
Leonard
> ir <- IRanges(1, 1)
> fl <- flank(ir, 1, c(a = FALSE))
> fl
IRanges of length 1
start end width
[1] 2 2 1
> start(fl)
a
2
> end(fl)
a
2
> findOverlaps(ir, fl)
Error in findOverlaps(query, IntervalTree(subject), maxgap = maxgap,
minoverlap = minoverlap, :
error in evaluating the argument 'subject' in selecting a method for
function 'findOverlaps': Error in validObject(from) :
invalid class "IRanges" object: 'start(x)' must be an unnamed
integer vector with no NAs
>
> sessionInfo()
R Under development (unstable) (2013-12-03 r64376)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] parallel stats graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] IRanges_1.21.36 BiocGenerics_0.9.3
loaded via a namespace (and not attached):
[1] stats4_3.1.0
>
More information about the Bioc-devel
mailing list