SAS例題及程序輸出

1、頁(yè)眉1/ 8地質(zhì)勘探中，在 A,B,C 三個(gè)地區(qū)采集了一些巖石，測(cè)量其部分化學(xué)成分，其數(shù)據(jù)見(jiàn)表 3.5。假定這三個(gè)地區(qū)掩飾的成分遵從N3()i(i 123)(0.05)(1)檢驗(yàn)不全 H:!2=3；出：!,2,3不全等;H0:12=3；H1:1，2，3不全等，在 Ho成立時(shí)，取近似檢驗(yàn)統(tǒng)計(jì)量為2(f)統(tǒng)計(jì)量：=1 d M 2 1 d l n4*。由樣本值計(jì)算三個(gè)總體的樣本協(xié)方差陣：S 士 A n1(X()X)(X()X)n 1n1110.243081=0.64264 9.28552，40.014060.02052 0.00452SiO2FeOK2O47.225.060.1047.454.350

2、.15A 地區(qū)47.526.850.1247.864.190.1747.317.570.1854.336.220.12B 地區(qū)56.173.310.1554.402.430.2252.625.920.1243.1210.330.05C 地區(qū)42.059.670.0842.509.620.0240.779.680.04(3)檢驗(yàn)Ho:(i)；比：存在 i(i)j,使(j)。表 3.5 巖石部分化學(xué)成分?jǐn)?shù)據(jù)解:(1)檢驗(yàn)假設(shè)(2)檢驗(yàn)Ho:(i)(2) (1) (2)；H1:頁(yè)眉2/ 8S2亠 A3ni(x(2)X(2)(X(2)X(2)n21n21i6.30461二-4.7567 10.6722

3、,30.05570.2388 0.0066751 1n1S3亠人亠(x()X)(x()X(3)n31n3112.97141二-0.63370.3421。40.0001 0.00295 0.001875進(jìn)一步計(jì)算可得0.0000942, S20.0011851, S30.0000417,M 24.52397, d 0.433333, f 12,(1 d)M =13.896916。對(duì)給定顯著性水平 =0.05，利用軟件 SAS9.3 進(jìn)行檢驗(yàn)時(shí)，首先計(jì)算 p 值:p=P詡 3.896916=0.3073394。因?yàn)?p 值=0.30733940.05,故接收 H。，即認(rèn)為方差陣之間無(wú)顯著性差異。p

4、roc iml ;n1= 5 ;n2= 4;n3= 4;n=n1+n 2+n3;k:=3;p= 3;x1= 47.225.060.1 ,47.454.350.15J47.526.850.12J47.864.190.17J47.31;7.570.18x2=54.336.220.12 ,56.173.310.15J54.42.430.22J52.62;5.920.12x3= 43.1210.330.05 ,42.059.670.08J42.59.620.02 ,40.779.680.041S A 0.0018318, S1頁(yè)眉3/ 8;xx=x1/x2/x3;/* 三組樣本縱向拼接 */頁(yè)眉4/

5、8mm1=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;pri nta1;a2=x2*mm2*x2;pri nta2;a3=x3*mm2*x3;pri nta3;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; /* 卡方分

6、布自由度 */kc=( 1 -dd)*m;/* 統(tǒng)計(jì)量值 */pr int da da1 da2 da3 m dd df;p0= 1 -probchi(kc,df);/* 顯著性概率 */ pri nt kc p0;quitSAS系統(tǒng)al0 243C6 0 042840 0KB-0鈿2別0.0102& JCM8-4. /Eft/CLOfl57-4 7567-6対88| 06 勺-0 2388oaean2 B7HQ S337-a CDor0 60D 34210 0C曲 u-0 0001 00 001汎tt312 46343 -132 BQS 0 W17O77-132. 50tt2與 98

7、23K!-1.別矽82 5157077 -1 =5 強(qiáng) B46Q Ml OTS9頁(yè)眉5/ 8-4 765610. O6&66-4 765S420 2曲號(hào)-0. 215330. 06966 -0. 215330. 0130?dodaldc2dal1dd df0 0018318Q 0CTQG吃 00118510 00001724 5E3870 H33333312kcPO13 BS691 6Ll. 3O73JH42nm2TD 134.81554,n mn m p 1FT232.09893&(n m 2) p對(duì)給定顯著性水平 =0.05，利用軟件 SAS9.3 進(jìn)行檢驗(yàn)時(shí)，首先計(jì)算 p

8、 值:p=P F 為 2.098939=0.0010831。因?yàn)?p 值=0.00108310.05,故否定 H0，即認(rèn)為 A，B 兩地巖石化學(xué)成分?jǐn)?shù) 據(jù)存在顯著性差異。在這種情況下，可能犯第一類(lèi)錯(cuò)誤，且犯第一類(lèi)錯(cuò)誤的概率為0.05OSAS 程序及結(jié)果如下：提出假設(shè)(1) H0:，H1:取檢驗(yàn)統(tǒng)計(jì)量為n+m p 1 丁2m 2)(n(p 3,n 6,m 9),由樣本值計(jì)算得:X=(47.472.5.604,0.144),X(2)=(54.38,4.47,0.1525)，0.24308Ai=0.642649.285520.014060.02052 0.004526.3046進(jìn)一步計(jì)算得:D2(n

9、A2=4.75670.05572)(X10.67220.2388 0.006675X)(A A2)1(X( (1)X( (2)=60.666995，頁(yè)眉6/ 8prociml ;n= 5;m= 4; p=3;x=47.225.060.1 ,47.454.350.15 ,47.526.850.12 ,47.864.19 0.17 ,47.317.570.18;In=5 1;x0=(I n*x)/n;pri nt x0;mx=i( n)-j( n,n.1)/n;a1=x *mx*x;print a1;y=54.336.22 0.12 ,56.173.31 0.15 ,54.42.430.22 ,5

10、2.625.92 0.12;lm=41;yO=(lm*y)/m;pri nty0;my=i(m)-j(m,m,1 )/m;a2=y*my*y; print a2; a=a1+a2; xy=xO-yO;ai=in v(a);pr int 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=finv( 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-

11、1); |pr intpp; quit ;頁(yè)眉7/ 8SAS系統(tǒng)xO4/ 4了25.伽a26.3CJ|6-4頁(yè)的仃0657-4. /56710 S722J皿皿;QC95?0 OOS375al0 2D8-0.目42S40 UH06-0 Q42G4fl 285520 OHW0. 020520 00452a wVO5d 384 470 I 525(3)檢驗(yàn)假設(shè)Ho：(i)因似然比統(tǒng)計(jì)量amia 54738 -5 39934D QSB70 0 1870B560 QITFiAR Yi迢閃辭-5 39934 19 D5772-D 21fi280 OJOB990. 0767294II? 059 Ti -0

12、21523J ?1H95-0陽(yáng)菲B21160145 11. 82S97d212fbo L10-60.05，故否定 H。，即認(rèn)為 A, B, C 三地巖石化學(xué)成 分?jǐn)?shù)據(jù)存在顯著性差異。在這種情況下，可能犯第一類(lèi)錯(cuò)誤，且犯第一類(lèi)錯(cuò)誤的 概率為 0.05。proc iml ;n1= 5 ;n2= 4;n3= 4;n=n1+n2+n 3;k=3 ;p= 3;x1= 47.225.060.1 ,47.454.350.15J47.526.850.12J47.864.190.17J47.31；7.570.18x2= 54.336.220.1256.173.310.15J54.42.430.22J52.62；

13、5.920.12x3= 43.1210.330.0542.059.670.08J42.59.620.02 ,40.779.680.04;xx=x1/x2/x3;/*三組樣本縱向拼接 */In=5 1 ;lnn4 1 ;lnnn=13 1;x10=(l n*x1)/n1;|x20=(l nn *x2)/n2;x30=(l nn *x3)/n3;xx0=(l nnn *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;頁(yè)眉10/ 8a2=x2*mm2*x2;a3=x

14、3*mm2*x3;tt=xx*mm*xx;pri nttt; /*總離差陣*/a=a1+a2+a3;pri nti; /*組內(nèi)離差陣*/da=det(a); /*合并樣本協(xié)差陣*/dt=det(tt);a0=da/dt;pri nt da dt a0;b=sqrt(aO);print b; |f=( n-k-p+1 )*( 1 -b)/(b*p);|df1= 2*p;df2=2*(n-k-p+1); |pO= 1-probf(f,df1,df2);/* 顯著性概率 */pri nt f p0;f1=(tt1, 1 -a1, 1 )*( n-k)/(k-1)*a1,1);p1= 1 -probf(f1,k-1,n-k);fa=finv( 0.95 ,k- 1 ,n-k);|pri nt fa f1 p1;quit ;SAS 系統(tǒng)tt312. 46343-132. 5062 2 54