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. Paul Gilbert -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._