統(tǒng)計(jì)分析系統(tǒng)－－SAS簡介3

1、統(tǒng)計(jì)分析系統(tǒng)SAS6、方差分析概述對(duì)于兩個(gè)平均數(shù)的假設(shè)測驗(yàn)，一般可以采用t測驗(yàn)的方法完成。對(duì)于多個(gè)平均數(shù)的假設(shè)測驗(yàn)，由于測驗(yàn)程序繁雜、誤差估計(jì)不精確，尤其是犯第一類錯(cuò)誤的概率增加等原因，一般不再利用t測驗(yàn)兩兩進(jìn)行，而采用一種新的統(tǒng)計(jì)方法方差分析法（ANALYSIS OF VARIANCE，簡計(jì)為ANOVA）。6.1 概述概述方差分析的基本思想是將所有觀察值的總變異分解成不同的變異來源，即對(duì)總變異的自由度和平方和進(jìn)行分解，進(jìn)而獲得不同變異來源的總體方差的估值。通過構(gòu)建適當(dāng)?shù)腇值，進(jìn)行F測驗(yàn)，完成多個(gè)樣本平均數(shù)之間差異顯著性測驗(yàn)。當(dāng)處理效應(yīng)為固定效應(yīng)時(shí)，尚可對(duì)各個(gè)平均數(shù)進(jìn)行多重比較。當(dāng)處理效應(yīng)為隨

2、機(jī)效應(yīng)時(shí)，進(jìn)行方差分量的估計(jì)。6.1 概述概述方差分析是生物領(lǐng)域中應(yīng)用最為廣泛的統(tǒng)計(jì)方法之一，在科學(xué)研究和生產(chǎn)實(shí)踐中有著極其重要的用途。在SAS中用于方差分析的主要過程有方差分析（ANOVA）和廣義線性模型（GLM）。對(duì)于平衡資料（數(shù)據(jù)均衡，沒有發(fā)生缺?。话銘?yīng)用ANOVA過程。對(duì)于非平衡資料，應(yīng)采用GLM過程。6.2 常用過程的使用方法1. ANOVA過程過程格式Proc anova 選項(xiàng);Class 變量表;Model 依變量效應(yīng)表/選項(xiàng);Means 效應(yīng)表/選項(xiàng);Freq 變量;Run;6.2 常用過程的使用方法語句說明：PROC ANOVA 選項(xiàng)DATA輸入數(shù)據(jù)集OUTSTAT輸出數(shù)

3、據(jù)集用于存儲(chǔ)方差分析結(jié)果CLASS語句用于指明分類變量。此語句一定要設(shè)定，并且應(yīng)出現(xiàn)在model語句之前。MODEL 語句定義分析所用的線性數(shù)學(xué)模型。6.2 常用過程的使用方法Means語句計(jì)算效應(yīng)平均數(shù)。選項(xiàng)用于設(shè)定多重比較方法等。多重比較方法選擇有T或LSD（配對(duì)t測驗(yàn)或Fisher氏最小顯著差數(shù)法），DUNCAN（Duncan氏新復(fù)極差測驗(yàn)法），TUCKY（Tukey氏固定極差測驗(yàn)法）和DUNNETT（Dunnett氏新復(fù)極差測驗(yàn)法）等。顯著水平的確定采用ALPHA設(shè)定，例如ALPHA0.01，將顯著水平設(shè)定為0.01，缺省時(shí)為0.05。6.2 常用過程的使用方法線性數(shù)學(xué)模型上述語句中，

4、關(guān)鍵語句在于定義線性數(shù)學(xué)模型。同一試驗(yàn)資料，分析結(jié)果依模型不同而異。因而應(yīng)依據(jù)試驗(yàn)設(shè)計(jì)，給出正確的線性數(shù)學(xué)模型。model y=a; 單因素方差分析。Model y=a b; 兩因素主效模型。Model y=a b a*b; 兩因素帶互作模型。Model y=a b(a)；嵌套（nested）模型，用于系統(tǒng)分組資料。模型定義中可用|和n簡化表達(dá)方法ANOVA模型的應(yīng)用實(shí)例6.2 常用過程的使用方法2、GLM過程GLM過程的應(yīng)用范圍較廣，除了用于方差分析外，還可用于回歸分析、協(xié)方差分析等。當(dāng)用于方差分析時(shí)，主要處理非平衡資料，可進(jìn)行多個(gè)變量的對(duì)比檢驗(yàn)。與ANOVA過程相比，GLM過程運(yùn)行時(shí)一般占

