第六章

多個樣本均數(shù)比較的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程1第一節(jié)完全隨機設(shè)計資料的方差分析（掌握）第二節(jié)隨機區(qū)組設(shè)計資料的方差分析（掌握）第五節(jié)析因設(shè)計資料的方差分析（了解）兩因素兩水平的析因分析能對SAS程序的輸出結(jié)果作出合理解釋2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程2教學(xué)內(nèi)容方差分析 SAS常用的過程ANOVA過程(AnalysisOfVariance)GLM過程(GeneralLinearModel)。。。2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程3anova過程的一般形式procanova<options>;classvariables;modeldependents=effects</options>;absorbvariables;byvariables;freqvariable;manova<test-options></detail-options>;meanseffects</options>;repeatedfactor-specification</options>;test<h=effects>e=effect;run;procanova語句procanova語句用以調(diào)用anova過程，語句中包含若干控制選項，詳見下表。選項功能與用法data=指定用于分析的輸入數(shù)據(jù)集。manova要求以多變量方式除去包含缺失值的觀測。如果反應(yīng)變量中的任何一個為缺失值，相應(yīng)觀測即被排除出分析過程。multipass要求anova過程在必要時重新讀入數(shù)據(jù)，而不是將數(shù)據(jù)寫入某個中間文件。其作用在于以增加程序運行時間為代價而減少對磁盤空間的占用，這在大樣本數(shù)據(jù)處理時非常有用。namelen=為效應(yīng)名稱指定字符串的長度。須設(shè)置為20~200的正整數(shù)。默認設(shè)置為“namelen=20”。noprint禁止在結(jié)果窗口顯示分析結(jié)果。此選項在僅需要anova過程創(chuàng)建輸出數(shù)據(jù)集時有用。order=指定anova過程對分組變量（class變量）各水平的排序方式，此處的排序方式?jīng)Q定了模型中分組變量各水平所對應(yīng)的參數(shù)。此選項對所有分組變量有效?？稍O(shè)置的值及其含義：“data”——與輸入數(shù)據(jù)集中各水平的排列順序（首次出現(xiàn)的順序）相同；“formatted”——以變量的格式化值作升序排列；“freq”——以包含觀測數(shù)的多少降序排列；“internal”——以變量的非格式化值作升序排列。默認設(shè)置為“order=formatted”。outstat=要求創(chuàng)建用于存儲方差分析表（包括平方和、自由度、f統(tǒng)計量以及p值等）的輸出數(shù)據(jù)集。class語句語句形式：classvariables;class語句用來指定作為分組因素的變量，class變量可為數(shù)值型或字符型。對于anova過程，class語句是必需的（即模型中應(yīng)至少包含一個分組變量），且必須位于model語句之前。class變量的水平由其格式化值所決定，因此可以通過格式化的方法確定class變量的水平。manova語句語句形式：manova<test-options></detail-options>;當model語句定義的模型中包含多個反應(yīng)變量時，使用manova語句要求anova過程執(zhí)行多元方差分析的操作。語句中可設(shè)置兩種形式的選項，即檢驗和細節(jié)選項檢驗選項用來定義所要檢驗的效應(yīng)。細節(jié)選項用來指定如何執(zhí)行檢驗過程以及給出哪些檢驗結(jié)果。model語句的模型表達式語句中以等號連接起來的變量（包括應(yīng)變量和自變量）列表或組合稱為模型表達式。其中“dependents”項代表應(yīng)變量，可為輸入數(shù)據(jù)集中的一個或多個數(shù)值型變量，多個應(yīng)變量之間以空格相分隔?！癳ffects”項為方差分析模型的效應(yīng)項，是由自變量（分組變量）以特定方式組合而成的表達式（自變量表達式），一個模型表達式中可以同時包含多個效應(yīng)項。means語句語句形式：meanseffects</options>;通過使用means語句，anova過程可以對model語句所定義的效應(yīng)計算其各水平下應(yīng)變量的均值及其標準差。還可通過設(shè)置必要的選項實現(xiàn)對指定主效應(yīng)的組間多重比較。同一過程步中可同時使用多條means語句，但均須位于model語句之后。means語句中可設(shè)置眾多的選項，均與指定效應(yīng)的組間多重比較有關(guān)。means語句選項及其功能（2）選項功能及用法lsd同“t”選項。snk對means語句中指定的全部主效應(yīng)執(zhí)行student-newman-keuls組間多重比較過程。regwq對means語句中指定的全部主效應(yīng)執(zhí)行ryan-einot-gabriel-welsch組間多重比較過程。scheffe對means語句中指定的全部主效應(yīng)執(zhí)行scheffe組間多重比較過程。sidak依據(jù)sidak不等式調(diào)整各組均數(shù)，然后進行組間比較的兩兩t檢驗過程。smm同“gt2”選項t在各設(shè)計單元樣本量相等的情況下，執(zhí)行組間均數(shù)的兩兩t檢驗過程，等同于fisher的最小顯著差異（lsd,leastsignificantdifference）檢驗。tukey對means語句中指定的全部主效應(yīng)執(zhí)行組間多重比較的turkey檢驗過程。waller對means語句中指定的全部主效應(yīng)執(zhí)行組間多重比較的waller-duncant檢驗過程。welch要求執(zhí)行welch方差加權(quán)的單因素方差分析過程。對于方差齊性的前提條件，此方法要比一般的方差分析方法更為穩(wěn)健。設(shè)置“welch”選項時model語句中定義的模型必須為單因素模型，否則該選項將被忽略。test語句語句形式：test<h=effects>e=effect;test語句用來進行其它類型的F檢驗，這種檢驗不同于通常方差分析中以誤差均方為分母的F檢驗，我們可以指定此F檢驗中所使用的分母項（誤差項）。當數(shù)據(jù)的誤差結(jié)構(gòu)為非常規(guī)狀態(tài)（如裂區(qū)設(shè)計）時，就必須使用test語句執(zhí)行特定的F檢驗過程。test語句須置于model語句之后?！癶=”指定需要檢驗的效應(yīng)項（作為F檢驗的分子），此效應(yīng)項必須為model語句中所包含的效應(yīng)。語句中可指定多個效應(yīng)項?！癳=”用來指定作為誤差項的效應(yīng)項（作為F檢驗的分母），語句中只能指定一個誤差項。1.兩兩比較常用的方法有SNK、Bonferonni、Dunnett等語句means分類變量/snkmeans分類變量/bonmeans分類變量/dunnett…2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程162.多個處理組和一個對照組的比較---dunnetttest

