seiseki <- read.csv("seiseki.csv", head=T)
result <- prcomp(seiseki, scale=T)
summary(result)
result
round(result$rotation, 3)
biplot(result)
plot(result$x)
plot(result$x[1:100,], col="blue", xlim=c(-6,6) ylim=c(-3,3))
plot(result$x[1:100,], col="blue", xlim=c(-6,6), ylim=c(-3,3))
par(new=T)
plot(result$x[101:166,], col="red", xlim=c(-6,6), ylim=c(-3,3))
dev.copy(png, file="seiseki-chart000.png")
dev.list()
dev.off()
q()
q()
