ggbiplot中的辅助变量

Question

我试着用ggbiplot绘制PCA结果，我如何绘制补充变量？我为MCA结果找到了this discussion，但是我也想要箭头.

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

​

​

此外，这段代码给了我一个错误：

 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]' 

Answer 1

我尝试过您的代码，没有复制您的错误，但是有其他问题。我搜索了PCA()，发现用来做PCA的包是FactoMineR。在查看了文档之后，我还将scale.更改为scale.unit，将quanti.sup更改为quali.sup，给出了分类变量所在的正确列。

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

这应该会给出正确的输出。

​

​