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")
R 語言 因子分析
因子分析 options digits 2 covariances ability.cov cov 將協方差矩陣轉化為相關係數矩陣 correlations cov2cor covariances 第一步 判斷需提取的公共因子數 library psych covariances ability.c...
R做因子分析
factanal x,factors,data null,covmat null,n.obs na,subset,na.actionstart null,scores c none regression bartlett rotation varimax control null,其中x是資料的公式...
因子分析 因子得分
因子分析的最後一步了，悲傷 在因子分析中，人們一般關注的重點是估計因子模型的引數 載荷矩陣 有時公共因子的估計 即所謂的因子得分 也是需要的，因子得分可以用於模型診斷，也可以作為下一步分析的原始資料。需要指出的是，因子的得分計算並不是通常意義下的引數估計，他是對不可觀測的隨機變數 取值的估計 通常可...