# R-alpha: bug in svd

**Peter Dalgaard BSA
**
p.dalgaard@kubism.ku.dk

*17 Oct 1997 14:42:22 +0200*

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
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-