Message-ID: <b19557450902060610j58c06e06yf87cc71b8ce08f0d@mail.gmail.com>
Date: 2009-02-06T14:10:27Z
From: Zhou Fang
Subject: Finding a basis in a set of vectors
Hi,
Okay, I have a n x p matrix X, which I know is not full rank. In
particular, there may be linear dependencies amongst the columns (but
not that many). What is a fast way of finding a linearly independent
subset of the columns of X that will span the column space of X, in R?
If it helps, I have the QR decomposition of the original X 'for free'.
I know that it's possible to do this directly by looping over the
columns and adding them, but at the very least, a solution without
horrible slow loops would be nice.
Any ideas welcome.
Zhou Fang