常用生物統(tǒng)計方法SAS程序

1、常用生物統(tǒng)計方法的SAS 一、SAS系統(tǒng)簡介二、SAS系統(tǒng)的啟動與關(guān)閉 三、SAS程序結(jié)構(gòu)、程序的輸入、修改 調(diào)試和運行 四、常用生物統(tǒng)計方法的SAS程序 t檢驗 方差分析 線性回歸與相關(guān)分析 1一、SAS系統(tǒng)簡介 SAS（Statistical Analysis System）是美國使用最為廣泛的三大著名統(tǒng)計分析軟件（SAS，SPSS，SYSTAT）之一； 1966年美國North Carolina 州立大學(xué)開始開發(fā) 目前國際上最流行的一種大型統(tǒng)計分析系統(tǒng) 統(tǒng)計分析的標(biāo)準(zhǔn)軟件 其基本部分是SAS/BASE軟件，在此基礎(chǔ)上可根據(jù)需要加上SAS系統(tǒng)的其他模塊2 SAS/STAT（統(tǒng)計分析）；SA

2、S/ETS（經(jīng)濟預(yù)測或時間序列分析）；SAS/QC（質(zhì)量控制管理）；SAS/OR（運籌學(xué)）；SAS/IML（矩陣運算）；SAS/GRAPH（繪圖）等 信息存儲簡單（能讀取任何形式的數(shù)據(jù)） 語言編程能力強 對數(shù)據(jù)連續(xù)處理 統(tǒng)計分析方法豐富，使用簡單 SAS過程傻瓜系統(tǒng)3二、SAS系統(tǒng)的啟動與關(guān)閉 （一）啟動 在WINDOWS操作系統(tǒng)中，直接用鼠標(biāo)雙擊桌面上SAS系統(tǒng)的快捷鍵圖標(biāo)，即進入SAS系統(tǒng)。 SAS系統(tǒng)軟件有時間租期限制，因此只有機器時間在軟件有效期內(nèi)才能運行 OUTPUT（SAS結(jié)果輸出窗口） LOG（SAS日志窗口，隨著SAS語句的執(zhí)行，顯示出SAS系統(tǒng)的信息和已執(zhí)行的語句） PGM（

3、SAS程序編輯窗口）（二）退出 4四、SAS程序結(jié)構(gòu)（一）程序結(jié)構(gòu) 在SAS系統(tǒng)中任何一個完整的處理過程均可分為兩大步數(shù)據(jù)步和過程步數(shù)據(jù)步將不同來源的數(shù)據(jù)讀入SAS系統(tǒng)建立起SAS數(shù)據(jù)集。由DATA語句開始，以RUN語句結(jié)束。 數(shù)據(jù)可以由外部或過程讀入；或者由CARDS語句開始，直接輸入數(shù)據(jù)，數(shù)據(jù)輸入完后提行單獨加“;”結(jié)束5 過程步調(diào)用SAS系統(tǒng)中已編好的各種過程來處理和分析數(shù)據(jù)集中的數(shù)據(jù)。 每一個過程步均以PROC語句開始，RUN語句結(jié)束，并且每個語句后均以“;”結(jié)束。（二）程序的輸入、修改調(diào)試和運行 SAS程序只能在PGM窗口輸入、修改，可在各行的任何位置開始SAS語句的書寫。一個語句可

4、以連續(xù)寫在幾行中，一行中也可以同時寫上幾個語句，但每個語句后面必須用“；”隔開。6 當(dāng)一個程序輸入完后，發(fā)送程序的命令為F10功能鍵或SUBMIT。當(dāng)程序發(fā)送到SAS系統(tǒng)后，PGM的程序語句全部自動清除，LOG窗口將逐步記下程序運行的過程和出現(xiàn)的錯誤信息（用紅色提示錯誤）。 如果過程步?jīng)]有錯誤，運行完成后，通常會在OUTPUT窗口打印出結(jié)果 如果程序運行出錯，則需要在PGM窗口用RECALL（或F9）命令調(diào)回已發(fā)送的程序進行修改。7四、常用生物統(tǒng)計方法的SAS程序 （一）t 檢驗 1、樣本平均數(shù)與總體平均數(shù)的差異顯著性檢驗(教材【例5.1】) DATA A; INPUT y; y=y-4.50

