matlab与典型相关分析

1、第十章典型相关分析安庆师范学院 胡云峰习题10.2 下表给出著名统计学家rao在1952年对25个家庭的成年长子的头长（x1）、头宽（x2）、与次子头长（y1）、头宽（y2）进行调查所得数据如下:长子次子长子次子头长头宽头长头宽头长头宽头长头宽1911551791451901591951571951492011521881511871581811481851491631371611301831531881491951551831581761441711421861531731482081571921521811451821461891501901491751401651371971591891

2、52192154185152188152197159174143178147192150187151176139176143179158186148197167200158183147174147190163187150174150185152x=191 155 179 145;195 149 201 152;181 148 185 149;. 183 153 188 149;176 144 171 142;208 157 192 152;. 189 150 190 149;197 159 189 152;188 152 197 159;. 192 150 187 151;179 158 18

3、6 148;183 147 174 147;. 174 150 185 152;190 159 195 157;188 151 187 158;. 163 137 161 130;195 155 183 158;186 153 173 148;. 181 145 182 146;175 140 165 137;192 154 185 152;. 174 143 178 147;176 139 176 143;197 167 200 158;.190 163 187 150第一步计算相关系数矩阵程序 r=corrcoef(x)输出结果r = 1.0000 0.7346 0.7108 0.7040

4、 0.7346 1.0000 0.6932 0.7086 0.7108 0.6932 1.0000 0.83930.7040 0.7086 0.8393 1.0000计算a、 b的特征值特征向量程序r11=r(1,2,1,2);r12=r(1,2,3,4);r21=r(3,4,1,2);r22=r(3,4,3,4);a=(r11(-1)*r12*(r22(-1)*r21;b=(r22(-1)*r21*(r11(-1)*r12;x1 b1=eig(a);x2 b2=eig(b);s=cov(x);s1=s(1 2,1 2);s2=s(3,4,3,4);s1(1,2)=0;s1(2,1)=0;s2(1,2)=0;s2(2,1)=0;b1=b1(1/2)l=(s1(-1)*x1b2=b2(1/2)m=(s2(-1)*x2输出结果b1 = 0.7885 0 0 0.0537l = 0.0076 -0.0074 0.0126 0.0131b2 = 0.0537 0 0 0.7885m = -0.0070 -0.0068 0.0157 -0.0162从而得出典型相关系数和典型变量第三步 典型相关系数的显著性检验这里由于数据比较少就不用程序计算了，直接手算（1）所以第一个典型相关系数为高度显著（2）所以第二个典型相关系数不显著，对第二个典