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

畜禽場(chǎng)育種數(shù)據(jù)統(tǒng)計(jì)分析

——SAS應(yīng)用初步

王楚端鄭江霞一、SAS系統(tǒng)簡(jiǎn)介

—SAS（StatisticalAnalysisSystem）是當(dāng)前使用最為廣泛的統(tǒng)計(jì)分析軟件YSTAT）之一；—1966年美國(guó)NorthCarolina州立大學(xué)開始開發(fā)

—目前國(guó)際上最流行的一種大型統(tǒng)計(jì)分析系統(tǒng)—統(tǒng)計(jì)分析的標(biāo)準(zhǔn)軟件—其基本部分是SAS/BASE軟件，在此基礎(chǔ)上可根據(jù)需要加上SAS系統(tǒng)的其他模塊

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)……SAS系統(tǒng)的組成系統(tǒng)管理程序

功能程序庫(kù)

SAS/BASE-SAS/STATSAS/GRAPHSAS/QCSAS/ORSAS/ETSSAS/IMLSAS/ASSISTSAS系統(tǒng)二、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程序編輯窗口）（二）退出

安裝：農(nóng)大資源僅供學(xué)習(xí)啟動(dòng)fileeditviewtoolssolutionswindowhelp退出出四、、SAS程程序序結(jié)結(jié)構(gòu)構(gòu)（一一））程程序序結(jié)結(jié)構(gòu)構(gòu)在SAS系系統(tǒng)統(tǒng)中中任任何何一一個(gè)個(gè)完完整整的的處處理理過(guò)過(guò)程程均均可可分分為為兩兩大大步步————數(shù)據(jù)據(jù)步步和和過(guò)過(guò)程程步步數(shù)據(jù)據(jù)步步———將將不不同同來(lái)來(lái)源源的的數(shù)數(shù)據(jù)據(jù)讀讀入入SAS系系統(tǒng)統(tǒng)建建立立起起SAS數(shù)數(shù)據(jù)據(jù)集集。。由由DATA語(yǔ)語(yǔ)句句開始始，，以以RUN語(yǔ)句句結(jié)結(jié)束束。?！獢?shù)據(jù)據(jù)可可以以由由外外部部或或過(guò)過(guò)程程讀讀入入；；或或者者由CARDS語(yǔ)句句開始，直接輸入入數(shù)據(jù)，數(shù)據(jù)輸輸入完后后提行單單獨(dú)加““;”結(jié)束簡(jiǎn)單運(yùn)行樣樣例假設(shè)我們有有一個(gè)班學(xué)學(xué)生的數(shù)學(xué)學(xué)成績(jī)和語(yǔ)語(yǔ)文成績(jī)，，數(shù)學(xué)滿分分為100，語(yǔ)文滿滿分為120，希望望計(jì)算學(xué)生生的平均分分?jǐn)?shù)（按百百分制）并并按此排名名，可以在在程序窗口口輸入此程程序：title'0401班學(xué)學(xué)生成成績(jī)排排名';datac0401;inputname$sex$mathchinese;avg=math*0.5+chinese/120*100*0.5;cards;李明男男9298張紅藝藝女女89106王思明明男男8690張聰男男98109劉潁女女80110;run;procprint;run;procsortdata=c0401;bydescendingavg;run;procprint;run;如果果數(shù)數(shù)據(jù)據(jù)要要永永久久保保存存……libnamea'd:\SAS82\SAS82\SAS\MYWORK\';/*mywork子子目目錄錄已已事事先先存存在在*/dataa.aaaa;inputx@@;cards;12345;procprint;varx;run;過(guò)程程步步———調(diào)調(diào)用用SAS系系統(tǒng)統(tǒng)中中已已編編好好的的各各種種過(guò)過(guò)程程來(lái)來(lái)處處理理和和分分析析數(shù)數(shù)據(jù)據(jù)集集中中的的數(shù)數(shù)據(jù)據(jù)。?！棵恳灰粋€(gè)個(gè)過(guò)過(guò)程程步步均均以以PROC語(yǔ)語(yǔ)句句開始始，，RUN語(yǔ)語(yǔ)句句結(jié)結(jié)束束，，并并且且每每個(gè)個(gè)語(yǔ)語(yǔ)句句后后均均以以““;””結(jié)結(jié)束束。。（二二））程程序序的的輸輸入入、、修修改改調(diào)調(diào)試試和和運(yùn)運(yùn)行行SAS程程序序只只能能在在PGM窗窗口口輸輸入入、、修修改改，，可可在在各各行行的的任任何何位位置置開開始始SAS語(yǔ)語(yǔ)句句的的書書寫寫。。一一個(gè)個(gè)語(yǔ)語(yǔ)句句可可以以連連續(xù)續(xù)寫寫在在幾幾行行中中，，一一行行中中也也可可以以同同時(shí)時(shí)寫寫上上幾幾個(gè)個(gè)語(yǔ)語(yǔ)句句，，但但每每個(gè)個(gè)語(yǔ)語(yǔ)句句后后面面必必須須用用““；；””隔隔開開。。當(dāng)一一個(gè)個(gè)程程序序輸輸入入完完后后，，發(fā)發(fā)送送程程序序的的命命令令為為F10功功能能鍵鍵或或SUBMIT。。當(dāng)當(dāng)程程序序發(fā)發(fā)送送到到SAS系系統(tǒng)統(tǒng)后后，，PGM的的程程序序語(yǔ)語(yǔ)句句全全部部自自動(dòng)動(dòng)清清除除，，LOG窗窗口口將將逐逐步步記記下下程程序序運(yùn)運(yùn)行行的的過(guò)過(guò)程程和和出出現(xiàn)現(xiàn)的的錯(cuò)錯(cuò)誤誤信信息息（（用用紅紅色色提提示示錯(cuò)錯(cuò)誤誤））。?！绻^(guò)程程步?jīng)]沒(méi)有錯(cuò)錯(cuò)誤，運(yùn)行行完成成后，，通常常會(huì)在在OUTPUT窗口口打印印出結(jié)結(jié)果—如果果程序序運(yùn)行行出錯(cuò)錯(cuò)，則需要要在PGM窗口口用RECALL（（或F9））命令令調(diào)回回已發(fā)發(fā)送的的程序序進(jìn)行行修改改。保存SAS程序序及其其他窗窗口內(nèi)內(nèi)容激活相相應(yīng)窗窗口，，然后后選擇擇：方式一一：鼠鼠標(biāo)點(diǎn)點(diǎn)擊““存盤盤”工工具按按鈕方式二二：執(zhí)執(zhí)行菜菜單命命令【【File】→→【Save】】/【【Saveas】】可保存存程序序、數(shù)數(shù)據(jù)文文件、、運(yùn)行行結(jié)果果、日日志等等。SAS保存存成純純文本本格式式文件件。打開已已保存存的文文件只能在在ProgramEdit窗口口打開開文件件，選選擇：：方式一一：鼠鼠標(biāo)點(diǎn)點(diǎn)擊““打開開”工工具按按鈕方式二二：執(zhí)執(zhí)行菜菜單命命令【【File】→→【Open】】可打開開以前前保存存的程程序、、數(shù)據(jù)據(jù)文件件、運(yùn)運(yùn)行結(jié)結(jié)果、、日志志等文文件及及其他他文本本格式式的文文件。。對(duì)于定定量資資料的的統(tǒng)計(jì)計(jì)描述述和簡(jiǎn)簡(jiǎn)單推推斷，，SAS提提供了了三個(gè)個(gè)強(qiáng)有有力的的程序序步，，它們們是：：MEANS過(guò)程程提提供單單個(gè)或或多個(gè)個(gè)變量量的簡(jiǎn)簡(jiǎn)單描描述，，對(duì)于于多個(gè)個(gè)變量量，它它的輸輸出格格式緊緊湊，，便于于閱讀讀。UNIVARIATE過(guò)過(guò)程提提供供單個(gè)個(gè)變量量的詳詳細(xì)描描述和和對(duì)其其分布布類型型的檢檢驗(yàn)。。TTEST過(guò)程對(duì)對(duì)變量進(jìn)行行t/u檢檢驗(yàn)。四、定量資料的的統(tǒng)計(jì)描述述和t檢驗(yàn)驗(yàn)PROCMEANS[選擇項(xiàng)項(xiàng)]計(jì)算描述統(tǒng)統(tǒng)計(jì)量的MEANS過(guò)程VAR變量表BY變量表CLASS變量WEIGHT變量OUTPUT選擇項(xiàng)FREQ變量1.N13.STDERR2.NMISS14.T3.MEAN15.PRT4.STD16.SUMWGT5.MIN17.SKEWNESS6.MAX18.KURTOSIS7.RANGE8.SUM9.VAR10.USS11.CSS12.CVPROCUNIVARIATE[選擇項(xiàng)項(xiàng)]描述統(tǒng)計(jì)計(jì)量的UNIVARIATE過(guò)程程VAR變量表BY變量表FREQ變量WEIGHT變量OUTPUTOUT=1.DATA=SAS數(shù)據(jù)集集2.PLOT做莖葉圖圖、盒盒式圖、、正態(tài)概概率圖3.FREQ做頻數(shù)分分布表4.ＮＯＯＲＭＡＡＬ進(jìn)行正態(tài)態(tài)性檢驗(yàn)驗(yàn)/*LT3*/DATAA;INPUTYIELD@@;CARDS;17721519797123159245119119131149152167104161214125175219118192176175951361991161652149515883137801381511871261961342061379897129143179174159194136108101141148168163176102158145173751301491501611551112051311899114214015415216312318714915513120918397119181149254131215111186118150155197116184239160172179151198124179135151168169173181188211197175122159171166175143190213192231163159158159177147194227141169124165;PROCMEANSMEANSTDCVSTDERRMAXDEC=2;PROCCHART;VARYIELD;PROCFREQ;PROCUNIVARIATENORMALPLOT;RUN;數(shù)值值保保留留兩兩位位小小數(shù)數(shù)t檢驗(yàn)驗(yàn)1、、樣本本平平均均數(shù)數(shù)與與總總體體平平均均數(shù)數(shù)的的差差異異顯顯著著性性檢檢驗(yàn)驗(yàn)DATAA;INPUTy@@;y=y-4.50;cards;4.334.623.894.144.784.644.524.484.554.26;PROCMEANSNMEANSTDERRTPRT;RUN;程序說(shuō)明：：樣本平均數(shù)數(shù)與總體平平均數(shù)的差差異顯著性性檢驗(yàn)可調(diào)調(diào)用MEANS過(guò)程。DATA語(yǔ)語(yǔ)句產(chǎn)生臨時(shí)數(shù)數(shù)據(jù)集A，，表明數(shù)據(jù)據(jù)步的開始始；INPUT語(yǔ)句指明讀取變變量y，@@表示讀入一一條觀測(cè)值值后不換行行，連續(xù)讀讀入數(shù)據(jù)，，減少數(shù)據(jù)據(jù)輸入行；；CARDS語(yǔ)句表明以下為為數(shù)據(jù)行，，數(shù)據(jù)行下下的“；””表示數(shù)據(jù)據(jù)行結(jié)束；；PROCMEANS語(yǔ)句指明調(diào)調(diào)用MEANS過(guò)程程對(duì)數(shù)據(jù)集集A進(jìn)行分分析，輸出出樣本含量量N、平均均數(shù)MEAN、平均均數(shù)的標(biāo)準(zhǔn)準(zhǔn)誤STDERR、、學(xué)生氏T值和t值概率PRT；RUN語(yǔ)句句表示過(guò)程步步結(jié)束，開開始運(yùn)行過(guò)過(guò)程步。2、配對(duì)試驗(yàn)資資料的t檢驗(yàn)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;3、非配對(duì)試試驗(yàn)資料料的t檢驗(yàn)DATAC;INPUTbreedy@@;CARDS;1572891120236110128211372501119239111723211042571732821532961682311118;PROCTTEST;CLASSbreed;VARy;RUN;五、、方方差差分分析析固定定模模型型隨機(jī)機(jī)模模型型混合合模模型型常用用的的方方差差分分析析模模型型表表達(dá)達(dá)主效效應(yīng)應(yīng)模模型型Y=ABC互作作效效應(yīng)應(yīng)模模型型Y=ABCA*BA*CB*CA*B*C嵌套套效效應(yīng)應(yīng)模模型型Y=AB(A)3、、方方差差分分析析SAS過(guò)過(guò)程程對(duì)于于一一般般的的方方差差分分析析（（平平衡衡資資料料，，即即各各處處理理重重復(fù)復(fù)數(shù)數(shù)相相等等））可可用用ANOVA過(guò)程；對(duì)對(duì)于非平平衡資料料（各處處理重復(fù)復(fù)數(shù)不等等）的方方差分析析可用GLM過(guò)程。ANOVA過(guò)程調(diào)用用格式PROCANOVA；CLASS效效應(yīng)名；；MODEL因因變量名名=效應(yīng)應(yīng)名；MEANS效效應(yīng)列表表/選項(xiàng)項(xiàng)；LSD最最小顯顯著差數(shù)數(shù)法SNK最最小顯顯著極差差法DUNCN鄧鄧肯氏均均數(shù)比較較檢驗(yàn)T兩兩兩T檢驗(yàn)驗(yàn)ALPHA=0.05檢驗(yàn)驗(yàn)顯著水水平GLM過(guò)程調(diào)用用格式PROCGLM；；CLASS效效應(yīng)名；；MODEL因因變量名名=效應(yīng)應(yīng)名；MEANS效效應(yīng)表/選項(xiàng);LSMEANS效應(yīng)應(yīng)表/選選項(xiàng)；RANDOM；；STDERR最最小二二乘均數(shù)數(shù)標(biāo)準(zhǔn)誤誤TDIFF最小小二乘均均數(shù)差異異顯著性性檢驗(yàn)單因素試試驗(yàn)（各各處理重重復(fù)數(shù)相相等）DATAE;INPUTTR$y@@;cards;A115A116A115A117A118A245A242A250A238A239A330A335A329A331A335A431A428A420A425A430A540A535A531A532A530;PROCanova;CLASSTR;modely=TR;MEANSTR/T;RUN;ClassLevelsValuesTR5A1A2A3A4A5Numberofobservations25DependentVariable:ySumofSourceDFSquaresMeanSquareFValuePr>FModel41905.440000476.36000034.32<.0001Error20277.60000013.880000CorrectedTotal242183.040000R-SquareCoeffVarRootMSEyMean0.87283812.303793.72558730.28000SourceFAnovaSSMeanSquareFValuePr>FTR41905.440000476.36000034.32<.0001tTests(LSD)foryNOTE:ThistestcontrolstheTypeIcomparisonwiseerrorrate,nottheexperimentwiseerrorrate.Alpha0.05ErrorDegreesofFreedom20ErrorMeanSquare13.88CriticalValueoft2.08596LeastSignificantDifference4.9151Meanswiththesameletterarenotsignificantlydifferent.GroupingMeanNTRA42.8005A2B33.6005A5BB32.0005A3C26.8005A4D16.2005A1DATAF;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;2.單因因素試驗(yàn)（（各處理重重復(fù)數(shù)不等等）ClassLevelsValuesBREED5B1B2B3B4B5Numberofobservations25DependentVariable:ySumofSourceDFSquaresMeanSquareFValuePr>FModel446.4983333311.624583335.990.0025Error2038.841666671.94208333CorrectedTotal2485.34000000R-SquareCoeffVarRootMSEyMean0.5448607.5656161.39358718.42000SourceDFTypeISSMeanSquareFValuePr>FBREED446.4983333311.624583335.990.0025SourceDFTypeIIISSMeanSquareFValuePr>FBREED446.4983333311.624583335.990.0025Comparisonssignificantatthe0.05levelareindicatedby***.DifferenceBREEDBetween95%ConfidenceComparisonMeansLimitsB1-B40.5417-1.33482.4181B1-B31.86670.10643.6269***B1-B23.00001.32174.6783***B1-B53.54171.66525.4181***B4-B1-0.5417-2.41811.3348B4-B31.3250-0.62513.2751B4-B22.45830.58194.3348***B4-B53.00000.94455.0555***B3-B1-1.8667-3.6269-0.1064***B3-B4-1.3250-3.27510.6251B3-B21.1333-0.62692.8936B3-B51.6750-0.27513.6251B2-B1-3.0000-4.6783-1.3217***B2-B4-2.4583-4.3348-0.5819***B2-B3-1.1333-2.89360.6269B2-B50.5417-1.33482.4181B5-B1-3.5417-5.4181-1.6652***B5-B4-3.0000-5.0555-0.9445***B5-B3-1.6750-3.62510.2751B5-B2-0.5417-2.41811.3348DATAF;INPUTA$B$y@@;cards;A1B13.5A1B22.3A1B32.0A2B12.6A2B22.0A2B31.9A3B12.0A3B21.5A3B31.2A4B11.4A4B20.8A4B30.3;PROCANOVA;CLASSAB;modely=AB;MEANSAB/DUNCAN;RUN;3.兩兩因素交交叉分組組單獨(dú)觀觀察值ClassLevelInformationClassLevelsValuesA4A1A2A3A4B3B1B2B3DependentVariable:ySumofSourceDFSquaresMeanSquareFValuePr>FModel57.510833331.5021666734.890.0002Error60.258333330.04305556CorrectedTotal117.76916667R-SquareCoeffVarRootMSEyMean0.96674911.581300.2074981.791667SourceDFAnovaSSMeanSquareFValuePr>FA35.289166671.7630555640.950.0002B22.221666671.1108333325.800.0011Duncan'sMultipleRangeTestforyNOTE:ThistestcontrolstheTypeIcomparisonwiseerrorrate,nottheexperimentwiseerrorrate.Meanswiththesameletterarenotsignificantlydifferent.DuncanGroupingMeanNAA2.60003A1B2.16673A2C1.56673A3D0.83333A4DuncanGroupingMeanNBA2.37504B1B1.65004B2BB1.35004B3六、相關(guān)分分析PROCCORR選項(xiàng)項(xiàng)；VAR變變量名；WITH變變量名名；PARTIAL變變量名表；；(PEARSON)SPEARMANNOSIMPLE指明偏相關(guān)關(guān)變量（即即保持固定定的變量）），不能出出現(xiàn)在VAR和和WITH語(yǔ)句中中。DATAG;INPUTxy@@;CARDS;802350862400982720902500120315010226809526308324001133080105292011029601002860;PROCCORR;RUN;TheCORRProcedure2Variables:xySimpleStatisticsvariableNMeanStdDevSumMinimumMaximumx1298.5000012.37666118280.00000120.00000y122721274.936633265023503150PearsonCorrelationCoefficients,N=12Prob>|r|underH0:Rho=0Xyx1.000000.97740<.0001y0.977401.00000<.0001七、線性回歸PROCREG選項(xiàng);MODEL因變量=自自變量/選選項(xiàng);BY變量量；PLOTy軸變變量