[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