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

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

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

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

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

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

SAS/STAT（統(tǒng)計分析）；SAS/ETS（經(jīng)濟(jì)預(yù)測或時間序列分析）；SAS/QC（質(zhì)量控制管理）；SAS/OR（運籌學(xué)）；SAS/IML（矩陣運算）；SAS/GRAPH（繪圖）等—信息存儲簡單（能讀取任何形式的數(shù)據(jù)）—語言編程能力強(qiáng)—對數(shù)據(jù)連續(xù)處理—統(tǒng)計分析方法豐富，使用簡單—SAS過程傻瓜系統(tǒng)……SAS系統(tǒng)的組成系統(tǒng)管理程序

功能程序庫

SAS/BASE-SAS/STATSAS/GRAPHSAS/QCSAS/ORSAS/ETSSAS/IMLSAS/ASSISTSAS系統(tǒng)二、SAS系統(tǒng)的啟動與關(guān)閉

（一）啟動

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

OUTPUT（SAS結(jié)果輸出窗口）LOG（SAS日志窗口，隨著SAS語句的執(zhí)行，顯示出SAS系統(tǒng)的信息和已執(zhí)行的語句）PGM（SAS程序編輯窗口）（二）退出

安裝：農(nóng)大資源僅供學(xué)習(xí)啟動fileeditviewtoolssolutionswindowhelp退出出四、、SAS程程序序結(jié)結(jié)構(gòu)構(gòu)（一一））程程序序結(jié)結(jié)構(gòu)構(gòu)在SAS系系統(tǒng)統(tǒng)中中任任何何一一個個完完整整的的處處理理過過程程均均可可分分為為兩兩大大步步————數(shù)據(jù)據(jù)步步和和過過程程步步數(shù)據(jù)據(jù)步步———將將不不同同來來源源的的數(shù)數(shù)據(jù)據(jù)讀讀入入SAS系系統(tǒng)統(tǒng)建建立立起起SAS數(shù)數(shù)據(jù)據(jù)集集。。由由DATA語語句句開始始，，以以RUN語句句結(jié)結(jié)束束。?！獢?shù)據(jù)據(jù)可可以以由由外外部部或或過過程程讀讀入入；；或或者者由CARDS語句句開始，直接輸入入數(shù)據(jù)，數(shù)據(jù)輸輸入完后后提行單單獨加““;”結(jié)束簡單運行樣樣例假設(shè)我們有有一個班學(xué)學(xué)生的數(shù)學(xué)學(xué)成績和語語文成績，，數(shù)學(xué)滿分分為100，語文滿滿分為120，希望望計算學(xué)生生的平均分分?jǐn)?shù)（按百百分制）并并按此排名名，可以在在程序窗口口輸入此程程序：title'0401班學(xué)學(xué)生成成績排排名';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;過程程步步———調(diào)調(diào)用用SAS系系統(tǒng)統(tǒng)中中已已編編好好的的各各種種過過程程來來處處理理和和分分析析數(shù)數(shù)據(jù)據(jù)集集中中的的數(shù)數(shù)據(jù)據(jù)。?！棵恳灰粋€個過過程程步步均均以以PROC語語句句開始始，，RUN語語句句結(jié)結(jié)束束，，并并且且每每個個語語句句后后均均以以““;””結(jié)結(jié)束束。。（二二））程程序序的的輸輸入入、、修修改改調(diào)調(diào)試試和和運運行行SAS程程序序只只能能在在PGM窗窗口口輸輸入入、、修修改改，，可可在在各各行行的的任任何何位位置置開開始始SAS語語句句的的書書寫寫。。一一個個語語句句可可以以連連續(xù)續(xù)寫寫在在幾幾行行中中，，一一行行中中也也可可以以同同時時寫寫上上幾幾個個語語句句，，但但每每個個語語句句后后面面必必須須用用““；；””隔隔開開。。當(dāng)一一個個程程序序輸輸入入完完后后，，發(fā)發(fā)送送程程序序的的命命令令為為F10功功能能鍵鍵或或SUBMIT。。當(dāng)當(dāng)程程序序發(fā)發(fā)送送到到SAS系系統(tǒng)統(tǒng)后后，，PGM的的程程序序語語句句全全部部自自動動清清除除，，LOG窗窗口口將將逐逐步步記記下下程程序序運運行行的的過過程程和和出出現(xiàn)現(xiàn)的的錯錯誤誤信信息息（（用用紅紅色色提提示示錯錯誤誤））。?！绻^程程步?jīng)]沒有錯錯誤，運行行完成成后，，通常常會在在OUTPUT窗口口打印印出結(jié)結(jié)果—如果果程序序運行行出錯錯，則需要要在PGM窗口口用RECALL（（或F9））命令令調(diào)回回已發(fā)發(fā)送的的程序序進(jìn)行行修改改。保存SAS程序序及其其他窗窗口內(nèi)內(nèi)容激活相相應(yīng)窗窗口，，然后后選擇擇：方式一一：鼠鼠標(biāo)點點擊““存盤盤”工工具按按鈕方式二二：執(zhí)執(zhí)行菜菜單命命令【【File】→→【Save】】/【【Saveas】】可保存存程序序、數(shù)數(shù)據(jù)文文件、、運行行結(jié)果果、日日志等等。SAS保存存成純純文本本格式式文件件。打開已已保存存的文文件只能在在ProgramEdit窗口口打開開文件件，選選擇：：方式一一：鼠鼠標(biāo)點點擊““打開開”工工具按按鈕方式二二：執(zhí)執(zhí)行菜菜單命命令【【File】→→【Open】】可打開開以前前保存存的程程序、、數(shù)據(jù)據(jù)文件件、運運行結(jié)結(jié)果、、日志志等文文件及及其他他文本本格式式的文文件。。對于定定量資資料的的統(tǒng)計計描述述和簡簡單推推斷，，SAS提提供了了三個個強(qiáng)有有力的的程序序步，，它們們是：：MEANS過程程提提供單單個或或多個個變量量的簡簡單描描述，，對于于多個個變量量，它它的輸輸出格格式緊緊湊，，便于于閱讀讀。UNIVARIATE過過程提提供供單個個變量量的詳詳細(xì)描描述和和對其其分布布類型型的檢檢驗。。TTEST過程對對變量進(jìn)行行t/u檢檢驗。四、定量資料的的統(tǒng)計描述述和t檢驗驗PROCMEANS[選擇項項]計算描述統(tǒng)統(tǒng)計量的MEANS過程VAR變量表BY變量表CLASS變量WEIGHT變量OUTPUT選擇項FREQ變量1.N13.STDERR2.NMISS14.T3.MEAN15.PRT4.STD16.SUMWGT5.MIN17.SKEWNESS6.MAX18.KURTOSIS7.RANGE8.SUM9.VAR10.USS11.CSS12.CVPROCUNIVARIATE[選擇項項]描述統(tǒng)計計量的UNIVARIATE過程程VAR變量表BY變量表FREQ變量WEIGHT變量OUTPUTOUT=1.DATA=SAS數(shù)據(jù)集集2.PLOT做莖葉圖圖、盒盒式圖、、正態(tài)概概率圖3.FREQ做頻數(shù)分分布表4.ＮＯＯＲＭＡＡＬ進(jìn)行正態(tài)態(tài)性檢驗驗/*LT3*/DATAA;INPUTYIELD@@;CARDS;17721519797123159245119119131149152167104161214125175219118192176175951361991161652149515883137801381511871261961342061379897129143179174159194136108101141148168163176102158145173751301491501611551112051311899114214015415216312318714915513120918397119181149254131215111186118150155197116184239160172179151198124179135151168169173181188211197175122159171166175143190213192231163159158159177147194227141169124165;PROCMEANSMEANSTDCVSTDERRMAXDEC=2;PROCCHART;VARYIELD;PROCFREQ;PROCUNIVARIATENORMALPLOT;RUN;數(shù)值值保保留留兩兩位位小小數(shù)數(shù)t檢驗驗1、、樣本本平平均均數(shù)數(shù)與與總總體體平平均均數(shù)數(shù)的的差差異異顯顯著著性性檢檢驗驗DATAA;INPUTy@@;y=y-4.50;cards;4.334.623.894.144.784.644.524.484.554.26;PROCMEANSNMEANSTDERRTPRT;RUN;程序說明：：樣本平均數(shù)數(shù)與總體平平均數(shù)的差差異顯著性性檢驗可調(diào)調(diào)用MEANS過程。DATA語語句產(chǎn)生臨時數(shù)數(shù)據(jù)集A，，表明數(shù)據(jù)據(jù)步的開始始；INPUT語句指明讀取變變量y，@@表示讀入一一條觀測值值后不換行行，連續(xù)讀讀入數(shù)據(jù)，，減少數(shù)據(jù)據(jù)輸入行；；CARDS語句表明以下為為數(shù)據(jù)行，，數(shù)據(jù)行下下的“；””表示數(shù)據(jù)據(jù)行結(jié)束；；PROCMEANS語句指明調(diào)調(diào)用MEANS過程程對數(shù)據(jù)集集A進(jìn)行分分析，輸出出樣本含量量N、平均均數(shù)MEAN、平均均數(shù)的標(biāo)準(zhǔn)準(zhǔn)誤STDERR、、學(xué)生氏T值和t值概率PRT；RUN語句句表示過程步步結(jié)束，開開始運行過過程步。2、配對試驗資資料的t檢驗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、非配對試試驗資料料的t檢驗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過過程程對于于一一般般的的方方差差分分析析（（平平衡衡資資料料，，即即各各處處理理重重復(fù)復(fù)數(shù)數(shù)相相等等））可可用用ANOVA過程；對對于非平平衡資料料（各處處理重復(fù)復(fù)數(shù)不等等）的方方差分析析可用GLM過程。ANOVA過程調(diào)用用格式PROCANOVA；CLASS效效應(yīng)名；；MODEL因因變量名名=效應(yīng)應(yīng)名；MEANS效效應(yīng)列表表/選項項；LSD最最小顯顯著差數(shù)數(shù)法SNK最最小顯顯著極差差法DUNCN鄧鄧肯氏均均數(shù)比較較檢驗T兩兩兩T檢驗驗ALPHA=0.05檢驗驗顯著水水平GLM過程調(diào)用用格式PROCGLM；；CLASS效效應(yīng)名；；MODEL因因變量名名=效應(yīng)應(yīng)名；MEANS效效應(yīng)表/選項;LSMEANS效應(yīng)應(yīng)表/選選項；RANDOM；；STDERR最最小二二乘均數(shù)數(shù)標(biāo)準(zhǔn)誤誤TDIFF最小小二乘均均數(shù)差異異顯著性性檢驗單因素試試驗（各各處理重重復(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.單因因素試驗（（各處理重重復(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.兩兩因素交交叉分組組單獨觀觀察值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選項項；VAR變變量名；WITH變變量名名；PARTIAL變變量名表；；(PEARSON)SPEARMANNOSIMPLE指明偏相關(guān)關(guān)變量（即即保持固定定的變量）），不能出出現(xiàn)在VAR和和WITH語句中中。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選項;MODEL因變量=自自變量/選選項;BY變量量；PLOTy軸變變量