转 R 主成分分析

# 1、数据准备

tt<-USJudgeRatings 

# 2、主成分个数判断

# 方法1、由于我们是通过特征值大小来找主成分的次序,所以KAISER-Harris准则保留特征值大于1的主成分—-碎石图
install.packages("psych")
library(psych)
fa.parallel(USJudgeRatings[,-1],fa="pc",n.iter = 100,show.legend = F)
abline(h=1,col="green")
# 方法2、信息保留比例，进行判断
pr <- princomp(tt, cor = T,score=T) 
summary(pr, loadings = F)

# 3、得出主成分结果

principal(tt,nfactors = 2,scores = T)
pc$scores  #主成分得分
pc$weights #主成分计算系数

备注：