常用生物統(tǒng)計(jì)方法的程

常用生物統(tǒng)計(jì)方法的程第一頁(yè)，共二十三頁(yè)，2022年，8月28日一、SAS系統(tǒng)簡(jiǎn)介

—SAS（StatisticalAnalysisSystem）是美國(guó)使用最為廣泛的三大著名統(tǒng)計(jì)分析軟件（SAS，SPSS，SYSTAT）之一；

—1966年美國(guó)NorthCarolina州立大學(xué)開(kāi)始開(kāi)發(fā)

—目前國(guó)際上最流行的一種大型統(tǒng)計(jì)分析系統(tǒng)

—統(tǒng)計(jì)分析的標(biāo)準(zhǔn)軟件

—其基本部分是SAS/BASE軟件，在此基礎(chǔ)上可根據(jù)需要加上SAS系統(tǒng)的其他模塊第二頁(yè)，共二十三頁(yè)，2022年，8月28日

SAS/STAT（統(tǒng)計(jì)分析）；SAS/ETS（經(jīng)濟(jì)預(yù)測(cè)或時(shí)間序列分析）；SAS/QC（質(zhì)量控制管理）；SAS/OR（運(yùn)籌學(xué)）；SAS/IML（矩陣運(yùn)算）；SAS/GRAPH（繪圖）等

—信息存儲(chǔ)簡(jiǎn)單（能讀取任何形式的數(shù)據(jù)）

—語(yǔ)言編程能力強(qiáng)

—對(duì)數(shù)據(jù)連續(xù)處理

—統(tǒng)計(jì)分析方法豐富，使用簡(jiǎn)單

—SAS過(guò)程傻瓜系統(tǒng)……第三頁(yè)，共二十三頁(yè)，2022年，8月28日二、SAS系統(tǒng)的啟動(dòng)與關(guān)閉

（一）啟動(dòng)

在WINDOWS操作系統(tǒng)中，直接用鼠標(biāo)雙擊桌面上SAS系統(tǒng)的快捷鍵圖標(biāo)，即進(jìn)入SAS系統(tǒng)。?SAS系統(tǒng)軟件有時(shí)間租期限制，因此只有機(jī)器時(shí)間在軟件有效期內(nèi)才能運(yùn)行

OUTPUT（SAS結(jié)果輸出窗口）

LOG（SAS日志窗口，隨著SAS語(yǔ)句的執(zhí)行，顯示出SAS系統(tǒng)的信息和已執(zhí)行的語(yǔ)句）

PGM（SAS程序編輯窗口）（二）退出

第四頁(yè)，共二十三頁(yè)，2022年，8月28日四、SAS程序結(jié)構(gòu)（一）程序結(jié)構(gòu)在SAS系統(tǒng)中任何一個(gè)完整的處理過(guò)程均可分為兩大步——數(shù)據(jù)步和過(guò)程步數(shù)據(jù)步——將不同來(lái)源的數(shù)據(jù)讀入SAS系統(tǒng)建立起SAS數(shù)據(jù)集。由DATA語(yǔ)句開(kāi)始，以RUN語(yǔ)句結(jié)束。

—數(shù)據(jù)可以由外部或過(guò)程讀入；或者由CARDS語(yǔ)句開(kāi)始，直接輸入數(shù)據(jù)，數(shù)據(jù)輸入完后提行單獨(dú)加“;”結(jié)束第五頁(yè)，共二十三頁(yè)，2022年，8月28日

過(guò)程步——調(diào)用SAS系統(tǒng)中已編好的各種過(guò)程來(lái)處理和分析數(shù)據(jù)集中的數(shù)據(jù)。

—每一個(gè)過(guò)程步均以PROC語(yǔ)句開(kāi)始，RUN語(yǔ)句結(jié)束，并且每個(gè)語(yǔ)句后均以“;”結(jié)束。（二）程序的輸入、修改調(diào)試和運(yùn)行

SAS程序只能在PGM窗口輸入、修改，可在各行的任何位置開(kāi)始SAS語(yǔ)句的書(shū)寫(xiě)。一個(gè)語(yǔ)句可以連續(xù)寫(xiě)在幾行中，一行中也可以同時(shí)寫(xiě)上幾個(gè)語(yǔ)句，但每個(gè)語(yǔ)句后面必須用“；”隔開(kāi)。第六頁(yè)，共二十三頁(yè)，2022年，8月28日

當(dāng)一個(gè)程序輸入完后，發(fā)送程序的命令為F10功能鍵或SUBMIT。當(dāng)程序發(fā)送到SAS系統(tǒng)后，PGM的程序語(yǔ)句全部自動(dòng)清除，LOG窗口將逐步記下程序運(yùn)行的過(guò)程和出現(xiàn)的錯(cuò)誤信息（用紅色提示錯(cuò)誤）。

—如果過(guò)程步?jīng)]有錯(cuò)誤，運(yùn)行完成后，通常會(huì)在OUTPUT窗口打印出結(jié)果

—如果程序運(yùn)行出錯(cuò)，則需要在PGM窗口用RECALL（或F9）命令調(diào)回已發(fā)送的程序進(jìn)行修改。第七頁(yè)，共二十三頁(yè)，2022年，8月28日四、常用生物統(tǒng)計(jì)方法的SAS程序

（一）t檢驗(yàn)

1、樣本平均數(shù)與總體平均數(shù)的差異顯著性檢驗(yàn)(教材【例5.1】)

DATAA;INPUTy@@;y=y-4.50;

cards;4.334.623.894.144.784.644.524.484.554.26;PROCMEANSNMEANSTDERRTPRT;RUN;第八頁(yè)，共二十三頁(yè)，2022年，8月28日程序說(shuō)明：樣本平均數(shù)與總體平均數(shù)的差異顯著性檢驗(yàn)可調(diào)用MEANS過(guò)程。DATA語(yǔ)句產(chǎn)生臨時(shí)數(shù)據(jù)集A，表明數(shù)據(jù)步的開(kāi)始；INPUT語(yǔ)句指明讀取變量y，@@表示讀入一條觀測(cè)值后不換行，連續(xù)讀入數(shù)據(jù)，減少數(shù)據(jù)輸入行；CARDS語(yǔ)句表明以下為數(shù)據(jù)行，數(shù)據(jù)行下的“；”表示數(shù)據(jù)行結(jié)束；PROCMEANS

語(yǔ)句指明調(diào)用MEANS過(guò)程對(duì)數(shù)據(jù)集A進(jìn)行分析，輸出樣本含量N、平均數(shù)MEAN、平均數(shù)的標(biāo)準(zhǔn)誤STDERR、學(xué)生氏T值和t值概率PRT；RUN語(yǔ)句表示過(guò)程步結(jié)束，開(kāi)始運(yùn)行過(guò)程步。第九頁(yè)，共二十三頁(yè)，2022年，8月28日2、配對(duì)試驗(yàn)資料的t檢驗(yàn)（教材【例5.5】）

DATAB;INPUTIDx1x2;d=x1-x2;CARDS;137.837.9238.239.0338.038.9437.638.4537.937.9638.139.0738.239.5837.538.6938.538.81037.939.0;PROCMEANSMEANSTDERRTPRT;VARd;RUN;第十頁(yè)，共二十三頁(yè)，2022年，8月28日3、非配對(duì)試驗(yàn)資料的t檢驗(yàn)（教材【例5.3】）

DATAC;INPUTbreedy@@;CARDS;1572891120236110128211372501119239111723211042571732821532961682311118;

PROCTTEST;

CLASSbreed;VARy;RUN;第十一頁(yè)，共二十三頁(yè)，2022年，8月28日（二）方差分析

對(duì)于一般的方差分析（平衡資料，即各處理重復(fù)數(shù)相等）可用ANOVA過(guò)程；對(duì)于非平衡資料（各處理重復(fù)數(shù)不等）的方差分析可用GLM過(guò)程。

第十二頁(yè)，共二十三頁(yè)，2022年，8月28日1.單因素試驗(yàn)（各處理重復(fù)數(shù)相等，教材【例6.1】

）2.單因素試驗(yàn)（各處理重復(fù)數(shù)不等，教材【例6.3】

）3.兩因素交叉分組單獨(dú)觀察值（教材【例6.4】

