From dweeks@watson.hgen.pitt.edu  Fri Feb 14 21:37:02 2003
From: dweeks@watson.hgen.pitt.edu (dweeks@watson.hgen.pitt.edu)
Date: Fri Feb 14 21:37:02 2003
Subject: [Rd] eigen() error: R Version 1.6.1 on Mac OS X (PR#2550)
Message-ID: <200302142037.h1EKb235007465@pubhealth.ku.dk>

Consider this matrix:

> sg
     X1    X2    X3    X4    X5
1 3.240 2.592 2.592 2.592 2.592
2 2.592 3.240 2.592 2.592 2.592
3 2.592 2.592 3.240 2.592 2.592
4 2.592 2.592 2.592 3.240 2.592
5 2.592 2.592 2.592 2.592 3.240

If I compute the eigenvalues of the 'sg' matrix using R Version 1.5.0
(2002-04-29) under Linux (or using Version 1.4.0  (2001-12-19) under
Solaris), I obtain:

> eigen(sg)
$values
[1] 13.608  0.648  0.648  0.648  0.648

I also obtain the same eigenvalues under R Version 1.5.1  (2002-06-17) under
Mac OS X 10.2.3.

However, if I compute the eigenvalues of the 'sg' matrix using R Version
1.6.1  (2002-11-01) under Mac OS X 10.2.3, I obtain a different answer:

> eigen(sg)
$values
[1] 13.123930  0.648000  0.648000  0.648000 -1.732930
> det(sg)
[1] 2.399354
> det(sg,method="eigenvalues")
[1] -7.312317


Thanks,
 -- Dan --



[Rd] eigen() error: R Version 1.6.1 on Mac OS X (PR#2550)
dweeks@watson.hgen.pitt.edu