我正在使用R功能 prcomp 运行PCA。这是功能:
d2.pca <- prcomp(sel.d2,center = TRUE,scale. = TRUE)
所以变量按比例缩放(这总是必须完成,对吧?)。
这是我原来的加载矩阵:
PC1 PC2 PC3 PC4
var1 0.551 -0.246 0.576 -0.551
var2 -0.545 -0.233 0.736 0.328
var3 -0.427 -0.704 -0.333 -0.460
var4 -0.467 0.625 0.126 -0.613
当我应用varmx rotation时:
varimax(d2.pca$rotation)
输出就是这个:
$loadings
Loadings:
PC1 PC2 PC3 PC4
var1 1
var2 1
var3 -1
var4 -1
PC1 PC2 PC3 PC4
SS loadings 1.00 1.00 1.00 1.00
Proportion Var 0.25 0.25 0.25 0.25
Cumulative Var 0.25 0.50 0.75 1.00
$rotmat
[,1] [,2] [,3] [,4]
[1,] 0.551 0.427 -0.545 0.466
[2,] -0.246 0.704 -0.232 -0.625
[3,] 0.576 0.333 0.736 -0.125
[4,] -0.551 0.461 0.328 0.613
对我来说这看起来很奇怪,如何在varimax旋转后解释加载(“1”和“-1”值)矩阵?任何帮助或建议将不胜感激,我可能会遗漏一些东西......
注意:相关矩阵的KMO为0.6。以防万一,这里是相关矩阵:
var1 var2 var3 var4
var1 1.000 -0.680 -0.491 -0.771
var2 -0.680 1.000 0.697 0.550
var3 -0.491 0.697 1.000 0.166
var4 -0.771 0.550 0.166 1.000