[R] Possible bug of QR decomposition in package Matrix

C6H5NO2 c6h5no2 at gmail.com
Thu Aug 4 10:03:54 CEST 2011


Thank you very much, Josh!

As you suggested, I will contact the developers of "Matrix".

PS, C6 are just initial characters of my email account :-)

Best wishes,

C6

2011/8/4 Joshua Wiley <jwiley.psych at gmail.com>:
> 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 <c6h5no2 at 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 at 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/
>



More information about the R-help mailing list