使用这些数据:
head(USArrests)
nrow(USArrests)
我可以这样做一个PCA:
plot(USArrests)
otherPCA <- princomp(USArrests)
我可以得到新的组件
otherPCA$scores
和分量解释的方差比例
summary(otherPCA)
但是,如果我想知道哪些变量主要由哪些主成分解释呢?反之亦然:例如 PC1 或 PC2 主要由murder
?我怎样才能做到这一点?
例如,我可以说 PC1 80% 由murder
or解释assault
吗?
我认为负载在这里对我有帮助,但是它们显示了方向性,而不是我所理解的解释的方差,例如
otherPCA$loadings
Loadings:
Comp.1 Comp.2 Comp.3 Comp.4
Murder 0.995
Assault -0.995
UrbanPop -0.977 -0.201
Rape -0.201 0.974