5、; cards; 4.33 4.62 3.89 4.14 4.78 4.64 4.52 4.48 4.55 4.26 ;PROC MEANS N MEAN STDERR T PRT;RUN;8程序說明：樣本平均數(shù)與總體平均數(shù)的差異顯著性檢驗可調(diào)用MEANS過程。DATA語句產(chǎn)生臨時數(shù)據(jù)集A，表明數(shù)據(jù)步的開始；INPUT語句指明讀取變量y，表示讀入一條觀測值后不換行，連續(xù)讀入數(shù)據(jù)，減少數(shù)據(jù)輸入行；CARDS語句表明以下為數(shù)據(jù)行，數(shù)據(jù)行下的“；”表示數(shù)據(jù)行結(jié)束；PROC MEANS 語句指明調(diào)用MEANS過程對數(shù)據(jù)集A進行分析，輸出樣本含量N、平均數(shù)MEAN、平均數(shù)的標(biāo)準(zhǔn)誤STDERR、學(xué)生氏T值

6、和t值概率PRT；RUN語句表示過程步結(jié)束，開始運行過程步。92、配對試驗資料的t檢驗（教材【例5.5】） DATA B;INPUT ID x1 x2;d=x1-x2;CARDS; 1 37.8 37.9 2 38.2 39.0 3 38.0 38.9 4 37.6 38.4 5 37.9 37.9 6 38.1 39.0 7 38.2 39.5 8 37.5 38.6 9 38.5 38.8 10 37.9 39.0;PROC MEANS MEAN STDERR T PRT; VAR d;RUN;103、非配對試驗資料的t檢驗（教材【例5.3】） DATA C;INPUT breed y;C

7、ARDS; 1 57 2 89 1 120 2 36 1 101 2 82 1 137 2 50 1 119 2 39 1 117 2 32 1 104 2 57 1 73 2 82 1 53 2 96 1 68 2 31 1 118; PROC TTEST; CLASS breed; VAR y; RUN;11（ 二）方差分析 對于一般的方差分析（平衡資料，即各處理重復(fù)數(shù)相等）可用ANOVA過程；對于非平衡資料（各處理重復(fù)數(shù)不等）的方差分析可用GLM過程。 121. 單因素試驗（各處理重復(fù)數(shù)相等， 教材【例6.1】 ）2. 單因素試驗（各處理重復(fù)數(shù)不等， 教材【例6.3】 ）3. 兩因素交叉

8、分組單獨觀察值（教材【例6.4】 ）4. 兩因素交叉分組重復(fù)觀察值（教材【例6.5】 ）5. 兩因素系統(tǒng)分組次級樣本含量相等（教材【例6.6】 ）6. 隨機單位組試驗資料（教材【例12.3】）7. 拉丁方試驗資料（教材【例12.4】）13教材【例6-1】 DATA E;INPUT TR$ y;cards; A1 15 A1 16 A1 15 A1 17 A1 18 A2 45 A2 42 A2 50 A2 38 A2 39 A3 30 A3 35 A3 29 A3 31 A3 35 A4 31 A4 28 A4 20 A4 25 A4 30 A5 40 A5 35 A5 31 A5 32 A5

9、 30; PROC anova; CLASS TR; model y=TR; MEANS TR/T;RUN;14DATA F;INPUT BREED$ y;cards; B1 21.5 B1 19.5 B1 20.0 B1 22.0 B1 18.0 B1 20.0 B2 16.0 B2 18.5 B2 17.0 B2 15.5 B2 20.0 B2 16.0 B3 19.0 B3 17.5 B3 20.0 B3 18.0 B3 17.0 B4 21.0 B4 18.5 B4 19.0 B4 20.0 B5 15.5 B5 18.0 B5 17.0 B5 16.0; PROC GLM; CLAS

