problem with eigen() ??

3 messages · Paul Rathouz, Don MacQueen, Thomas Lumley

Original

1

3

Paul Rathouz

Fri, Oct 17, 2003 10:11 AM #

Hi All -- I am running a dual processor G4 with R v 1.6.2 for X11 which I
installed from fink.  Look what the eigen() function produces.  Am I
thinking about this wrong, or should  all of the eigen values be 1? -- pr

-----------------------------

[,1] [,2] [,3]
[1,]    1    0    0
[2,]    0    1    0
[3,]    0    0    1

[,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    1    0    0
[3,]    0    0    1    0
[4,]    0    0    0    1

[,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

[1] 1 1 1 1 0

[1] 1 1 1 1 1 1 1 1 1 1 1 1 0 0

[1] 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0

[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0

[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0

-----------------------------

==========================================================================
Paul Rathouz, Asst. Professor        ph   773-834-1970
Dept. of Health Studies, Rm. W-264   fax  773-702-1979
University of Chicago                prathouz@health.bsd.uchicago.edu
5841 S. Maryland Ave. MC 2007
Chicago, IL  60637

Don MacQueen

Fri, Oct 17, 2003 10:33 AM #

Whether or not this is a Mac issue, I would suggest upgrading R to 1.8.0.

-Don

[1] 1 1 1 1 1

_                       
platform powerpc-apple-darwin6.7.5
arch     powerpc                 
os       darwin6.7.5             
system   powerpc, darwin6.7.5    
status   Patched                 
major    1                       
minor    8.0                     
year     2003                    
month    10                      
day      13                      
language R

At 10:12 AM -0500 10/17/03, Paul Rathouz wrote:

Hi All -- I am running a dual processor G4 with R v 1.6.2 for X11 which I
installed from fink.  Look what the eigen() function produces.  Am I
thinking about this wrong, or should  all of the eigen values be 1? -- pr

-----------------------------

 diag(rep(1,3))

     [,1] [,2] [,3]
[1,]    1    0    0
[2,]    0    1    0
[3,]    0    0    1

 diag(rep(1,4))

     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    1    0    0
[3,]    0    0    1    0
[4,]    0    0    0    1

 diag(rep(1,5))

     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

 > eigen(diag(rep(1,5)))$val

[1] 1 1 1 1 0

 eigen(diag(rep(1,14)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 0 0

 eigen(diag(rep(1,15)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0

 eigen(diag(rep(1,16)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

 eigen(diag(rep(1,17)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0

 eigen(diag(rep(1,20)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

 eigen(diag(rep(1,22)))$val

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0

 q("no")

-----------------------------

==========================================================================
Paul Rathouz, Asst. Professor        ph   773-834-1970
Dept. of Health Studies, Rm. W-264   fax  773-702-1979
University of Chicago                prathouz@health.bsd.uchicago.edu
5841 S. Maryland Ave. MC 2007
Chicago, IL  60637

_______________________________________________
R-SIG-Mac mailing list
R-SIG-Mac@stat.math.ethz.ch
https://www.stat.math.ethz.ch/mailman/listinfo/r-sig-mac

--------------------------------------
Don MacQueen
Environmental Protection Department
Lawrence Livermore National Laboratory
Livermore, CA, USA

Fri, Oct 17, 2003 1:49 PM #

On Fri, 17 Oct 2003, Paul Rathouz wrote:

They should be 1, and if you were using more recent version of R they
would be. I seem to recall that there were a number of problems with
getting linear algebra code linked in properly back then.

	-thomas