[R] Finding solution set of system of linear equations.

[dslowik]

Thanks Robert. That all seems to work. I also found the MASS::Null() function
that gives the null space for the matrix(transpose) given as argument. I am
still trying to appreciate the math behind the Moore-Penrose inverse matrix.
If you have any suggestions for understanding how to use R to solve these
kinds of systems, please send me the pointer.

 I think the R documentation for solve() should state right up front that it
only solves non-singular systems, and it should point me to MASS::ginv() and
how to use it to solve this obvious kind of problem. Better, there could be
a generalized solve() which produces a particular solution, the null space,
image space (basis thereof). This other solution (using ginv() and giA%*%A -
I for the kernel) seems deeply embedded in a particular solution technique
(and should be available), but the generalized solve_lin_sys(), as I
suggested, seems generally quite useful..

Don

--
View this message in context: http://r.789695.n4.nabble.com/Finding-solution-set-of-system-of-linear-equations-tp3541490p3542930.html
Sent from the R help mailing list archive at Nabble.com.