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(eig