R因子、聚类

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")