关于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$values[i])*eig$vectors[,i]
h<-diag(A%*%t(A))
rowname<-c("SS loadings","Proportion Var","Cumulative Var")
B<-matrix(0, nrow=3, ncol=m,
dimnames=list(rowname, colname))
for (i in 1:m){
B[1,i]<-sum(A[,i]^2)
B[2,i]<-B[1,i]/sum_rank
B[3,i]<-sum(B[1,1:i])/sum_rank
}
method<-c("Principal Component Method")
list(mehod=method, loadings=A,
var=cbind(common=h, spcific=diag_S-h), B=B)
}
这个只做到得到载荷矩阵 求大神们帮忙一下啊 用载荷矩阵旋转正交得到新的因子载荷矩阵然后计算得分。

求解

郭晓宇 发表于 2012-12-30 03:48
求解

直接varimax(A) A是矩阵嗯。。。我惊奇发现如此简单