[R] r: eviews and r // eigen analysis
    Clark Allan 
    Allan at STATS.uct.ac.za
       
    Mon Mar 14 11:59:45 CET 2005
    
    
  
hi all
i have a question that about the eigen analysis found in R and in
eviews.
i used the same data set in the two packages and found different
answers. which is incorrect?
the data is:
aa ( a correlation matrix)
1	0.9801	0.9801	0.9801	0.9801
0.9801	1	0.9801	0.9801	0.9801
0.9801	0.9801	1	0.9801	0.9801
0.9801	0.9801	0.9801	1	0.9801
0.9801	0.9801	0.9801	0.9801	1
now
> svd(aa)
$d
[1] 4.9204 0.0199 0.0199 0.0199 0.0199
$u
           [,1]          [,2]          [,3]          [,4]       [,5]
[1,] -0.4472136  9.283999e-18  1.939587e-17 -2.101554e-15  0.8944272
[2,] -0.4472136  8.089763e-01  7.115435e-17  3.091235e-01 -0.2236068
[3,] -0.4472136  2.178563e-02  8.226578e-18 -8.657513e-01 -0.2236068
[4,] -0.4472136 -4.153810e-01 -7.071068e-01  2.783139e-01 -0.2236068
[5,] -0.4472136 -4.153810e-01  7.071068e-01  2.783139e-01 -0.2236068
$v
           [,1]        [,2]          [,3]       [,4]       [,5]
[1,] -0.4472136  0.00000000  0.000000e+00  0.0000000  0.8944272
[2,] -0.4472136  0.80897632 -4.976488e-17  0.3091235 -0.2236068
[3,] -0.4472136  0.02178563  1.077421e-17 -0.8657513 -0.2236068
[4,] -0.4472136 -0.41538097 -7.071068e-01  0.2783139 -0.2236068
[5,] -0.4472136 -0.41538097  7.071068e-01  0.2783139 -0.2236068
the results from Eviews is:
eigenvectors = (note that Eviews arranges the eigen vectors in ascending
order of the size of the eigen values)
 0.305963	-0.266024	-0.219066	-0.766569	 0.447214
 0.315257	-0.603871	 0.075772	 0.574640	 0.447214
 0.461958	 0.726861	 0.199385	 0.136062	 0.447214
-0.482609	 0.185567	-0.700705	 0.204124	 0.447214
-0.600569	-0.042534	 0.644614	-0.148258	 0.447214
eigen values =
 0.019900
 0.019900
 0.019900
 0.019900
 4.920400
NOTE THAT THE EIGEN VALUES ARE THE SAME BUT THE EIGEN VECTORS ARE NOT!!!
why is this so?
if one sets aa=
 1.000000	 0.500000	 0.250000
 0.500000	 1.000000	 0.350000
 0.250000	 0.350000	 1.000000
then both packages give the same answers.
    
    
More information about the R-help
mailing list