r - 相关矩阵的特征值分解

Question

我有一个相关矩阵：

cor.table <- matrix( sample( c(0.9,-0.9) , 2500 , prob = c( 0.8 , 0.2 ) , repl = TRUE ) , 50 , 50 )
diag(cor.table) <- 1

我尝试做特征值分解：

library(psych)
fit<-principal(cor.table, nfactors=50,rotate="none")

或者

stopifnot( eigen( cor.table )$values > 0 )

在这两种情况下，我都会收到错误：

Error in eigens$values < .Machine$double.eps : 
  invalid comparison with complex values

我究竟做错了什么？

Answer 1

这与您之前提出的问题相同。你需要一个对称矩阵。

set.seed(1)
cor.table <- matrix(sample(c(0.9,-0.9),2500,prob=c(0.8,0.2),repl=TRUE),50,50)
ind <- lower.tri(cor.table)
cor.table[ind] <- t(cor.table)[ind]
diag(cor.table) <- 1

现在，当您尝试使用时，eigen您不会收到错误消息。

your.eigen <- eigen(cor.table)
> summary(your.eigen)
Length Class  Mode   
values    50   -none- numeric
vectors 2500   -none- numeric