r - PCA のすべてのコンポーネントをプロットする

Question

次のように、個人を第 1 および第 2 成分以上にプロットする方法を知りたいです。

Answer 1

これはうまくいくかもしれません：

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

Answer 2

別のオプションは、パッケージsplomから使用することです。latticeここで私は使用します：

1. splom行列散布図の場合
2. prcompPCAを計算する
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='*')