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
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
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_