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是數據的公式、矩陣或數據框，factors是因子的個數，scores是因子得分的方法，scores ="regression"表示迴歸方法計算因子得分，scores ="Bartlett"表示用Bartlett方法計算因子得分，缺省時不計算因子得分，rotation表示旋轉，缺省時爲方差最大旋轉。
例子：
主成分法的R程序（程序名：factor.analy1.R）

factor.analy1<-function(S, m){
p<-nrow(S); diag_S<-diag(S); sum_rank<-sum(diag_S)
rowname<-paste(“X”, 1:p, sep="")
colname<-paste(“Factor”, 1:m, sep="")
A<-matrix(0, nrow=p, ncol=m,
dimnames=list(rowname, colname))
eig<-eigen(S)
for (i in 1:m)
A[,i]<-sqrt(eigvalues[i])∗