Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Feng Zhang
Thank, Jerome The question is if this generalized inverse can make their product to be identity matrix? - Original Message - From: Jerome Asselin [EMAIL PROTECTED] To: Feng Zhang [EMAIL PROTECTED]; R-Help [EMAIL PROTECTED] Sent: Thursday, August 14, 2003 11:52 AM Subject: Re: [R] How to

Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Gabor Grothendieck
As already pointed out ZZ' is not invertible. That is, it is not one-to-one and onto. What we can do is restrict the domain and range of ZZ' to the range of ZZ' or equivalently to the range of Z. In operational terms this means if y = ZZ'x then we only consider y's expressable as y = Zu for

Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Jerome Asselin
No. If B is singular, it's impossible to find a matrix A such that A%*%B be the identity matrix (unless you can find a number x such that x*0=1). Cheers, Jerome On August 14, 2003 10:02 am, Feng Zhang wrote: Thank, Jerome The question is if this generalized inverse can make their product

Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Jerome Asselin
Singular matrices are not invertible. However you can calculate the generalized inverse with the function ginv() from package MASS. HTH, Jerome On August 14, 2003 09:24 am, Feng Zhang wrote: Dear R-listers, I have a dxr matrix Z, where d r. And the product Z*Z' is a singular square