# 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 #主成分计算系数
备注: