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)

这应该给出正确的输出。

ggbiplot PCA 中的补充变量

Supplementary variables in ggbiplot PCA

r

pca

ggbiplot