armin@approximity.com
2002-Jan-05 00:48 UTC
[Rd] possible bug in det using method="qr" (PR#1244)
Full_Name: Armin Roehrl Version: 1.4.0(2001-12-19) OS: Linux; Suse 7.3 Submission from: (NULL) (217.84.18.204) Hello, On a given matrix det returns different results whether I use the method "qr" or "eivenvalues". The matrix:> m2[,1] [,2] [,3] [1,] 9.822616e+09 3.841723e+09 79790.09 [2,] 3.841723e+09 1.502536e+09 31251.82 [3,] 7.979009e+04 3.125182e+04 64156419.36> det(m2, method="eigenvalues")[1] -9.331893e+19> det(m2, method="qr")[1] 9.331893e+19 I think the result based on the eigenvalues is correct. Good night. Armin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Peter Dalgaard BSA
2002-Jan-05 02:10 UTC
[Rd] possible bug in det using method="qr" (PR#1244)
armin@approximity.com writes:> Full_Name: Armin Roehrl > Version: 1.4.0(2001-12-19) > OS: Linux; Suse 7.3 > Submission from: (NULL) (217.84.18.204) > > > Hello, > > On a given matrix det returns different results whether > I use the method "qr" or "eivenvalues". > > The matrix: > > m2 > [,1] [,2] [,3] > [1,] 9.822616e+09 3.841723e+09 79790.09 > [2,] 3.841723e+09 1.502536e+09 31251.82 > [3,] 7.979009e+04 3.125182e+04 64156419.36That format is a bit awkward to key back in. Could you do a dput(m2) so that it becomes easier for us to reproduce the effect?> > det(m2, method="eigenvalues") > [1] -9.331893e+19 > > > > det(m2, method="qr") > [1] 9.331893e+19 > > I think the result based on the eigenvalues is correct. > > Good night. > Armin > > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- > 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 > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._ >-- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._