第六章 因子分析

load("D:\\文件\\2022下（大三上）\\classR\\pollu_ch6.RData")

pollu=pollu.ch6

head(pollu)

#----------第六章 因子分析------------

# 主成分分析-载荷阵不具有良好可分性，不易找出实际意义

# 因子分析-初始载荷矩阵

fa=factanal(~.,factors=2,data=pollu) #1.因子分析

fa

fa=factanal(~.,factors=2,data=pollu,scores='regression')

# 2.因子分析的有效性：贡献率、累计贡献率、共同度

# 3.因子命名

Factor=fa$scores # 4.公因子得分

Score=-0.490*Factor[,1]+0.397*Factor[,2] # 5.公因子的应用 ①综合评价公式

Score

sort(Score,decreasing=T)

rank(-Score)

library(factoextra) # ②利用公因子进行聚类

fviz_nbclust(Factor,kmeans,method="wss")

km=kmeans(Factor,3)

km

fviz_cluster(km,data=Factor)