应用多元分析第四章 多元正态总体1.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 或等价于 其中 例1 proc iml; n1=20;n2=20;n3=20; n=n1+n2+n3;k=3; p=4; x1={ 260 75 40 18 , 200 72 34 17 , 240 87 45 18 , 170 65 39 17 , 270 110 39 24 , 205 130 34 23 , 190 69 27 15 , 200 46 45 15 , 250 117 21 20 , 200 107 28 20 , 225 130 36 11 , 210 125 26 17 , 170 64 31 14 , 270 76 33 13 , 190 60 34 16 , 280 81 20 18 , 310 119 25 15 , 270 57 31 8 , 250 67 31 14 , 260 135 39 29 } ; x2={ 310 122 30 21 ,310 60 35 18 ,190 40 27 15 ,225 65 34 16 , 170 65 37 16 ,210 82 31 17 ,280 67 37 18 ,210 38 36 17 , 280 65 30 23 ,200 76 40 17 ,200 76 39 20 ,280 94 26 11 , 190 60 33 17 ,295 55 30 16 ,270 125 24 21 ,280 120 32 18 , 240 62 32 20 ,280 69 29 20 ,370 70 30 20 ,280 40 37 17 }; x3={ 320 64 39 17 ,260 59 37 11 ,360 88 28 26 ,295 100 36 12 , 270 65 32 21 ,380 114 36 21 ,240 55 42 10 ,260 55 34 20 , 260 110 29 20 ,295 73 33 21 ,240 114 38 18 ,310 103 32 18 , 330 112 21 11 ,345 127 24 20 ,250 62 22 16 ,260 59 21 19 , 225 100 34 30 ,345 120 36 18 ,360 107 25 23 ,250 117 36 16 }; xx=x1//x2//x3; ln={[20] 1} ; x10=(ln*x1)/n1; print x10; mm1=i(n1)-j(n1,n1,1)/n1; mm=i(n)-j(n,n,1)/n; s1=x1`*mm1*x1/(n1-1); s2=x2`*mm1*x2/(n2-1); s3=x3`*mm1*x3/(n3-1); tt=xx`*mm*xx/(n-1); s=(s1+s2+s3)*(n1-1)/(n-3); ds1=det(s1); ds2=det(s2); ds3=det(s3); ds=det(s); d=(2*p**2+3*p-1)*(k+1)/((6*p+1)*(n-k)); f=p*(p+1)*(k-1)/2; m=(n-k)*log(det(s))-(n1- 1)*(log(det(s1))+ log(det(s2))+log(det(s3))); kc=(1-d)*m; p0=1-probchi(kc,df) ; print s1 s2 s3 s f d m kc; 记 其中 对例1 中的数据，判断三个组(即三个总体)的均值向量和协方差矩阵是否全相等(α=0.05) proc iml; n1=20;n2=20;n3=20