Re: [R] detecting singular matrices

2003-10-09 Thread Thomas Lumley
On Fri, 10 Oct 2003, Murray Jorgensen wrote: > > > > with a condition number of 1.236260e+16 I think I'd like it to tell me > > it's singular like it used to :) and I know it is singular This was found during code freeze (and I suggested to Martin that he file a bug report on it then, so that we

[R] detecting singular matrices

2003-10-09 Thread Murray Jorgensen
My colleague runs R 1.7.1 under Windows XP. He remarks: A [,1] [,2] [,3] [1,]123 [2,]456 [3,]789 b [1] 1 2 3 solve(A,b) [1] -0.333 0.667 0.000 solve(A) [,1] [,2][,3] [1,] -4.5036e+15 9.00720e+15 -4.5036e+15 [2,] 9.0