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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
