如题  数据：
程序如下
data<-read.csv("年度数据.csv",header = TRUE)
cor(data[,-1])
eigen(cor(data[,-1]))
PCA <- princomp(data[,-1],cor=T)
summary(PCA,loadings = TRUE)
PCA$loadings
predict(PCA)
biplot(PCA)
biplot(PCA, choices = 2:3)
screeplot(PCA,type="lines")
eigen(cor(data[,-1]))$values

PCA$scores[,1:3]  #每个样本的主成分得分      
PCs1 <- as.matrix(PCA$scores[,1:2])
PCs2 <- as.matrix(PCA$scores[,2:3])
plot(PCs1)
abline(h=0,v=0,lty=3)
text(PCs1,label=data[,1],pos=1.8,adj=0.5,cex=0.85)
plot(PCs2)
abline(h=0,v=0,lty=3)
text(PCs2,label=data[,1],pos=1.8,adj=0.5,cex=0.85)