R语言 主成分分析 代码

代码

#主成分分析
getwd()
library(factoextra)
install.packages("psych")
library(psych)
df<-read.csv("1.csv",header = T,row.names = 1)
head(df)
#数据标准化
d <- scale(df)
d
#确定提取主城分数，绘制碎石图
fa.parallel(d,fa = 'pc')
#提取主成分
p <- principal(d, nfactors = 2, rotate = "none")
#看看主成分结果有什么
names(p)
p$values#输出特征根
#前两个特征根的方差贡献率
sum(p$values[1:2])/sum(p$values)

#主成分载荷矩阵
p$loadings

#主成分得分
p$scores