SAS例题及程序输出6.docx

地质勘探中,在A,B,C三个地区采集了一些岩石,测量其部分化学成分,。假定这三个地区掩饰的成分遵从。(1)检验不全(2)检验;(3)检验。:(1)检验假设,在H0成立时,取近似检验统计量为统计量:。由样本值计算三个总体的样本协方差阵:进一步计算可得对给定显著性水平,,首先计算p值:p=P{ξ≥}=。因为p值=>,故接收,即认为方差阵之间无显著性差异。prociml;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;printa1;a2=x2`*mm2*x2;printa2;a3=x3`*mm2*x3;printa3;tt=xx`*mm*xx;printtt;/*总离差阵*/a=a1+a2+a3;printa;/*组内离差阵*/da=det(a/(n-k));/*合并样本协差阵*/da1=det(a1/(n1-1));/*每个总体的样本协差阵阵*/da2=det(a2/(n2-1));da3=det(a3/(n3-1));m=(n-k)*log(da)-(4*log(da1)+3*log(da2)+3*log(da3));dd=(2*p*p+3*p-1)*(k+1)/(6*(p+1)*(n-k));df=p*(p+1)*(k-1)/2;/*卡方分布自由度*/kc=(1-dd)*m;/*统计量值*/printdada1da2da3mdddf;p0=1-probchi(kc,df);/*显著性概率*/printkcp0;quit;(2)提出假设。取检验统计量为,由样本值计算得:进一步计算得:对给定显著性水平,,首先计算p值:p=P{F≥}=。因为p值=<,故否定,即认为A,B两地岩石化学成分数据存在显著性差异。在这种情况下,可能犯第一类错误,。SAS程序及结果如下:prociml;n=5;m=4;p=3;x={