This is a test version of Biostars. For the public version, visit https://www.biostars.org.
R long vectors not supported yet on R 3.3.2

0 down vote favorite

I am trying to simulate hilbertcurve (of Bioconductor package) of level 16 in R. It takes about 4^16=4 Billion points. I want to generate the hilbert curve of human genome (size about 3 billion).

But I am getting this error

long vectors not supported yet: memory.c:1668 I am using 64 bit version (R 3.3.2) so my guess is I can use vectors of length > 2^31. Also, my RAM is about 350GB.

The command I am using is

itr=4^16 hc = HilbertCurve(1, itr, 16, mode = "pixel", title = "pixel mode",start_from = "topleft")
Even when I am reading the whole genome sometimes R is crashing in the process.

I have read the other similar questions on this topic but could not find the solution. Please help me what should I use for this problem.

r genome sequence

You may have better luck asking on the R mailing list since this is about R and not a bioinformatics question.

EDIT: As stated in the docs:

Whether matrix algebra (including %*%, crossprod, svd, qr, solve and eigen) will actually work is somewhat implementation dependent, including the Fortran compiler used and if an external BLAS or LAPACK is used.

Hello user12345!

We believe that this post does not fit the main topic of this site.

For questions very specific to the implementation of Bioconductor packages, please use: https://support.bioconductor.org/

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

0 answers

No answers yet.

Log in to answer this question.