SAS例題及程序輸出6[學(xué)習(xí)相關(guān)]

1、地質(zhì)勘探中，在A,B,C三個地區(qū)采集了一些巖石，測量其部分化學(xué)成分，其數(shù)據(jù)見表3.5。假定這三個地區(qū)掩飾的成分遵從。（1）檢驗不全（2）檢驗；（3）檢驗。表3.5 巖石部分化學(xué)成分數(shù)據(jù)SiO2FeOK2OA地區(qū)47.225.060.1047.454.350.1547.526.850.1247.864.190.1747.317.570.18B地區(qū)54.336.220.1256.173.310.1554.402.430.2252.625.920.12C地區(qū)43.1210.330.0542.059.670.0842.509.620.0240.779.680.04解：(1)檢驗假設(shè)，在H0成立時，取近

2、似檢驗統(tǒng)計量為 統(tǒng)計量：。由樣本值計算三個總體的樣本協(xié)方差陣：進一步計算可得對給定顯著性水平，利用軟件SAS9.3進行檢驗時，首先計算p值：p=P13.896916=0.3073394。因為p值=0.30733940.05，故接收，即認為方差陣之間無顯著性差異。proc iml;n1=5;n2=4;n3=4;n=n1+n2+n3;k=3;p=3;x1=47.22 5.06 0.1,47.45 4.35 0.15,47.52 6.85 0.12,47.86 4.190.17,47.31 7.570.18;x2=54.33 6.22 0.12,56.17 3.310.15,54.4 2.430.2

3、2,52.62 5.920.12;x3=43.12 10.33 0.05,42.05 9.670.08,42.5 9.62 0.02,40.77 9.680.04;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;tt=xx*mm*xx;print tt;/*總離差陣*/a=a1+a2+a3;print a;/*組內(nèi)離差陣*/da=det(a/

4、(n-k);/*合并樣本協(xié)差陣*/da1=det(a1/(n1-1);/*每個總體的樣本協(xié)差陣陣*/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; /*統(tǒng)計量值*/print da da1 da2 da3 m dd df;p0=1-probchi(kc,df); /*顯著性概率*/print kc p0;qui

5、t;(2) 提出假設(shè)。取檢驗統(tǒng)計量為，由樣本值計算得：進一步計算得：對給定顯著性水平，利用軟件SAS9.3進行檢驗時，首先計算p值：p=PF32.098939=0.0010831。因為p值=0.00108310.05，故否定，即認為A，B兩地巖石化學(xué)成分數(shù)據(jù)存在顯著性差異。在這種情況下，可能犯第一類錯誤，且犯第一類錯誤的概率為0.05。SAS程序及結(jié)果如下：proc iml;n=5;m=4; p=3;x= 47.22 5.06 0.1,47.45 4.35 0.15,47.52 6.85 0.12,47.86 4.190.17,47.31 7.570.18 ;ln=5 1 ;x0=(ln*x)

6、/n; print x0;mx=i(n)-j(n,n,1)/n;a1=x*mx*x; print a1;y= 54.33 6.22 0.12,56.17 3.310.15,54.4 2.430.22,52.62 5.920.12 ;lm=4 1 ;y0=(lm*y)/m; print y0;my=i(m)-j(m,m,1)/m;a2=y*my*y; print a2;a=a1+a2; xy=x0-y0;ai=inv(a); print a ai;dd=xy*ai*xy; d2=(m+n-2)*dd;t2=n*m*d2/(n+m) ;f=(n+m-1-p)*t2/(n+m-2)*p);fa=fi

7、nv(0.95,p,m+n-p-1);beta=probf(f,p,m+n-p-1,t2);print d2 t2 f beta;pp=1-probf(f,p,m+n-p-1);print pp; quit;(3) 檢驗假設(shè)；因似然比統(tǒng)計量 ，本題中k-1=2，可以利用統(tǒng)計量與F統(tǒng)計量的關(guān)系，去檢驗統(tǒng)計量為F統(tǒng)計量：由樣本值計算得：及 ，進一步計算得：對給定顯著性水平，利用軟件SAS9.3進行檢驗時，首先計算p值：p=PF18.390234=2.345110-6。因為p值=2.345110-60.05，故否定，即認為A，B，C三地巖石化學(xué)成分數(shù)據(jù)存在顯著性差異。在這種情況下，可能犯第一類錯誤，

8、且犯第一類錯誤的概率為0.05。proc iml;n1=5;n2=4;n3=4;n=n1+n2+n3;k=3;p=3;x1=47.22 5.06 0.1,47.45 4.35 0.15,47.52 6.85 0.12,47.86 4.190.17,47.31 7.570.18;x2=54.33 6.22 0.12,56.17 3.310.15,54.4 2.430.22,52.62 5.920.12;x3=43.12 10.33 0.05,42.05 9.670.08,42.5 9.62 0.02,40.77 9.680.04;xx=x1/x2/x3; /*三組樣本縱向拼接*/ln=51;ln

9、n41;lnnn=131;x10=(ln*x1)/n1;x20=(lnn*x2)/n2;x30=(lnn*x3)/n3;xx0=(lnnn*x1)/n1;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;a2=x2*mm2*x2;a3=x3*mm2*x3;tt=xx*mm*xx;print tt;/*總離差陣*/a=a1+a2+a3; print a;/*組內(nèi)離差陣*/da=det(a);/*合并樣本協(xié)差陣*/dt=det(tt);a0=da/dt;print da dt a0;b=sqrt(a0); print b;f=(n-k-p+1)*(1-b)/(b*p);df1=2*p;df2=2*(n-k-p+1);p0=