第五章 主成分分析

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

#----------第五章 主成分分析---------

cor(fashion)      #1

res.pc=princomp(~sale+price+collect+Baidu+review+reviewp,data=fashion,cor=T)  #2.主成分分析的实现
summary(res.pc,loading=T)

score=predict(res.pc)           #3.每条观测对应的主成分得分
head(score)
score=predict(res.pc)[,1:2]

screeplot(res.pc)              #碎石图     4.保留主成分的合理性：累计贡献率、共同度、碎石图
#or
screeplot(res.pc,type="line")

load=res.pc$loadings[,1:2]         # 共同度
c<-1.8852471^2*load[,1]^2+1.735260^2*load[,2]^2
c

plot(score,type="n")                           # 5.主成分的应用1--利用第一、二个主成分进行分类(对各样本)
text(score,labels = row.names(fashion),cex=0.7)

load=res.pc$loadings[,1:2]                  # 主成分的应用2--利用主成分载荷和对指标进行分组
plot(load,type="n")
text(load,labels=names(fashion)[-c(2,8)])
#另
biplot(res.pc)  # (绘制前两个主成分散点图)