[R] eigen bug?
Christian Hennig
fm3a004 at math.uni-hamburg.de
Wed Apr 12 17:36:56 CEST 2000
Sorry if this bug is fixed in the meantime, but what is this??
> m
[,1] [,2] [,3]
[1,] 1 2 1
[2,] 0 2 0
[3,] 1 2 1
> eigen(m)
$values
[1] 2 2 0
$vectors
[,1] [,2] [,3]
[1,] -4.194304e+06 0.7071068 -0.7071068
[2,] -4.656613e-10 0.0000000 0.0000000
[3,] -4.194304e+06 0.7071068 0.7071068
> eigen(m, symmetric=T)
$values
[1] 3.709275 1.193937 -0.903212
$vectors
[,1] [,2] [,3]
[1,] 0.2331920 0.8876503 -0.3971125
[2,] 0.7392387 -0.4271323 -0.5206574
[3,] 0.6317813 0.1721479 0.7557893
As far as I can see, m is symmetric. Why do these two not coincide?
(The first is correct.)
Best regards,
Christian
***
Christian Hennig
Fachbereich Mathematik
(Schwerpunkt Mathematische Statistik und Stochastische Prozesse,
Zentrum fuer Modellierung und Simulation)
Bundesstrasse 55
D-20146 Hamburg
Tel: x40/42838 4907, privat x40/631 62 79
hennig at math.uni-hamburg.de
http://www.math.uni-hamburg.de/home/hennig/
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list