Re: [R] Possible bug of QR decomposition in package Matrix

2011-08-04 Thread C6H5NO2
Thank you, Martin. Yes, on my computer (48Gb memory) it also gives a message "caught segfault" as you got. So it is because the matrix size is too large. I apologize for the offence I've caused. 2011/8/4 Martin Maechler : >>>>>> C6H5NO2   >>>>&g

Re: [R] Possible bug of QR decomposition in package Matrix

2011-08-04 Thread C6H5NO2
ix_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 wrote: >> Hello R users, >> >> I am trying to give the QR decomposition for a large

[R] Possible bug of QR decomposition in package Matrix

2011-08-03 Thread C6H5NO2
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 window

[R] columnwise nnzero for dgCMatrix

2011-06-26 Thread C6H5NO2
Hi R users, I want to know whether there is a fast method to compute the nonzeroes for each column of dgCMatrix. For summation or average I can use colSums or colMeans. To count the non-zeroes I write a function colCounts <- function(Mat) { M1 <- apply(Mat, MARGIN=2, FUN=nnzero) } But it se