Entering edit mode
Dear all,
I observed this problem regarding the maximal length of a Rle vector:
> rle = Rle(rep(0, 1000000000))
> length(rle)
[1] 1000000000
> length(c(rle, rle, rle))
[1] -1294967296
Probably, it is caused by the maximum positive number (~2.1E9) that
can
be represented by an integer variable. However, there is no warning
message.
I noticed this problem when I wanted to calculate the average coverage
of a sequencing project accross the human genome. I used the
coverage()
method and then concatenated all chromosomes. This should give me an
Rle
vector of length ~3*109, but mean() does not work on that vector.
Best,
Hans-Ulrich
> sessionInfo()
R version 2.14.0 (2011-10-31)
Platform: x86_64-pc-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=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] IRanges_1.12.3
[[alternative HTML version deleted]]