R程序求正交旋转矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

y1 y2 y3 y4  y5 之间的相关系数矩阵为
          y1        y2        y3        y4        y5
y1 1.0000000 0.8680120 0.9738880 0.9485825 0.9683039
y2 0.8680120 1.0000000 0.9494782 0.9714577 0.9348399
y3 0.9738880 0.9494782 1.0000000 0.9866021 0.9827903
y4 0.9485825 0.9714577 0.9866021 1.0000000 0.9883761
y5 0.9683039 0.9348399 0.9827903 0.9883761 1.0000000

经主成分分析后，用两个因子来描述他们，现要求他们之间的最大方差旋转矩阵，我下面的编程错在哪啊
ot=c(1.00,0.868,0.974,0.949,0.968,
0.868,1.000,0.949,0.971,0.935,
0.974,0.949,1.00,0.987,0.983,
0.949,0.971,0.987,1.00,0.988,
0.968,0.935,0.983,0.988,1.00)
R=matrix(ot,nrow=5)
fa=factanal(covmat=R,factors=2)
fa
另若要求这两个因子在分别在 y1  y2  y3  y4  y5 上的得分系数，R语言怎么写呢？？
谢谢各位~~~~
附：
y1=c(7098.00,7498.00,7848.00,8254.00,8761.00,8801.12,8951.32,9325.03,9683.90,10200.00,11000.00,12360.74),
y2=c(152.7,165.68,181.4,204.26,272.53,319.87,491.92,571.35,610.23,652.36,640.37,628.92),
y3=c(646,719.33,775.00,805.78,993.82,1079.80,1216.34,1308.41,1420.00,1611.30,1612.41,1913.80),
y4=c(12.45,13.64,16.13,18.25,22.75,24.65,30.89,36.08,40.36,47.40,50.27,51.32),
y5=c(3688.20,4034.96,4551.15,4950.84,5408.76,6250.81,7450.27,9143.95,10297.00,12001.16,13698.15,14900.93)

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：旋转矩阵 R程序 factors matrix factor 旋转矩阵 程序

yr<-cor(y)
ans<-factor.analy1(yr,2)
vm<-varimax(ans$loadings,normalize=F)
vm

计算因子得分
fa<-factanal(~.,factor=2,data=y,scores="regression")
试试吧