代码：
sports.cor<-read.table("sports.dat",head=T)
sports.cor
x<-as.matrix(sports.cor)
princomp.x<-princomp(covmat=x,cor=T)
summary(princomp.x)
screeplot(princomp.x,type="lines")
fa.x<-factanal(factors=2,covmat=x,rotation="varimax")
fa.x
apply((fa.x$loadings)^2,1,sum)
结果：
Factor1
Factor2
Communalities
100 米
0.701
0.22
0.539
跳远
0.733
0.294
0.624
铅球
0.197
0.939
0.920
跳高
0.508
0.297
0.347
400 米
0.723
0.525
110 米跨栏
0.539
0.281
0.369
铁饼
0.164
0.741
0.576
撑竿跳高
0.378
0.192
0.180
标枪
0.108
0.444
0.209
1500 米
0.24
-0.135
0.076
Cumulative Var
0.238
0.436
评论13