主成分分析 PCA

各パラメータは関係があったりすることがあるので、無関係な主成分にしてみる。

pca_score <- scale(data0) %*% eigen(cor(data0))$vectors *sqrt(nrow(data0)/(nrow(data0) - 1))
library(rgl)
plot3d(pca_score, type="n", xlab="PC1", ylab="PC2", zlab="PC3", col=cols)
text3d(pca_score, texts=seq(nrow(data0)))

plot(as.data.frame(pca_score), pch=16, col=cols)

それっぽく突出してはみ出ているやつが現れた。