R-alpha: bug in svd
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 =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-