全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
7094 6
2012-12-23
我找到的函数是主成分估计方法:
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)
}
这个只做到得到载荷矩阵 求大神们帮忙一下啊 用载荷矩阵旋转正交得到新的因子载荷矩阵然后计算得分。
二维码

扫码加我 拉你入群

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

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

全部回复
2012-12-23 19:24:08
psych package
二维码

扫码加我 拉你入群

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

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

2012-12-23 23:07:29
zyz0329 发表于 2012-12-23 19:24
psych package
里面有极大似然法ml、最小残差法minres、主轴因子pa、加权最小二乘wls、广义最小二乘gls,然后我没有找到主成分的啊
二维码

扫码加我 拉你入群

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

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

2012-12-24 09:14:56
公子小七 发表于 2012-12-23 23:07
里面有极大似然法ml、最小残差法minres、主轴因子pa、加权最小二乘wls、广义最小二乘gls,然后我没有找到 ...
有啊 你去论坛下那本 R in action data analysis  上面有
二维码

扫码加我 拉你入群

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

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

2012-12-24 14:56:59
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).
嗯。。大概就是这样子
二维码

扫码加我 拉你入群

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

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

2012-12-30 03:47:05
同问
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群