1.因子分析
a=read.csv("F:/2014.csv");head(a)
x=scale(a[,-1])
library(mvstats)
(Fac1=factpc(x,4,rot="varimax"))
(Fac0=factpc(x,4))#主成分法,因子分析
Fac0$Vars#求方差及贡献率
(Fac1=factpc(x,4,rot="varimax"))
Fac1$Vars#求旋转后方差及贡献率
Fac1$Rank#排名
2.聚类
a=read.csv("F:/聚类.csv",header=T);head(a);attach(a)
x=scale(a[,c('偿债能力均值','盈利能力均值','成长能力均值','营运能力均值')])
row.names(x)=a$公司名称
hc=hclust(dist(x),'comp')
plot(hc)
re1<-rect.hclust(hc, k=4, border="red")