）4.兩因素交叉分組重復(fù)觀察值（教材【例6.5】

）5.兩因素系統(tǒng)分組次級(jí)樣本含量相等（教材【例6.6】

）6.隨機(jī)單位組試驗(yàn)資料（教材【例12.3】）7.拉丁方試驗(yàn)資料（教材【例12.4】）第十三頁(yè)，共二十三頁(yè)，2022年，8月28日教材【例6-1】

DATAE;INPUTTR$y@@;cards;A115A116A115A117A118A245A242A250A238A239A330A335A329A331A335A431A428A420A425A430A540A535A531A532A530;

PROCanova;

CLASSTR;

modely=TR;

MEANSTR/T;RUN;第十四頁(yè)，共二十三頁(yè)，2022年，8月28日DATAF;INPUTBREED$y@@;cards;

B121.5B119.5B120.0B122.0B118.0B120.0B216.0B218.5B217.0B215.5B220.0B216.0B319.0B317.5B320.0B318.0B317.0B421.0B418.5B419.0B420.0B515.5B518.0B517.0B516.0;

PROCGLM;

CLASSBREED;

modely=BREED;MEANSBREED/T;RUN;第十五頁(yè)，共二十三頁(yè)，2022年，8月28日DATAF;INPUTA$B$y@@;cards;A1B13.5A1B22.3A1B32.0A2B12.6A2B22.0A2B31.9A3B12.0A3B21.5A3B31.2A4B11.4A4B20.8A4B30.3;

PROCANOVA;CLASSAB;

modely=AB;MEANSAB/DUNCAN;RUN;第十六頁(yè)，共二十三頁(yè)，2022年，8月28日

DATACaP;DOa=1to4;DOb=1to4;DOn=1to3;inputy@@;output;end;end;end;DROPn;CARDS;22.026.524.430.027.526.032.426.527.030.527.025.123.525.827.033.228.530.138.035.533.926.524.025.030.526.825.536.534.033.528.030.524.620.522.519.534.531.429.329.027.528.027.526.328.518.520.019.0;

PROCANOVA;CLASSAB;MODELy=ABA*B;RUN;第十七頁(yè)，共二十三頁(yè)，2022年，8月28日DATAG;DOa=1to3;DOb=1to3;DOn=1to2;inputy@@;output;end;end;end;DROPn;CARDS;82.582.487.186.584.083.986.685.886.285.787.087.682.081.580.080.579.580.3;

PROCANOVA;CLASSAB;

MODELy=AB(A);RUN;第十八頁(yè)，共二十三頁(yè)，2022年，8月28日DATAR;INPUTT$B$X@@;CARDS;A1B1205A1B2168A1B3222A1B4230A2B1230A2B2198A2B3242A2B4255A3B1252A3B2248A3B3305A3B4260A4B1200A4B2158A4B3183A4B4196A5B1265A5B2275A5B3315A5B4282;PROCGLM;CLASSTB;MODELX=TB;MEANST/DUNCAN;RUN;第十九頁(yè)，共二十三頁(yè)，2022年，8月28日DATAL;INPUTT$ROWCOLX@@;CARDS;D1123E1221A1324B1421C1519A2122C2220E2320D2421B2522E3120A3225B3326C3422D3523B4125D4222C4325E4421A4523C5119B5220D5324A5422E5519;PROCGLM;CLASSTROWCOL;MODELX=TROWCOLl;MEANST/DUNCAN;RUN;第二十頁(yè)，共二十三頁(yè)，2022年，8月28日（三）一元線性回歸與相關(guān)分析（教材【例8.1】）

DATAG；INPUTx

y@@；CARDS；

802350862400982720902500120315010226809526308324001133080105292011029601002860；PROC

REG

CORR；MODELy=x

；RUN；第二十一頁(yè)，共二十三頁(yè)，2022年，8月28日（四）協(xié)方差分析（教材【例10.1】）DATAK;INPUTT$XY@@;CARDS;CK1.5012.40CK1.8512.00CK1.3510.80CK1.4510.00CK1.4011.00CK1.4511.80CK1.5012.50CK1.5513.40CK1.4011.20CK1.5011.60CK1.6012.60CK1.7012.5011.351