sas例题及程序输出6.doc

SAS例题及程序输出6
地质勘探中,在A,B,C三个地区采集了一些岩石,测量其部分化学成分,其
(i),?i?(i?1,2,3)(??)。。假定这三个地区掩饰的成分遵从N3??
(1)检验不全H0:?1??2=?3;H1:?1,?2,?3不全等;
(2)检验H0:?(1)??(2);H1:?(1)??(2);
(3)检验H0:?(1)??(2)??(3);H1:存在i?j,使?(i)??(j)。
(1)检验假设
H0:?1??2=?3;H1:?1,?2,?3不全等,
在H0成立时,取近似检验统计量为?2(f) 统计量:
?=?1?d?M??2?1?d?ln?4*。
由样本值计算三个总体的样本协方差阵:
11n1(1)(1)(1)S1?A1?(X((1)?)(X?)???)(?)n1?1n1?1??1
??1??=???,4????
11n1(2)(2)S2?A3?(X((2))(X((2))???)??)?n2?1n2?1??1
??1??=??,?3?????
11n1(3)(3)S3?A3?(X((3))(X((3))???)??)?n3?1n3?1??1
??1??=?。?4?????
进一步计算可得
S? 1A?,S1?,S2?,S3?, 10
M?,d?,f?12,
??(1?d)M=。
对给定显著性水平?=,,首先计算p值:
p=P{ξ≥}=。
因为p值=>,故接收H0,即认为方差阵之间无显著性差异。 proc iml;
n1=5;n2=4;n3=4;
n=n1+n2+n3;k=3;p=3;
x1={ ,
,
,
,
};
x2={ ,
,
,
};
x3={ ,
,
,
};
xx=x1//x2//x3; /*三组样本纵向拼接*/
mm1=i(5)-j(5,5,1)/n1;
mm2=i(4)-j(4,4,1)/n2;
mm=i(n)-j(n,n,1)/n;
a1=x1`*mm1*x1;print a1;
a2=x2`*mm2*x2;print a2;
a3=x3`*mm2*x3;print a3; t