[R] QR decomposition and rank of a matrix

From: Alex Nu <notulei_at_yahoo.com>
Date: Fri 02 Jul 2004 - 00:59:41 EST

 In summary.manova the qr decomposition of a NxN matrix
 is calculated and for some cases is giving me  a rank < N.

 However, following suggestions of professor Ripley to

 calculate the rank of a Matrix

On 7 Jun 2002, Brian Ripley wrote:
> For a more reliable answer, look at the SVD
> (function svd) and look at the
> singular values. For example (from lda.default)

    X.s <- svd(X, nu = 0)
    rank <- sum(X.s$d > tol * X.s$d[1])

I'm getting rrank = N,

 I wonder wether use this new method for  calculating the rank inside the manova.R

 I'd appreciate any suggestions.

 Thanks

 Alex



R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Jul 02 01:05:36 2004

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 08:55:49 EST