r - 绘制 PCA 的所有组件

Question

我想知道我们如何才能将个人绘制在第一个和第二个组件之外，如下所示：

Answer 1

这可能有效：

pc.cr <- princomp(USArrests, cor = TRUE)
pairs(pc.cr$loadings, col=c("red", "green", "blue"))

Answer 2

另一种选择是splom从lattice包中使用。我在这里使用：

1. splom用于矩阵散点图
2. prcomp计算 PCA
3. mclustBIC对数据进行聚类。因为您需要按某些标准对点进行分组。请注意，仅提供颜色只会回收 col 参数。

你可以写例如：

library(lattice)
library(mclust)
dat <- iris[,1:4]
dat.pca <- prcomp(iris[,1:4])
dat.em <- mclustBIC(dat) 
splom(as.data.frame(dat.pca$x), 
      col=summary(dat.em,data=dat)$classification, cex=2,pch='*')