package psychのプリンシパルを使用してPCA出力から "Proportion Var"ラインを取得する方法は？

Question

以下の "pca $ loadings"のオブジェクトまたはベクトルのような "Proportion Var"行を取得して、PCAグラフィックの値を使用したいとします。

私は次のようでした：

library(psych) 

data(iris) 

pca <- principal(iris[1:4],nfactors = 2,rotate="varimax",scores=TRUE) 

pca$loadings 

どのように私は、比例VARを取得する必要がありますか？

Answer 1

もう一つの方法は、手動で計算するのですが、最初にあなたはあなたの例では4であるあなたが持っている変数の数、などnfactorsを指定する必要があり

すべての要因（例えば、すべての軸）を抽出する必要があります。 、そして、

pca <- principal(iris[1:4],nfactors = 4,rotate="varimax",scores=TRUE) 

あなたはコンポーネントごとの二乗負荷（この場合はRC）の合計を取得することにより、比例分散を計算し、合計することによって、それらを分割することができた後pca$loadingsを抽出します。だから、このようなようにする必要がありますあなたがこれで行うことができる二乗荷重：

colSums(pca$loadings[ , ]^2)/sum(pca$loadings[ , ]^2) 

これは、pca$loadingsの割合の分散と同じ情報を提供するはずですが、すべてのコンポーネント（この場合はRC1〜RC4）に適用されます。