5、用較多的計(jì)算機(jī)內(nèi)存和時(shí)間。2、GLM過程GLM過程格式Proc glm 選項(xiàng);Class 變量表;Model 依變量效應(yīng)表/選項(xiàng);Means 效應(yīng)表/選項(xiàng);Random 效應(yīng)表/選項(xiàng);Run;2、GLM過程語句說明Proc glm語句設(shè)定分析數(shù)據(jù)集和輸出數(shù)據(jù)集等。Class語句指定分類變量，此語句必須設(shè)定，并且應(yīng)出現(xiàn)在model語句之前。Model語句用于定義數(shù)學(xué)模型和結(jié)果輸出項(xiàng)。2、GLM過程語句說明：Means語句計(jì)算平均數(shù)，并可選用多種方法進(jìn)行多重比較。Random語句指定模型中的隨機(jī)效應(yīng)。在SAS中，數(shù)據(jù)行中的小數(shù)點(diǎn)“.”代表缺失數(shù)據(jù)。GLM過程應(yīng)用實(shí)例6.3應(yīng)用舉例兩向分組資料的方差

6、分析2.4.2 兩向分組資料的方差分析6.3 應(yīng)用舉例兩向分組資料的方差分析（ANOVA過程）二因素又重復(fù)資料的方差分析二因素又重復(fù)資料的方差分析二因素有重復(fù)資料的方差分析（ANOVA過程）單因素隨機(jī)區(qū)組試驗(yàn)二因素隨機(jī)區(qū)組試驗(yàn)資料的方差分析協(xié)方差分析協(xié)方差分析解決的典型問題如不同施肥量對(duì)果樹增產(chǎn)效果的影響的研究。通常情況下，果樹的基礎(chǔ)生產(chǎn)力（x，可用上一年的產(chǎn)量表示）不同，因而試驗(yàn)結(jié)果（y，果樹的產(chǎn)量）除了有肥料效應(yīng)之外，一般還包含有每棵果樹的基礎(chǔ)生產(chǎn)力x不同造成的差異。要真正反映施肥量的效應(yīng)，應(yīng)先將變量x對(duì)變量y的影響剔除，再進(jìn)行方差分析，這種分析即是協(xié)方差分析。協(xié)方差分析類似于方差分析中總

7、變異自由度和平方和的分解，協(xié)方差分析將總協(xié)方同樣分解，即分解為不同變異來源的乘積和與自由度，利用回歸分析的原理分析變量y與變量x之間的關(guān)系，從而將方差分析和相關(guān)、回歸分析結(jié)合起來。協(xié)方差分析的SAS過程為廣義線性模型（GLM）。協(xié)方差分析協(xié)方差分析分析方法先對(duì)x、y進(jìn)行方差分析，檢驗(yàn)其差異顯著性，然后對(duì)處理內(nèi)項(xiàng)（誤差）做回歸分析，測驗(yàn)x與y是否存在線性關(guān)系。若無，表明x對(duì)y無影響，對(duì)y作方差分析的結(jié)果即能說明3中肥料對(duì)蘋果樹茶兩的實(shí)際效應(yīng)。若x和y之間存在線性關(guān)系，說明基礎(chǔ)生產(chǎn)力對(duì)實(shí)驗(yàn)結(jié)果有影響，因此不能用原有的y值進(jìn)行方差分析，必須消去x的不同對(duì)y帶來的影響，即通過求y依x的線性回歸方程，將

8、各處理的yi都矯正到x在同意水平時(shí)的值。最后對(duì)矯正平均數(shù)做方差分析，比較3種肥料對(duì)蘋果樹產(chǎn)量的影響有無顯著差異。協(xié)方差分析GLM過程過程格式PROC GLM 選項(xiàng)；CLASS 變量表；MODEL 依變量=效應(yīng)/選項(xiàng)；MEANS 效應(yīng)/選項(xiàng)；LSMEANS 效應(yīng)/選項(xiàng)；RUN；協(xié)方差分析語句說明PROC GLM語句選項(xiàng)為可設(shè)定分析數(shù)據(jù)集等。CLASS 語句指明分類變量，協(xié)方差分析時(shí)必須設(shè)立，且必須出現(xiàn)在MODEL語句之前。MODEL 語句定義協(xié)方差分析的線性數(shù)學(xué)模型。例如：MODEL y=a t；選項(xiàng)SOLUTION給出參數(shù)的估計(jì)值。MEANS 語句用于計(jì)算依變量的平均數(shù)。選項(xiàng)用于多重比較。LSMEANS 語句計(jì)算效應(yīng)的最小二乘估