主成分分析R语言

1、第一题data=read.table(q1.txt,head=T)head(data)#前六行数据x3K.S10100590350340启Q40加011-007O5m10Q口吐三Q51oqmQ52n三丄0490.350.421.000.330.190.360.730.240.4-0.030日4051。霁10002吕Q4石0沖0ivais0630理9Q丄0殆丄DO0孑4013。知05203咅04咗0M10003；30150001attach(data)options(digits=2)data=as.matrix(data)#转化为矩阵v1-eigen(data)#计算相关系数矩阵特征值和特征向

2、量v113.791.521.110.310.720.590.530.330.20.21vector3f3【沁【3fe;【八【4【川【Vfl：山1,0.355-0.176-0.4330216-0丄590.1310.1540233-0.0760.625540.405-0.149-0.182-O.O97-0.065-0.0960.181-0.833-0.129-0.120163f03610.4340.0450.1S30.2550.0340.109-0.0530仙0.023700.346-0.1090.131-0.4140.274-0.3-510.5320.3S2-0.133-0.112095门0.3

3、18-0-447一0030.438-0.130-0.10-S0.255O.OSS-0.625S36r0.353-0.034-0.153-0.3610.054-0.321-0.7690.0950.05-50.102620.3190.4230.丄040.2450.4200.293-0.1530.033-0.5DS-0.0569S弘0.277-0.0710.389-CL吐60-D.3220.665-0.0200.070.0-30.000933f0.2230.3560.3520.：45-.669-0.4010.0270.056-0.1480.0717-910.0.075-0.4340.6230.232

4、D.255-0.105-0.101-0.130.030.40261zhuchengfenx-zhuchengfen$loadings；x#主成分载荷Loadings:Comp.1Ccnp.2Conp.3Comp.4Comp.5Comp.6Comp.Camp.5Comp.?Con?).10X1-0.355-0.1760.493-0.216-0.159O.lrl-0.13-0.2330.623注-0.405-C.140.15-0.1S10.533-0.120 x3-0.361C.34-0.1330.253-0.1030他-0.346-C.109-0.1310.4140.274-0.391-0.53

5、2-D.3-32-0.1S3-0.112k5-0.318-a.447-0.43B-0.130O.lOE-D.255-0.626x6-0.3530.1530.361-0.3210.76?0.103x7-0.319C.423-0.104-0.2450.400.2930.15c-0.595君-0応-0.3890珏0-0.3220&5-O.；231.3龙-0.352-0.245-0.E69-0.4C1-0.143xlD-a.484-0.623-0.2920.255-0.1050.101D.1430.403Comp.1CQinp.2Comp.3Zonp.4Comp.5Ocmp.6Corp.7Comp.=

6、匚omp.9Comp.1055loadings1.01.01.01.01.01,01.31.01ProportionVar0.LC.l0.10.10.10,10.1hl0CmmlativeVar0.1C.20.30.40.50.60.79.E0y-zhuchengfen$scores;y#主成分得分NULL因为只有相关系数矩阵，所以没有主成分的得分screeplot(zhuchengfen,type=lines)筈LIE匸ER根据累计方差贡献率和Kaiser准则，应该选择三个主成分，即PC1、PC2、PC3。累计方差贡献率为64%。PC1、PC2、PC3对方差贡献率分别是38%、15%、11%

7、。结果如下：运动项目第一主成分第二主成分第三主成分100米-0.355-0.1760.493跳远-0.405-0.1490.182铅球-0.3610.434跳咼-0.346-0.109-0.131400米-0.318-0.447110米跨栏-0.3530.153铁饼-0.3190.423-0.104撑杆跳高-0.277-0.389标枪-0.2230.356-0.3521500米-0.484-0.623方差贝献率0.380.150.11累积方差贡献率0.380.530.64解释：前三个主成分解释了大部分信息；第一次主成分测量的是综合运动能力，第二个主成分代表了爆发力臂力，第三主成分测量的是爆发性

8、腿力。附录：library(psychjzhu|这是psych软件包中计算主成分的结果，主成分载荷与princomp得到的结果有很大的差异，并且这个结果似乎更容易解释第二题library(RODBC)data=odbcConnectExcel(G:/R/q2.xls)data1=sqlFetch(data,Sheet1)close(data)#关闭通道attach(data1)head(datal)#查看前六行数据data2=data1,-1#去掉第一列数据options(digits=2)#数据均保留小数点后两位q2x-q2$loadings;x#主成分载荷Loadings:Comp.1Co

9、mp.2Comp.3Comp.4Ccnip.SCoiuip石Ccrnp.：100rr.口-0.3650.4900.2560.319-0.23丄0.20200rr.口-0.36S0.5370.230-0.711-0.丄09400rr.口-0.3520.247-0.515-0.3470.5720.1-910.203BOOrr.(min)-0.3B5-0.155-0.SEE-0.620-0.3151500rr.(mi对-0.3B3-o.seoC.430-0.2310.6D33000rr.(mi对-0.3B3-0.3450.153C.330.4-63-0.5D-BMarathon(rain)-0.367-o.ses0.454-0.-0.1310.142Comp工匚口啤2匚削野.3CQiap呛匚sop5Coiop655loadingsJ.曲:L3:L3:La1Q01皿1Q0ProportionV己工Q：l理Q工呛Q：1.呛Q工呛Q.理0Q.理CumulativeV己工：！理Q4Q,弓2yorder(y,1)#排序TOC o 1-5 h z155123641231335162646517425234322522109133511253037623444731501527433244043333451B45223173434321