Paul Gilbert <pgilbert@bank-banque-canada.ca> writes:
> Eigen seems to be giving different results in 0.64 :
>
> eigen(matrix(
> c(0.00000000, 0.0000000, 0.00000000, 1.000000000,
> 0.01262442, -0.1141633, 0.04989433, 0.006112607,
> 0.04781087, 0.4177869, -0.15569656, 0.074242237,
> 0.20167331, 0.2531109, -0.04973573, 0.034046467) ,4,4))$values
>
> R:
> [1] 0.35519256 -0.29294606 -0.26637369 -0.03168619
>
> S:
> [1] 0.4689288205 -0.4246063863 -0.2810596564 0.0009238293
>
> I am fairly sure that S and R previously gave the same results.
Indeedy. At least the pre 63.3 binary I had lying about did give the S
results. The eigenvectors are also bogus as far as I can tell.
The weird thing is that eigen() goes via .Fortran to eigen.c
(which was f2c'ed at some point) and neither eigen.c or eigen.R has
been changed in 1999 according to the cvs logs.
--
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 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._