> -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Rosario Garcia Gil > Sent: Monday, March 28, 2011 7:51 AM > To: r-help@r-project.org > Subject: [R] matrix inverstion > > Hello > > I have this matrix which I am trying to invert. I get a message about > reciprocal condition number, what that does mean? > > > XT_X > [,1] [,2] [,3] [,4] [,5] > [1,] 3 0 0 2 1 > [2,] 0 2 0 1 1 > [3,] 0 0 2 1 1 > [4,] 2 1 1 4 0 > [5,] 1 1 1 0 3 > > iXT_X <- solve(XT_X) > Error in solve.default(XT_X) : > system is computationally singular: reciprocal condition number = > 1.11022e-17 > >
Well, it means exactly what the message says. Within the precision of your computer, the matrix is singular and has no inverse. If you try the following code you will see that column 1 is perfectly correlated with the remaining columns. summary(lm(XT_X [,1] ~ XT_X[,-1])) Hope this is helpful, Dan Daniel Nordlund Bothell, WA USA ______________________________________________ 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.
