R-alpha: bug in svd

[Peter Dalgaard BSA]

Andreas Weingessel <Andreas.Weingessel@ci.tuwien.ac.at> writes:

> I use R Version 0.60 Alpha (September 18, 1997) on a Linux Pentium
                  ^^^^^^^^^^
> (Debian 1.3) and on a Sparc-Sun-Solaris 2.5. 

As a matter of principle, that version is not released yet, so it may
contain bugs that shouldn't be seen outside the core team. However,
I seem to get the same errors with 0.50.x (with x==2 because I was to
laz^H^H^Hbusy to upgrade it...), so I suppose it can be discussed here
anyway.

> R> svd(matrix(1:16,4,4))
> 	yields on both machines
> Error: error  4  in dsvdc
> 
> R> svd(matrix(1:20,4,5))
> 	gives a result on the Linux computer

Both examples are rank-2 matrices, the problems disappear with random
matrices, so I suspect that the code is unhappy about dealing with
matrices of less than full rank?

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-