2016-06-01 9 views
2

ggbiplotでPCA結果を描画しようとしていますが、どのように補助変数を描画できますか? 私はMCA結果のthis discussionを見つけましたが、私は同様に矢印を持っているしたいと思います...ggbiplotの補助変数PCA

data(wine) 
wine.pca <- PCA(wine, scale. = TRUE, quanti.sup = c(4,5)) 
plot(wine.pca) 
ggbiplot(wine.pca) 

With plot.PCA With gbiplot

また、このコードは私にエラーを与える:

1: In sweep(pcobj$ind$coord, 2, 1/(d * nobs.factor), FUN = "*") : 
    STATS is longer than the extent of 'dim(x)[MARGIN]' 
2: In sweep(v, 2, d^var.scale, FUN = "*") : 
    STATS is longer than the extent of 'dim(x)[MARGIN]' 

答えて

1

私はあなたのコードを試したが、あなたのエラーを再現しなかったが、他の問題があった。私はgoogled PCA()と、PCAを行うために使用されるパッケージがFactoMineRだったことがわかった。ドキュメントを見た後、私はまた、カテゴリ変数はである正しい列を与え、quali.supscale.unitquanti.supscale.を変更しました。正しい出力が得られます

library(FactoMineR) 
data(wine) 
wine.pca <- PCA(wine, scale.unit = TRUE, quali.sup = c(1,2)) 
plot(wine.pca) 
ggbiplot(wine.pca) 

を。

enter image description here