SAS軟件的統(tǒng)計系統(tǒng)

SAS—統(tǒng)計分析系統(tǒng)StatisticalAnalysisSystem公共衛(wèi)生學(xué)院信息數(shù)據(jù)處理教學(xué)實驗室SAS軟件簡介1966年美國NorthCarolina大學(xué)的兩位生物統(tǒng)計學(xué)研究生開始研制SAS系統(tǒng)，直至1976年成立了SAS軟件研究所，正式推出了SAS軟件，1985年推出微機版SAS。SAS是用于決策支持的大型集成信息系統(tǒng)，軟件系統(tǒng)最早的功能限于統(tǒng)計分析，現(xiàn)在統(tǒng)計分析功能也仍是它的重要組成部分和核心功能。SAS現(xiàn)在的最新版本為9.0版，大小約為1GB。

在數(shù)據(jù)處理和統(tǒng)計分析領(lǐng)域，SAS系統(tǒng)被譽為國際上的標(biāo)準(zhǔn)軟件系統(tǒng)，并在96～97年度被評選為建立數(shù)據(jù)庫的首選產(chǎn)品，堪稱統(tǒng)計軟件界的巨無霸。經(jīng)過多年的發(fā)展，SAS已經(jīng)被全世界120多個國家和地區(qū)的近三萬家機構(gòu)所采用，直接用戶超過三百萬人，遍及金融、醫(yī)藥衛(wèi)生、生產(chǎn)、運輸、通訊、政府和教育科研等領(lǐng)域。由于SAS系統(tǒng)是從大型機上的系統(tǒng)發(fā)展而來，在設(shè)計上也是完全針對專業(yè)用戶進行設(shè)計，因此操作以SAS語言編程為主，另外也有非編程方式的菜單系統(tǒng)，但人機對話的界面不太友好。學(xué)習(xí)SAS軟件時，需要用戶有基本的統(tǒng)計學(xué)知識，對所選用的統(tǒng)計方法有較清楚的了解，沒有統(tǒng)計學(xué)基礎(chǔ)知識要掌握SAS軟件較為困難。第一節(jié)SAS軟件的使用基礎(chǔ)一、SAS的啟動和退出

啟動SAS雙擊桌面SAS圖標(biāo)開始菜單程序TheSasSystemTheSasSystemforWindowsV*.*直接雙擊SAS應(yīng)用程序退出SAS標(biāo)題欄關(guān)閉按鈕菜單欄FileExit在命令框鍵入“Bye”或“Endsas”標(biāo)題欄、菜單欄、工具欄PROGRAMEDITOR程序編輯窗口LOG 運行過程日志窗口 OUTPUT結(jié)果輸出窗口SAS常用功能鍵

三個窗口間切換F4RECALL(調(diào)回程序)F5PGM(程序)F8SUBMIT(發(fā)送程序)F6LOG(日志)F1HELP(幫助)F7OUTPUT(結(jié)果)F9KEYS(功能鍵列表)二、SAS的工作窗口例：一組體檢資料，包括姓名、性別、年齡、身高、體重，要求分性別計算身高、體重的均數(shù)、標(biāo)準(zhǔn)差、最大值、最小值。姓名性別年齡身高體重ZHANGLINF47156.347.1ZHAOHUAM38172.461.5WANGQANGM41169.264.5LIULIF52158.253.6SHIDONGF39160.148.0KONGYINGM29174.164.6LILINGF37152.342.2GUANFENM32166.260.2MIAOQINGM26180.374.3NIUHUAM31178.280.2TANGJIEF35169.260.8GAOJUNM24176.073.3SUNHONGF27158.349.9三、SAS程序構(gòu)成DATAAA1;INPUTNAME$SEXAGEHEIGHTWEIGHT;CARDS;

ZHANGLIN247156.347.1ZHAOHUA138172.461.5WANGQANG141169.264.5LIULI252158.253.6SHIDONG239160.148.0KONGYING129174.164.6LILING237152.342.2GUANFEN132166.260.2MIAOQING126180.374.3NIUHUA131178.280.2TANGJIE235169.260.8GAOJUN224176.073.3SUNHONG227158.349.9

;PROCMEANS;CLASSSEX;VARHEIGHTWEIGHT;RUN;數(shù)據(jù)步過程步DATASTEP(數(shù)據(jù)步)以DATA語句開始將數(shù)據(jù)讀入SAS系統(tǒng)，建立SAS數(shù)據(jù)集PROCEDURESTEP(過程步)以PROC語句開始調(diào)用各種已編好的過程處理和分析數(shù)據(jù)RUN

(結(jié)束)

整個程序的最后必須要有RUN語句結(jié)束，表示要執(zhí)行以上全部任務(wù)。1.必須在PGM窗口。SAS語句書寫格式自由，每個語句以“；”結(jié)束，可以在一行中寫幾個語句或一個語句寫在幾行中，不區(qū)分字母大小寫。

2.運行SAS程序工具欄執(zhí)行按鈕功能鍵“F8”（Submit

