典型相关分析.pdf

《典型相关分析.pdf》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

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

2、组变量的第一对线性组合为：uXvY1111(a,a,,a)11121p1(,,,)11121q1Var(u)Var(X)111111Var(v)Var(Y)11111221典型相关分析就是求和，使二者的相关系数达到最大。11Cov(u,v)Cov(X,Y)u,v11111121111欢迎下载。精品文档典型相关分析希望寻求a和b使得ρ达到最大，但是由于随机变量乘以常数时不改变它们的相关系数，为了防止不必要的结果重复出现，最好的限制是令Var（U）

3、=1和Var（V）=1。1.实测变量标准化；2.求实测变量的相关阵R；XX,...,XYY,...,Y1p1qrrrr111p111qRRXXXYrrrrRxxxyp1ppp1pqRRYXYYrrrryxyy111p111qrrrrq1qpq1qq(pq)(pq)3.求A和B；AXX1XYYY1YXBYY1YXXX1XY4、求A和B的特征根及特征向量；...12pA关于的特征向量(a,a

4、,…,a)，求B关于iii1i2ip的特征向量(bi,b,…,bi)1i2p5、计算Vi和Wi；VbXbX...bXii11i22ippWaYaY...aYii11i22iqq2欢迎下载。精品文档6、Vi和Wi的第i对典型相关系数rii应用典型相关分析的场合是：可以使用回归方法，但有两个或两个以上的因变量；特别是因变量或准则变量相互间有一定的相关性，无视它们之间相互依赖的关系而分开处理，研究就毫无意义。另一种有效用法是检验X变量集合和Y变量集合间的独立性。四、典型相关系数的检验典型相关分析是否恰当，应该取决于两组原变量之间是

5、否相关，如果两组变量之间毫无相关性而言，则不应该作典型相关分析。用样本来估计总体的典型相关系数是否有误，需要进行检验。在原假设为真的情况下，检验的统计量为：1Qn-1(pq1)ln020近似服从自由度为pq的2分布。在给定的显著性水平下，如果22(pq)，则拒绝原假设，认为至少第一对典型变量之间的相关性显著。相应的R编程如下：setwd("D:/data")ex1=read.table("9-1.txt",head=T)3欢迎下载。精品文档ex1x=ex1[,1:3];xy=ex1[,4:6];yx=as.

6、matrix(x)y=as.matrix(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[,1]a=sqrt(eigen(A)$values)#求典型相关系数=sqrt(特征值)axt(a)4欢迎下载。精品文档t(t(a)

7、)%*%xB=solve(s22)%*%s21%*%solve(s11)%*%s12Beigen(B)sqrt(eigen(B)$values)A0=prod(1-eigen(A)$values)A0Q0=-15.5*log(A0);Q0#求检验统计量pr=1-pchisq(Q0,9)#求P值prm1=cancor(x,y)#典型相关分析m1#相关系数的假设检验corcoef.test<-function(r,n,p,q,alpha=0.1){#r为相关系数n为样本个数且n>p+qm<-length(r);Q<-rep(0,m);lambda<

8、-1for(kinm:1){lambda<-lambda*(1-r[k]^2);#检验统计量5欢迎下载。精品文档Q[k]<--log(lambda)#