[BioC] unexpected result of rep(Rle(), times=<vector or Rle>)
Charles C. Berry
cberry at tajo.ucsd.edu
Sun Jan 31 21:40:52 CET 2010
Using IRanges_1.4.0 here (checked elsewhere with 1.4.8):
I would expect
rep( Rle(x), times=Rle(x) ) == Rle( rep(x, x) )
but
> x <- rep(1:4,c(1,2,1,2))
> Rle(rep(x,x))
'integer' Rle of length 16 with 4 runs
Lengths: 1 4 3 8
Values : 1 2 3 4
> rep(Rle(x),times=Rle(x))
'integer' Rle of length 18 with 4 runs
Lengths: 1 5 3 9
Values : 1 2 3 4
> rep(Rle(x),times=x)
'integer' Rle of length 18 with 4 runs
Lengths: 1 5 3 9
Values : 1 2 3 4
>
Is this the intended behavior?
If so, why?
Chuck
Charles C. Berry (858) 534-2098
Dept of Family/Preventive Medicine
E mailto:cberry at tajo.ucsd.edu UC San Diego
http://famprevmed.ucsd.edu/faculty/cberry/ La Jolla, San Diego 92093-0901
More information about the Bioconductor
mailing list