[Rscript]探索性因子分析

setwd("E:/天睿TERADATA/data_analysis")
dta<-read.csv("newdata.csv")

#判断需要提取的公共因子数
fa.parallel(dta,fa="both",n.iter=100,main="Scree plots with parallel analysis")
 #代码中使用了fa="both",从而会同时展示主成分和公共银子分析的结果
 #结果中最大拐角处前有八个变量,从而我们选择八个公共因子

 #若是用相关系数矩阵做
 #cov_dta<-cov(dta)
 #cor_dta<-cov2cor(cov_dta)
 #fa.parallel(cor_dta,n.obs=nrow(dta),fa="both",n.iter=100,main="Scree plots with parallel analysis")
 
#提取公共因子
fa<-fa(dta,nfactors=4,rotate="none",fm="pa")
fa
 
#因子旋转有利于因子解释
#正交旋转
fa.varimax<-fa(dta,nfactors=4,rotate="varimax",fm="pa")
fa.varimax
#斜交旋转允许因子之间相关

#绘制正交旋转的结果
factor.plot(fa.varimax,lables=rownames(fa.varimax$loadings))
fa.diagram(fa.varimax,simple=T)

#因子得分:相比于主成分分析,探索性因子分析不那么关注计算因子得分
fa.varimax$weights

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值