Hi C6 (were C1 - 5 already taken in your family?), I downloaded your data and can replicate your problem. R ceases responding and terminates. This does not occur with all uses of qr on a dgCMatrix object. I know nothing about sparse matrices, but if you believe this should not be occurring, you should contact the package maintainers. Here is my sessionInfo() (FYI, it would probably be helpful to report yours also in case the issue is version dependent):
R Under development (unstable) (2011-07-30 r56564) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] Matrix_0.999375-50 lattice_0.19-30 loaded via a namespace (and not attached): [1] grid_2.14.0 tools_2.14.0 Cheers, Josh On Wed, Aug 3, 2011 at 4:26 PM, C6H5NO2 <c6h5...@gmail.com> wrote: > 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 > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joshua Wiley Ph.D. Student, Health Psychology Programmer Analyst II, ATS Statistical Consulting Group University of California, Los Angeles https://joshuawiley.com/ ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.