Message-ID: <x24t6g8ssx.fsf@blueberry.kubism.ku.dk>
Date: 1997-10-17T12:42:22Z
From: Peter Dalgaard
Subject: R-alpha: bug in svd
In-Reply-To: Andreas Weingessel's message of Fri, 17 Oct 1997 16:42:49 +1000
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
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-