相應(yīng)的選項為dunnett

語句meansc/dunnett;2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程173.多重比較語句的其他設(shè)置設(shè)置alpha水平默認為0.05如meansc/dunnettalpha=0.01;2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程18第一節(jié)完全隨機設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程20第一節(jié)完全隨機設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程21第一節(jié)完全隨機設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程23完全隨機設(shè)計的方差分析表變異來源SSdfMSF組間(處理組間)

SS組間k-1SS組間/v組間MS組間/

MS組內(nèi)組內(nèi)(誤差)

SS組內(nèi)N-kSS組內(nèi)/v組內(nèi)總SS總N-1第一節(jié)完全隨機設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程26第一節(jié)完全隨機設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程27第二節(jié)隨機區(qū)組設(shè)計資料的方差分析隨機區(qū)組設(shè)計：處理因素（處理組）區(qū)組因素（區(qū)組）步驟除處理因素外其它條件相似的對象歸入一個區(qū)組將區(qū)組內(nèi)的受試對象隨機分配到不同處理組優(yōu)點保證同一區(qū)組內(nèi)的受試對象接受的處理不同實驗效應(yīng)差異主要由處理因素引起2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程28第二節(jié)隨機區(qū)組設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程30第二節(jié)隨機區(qū)組設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程33隨機區(qū)組設(shè)計的方差分析表變異來源SSdfMSF處理組

SS處理k-1SS處理/k-1MS處理/

MS誤差區(qū)組

SS區(qū)組b-1SS區(qū)組/b-1MS區(qū)組/

MS誤差誤差

SS誤差N-k-b+1SS誤差/v誤差總SS總N-1第二節(jié)隨機區(qū)組設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程35第二節(jié)隨機區(qū)組設(shè)計資料的方差分析2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程36第五節(jié)析因設(shè)計資料的方差分析例6-5將20只家兔隨機等分4組，每組5只，進行神經(jīng)損傷后的縫合試驗。處理由兩個因素組合而成，A因素為縫合方法，有兩水平，一水平為外膜縫合，另一水平為束膜縫合

；B因素為縫合后的時間，有兩水平，一水平為縫合后1月

，另一水平為縫合后2月

。試驗結(jié)果為家兔神經(jīng)縫合后的軸突通過率(%)（注：測量指標，視為計量資料）。欲比較不同縫合方法及縫合后時間對軸突通過率的影響，試做析因方差分析。2022/10/24SAS9.1統(tǒng)計軟件應(yīng)用教程37第五節(jié)析因設(shè)計資料的方差分析2022/10/24SAS9.1