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