构造主成分的表达式(综合指标),并做分析test_test <- read.table("D:/R/R Code/5/Chap7/test_score.csv", sep=",", header=T)R <- round(cor(test), 3) <- princomp(test, cor=T)summary(test_PCA, loadings=T)# 找出第一个主成分(PC1)上得分较高的学生test[c(6,7,45,30,49),]# 找出第二个主成分(PC2)上得分较高的学生test[c(26,33,8),]# 获取这些学生在主成分上的得分samplePC <- (round(test_PCA$scores,3))[c(6,7,45,30,49,26,33,8),]rownames(samplePC) <- c(6,7,45,30,49,26,33,8)samplePC# 也可以使用另一种方式获取这些学生在主成分上的得分samplePC2 <- round(predict(test_PCA),3)[c(6,7,45,30,49,26,33,8),]rownames(samplePC2) <- c(6,7,45,30,49,26,33,8)samplePC2
最新发布