R语言主成分分析步骤
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包#fapc显示特征值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
页:
[1]