[Rd] possible bug in det using method="qr" (PR#1244)

Peter Dalgaard BSA p.dalgaard@biostat.ku.dk
05 Jan 2002 03:10:15 +0100


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.36

That 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._