10、S BREED; model y=BREED; MEANS BREED/T;RUN;15DATA F;INPUT A$ B$ y;cards; A1 B1 3.5 A1 B2 2.3 A1 B3 2.0 A2 B1 2.6 A2 B2 2.0 A2 B3 1.9 A3 B1 2.0 A3 B2 1.5 A3 B3 1.2 A4 B1 1.4 A4 B2 0.8 A4 B3 0.3; PROC ANOVA; CLASS A B; model y=A B; MEANS A B/DUNCAN;RUN;16 DATA CaP; DO a=1 to 4; DO b=1 to 4; DO n=1 to 3

11、; input y; output; end; end; end; DROP n; CARDS; 22.0 26.5 24.4 30.0 27.5 26.0 32.4 26.5 27.0 30.5 27.0 25.1 23.5 25.8 27.0 33.2 28.5 30.1 38.0 35.5 33.9 26.5 24.0 25.0 30.5 26.8 25.5 36.5 34.0 33.5 28.0 30.5 24.6 20.5 22.5 19.5 34.5 31.4 29.3 29.0 27.5 28.0 27.5 26.3 28.5 18.5 20.0 19.0; PROC ANOVA

12、; CLASS A B; MODEL y=A B A*B; RUN;17DATA G; DO a=1 to 3; DO b=1 to 3; DO n=1 to 2; input y; output; end; end; end;DROP n;CARDS; 82.5 82.4 87.1 86.5 84.0 83.9 86.6 85.8 86.2 85.7 87.0 87.6 82.0 81.5 80.0 80.5 79.5 80.3 ; PROC ANOVA; CLASS A B; MODEL y=A B(A); RUN;18DATA R;INPUT T$ B$ X;CARDS;A1 B1 20

13、5 A1 B2 168 A1 B3 222 A1 B4 230A2 B1 230 A2 B2 198 A2 B3 242 A2 B4 255A3 B1 252 A3 B2 248 A3 B3 305 A3 B4 260A4 B1 200 A4 B2 158 A4 B3 183 A4 B4 196A5 B1 265 A5 B2 275 A5 B3 315 A5 B4 282;PROC GLM;CLASS T B;MODEL X=T B;MEANS T/DUNCAN;RUN;19DATA L;INPUT T$ ROW COL X;CARDS;D 1 1 23 E 1 2 21 A 1 3 24 B

14、 1 4 21 C 1 5 19A 2 1 22 C 2 2 20 E 2 3 20 D 2 4 21 B 2 5 22E 3 1 20 A 3 2 25 B 3 3 26 C 3 4 22 D 3 5 23B 4 1 25 D 4 2 22 C 4 3 25 E 4 4 21 A 4 5 23C 5 1 19 B 5 2 20 D 5 3 24 A 5 4 22 E 5 5 19;PROC GLM;CLASS T ROW COL;MODEL X=T ROW COLl;MEANS T/DUNCAN;RUN;20（三）一元線性回歸與相關(guān)分析（教材【例8.1】） DATA G；INPUT x y；

15、CARDS； 80 2350 86 2400 98 2720 90 2500 120 3150102 2680 95 2630 83 2400 113 3080 105 2920110 2960 100 2860；PROC REG CORR；MODEL y=x ；RUN；21（四）協(xié)方差分析（教材【例10.1】）DATA K;INPUT T$ X Y;CARDS;CK 1.50 12.40 CK 1.85 12.00 CK 1.35 10.80CK 1.45 10.00 CK 1.40 11.00 CK 1.45 11.80CK 1.50 12.50 CK 1.55 13.40 CK 1.40

16、 11.20CK 1.50 11.60 CK 1.60 12.60 CK 1.70 12.50 1 1.35 10.20 1 1.20 9.40 1 1.45 12.20 1 1.20 10.30 1 1.40 11.30 1 1.30 11.40 1 1.15 12.80 1 1.30 10.90 1 1.35 11.60 1 1.15 8.50 1 1.35 12.20 1 1.20 9.30 2 1.15 10.00 2 1.10 10.60 2 1.10 10.40 2 1.05 9.20 2 1.40 13.00 2 1.45 13.50 2 1.30 13.00 2 1.70 14.80 2 1.40 12.30 2 1.45 13.20 2 1.25 12.00