）菜單欄LocalsSubmit命令框鍵入“Submit”命令自動激活OUTPUT窗口，出現(xiàn)程序運行結(jié)果四、SAS程序編寫、運行、調(diào)試、保存3.調(diào)試試SAS程程序沒出現(xiàn)結(jié)果果窗口，或或結(jié)果不理理想時。激活LOG窗口（F6）查看運行過過程及出錯錯信息黑色文字：：源程序及及運行情況況藍(lán)色文字：：程序運行行情況說明明信息紅色文字：：出錯信息息綠色文字：：警告信息息激活PGM窗口（F5）調(diào)回SAS程序（（F4）修改錯誤后后重新發(fā)送送4.保存存SAS程程序或其他他窗口內(nèi)容容激活相應(yīng)窗窗口常用工具欄欄存存盤盤按鈕菜單欄FileSave（或或Saveas））命令框輸入入File命令可以分別保保存程序、、運行結(jié)果果、日志5.打開開已保存的的程序、結(jié)結(jié)果等只能在PGM窗口常用工具欄欄打打開開按鈕菜單欄FileOpen命令框輸入入Include命命令可打開以前前保存的程程序、運行行結(jié)果、日日志等文件件及其他文文本格式的的文件。第二節(jié)SAS數(shù)數(shù)據(jù)集SAS數(shù)據(jù)據(jù)集(SASDataset)數(shù)據(jù)值被安安排在一個個矩陣式的的表格結(jié)構(gòu)構(gòu)中每一行數(shù)據(jù)據(jù)表示一個個觀測（Observation））每一列數(shù)據(jù)據(jù)表示一個個變量（Variable）變量命名：：1~8個個字符，第第1個字符符必須是字字母，后面面可以是是字母、數(shù)數(shù)字和下劃劃線變量類型：：數(shù)值型、、字符型和和日期時間間型字符型變量量名后加““$”表示，日期型:'16OCT98'D，時間型:'9:25:30'T變量缺失值值：用“．”表示一、SAS數(shù)據(jù)集的的種類臨時SAS數(shù)據(jù)集不保存，一一旦退出SAS就消消失。數(shù)據(jù)集名稱稱為1個單單詞，如DATAAA；永久SAS數(shù)據(jù)集保存在用戶戶指定的文文件夾里。。數(shù)據(jù)集名稱稱為2個單單詞，中間間用“.””分開，也也稱二級名名，如DATANEW.AA；LIBNAMENEW‘‘D:\MYDATA’;DATANEW.AA;在數(shù)據(jù)集的的二級名稱稱NEW.AA中；第1個詞NEW是庫標(biāo)記，，代表用戶戶指定的保保存SAS數(shù)據(jù)文件件的位置，，前面要用用LIBNAME語句指定它它代表的文文件夾和位位置。第2個詞AA是保存的SAS數(shù)據(jù)據(jù)集的文件件名。執(zhí)行這個過過程的結(jié)果果將會在D:\MYDATA文件夾夾中產(chǎn)生1個永久SAS數(shù)據(jù)據(jù)文件AA.SD2SD2是自自動產(chǎn)生的的SAS數(shù)數(shù)據(jù)文件默默認(rèn)擴展名名二、創(chuàng)建SAS數(shù)據(jù)據(jù)集1.在PGM窗口創(chuàng)創(chuàng)建數(shù)據(jù)集集DATAAA;命名將要創(chuàng)創(chuàng)建的數(shù)據(jù)據(jù)集名INPUTX1-X4；描述數(shù)據(jù)集集中的變量量名CARDS;以下開始是是數(shù)據(jù)247156.347.1138172.461.5……；數(shù)據(jù)結(jié)束INPUTXY＠＠＠；連連續(xù)讀數(shù)數(shù)據(jù)不換行行INPUTNAME$AGE;字字符型變變量名加＄＄2.在VIEWTABLE窗窗口創(chuàng)建數(shù)數(shù)據(jù)集菜單欄GlobalsManageOpentables打開VIEWTABLE窗口口直接輸入入數(shù)據(jù)保存數(shù)據(jù)：：FileSave(或或Saveas)在Tables右邊邊框中輸入入數(shù)據(jù)集名名稱在Path右邊框中中選默認(rèn)位置置*:\sas\saswork\建建立臨時數(shù)數(shù)據(jù)集用戶指定位位置可以建建立永久數(shù)數(shù)據(jù)集3.用分析析員模塊創(chuàng)創(chuàng)建數(shù)據(jù)集集菜單欄GlobalsAnalyzeSta-tisticalAnalysis打開Analyze(分析員)窗口直接輸輸入數(shù)據(jù)保存數(shù)據(jù)：：FileSaveLibrary：庫庫標(biāo)記名（（代表要保保存文件的的位置置，選Work可以以建臨時數(shù)數(shù)據(jù)集）TableName：數(shù)據(jù)據(jù)集文件名名Description：可可以對文件件進行描述述4.由其它它類型數(shù)據(jù)據(jù)文件轉(zhuǎn)換換成SAS數(shù)據(jù)集菜單欄FileimportSAS可以以直接與下下列格式文文件互相轉(zhuǎn)轉(zhuǎn)換dBASE/FoxBASEFile(*.dbf)Excel4.0/5.0File(*.xls)Lotus1/3/4File(*.wk1/3/4)DelimitedFile(*.*)空空格分隔TabDelimitedFile(*.txt)制制表符分分隔CommaSeparatedVelues(*.csv)逗號分隔隔注意：VisualFoxpro的數(shù)據(jù)表（（*.dbf)在VFP系系統(tǒng)中通過過菜單欄文文件導(dǎo)導(dǎo)出類型選dBASE、FoxBASE、Excel5.0、DelimitedFile均均可以Excel97/2000/2002的文件（（*.xls)在Excel系系統(tǒng)中通通過菜單單欄文文件另另存存為保存類型型選Excel5.0/4.0(*.xls)或或逗號分分隔(*.csv)制表符分分隔(*.txt)均均可以以Excel文件件中要求求第1行行是變量量名另外：在寫字板板、記事事簿、Word或Excel中中的數(shù)數(shù)據(jù)也可可以直接接采用復(fù)復(fù)制、粘粘貼的方方式，通通過PGM窗口口創(chuàng)建SAS數(shù)數(shù)據(jù)集。。通過菜單欄FileExport也可以把把SAS數(shù)據(jù)集集轉(zhuǎn)換成成相應(yīng)的的其它類類型數(shù)據(jù)據(jù)文件5.從從已有的的SAS數(shù)據(jù)集集中創(chuàng)建建新數(shù)據(jù)據(jù)集DATAA1(子子集)；；SETAA(原原數(shù)據(jù)集集)；(1)選擇部分分變量（（列）保留部分分變量：KEEPDATAA1；SETAA；；KEEPX1X2；；去掉部分分變量：DROPDATAA1；SETAA；；DROPX3X4；；賦值新的的變量Y=X3+X4；LY=LOG(Y)；常用SAS函數(shù)數(shù)：LOG(x)以以e為為底的自自然對數(shù)數(shù)LOG10(x)以以10為底底的常用用對數(shù)EXP(x)以以e為底的的反對數(shù)數(shù)SQRT(x)平平方根ABS(x)絕絕對值INT(x)取取整數(shù)(不不四舍五入)MOD(x,y)取取x/y的余數(shù)SIN(x)正正弦COS(x)余余弦三三角函函數(shù)(x是弧弧度)TAN(x)正正切(2)選擇部分觀測測（行）建子集語句：：IFDATAA1；SETAA；；IFSEX=1；；DATAA1；SETAA；；IFAGE>40；IFX1^=3(不等于)IFSEX=1ANDAGE>60IFAGE>20ANDAGE<=40IFAGE<21ORAGE>40(3)數(shù)據(jù)集連接縱向連接A1DATAAA；SETA1A2；A2要求A1和和A2中的的變量相同A1中有10行，A2中中有15行，，則AA有25行。橫向合并A1A2DATAAA；MERGEA1A2；；BYNum；匹配合并時，，要求A1和A2中中必須有1個個相同的變量量，且排列順順序也相同，，作為標(biāo)識變變量用BY語句指定定。(4)數(shù)據(jù)集轉(zhuǎn)置（（行列互換））PROCTRANSPOSE；；BYG；VARX1-X4；OBSX1X2X3123456789OBS_Name_Col1Col2Col3X1147X2258X3369選擇項項有：：OUT=新新數(shù)據(jù)據(jù)集名名（缺缺省按按DATAn自動動命名名）PREFIX=新變變量名名前綴綴（缺缺省按按Coln自自動命命名））NAME=為原原變量量名列列命名名（缺缺省自自動命命名為為_Name_）第三節(jié)節(jié)統(tǒng)統(tǒng)計描描述一、定定量資資料的的統(tǒng)計計描述述簡單描描述：：MEANS過程輸出樣樣本含含量、、均數(shù)數(shù)、標(biāo)標(biāo)準(zhǔn)差差、最最小值值、最最大值值(N、Mean、、StdDev、、Minimum、、Maximum)詳細(xì)細(xì)描描述述：：UNIVARIATE過程程輸出出更更多多指指標(biāo)標(biāo)：：樣樣本本含含量量、、均均數(shù)數(shù)、、標(biāo)標(biāo)準(zhǔn)準(zhǔn)差差、、最最小小值值、、最最大大值值、、方方差差、、標(biāo)標(biāo)準(zhǔn)準(zhǔn)誤誤、、中中位位數(shù)數(shù)、、四四分分位位數(shù)數(shù)間間距距、、常常用用百百分分位位數(shù)數(shù)、、眾眾數(shù)數(shù)、、平平方方和和、、離離均均差差平平方方和和等等，，可可以以作作正正態(tài)態(tài)性性檢檢驗驗另外外還還顯顯示示原原數(shù)數(shù)據(jù)據(jù)中中的的部部分分極極端端值值PROCMEANS;CLASSSEX;VARHEIGHTWEIGHT;RUN;CLASS指定定分分組組變變量量（（可可以以不不寫寫））VAR指定定需需要要分分析析的的變變量量(不不寫寫即即分分析析全全部部變變量量））MAXDEC=n指定定計計算算結(jié)結(jié)果果保保留留小小數(shù)數(shù)位位數(shù)數(shù)(默默認(rèn)認(rèn)2位位)可加加的的選選項項有有：：標(biāo)標(biāo)準(zhǔn)準(zhǔn)誤誤(StdErr)、、方方差差(Var)、、總總和和(Sum)、、平平方方和和(USS)、、離離均均差差平平方方和和(CSS)、、全全距距(Range)、、變變異異系系數(shù)數(shù)(CV)、、偏偏度度系系數(shù)數(shù)(Skewness)、、峰峰度度系系數(shù)數(shù)(Kurtosis)、、可可信信區(qū)區(qū)間間(CLM)等等還可可以以計計算算總總體體均均數(shù)數(shù)=0的的t檢檢驗驗(T、、Prt)不加加任任何何選選項項時時，，自自動動輸輸出出5項項結(jié)結(jié)果果：：N、、Mean、、StdDev、、Minimum、、Maximum頻數(shù)表資資料：X-f2.4-52.6-102.8-123.0-233.2-243.4-223.6-173.8-74.0-3DATAAA2;INPUTXF@@;CARDS;2.552.7102.9123.1233.3243.5223.7173.974.13；PROCMEANSNMEANSTDSTDERRCLM；VARX；FREQF；RUN；數(shù)據(jù)步中中：X是是組中中值，F(xiàn)是頻數(shù)數(shù)過程步中中：VARX；指指定分析析變量是是XFREQF；指定定F是是頻數(shù)數(shù)計算幾何何均數(shù)滴度人數(shù)1:411:851:1661:3221:6471:128101:25641:5125合計40DATAAA3;INPUTXF@@;Y=LOG10(X);CARDS;41851663226471281025645125;PROCMEANS;VARY;FREQF;OUTPUTOUT=CMEAN=M;&DATAB2;SETC;G=10**M;#PROCPRINTDATA=B2;VARG;*RUN;&輸出臨時時數(shù)據(jù)集集C，變變量M=均數(shù)，，#新建數(shù)據(jù)據(jù)集B2，取自自C，求求反對數(shù)數(shù)G=10M*輸出過程程步，輸輸出幾何何均數(shù)GPROCSORT;BYSEX;PROCUNIVARIATE;BYSEX;VARHEIGHTWEIGHT;RUN;BY指定分組組變量（（可以不不寫）SORT用BY之之前必須須先按分分組變量量排序加選項NORMAL可以作正正態(tài)性檢檢驗（W）在已知資資料不服服從正態(tài)態(tài)分布時時或?qū)Y資料的分分布類型型不確定定時可以以用UNIVARIATE過過程作探探索性分分析。UNIVARIATE過程Moments樣本含量量N7頻數(shù)之和和SumWgts7均數(shù)Mean53.55714ΣXSum374.9標(biāo)準(zhǔn)差StdDev10.46564方差Variance109.5295偏度系數(shù)Skewness1.246峰度系數(shù)Kurtosis1.376ΣX2USS20735.75CSS657.1771變異系數(shù)CV19.54106標(biāo)準(zhǔn)誤StdMean3.955638μ=0的t檢驗T:Mean=013.53944Pr>|T|0.0001非0變量個個數(shù)Num^=07>0變量個個數(shù)Num>07符號檢驗M(Sign)3.5Pr>=|M|0.0156符號秩和檢檢驗SgnRank14Pr>=|S|0.0156正態(tài)性檢驗驗W:Normal0.900687Pr<W0.3488Quantiles(Def=5)P100100%Max73.3P9999%73.3P7575%Q360.8P9595%73.3M50%Med49.9P9090%73.3P2525%Q147.1P1010%42.2P00%Min42.2P55%42.2P11%42.2全距Range31.1四分位數(shù)間間距Q3-Q113.7眾數(shù)Mode42.2極端值ExtremesLowestObsHighestObs最小序序號最最大序序號42.2(4)48(3)47.1(1)49.9(7)48(3)53.6(2)49.9(7)60.8(5)53.6(2)73.3(6)定量資料統(tǒng)統(tǒng)計描述———非編程程方式（分分析員模塊塊）(1)打打開分析員員窗口菜單欄GlobalsAnalyzeStatisticalAnalysis(2)打打開SAS數(shù)據(jù)集菜單欄FileOpenbySASname(3)作作簡單統(tǒng)計計分析菜單欄StatisticsDescriptiveSummaryStatistics在Analysis窗口口選選入入要要分分析析的的變變量量在Classification窗口口選選入入分分組組變變量量用Statistics可以以選選更更多多的的統(tǒng)統(tǒng)計計分分析析結(jié)結(jié)果果用Variables可以以對對變變量量有有更更多多的的選選擇擇方方式式選OK完成成二、、定定性性資資料料統(tǒng)統(tǒng)計計描描述述FREQ過過程程可以以產(chǎn)產(chǎn)生生頻頻數(shù)數(shù)表表和和n維維列列聯(lián)聯(lián)表表PROCFREQ；；TABLESA；；RUN；；TABLES指指定定要要分分析析的的定定性性變變量量TABLESA產(chǎn)產(chǎn)生生一一維維頻頻數(shù)數(shù)表表TABLESA*B產(chǎn)產(chǎn)生生二二維維列列聯(lián)聯(lián)表表(A是是行行B是是列列)TABLESA*B*C產(chǎn)產(chǎn)生生三三維維列列聯(lián)聯(lián)表表(A是是分分層層，，B是是行行，，C是是列列)TABLESA*（（BC））產(chǎn)產(chǎn)生生A*B、、A*C兩兩個二二維維表表TABLESX1-X3產(chǎn)產(chǎn)生生X1、、X2、、X3三三個個一一維維表表TABLES(X1-X4)*Y產(chǎn)產(chǎn)生生X1*YX2*YX3*YX4*Y四個個二二維維表表如描描述述以以下下三三個個定定性性變變量量（（258人人））A文文化化程程度度（（1=文文盲盲2=初初中中以以下下3=高高中中4=大大學(xué)學(xué)））B工工作作強強度度（（1=輕輕2=中中3=重重））C月月收收入入（（1=<3002=300-3=500-4=1000-)DATAAA4;INPUTABC@@;CARDS;323414413………;PROCFREQ;TABLESABC;PROCFREQ;TABLESA*B;PROCFREQ;TABLESA*B*C;PROCFREQ;TABLESA*(BC);RUN;一維頻數(shù)數(shù)表CumulativeCumulativeAFrequencyPercentFrequencyPercent頻數(shù)構(gòu)構(gòu)成比累累計頻數(shù)數(shù)累累計計構(gòu)成比比27830.27830.238532.916363.249536.8258100.0二維頻數(shù)數(shù)表A=4Frequency|Percent|CRowPct|ColPct|1|2|3|4|Total---------+--------+---------+---------+--------+1|0|6|37|26|69B|0.00|6.32|38.95|27.37|72.63|0.00|8.70|53.62|37.68||.|75.00|67.27|81.25|---------+--------+---------+---------+--------+2|0|2|18|6|26|0.00|2.11|18.95|6.32|27.37|0.00|7.69|69.23|23.08||.|25.00|32.73|18.75|---------+--------+---------+---------+--------+Total085532950.008.4257.8933.68100.00頻數(shù)構(gòu)成比行構(gòu)成比比列構(gòu)成比比定性資料料統(tǒng)計描描述———非編程程方式（（分析員員模塊））菜單欄GlobalsAnalyzeStatisticalAnalysis菜單欄FileOpenbySASname作一維頻頻數(shù)表菜單欄StatisticsDescriptiveFreguencycounts把變量選選入Frequencies窗口，OK作二維頻頻數(shù)表菜單欄StatisticsTableAnalysis分別把變變量選入入Row、、Column、Strata窗口，OK用Tables可以選每每個格中中出現(xiàn)更更多的結(jié)結(jié)果三、定量量資料頻頻數(shù)表先用MEANS過程找找出最大大值和最最小值MIN=2.5MAX=5.71RANGE=3.21確定組距距(I=0.3)和第第一組段段的下限限(L=2.5)然后用數(shù)數(shù)據(jù)步分分組賦值新的的變量：：Y=INT((X-L)/I)*I+L用FREQ過程程作頻頻數(shù)表表：DATAAA;SETAA5;Y=INT((X-2.5)/0.3)*0.3+2.5;PROCFREQ;TABLESY;RUN;定量資資料作作直方方圖PROCGCHART;VBARX/LEVELS=11SPACE=0;RUN;GCHART過程作作圖VBAR作縱向向條圖圖；；HBAR作橫向向條圖圖常用選選項有有：LEVELS分組的的數(shù)目目WIDTH指定豎豎條的的寬度度，以以字符符寬度度為單單位SPACE指定豎豎條之之間的的距離離，直直方圖圖=0MIDPOINTS=2.65TO5.65BY0.3指指定組組中值值直方圖圖———非編編程方方式（（分析析員模模塊））菜單欄欄GlobalsAnalyzeStatisticalAnalysis菜單欄欄FileOpenbySASname菜單欄欄GraphsHistogram把變量量選入入Analysis窗口，，OK用Display可以增增加更更多的的選項項第四節(jié)節(jié)定定量量資料料假設(shè)設(shè)檢驗驗一、樣樣本均均數(shù)與與總體體均數(shù)數(shù)比較較以往新新生兒兒出生生體重重平均均為3000g，某某醫(yī)院院對產(chǎn)產(chǎn)婦開開展孕孕期保保健指指導(dǎo)后后，測測得15名名新生生兒出出生體體重如如下，，問開開展孕孕期保保健指指導(dǎo)后后對新新生兒兒體重重是否否有影影響？？280034002900275030003100460039002350350026003700380029503200樣本本均均數(shù)數(shù)與與總總體體均均數(shù)數(shù)比比較較用用MEANS過程程DATAAA6;INPUTX@@;Y=X-3000;CARDS;280034002900275030003100460039002350350026003700380029503200;PROCMEANSMEANSTDSTDERRTPRT;VARY;RUN;結(jié)果果MeanStdDevStdErrorTProb>|T|236.67586.60151.461.560.1405P>0.05，，不不拒拒絕絕H0，尚尚不不能能認(rèn)認(rèn)為為開開展展孕孕期期保保健健指指導(dǎo)導(dǎo)對對新新生生兒兒體體重重有有影影響響。二、、配配對對設(shè)設(shè)計計t檢檢驗驗為了了比比較較新新舊舊兩兩個個測測聲聲計計的的測測定定效效果果，，隨隨機機選選定定10個個場場地地，，在在每每地地同同時時用用新新舊舊兩兩個個測測聲聲計計測測量量數(shù)數(shù)據(jù)據(jù)如如下下，，問問新新舊舊兩兩個個測測聲聲計計的的測測量量結(jié)結(jié)果果是是否否不不同同？？場地地12345678910舊：：87657495655563886154新：：86667795605362855955配對設(shè)計t檢驗用MEANS過程DATAAA7;INPUTX1X2@@;Y=X1-X2;CARDS;8786656674779595656055536362888561595455;PROCMEANSMEANSTDSTDERRTPRT;VARY;RUN;結(jié)果MeanStdDevStdErrorTProb>|T|0.902.280.721.250.2440P>0.05，不拒絕H0，尚不能認(rèn)為為新舊兩個測測聲計的測定定結(jié)果有差別別。三、兩樣本均均數(shù)比較試比較克山病病人與健康人人的血磷值(mmol/L)是否不不同？克山病人0.841.051.201.201.391.531.671.801.872.072.11健康人0.540.640.640.750.760.811.161.201.341.351.481.561.87兩樣本均數(shù)比比較用TTEST過程DATAAA8;INPUTX@@;IF_N_>11THENG=1;ELSEG=0;CARDS;0.841.051.201.201.391.531.671.801.872.072.110.540.640.640.750.760.811.161.201.341.351.481.561.87;PROCTTEST;CLASSG;VARX;RUN;在TTEST過程中，，分類變變量G只能是兩兩分類。。條件語句句IF––THEN/ELSEIF后面是條條件的表表達(dá)式，，如果滿滿足條件件就執(zhí)行行THEN后面操作作，否則則(不滿滿足條件件)，就就執(zhí)行ELSE后面的操操作(ELSE可選)。IFX4=0THENX4=.;IFCJ<60THENS+1;IFX4>2ANDX5=1THENY=1;ELSEY=2;_N_是SAS數(shù)據(jù)步步中自定定義的變變量，它它的值表表示觀測測的序號號。IF_N_>11THENG=1;ELSEG=0;也可用G=(_N_>11);賦值分組組變量G：前11個觀觀測G=0，超超過11時G=1。t檢驗驗結(jié)果(1)GNMeanStdDevStdError1111.520.42180.12722131.080.42210.1171(2)Variances方差TDFProb>|T|Unequal不齊2.523921.40.0196Equal齊2.523722.00.0193(3)ForH0:Variancesareequal,方差齊性性檢驗F'=1.00DF=(12,10)Prob>F'=1.00結(jié)論：P<0.05，，拒絕H0，接受H1，可以認(rèn)認(rèn)為兩組組人群血血磷值不不同，克克山病人人血磷值值較高。。兩大樣本均數(shù)數(shù)比較時，如如只有均數(shù)標(biāo)標(biāo)準(zhǔn)差，而無無原始數(shù)據(jù)時時，無法直接接調(diào)用SAS軟件中的某某個過程進行行分析，只能能在數(shù)據(jù)步中中按照計算公公式編寫程序序完成。調(diào)查健康人群群紅細(xì)胞數(shù)(1012/L)資料如如下：男性360人人，均數(shù)466.02，，標(biāo)準(zhǔn)差57.46；女性255人人，均數(shù)417.80，，標(biāo)準(zhǔn)差29.10。試比較男女紅紅細(xì)胞數(shù)有無無差別？DATABB;INPUTN1X1S1N2X2S2;SE=SQRT(S1*S1/N1+S2*S2/N2);U=ABS(X1-X2)/SE;P1=1-PROBNORM(U);P=2*P1;U=ROUND(U,0.001);P=ROUND(P,0.00001);IFP<0.0001THENP=0.0001;FILEPRINT;PUT#2@10U=@25P=;CARDS;360466.0257.46255417.8029.10;RUN;定量資料假設(shè)設(shè)檢驗——非非編程方式（（分析員模塊塊）菜單欄GlobalsAnalyzeStatisticalAnalysis菜單欄FileOpenbySASname(1)樣本均均數(shù)與總體均均數(shù)比較菜單欄StatisticsHypothesisTastsOne-Samplet-testforMeans往Variable窗口選入變量量，Mean窗口輸入μ,OKμ和σ己知時時可選One-Samplez-testforMeansMean窗口輸入μ,StdDev窗口輸入σ。(2)配對t檢驗菜單欄StatisticsHypothesisTastsTwo-SamplePairedt-testforMeans分別選入第一一組和第二組組變量，OK(3)兩樣本本均數(shù)比較t檢驗菜單欄StatisticsHypothesisTastsTwo-Samplet-testforMeans往Dependent窗口選入要分分析的變量往Group窗口選入分類類變量，OK四、單因素方方差分析用于完全隨機機設(shè)計（成組組設(shè)計）資料料例：比較三組組人群載脂蛋蛋白測定值有有無差別？糖尿病IGT異常健康人85.796.0144.0105.2124.5117.0109.5105.1110.096.076.4109.0115.295.3103.095.3110.0123.0110.095.2127.0106.599.0121.0125.6120.0115.0100.0159.0111.0單因素方差分分析用ANOVA過程DATAAA9;INPUTX@@;IF_N_<12THENG=1;ELSEIF_N_>20THENG=3;ELSEG=2;CARDS;85.7105.2109.596.0115.295.3……；PROCMEANS;CLASSG;VARX;PROCANOVA;CLASSG;MODELX=G;MEANSG/SNK;RUN;在SNK后后面面加加HOVTEST可以以作作方方差差齊齊性性檢檢驗驗AnalysisofVarianceProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>F變異異來來源源自自由由度度SSMSFPModel模型型22384.025505051192.012752535.850.0077Error誤差275497.83616162203.62356154Total總297881.86166667R-SquareR2C.V.RootMSESy.xXMean0.30247012.9351914.26967279110.31666667SourceDFAnovaSSMeanSquareFValuePr>F變異來源源自自由度SSMSFPG組間22384.025505051192.012752535.850.0077P<0.01，，拒絕H0，接受H1，可以認(rèn)認(rèn)為三組組人群載載脂蛋白白不同。。Student-Newman-Keulstestforvariable:X均數(shù)間兩兩兩比較較q檢驗驗（SNK法））Meanswiththesameletterarenotsignificantlydifferent.相同字母母表示的的均數(shù)無無顯著性性差別SNKGroupingMeanNGA122.800103B105.455111BB102.38992第三組高高于前兩兩組，即即健康人人的載脂脂蛋白高高于兩組組病人。。單因素方方差分析析時，如如只有各各組的均均數(shù)、標(biāo)標(biāo)準(zhǔn)差，，而無原原始數(shù)據(jù)據(jù)時，無無法直接接調(diào)用SAS軟軟件中的的某個過過程進行行分析，，只能在在數(shù)據(jù)步步中按照照計算公公式編寫寫程序完完成。例：研究不不同季節(jié)人人體腎上腺腺皮質(zhì)3ββ-HSD活性變化化季節(jié)NMeansStd春季420.780.13夏季400.690.22秋季320.680.14冬季360.580.20DATAAA;INPUTKXX@;DOI=1TOK;INPUTNXS@@;OUTPUT;SS1+(XX-X)*(XX-X)*N;SS2+S*S*(N-1);NN+N;END;SS1=SS1/(K-1);SS2=SS2/(NN-K);F=SS1/SS2;P=1-PROBF(F,K-1,NN-K);F=ROUND(F,0.001);P=ROUND(p,0.00001);IFP<0.0001THENP=0.0001;FILEPRINT;PUT#2@10F=@25P=;CARDS;40.6867(總均均數(shù)）420.780.13400.690.22320.680.14360.580.20;RUN;五、兩因素素方差分析析用于隨機區(qū)區(qū)組設(shè)計(配伍組設(shè)設(shè)計)資料料例：取同品品系同體重重24只小小鼠分為8組，每組組3只小鼠鼠來自同一一窩，分別別用A、B、C三種種營養(yǎng)素喂喂養(yǎng)，問三三周后所增增體重有無無差別。區(qū)組ABC150.158.264.5247.848.562.4353.153.858.6463.564.272.5571.268.479.3641.445.738.4761.953.051.2842.239.846.2兩因素方差差分析用ANOVA過程DATAAA10;DOI=1TO3；DOJ=1TO8；INPUTX@@;OUTPUT；END；END；CARDS;50.147.853.163.571.241.461.942.258.248.553.864.268.445.753.039.864.562.458.672.579.338.451.246.2；PROCMEANS;CLASSI;VARX;PROCANOVA;CLASSIJ;MODELX=IJ;RUN;AnalysisofVarianceProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel模型92521.29375280.1437511.520.0001Error誤差14340.5425024.32446Total總232861.83625R-SquareC.V.RootMSEXMean0.8810068.8605144.9319838155.6625SourceDFAnovaSSMeanSquareFValuePr>FI處理間2144.917572.458752.980.0836J區(qū)組間72376.37625339.4823213.960.0001處理間P>0.05，不拒絕絕H0，尚不能認(rèn)認(rèn)為三種營營養(yǎng)素喂養(yǎng)養(yǎng)小鼠增重重有差別；；區(qū)組間P<0.01，拒絕絕H0，接受H1，可以認(rèn)為為8組小鼠鼠體重增加加有差別。。循環(huán)語句在數(shù)據(jù)步的的執(zhí)行中反反復(fù)地執(zhí)行行某一部分分語句。DO循環(huán)變量=初值TO終值BY增值；END；增量=1時時可省略。。DATAB1;DATAB2;DOI=1TO10;DOI=1TO10BY2;S+I;S+I;OUTPUT;OUTPUT;END;END;PROCPRINT;PROCPRINT；RUN;RUN；循環(huán)語句可可以嵌套（（多重循環(huán)環(huán)）如打印乘法法表：1*1=1DATABB;2*1=22*2=4DOI=1TO9;3*1=33*2=63*3=9DOJ=1TOI;4*1=44*2=84*3=124*4=16X=I*J;5*1=55*2=105*3=15………OUTPUT;6*1=66*2=126*3=18………END;7*1=77*2=147*3=21………END;8*1=88*2=168*3=24………PROCPRINT;9*1=99*2=189*3=27………RUN;單因素方差差分析的數(shù)數(shù)據(jù)也可以以用雙重循循環(huán)語句編編寫，在各各組樣本含含量不等時時，加缺失失值“.”補滿。。糖尿病IGT異常健康人85.796.0144.0105.2124.5117.0109.5105.1110.096.076.4109.0115.295.3103.095.3110.0123.0110.095.2127.0106.599.0121.0125.6120.0115.0100.0159.0111.0DATABB;DOI=1TO11;DOJ=1TO3;INPUTX@@;OUTPUT;END;END;CARDS;85.796144105.2124.5117109.5105.11109676.4109115.295.310395.311012311095.2127106.599121125.6120115100.159111..;PROCANOVA;CLASSJ;MODELX=J;RUN;六、三因素素方差分析析用于拉丁方方設(shè)計資料料例：比較動動物對不同同藥物和用用藥次序的的反應(yīng)拉丁方設(shè)計計要求三因因素水平數(shù)數(shù)相同，因因有四種藥藥物A、B、C、D，故取4只動物，，反復(fù)用藥藥4次。動物編號用藥次序1234ⅠA75B29C27D42ⅡB45D71A81C53ⅢC25A71D80B23ⅣD47C31B48A87DATAAA11;DOCX=1TO4;DODW=1TO4;INPUTYW$X@@;OUTPUT;END;END;CARDS;A75B29C27D42B45D71A81C53C25A71D80B23D47C31B48A87;PROCANOVA;CLASSCXDWYW;MODELX=CXDWYW;MEANSYW/SNK;RUN;DependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel96396.5625710.72924.020.0521Error61059.8750176.6458Total157456.4375R-SquareC.V.RootMSEXMean0.85785825.4674313.2908176352.1875SourceDFAnovaSSMeanSquareFValuePr>FCX3773.1875257.72921.460.3167DW3270.6875090.22920.510.6894YW35352.687501784.229210.100.0092四種藥物之間間P<0.01，拒絕H0，接受H1，可以認(rèn)認(rèn)為動物對四四種藥物的反反應(yīng)不同，A>B、C。。而動物之間間和用藥次序序均無影響。。七、析因試驗驗設(shè)計的方差差分析可以分析各因因素之間的交交互作用。例：用新、舊舊兩法分別提提取某食品中中甲、乙兩種種化合物，觀觀察回收率有有無差別（包包括方法間、、化合物間及及方法與化合合物間交互作作用）。新法(A1)舊法(A2)甲物(B1)乙物(B2)甲物(B1)乙物(B2)52845247488844644490405244882645合計188350162208DATAAA12;DOA=1TO2;DOB=1TO2;DOI=1TO4;INPUTX@@;OUTPUT;END;END;END;CARDS;52484444848890885244402647645245;PROCANOVA;CLASSAB;MODELX=ABA*B;RUN;DependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel35309.001769.6733.390.0001Error12636.0053.00Total155945.00R-SquareC.V.RootMSEXMean0.89301912.828397.2856.75SourceDFAnovaSSMeanSquareFValuePr>FA11764.001764.0033.280.0001B12704.002704.0051.020.0001A*B1841.00841.0015.870.0018不同方方法、、不同同化合合物間間均有有差別別，且且交互互作用用也有有統(tǒng)計計學(xué)意意義。。新法高高于舊舊法，，乙物物高于于甲物物，以以新法法+乙乙物為為最佳佳組合合。八、交交叉試試驗設(shè)設(shè)計的的方差差分析析研究新新藥催催醒寧寧對麻麻醉劑劑的作作用，，將6只大大鼠配配成3對，，隨機機決定定每對對中1只接接受處處理的的順序序（催催醒寧寧或生生理鹽鹽水）），等等動物物完全全清醒醒3天天后將將實驗驗組和和對照照組動動物互互換，，同法法重復(fù)復(fù)試驗驗一次次。結(jié)結(jié)果如如下：：（Y藥物物組D對照照組）動物編號123456第1次Y15D31Y17D30D25Y5第2次D27Y25D28Y14Y18D26DATAAA13;DOCX=1TO2;DODW=1TO6;INPUTYW$X@@;OUTPUT;END;END;CARDS;Y15D31Y17D30D25Y5D27Y25D28Y14Y18D26;PROCANOVA;CLASSCXDWYW;MODELX=CXDWYW;RUN;DependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel7621.583388.79765.850.0535Error460.666715.1667Total11682.2500R-SquareC.V.RootMSEXMean0.91107917.905473.894421.75SourceDFAnovaSSMeanSquareFValuePr>FCX118.7518.751.240.3285DW5158.7531.752.090.2470YW1444.08444.0829.280.0057藥物間P<0.01，拒拒絕H0，接接受H1，可可以認(rèn)為實驗驗組與對照組組不同，催醒醒寧有效。方差分析———非編程方式式（分析員模模塊）菜單欄GlobalsAnalyzeStatisticalAnalysis菜單欄FileOpenbySASname(1)單因素素方差分析菜單欄StatisticsANOVAOne-wayANOVA變量X選入Dependent窗口，分組G選入Independent窗口。(2)兩因素素、三因素及及交叉設(shè)計的的方差分析菜單欄StatisticsANOVAFactorialANOVA變量X選入Dependent窗口，分組選選入Independent窗口。需要作均數(shù)間間兩兩比較時時選Means，，在ComparisonMathod窗口中選兩兩兩比較的方法法，在MainEffacts窗口選中分組組變量，AddOKOK析因設(shè)計的方方差分析：交交互作用分析析選Model,在StandardModels窗口中選Two-way(兩因素)或或Three-way(三因素)第五節(jié)定性性資料假設(shè)檢檢驗一、樣本率與與總體率比較較沒有SAS過過程可以直接接調(diào)用。可以按照二項項分布或Poisson分布的概率率函數(shù)在數(shù)據(jù)據(jù)步中編程計計算。己知某藥治愈愈率為60%。某醫(yī)生隨隨機抽取10名病人試用用此藥加VC同服，結(jié)果果有8人治愈愈，問加服VC是否有增增效作用？該題己知π=0.6,n=10,x=8,求n≥8的概率?？梢杂枚椃址植己瘮?shù)PROBBNML(π，，n,x)DATABB;Q=PROBBNML(0.6,10,7);P=1-Q;PROCPRINT;RUN;結(jié)果：OBSQP10.832710.16729P>0.05，不拒絕H0，尚不能認(rèn)為為VC有增效效作用(1)二項分分布直接概率率法(2)Poisson分分布直接概率率法國家規(guī)定飲用用水中大腸菌菌數(shù)平均不得得超過3個/L。現(xiàn)從某某水源地抽取取水樣1L測測得大腸菌4個，問該水水源是否符合合飲用水條件件？該資料服從Poisson分布，己己知μ=3，，求X≥3的的概率?？梢砸杂肞oisson分布布函數(shù)POISSON(μ,x)求得。DATABB;Q=POISSON(3,3);P=1-Q;PROCPRINT;RUN;結(jié)果：OBSQP10.647230.35277P>0.05，不拒絕H0，尚不能認(rèn)為為細(xì)菌數(shù)超標(biāo)標(biāo)。(3)大樣樣本正態(tài)近似似法—u檢驗驗一般胃潰瘍患患者中有20%會發(fā)生胃胃出血癥狀。。某醫(yī)醫(yī)院觀察65歲以上潰瘍瘍病人304例，有96人發(fā)生胃出出血。問老年年胃潰瘍患者者是否容易出出血？己知π=0.2，n=304，x=96，可以以用u檢驗。。DATABB;INPUTPNX;S=SQRT(P*(1-P)/N);U=ABS(X/N-P)/S;P1=1-PROBNORM(U);P=2*P1;U=ROUND(U,0.001);P=ROUND(P,0.00001);IFP<0.0001THENP=0.0001;FILEPRINT;PUT#2@10u=@25P=;CARDS;0.230496;RUN;結(jié)結(jié)果果：u=5.047P=0.0001二、樣本率之之間比較—χχ2檢驗用兩種藥物治治療十二腸潰潰瘍療效比較較藥物愈合未愈合合計愈合率呋喃硝胺5486287.10甲氰咪胍44206468.75合計982812677.78DATABB;DOI=1TO2;DOJ=1TO2;INPUTF@;OUTPUT;END;END;CARDS;5484420;PROCFREQ;TABLESI*J/CHISQ;WEIGHTF;RUN;WEIGHT語句指明變變量F是是頻數(shù)加選項EXACT可以對>2×2表作作Fisher確切切概率EXPECTED可以顯示每每個格子的的理論值CMH可以計算隊隊列研究和和病例對照照研究的RR值NOROW、NOCUL、NOPERCENT不顯示相應(yīng)應(yīng)百分?jǐn)?shù)NOPRINT不顯示表格格，只顯示示計算的統(tǒng)統(tǒng)計量χ2檢驗用FREQ過程程加CHISQ選項STATISTICSFORTABLEOFIBYJStatistic統(tǒng)計量量DFValueProbChi-Squareχ216.1330.013LikelihoodRatioChi-Square似然比χ216.3040.012ContinuityAdj.Chi-Square校正正χ215.1180.024Mantel-HaenszelChi-SquareMHχ216.0840.014Fisher‘‘sExactTest確切切概概率率(Left)0.997(Right)0.011(2-Tail)0.018PhiCoefficient(Phi系系數(shù)數(shù))0.221ContingencyCoefficient(列列聯(lián)聯(lián)系系數(shù)數(shù))0.215Cramer‘‘sV(克克萊萊姆姆系系數(shù)數(shù))0.221SampleSize=126由原原始始數(shù)數(shù)據(jù)據(jù)計計算算：：如258人人資料料，三三個定定性變變量A文化化程度度（1=文文盲2=初初中以以下3=高高中4=大大學(xué)））B工作作強度度（1=輕輕2=中3=重重）C月收收入（（1=<3002=300-3=500-4=1000-)試分析析文化化程度度與月月收入入之間間的關(guān)關(guān)系DATAAA4;INPUTABC@@;CARDS;323414413………;PROCFREQ;TABLESA*C/CHISQ；RUN；三、配配對定定性資資料的的χ2檢驗用甲乙乙兩種種血清清學(xué)方方法檢檢查鼻鼻咽癌癌患者者血清清410份份。問問兩法法的檢檢出率率有無無差別別？甲法乙法+-合計+26187348-313162合計292118410DATABB;DOI=1TO2;DOJ=1TO2;INPUTF@;OUTPUT;END;END;CARDS;261873131;PROCFREQ;TABLESI*J/AGREE;WEIGHTF;RUN;AGREE選項作定性資資料的一致性性檢驗，計算算χ2和Kappa統(tǒng)計量，Kappa<0.4不滿滿意，>0.75很滿意意。配對四格表資資料計算McNemar‘sχχ2統(tǒng)計量配對R×C表表(方表)計計算Bowker對稱性性檢驗QB統(tǒng)計量，QB近似服從DF=R(R-1)/2的的χ2分布。當(dāng)P<0.05，拒絕H0時，結(jié)論為不不服從對稱分分布。結(jié)果：McNemar'sTestStatistic=26.576DF=1Prob=0.001SimpleKappaCoefficient95%ConfidenceBoundsKappa=0.182ASE=0.0510.0830.281SampleSize=410P<0.01，拒絕H0，接受H1，可以認(rèn)為甲甲、乙兩法的的檢出率不同同，甲法檢出出率高于乙法法。Kappa值值的假設(shè)檢驗驗：H0：Kappa=0u=Kappa/ASE=0.182/0.051=3.57P<0.01拒絕H0，接受H1，雖然可以認(rèn)認(rèn)為甲、乙兩兩法有一致性性，但一致的的程度很差。。兩種方法檢查查冠心病人心心室壁收縮運運動的符合情情況核素法對比法正常減弱異常正常5823減弱1427異常8917結(jié)果：TestofSymmetry(對稱性檢驗驗)Statistic=2.856DF=3Prob=0.414KappaCoefficients(一致性檢驗)StatisticValueASE95%ConfidenceBoundsSimpleKappa0.6810.0500.5830.779WeightedKappa0.6610.0560.5520.771SampleSize=147對稱性檢驗：P>0.05，不拒絕絕H0，認(rèn)為為該資料服從從對稱分布。。一致性檢驗：u=Kappa/Ase=13.6，P<0.01，可可以認(rèn)為兩種種方法檢查結(jié)結(jié)果具有一致致性，且一致致性程度比較較高。χ2檢驗——非編程方式（（分析員模塊塊）菜單欄GlobalsAnalyzeStatisticalAnalysis在右面窗口內(nèi)內(nèi)輸入3列數(shù)數(shù)據(jù)：A列輸輸入行變量；；B列輸入列變變量代碼；C列輸入相應(yīng)應(yīng)的頻數(shù)。或打開己有的的SAS數(shù)據(jù)據(jù)集菜單欄StatisticsTableAnalysis把變量A選入Row，B選入Column，C選入Cellcounts，打開Statistics窗口Chi-sqareStatistics計算χ2統(tǒng)計量Measuresofassociation計算行列變量量的關(guān)聯(lián)性統(tǒng)統(tǒng)計量Mantel-HaenszelStatistics隊列研究和病病例對照研究究RR值MeasuresofAgreement一致性檢驗選中所要的項項目OKOKABC115412821442220第六節(jié)秩秩和檢驗(1)配對對設(shè)計資料UNIVARIATE過過程可以作符號檢檢驗和符號秩秩和檢驗例：8名男性性服某藥前后后精液中精子子濃度(萬/ml)變化化病人服藥前服藥后1600066022200056003590037004440050005600063006650012007260001800858002200DATAAA14;INPUTX1X2@@;D=X1-X2;CARDS;60006602200056005900370044005000600063006500120026000180058002200;PROCUNIVARIATE;VARD;RUN;MomentsN8SumWgts8Mean7017.5Sum56140StdDev8746.372Variance76499021Skewness1.39594Kurtosis1.066563USS9.2946E8CSS5.3549E8CV124.6366StdMean3092.309T:Mean=02.26934Pr>|T|0.0575Num^=08Num>06M(Sign)符號檢驗2Pr>=|M|0.2891SgnRank符號秩和檢驗驗15Pr>=|S|0.0391(2)成組設(shè)設(shè)計資料成組設(shè)計資料料秩和檢驗用用NPAR1WAY過程兩組比較時，，可作Wilcoxon法計算u值值多組比較時，，可作Kruskal-Walli