[R] Possible bug of QR decomposition in package Matrix

C6H5NO2 c6h5no2 at gmail.com
Thu Aug 4 01:26:09 CEST 2011


Hello R users,

I am trying to give the QR decomposition for a large sparse matrix in
the format of dgCMatrix. When I run qr function for this matrix, the R
session simply stops and exits to the shell.
The matrix is of size 108595x108595, and it has 4866885 non-zeros. I
did the experiment on windows 7 and linux mint 11 (both 64 bit), and
the results are the same.

I have uploaded my data file to http://ifile.it/elf2p6z/A.RData . The
file is 10.681 MB and I hope someone could kindly download it.
The code to see my problem is:
library(Matrix)
load("A.RData")
B <- qr(A)

Best wishes,
C6



More information about the R-help mailing list