[R] lm handling of ill-conditioned systems

[Tomislav Goles]

The lm() function in R seems to handle the inversion of singular X'X matrices
(where there is collinearity between regression inputs) in a way
where one of the inputs is dropped and this also seems to be the
default behavior in SAS (please let me know if i'm wrong about this).
In some other packages (i.e. octave ols() function) the pseudo 
inverse is computed where singular values less then some small threshold are 
not included in computation of inverse but the data columns themselves are not 
dropped).
I am wandering if few experts could comment on what is a better (best?) way to 
handle singular matrices in this context or to perhaps point me to some 
literature on this.
Thanks,
Tomislav Goles
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._