典型相关分析Word下载.doc
《典型相关分析Word下载.doc》由会员分享,可在线阅读,更多相关《典型相关分析Word下载.doc(6页珍藏版)》请在冰点文库上搜索。
d11.1生理指标与训练指标之间的典型相关性。
生理指标:
体重(x1)、腰围(x2)、脉搏(x3);
训练指标:
引体向上次数(y1)、起坐次数(y2)、跳跃次数(y3)。
>
X<
-read.table("
clipboard"
header=T)
R<
-cor(X)
R
x1x2x3y1y2y3
x11.00000.8702-0.36576-0.3897-0.4931-0.22630
x20.87021.0000-0.35289-0.5522-0.6456-0.19150
x3-0.3658-0.35291.000000.15060.22500.03493
y1-0.3897-0.55220.150651.00000.69570.49576
y2-0.4931-0.64560.225040.69571.00000.66921
y3-0.2263-0.19150.034930.49580.66921.00000
R11<
-R[1:
3,1:
3];
R12<
3,4:
6];
R21<
-R[4:
6,1:
R22<
6,4:
6]
A<
-solve(R11)%*%R12%*%solve(R22)%*%R21#A=(R11)-1R12(R22)-1R21
ev<
-eigen(A)$values#特征值
sqrt(ev)#典型相关系数
[1]0.795610.200560.07257
以上过程是一步一步计算的,接下来我们使用R自带的典型相关函数:
xy<
-scale(X)#数据标准化
ca<
-cancor(xy[,1:
3],xy[,4:
6])#典型相关分析
ca$cor#典型相关系数
ca$xcoef#x的典则载荷
[,1][,2][,3]
x1-0.17789-0.432300.04381
x20.362330.27086-0.11609
x3-0.01356-0.05302-0.24107
ca$ycoef#y的典则载荷
y1-0.08018-0.086160.29746
y2-0.241810.02833-0.28374
y30.164360.243680.09608
典型变量的系数载荷并不唯一,只要是它的任意倍数即可,所以每个软件得出的结果并不一样,而是相差一个倍数。
R自带的典型分析函数cancor()并不包括对典则相关系数的假设检验,为了方便,使用自编典型相关检验函数cancor.test()。
cancor.test(xy[,1:
6],plot=T)
$cor
$xcoef
$ycoef
$xcenter
x1x2x3
-5.551e-18-1.943e-171.821e-17
$ycenter
y1y2y3
-2.776e-173.331e-173.365e-17
cancortest:
rQP
[1,]0.7956116.254960.06174
[2,]0.200560.671850.95475
[3,]0.072570.071280.78948
经检验不拒绝原假设,即认为在0.05的水平上没有一个典型相关是显著的。
从典型相关图上也可以看出效果不是很理想,所以就不需要做进一步的典型相关分析了。
例2:
d11.2广东省能源消费量与经济增长之间的典型相关分析。
library(mvstats)
cancor.test(X[,1:
4],X[,5:
10],plot=T)
[1]0.99900.95490.73730.4267
[,1][,2][,3][,4]
x1-0.013980.2627-0.1634-0.05500
x20.118870.43591.5137-0.02025
x30.09036-0.7627-1.6045-0.96536
x40.036870.17240.19851.04168
[,1][,2][,3][,4][,5][,6]
y1-0.059011.77121.057001.5687-2.690-1.4398
y2-0.22982-5.07973.721011.568912.6366.4680
y30.058151.21930.091710.9209-1.509-3.9680
y40.327234.4700-2.21300-0.2103-12.580-6.8854
y50.08586-0.37870.95434-0.85822.1510.7645
y60.05439-1.9202-3.59565-2.93451.8844.9775
x1x2x3x4
3.418e-16-4.382e-17-8.035e-17-9.641e-17
y1y2y3y4y5y6
-5.551e-173.871e-17-9.641e-177.888e-173.494e-168.327e-17
rQP
[1,]0.9990120.6487.438e-15
[2,]0.954939.2645.851e-04
[3,]0.737310.3452.416e-01
[4,]0.42671.9095.914e-01
经检验,在0.05水平下有三个典型相关是显著的,即需要三个典型变量,于是可得出前三对典型变量的线性组合是:
对结果进行分析:
1)由于r1=0.999,说明u1、v1之间具有高度的相关关系(尤其是绝对值较大的权系数),而各自的线性组合中变量的系数大部分都为正,因此一般来说,能源消费越多,经济增长就越快。
2)在第一对典型变量u1、v1中,u1为能源消费指标的线性组合,其中x2(油品消费量)和x3(电力消费量)比其他变量有较大的载荷,说明油品、电力是能源消费量的主要指标,它们在能源消费中占主导地位。
x4(能源进口量)比x1(原煤消费量)有较大载荷,说明随着经济的逐渐发展,本地的能源逐渐不能满足经济发展的需要,进口能源逐渐展示其重要性。
。
从图中可以看出,散点在近似的一条直线上分布,两者之间呈线性相关关系。
这说明用典型相关分析的方法能较好的说明能源消费与经济增长之间的相关关系。
散点图上几乎没有离开群体的差异值,这表明能源消费量和经济增长之间的关系很稳定,波动也非常平稳。
练习题1:
e11.6对我国工农业产业系统的典型相关分析。
练习题2:
e11.7对各类投资资金与三大产业的典型相关分析。