d<-read.table("D:/R/Test/主成分.txt") library(psych) KMO(cor(d))#KMO检验值为0.35,不能作成分分析 cortest.bartlett(d)#bartlett球形度检验 #两种计算主成份1 princomp 基础包 test<-princomp(d,cor=TRUE) test$loadings pca_data <- predict(test)#主成份得分 screeplot(test,type="lines") #两种计算主成份2,这组与Spss 的计算是一致的 psych包 #fa pc显示特征值 fa.parallel(d,fa="pc",n.iter=100, show.legend=FALSE,main="Screen plot with parallelanalysis") test<-principal(d,scores = T,nfactors =2,rotate="none") test$loadings test$scores
|