r - ggbiplot PCA 中的补充变量

Question

我正在尝试用 ggbiplot 绘制 PCA 结果，如何绘制补充变量？我找到了关于 MCA 结果的讨论，但我也想要箭头......

 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)

那应该给出正确的输出。