典型相关分析

1、典型相关分析一、典型相关分析的概念典型相关分析（canonical correlation analysis）就是利用综合变量对之间的相关关系来反映两组指标之间的整体相关性的多元统计分析方法。它的基本原理是：为了从总体上把握两组指标之间的相关关系，分别在两组变量中提取有代表性的两个综合变量u1和v1（分别为两个变量组中各变量的线性组合），利用这两个综合变量之间的相关关系来反映两组指标之间的整体相关性。二、条件：典型相关分析有助于综合地描述两组变量之间的典型的相关关系。其条件是,两组变量都是连续变量,其资料都必须服从多元正态分布。三、相关计算如果我们记两组变量的第一对线性组合为： 典型相关分析就

2、是求a1和b1，使二者的相关系数r达到最大。1. 实测变量标准化；2. 求实测变量的相关阵r； 3. 求a和b；4、求a和b的特征根及特征向量； 典型相关分析希望寻求 a 和 b 使得 达到最大，但是由于随机变量乘以常数时不改变它们的相关系数，为了防止不必要的结果重复出现，最好的限制是令 var （u） =1 和 var（v）= 1。a关于 的特征向量(ai1,ai2,aip)，求b关于 的特征向量(bi1,bi2,bip) 5、计算vi和wi ；6、vi和wi 的第i对典型相关系数 应用典型相关分析的场合是：可以使用回归方法，但有两个或两个以上的因变量；特别是因变量或准则变量相互间有一定的相

3、关性，无视它们之间相互依赖的关系而分开处理，研究就毫无意义。另一种有效用法是检验x变量集合和y变量集合间的独立性。四、典型相关系数的检验典型相关分析是否恰当，应该取决于两组原变量之间是否相关，如果两组变量之间毫无相关性而言，则不应该作典型相关分析。用样本来估计总体的典型相关系数是否有误，需要进行检验。在原假设为真的情况下，检验的统计量为：近似服从自由度为pq的c2分布。在给定的显著性水平a下，如果c2³c2 (pq)，则拒绝原假设，认为至少第一对典型变量之间的相关性显著。相应的r编程如下：setwd("d:/data")ex1=("",head

4、=t)ex1x=ex1,1:3;xy=ex1,4:6;yx=(x)y=(y)x;ys11=cov(x);s11s22=cov(y);s22s12=cov(ex1)1:3,4:6;s12s21=cov(ex1)4:6,1:3;s21#求协方差矩阵a=solve(s11)%*%s12%*%solve(s22)%*%s21#矩阵相乘用%*%，solve:求逆矩阵aeigen(a)#求特征值及其对应的特征向量，eigen(a)$vectors,1a=sqrt(eigen(a)$values)#求典型相关系数=sqrt(特征值)axt(a)t(t(a)%*%xb=solve(s22)%*%s21%*%s

5、olve(s11)%*%s12beigen(b)sqrt(eigen(b)$values)a0=prod(1-eigen(a)$values)a0q0=*log(a0);q0#求检验统计量pr=1-pchisq(q0,9)#求p值prm1=cancor(x,y)#典型相关分析m1#相关系数的假设检验<-function(r, n, p, q, alpha= #r为相关系数 n为样本个数 且n>p+q m<-length(r); q<-rep(0, m); lambda <- 1 for (k in m:1) lambda<-lambda*(1-rk2); #检验统计量 qk<- -log(lambda) #检验统计量取对数 s<-0; i<-m for (k in 1:m) qk<- (n-k+1-1/2*(p+q+3)+s)*qk #统计量 chi<-1-pchisq(qk, (p-k+1)*(q-k+1) if (chi>alpha)