关于R语言主成分估计的因子分析有函数可以求方差最大的正交旋转和得分吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我找到的函数是主成分估计方法：
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)*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)
}
这个只做到得到载荷矩阵 求大神们帮忙一下啊 用载荷矩阵旋转正交得到新的因子载荷矩阵然后计算得分。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：因子分析 R语言 主成分 最大的 factor 旋转 语言 最大的

psych package

zyz0329 发表于 2012-12-23 19:24
psych package

里面有极大似然法ml、最小残差法minres、主轴因子pa、加权最小二乘wls、广义最小二乘gls，然后我没有找到主成分的啊

公子小七 发表于 2012-12-23 23:07
里面有极大似然法ml、最小残差法minres、主轴因子pa、加权最小二乘wls、广义最小二乘gls，然后我没有找到 ...

有啊 你去论坛下那本 R in action data analysis  上面有

zyz0329 发表于 2012-12-24 09:14
有啊 你去论坛下那本 R in action data analysis  上面有

嗯。。那个包里是有主成分分析和因子分析，然后因子分析里面的不同方法也只有上面我说的那几个，就我书上在因子分析里面还有一种主成分估计的方法 那个包里面木有 我找了 sigh
Unlike  PCA, there are many methods of extracting common factors. They include maximum likelihood (ml), iterated principal axis (pa), weighted least square (wls),generalized weighted least squares (gls), and minimum residual (minres).
嗯。。大概就是这样子

同问

我最近也在做与因子分析相关的东西   对于正交旋转   可以用varimax（）函数找到使方差最大的正交矩阵   详见附件图