[R] Matrix decomposition: orthogonal complement

2004-11-02 Thread Florin G. Maican
Hello, How I can compute in R the orthogonal complement of one matrix? If A (n x m ) matrix of full column rank (nm), its orthogonal complement is denoted by A_ . A_ is n X (n-m) matrix of full column rank and such that A'A_=0. I need to compute A_. How I can compute A_ in R? Best

Re: [R] Matrix decomposition: orthogonal complement

2004-11-02 Thread Prof Brian Ripley
On Tue, 2 Nov 2004, Florin G. Maican wrote: How I can compute in R the orthogonal complement of one matrix? If A (n x m ) matrix of full column rank (nm), its orthogonal complement is denoted by A_ . A_ is n X (n-m) matrix of full column rank and such that A'A_=0. It is far from

Re: [R] Matrix decomposition: orthogonal complement

2004-11-02 Thread Simon Wood
How I can compute in R the orthogonal complement of one matrix? use the qr decomposition. For example: A-matrix(rnorm(40),10,4) B - t(qr.Q(qr(A),complete=TRUE)[,5:10]) B%*%A best, Simon If A (n x m ) matrix of full column rank (nm), its orthogonal complement is denoted by A_ . A_