統(tǒng)計(jì)分析系統(tǒng)SAS軟件 第五章 方差分析和協(xié)方差分析

1、第五章 方差分析和協(xié)方差分析ANOVA過(guò)程 -速度快，功能有限,平衡資料GLM過(guò)程 -速度慢，但功能強(qiáng)大（非平衡資料 協(xié)方差分析）其實(shí)，這里的速度快慢只是相對(duì)而言，SAS的處理速度是首屈一指的。舉個(gè)例子，一個(gè)有6600條記錄的，7因素的，交叉設(shè)計(jì)的方差分析（是不是已經(jīng)有人喊頭痛了？），我先是用SPSS FOR WIN95 7.5來(lái)做，運(yùn)行了大約10分鐘才出結(jié)果。我又換用SAS FOR WIN95 6.12來(lái)做，結(jié)果用了2.47秒！5.1 方差分析數(shù)據(jù)集的建立技巧5.1.1 方差分析的數(shù)據(jù)集格式 統(tǒng)計(jì)分析所用的數(shù)據(jù)格式和我們?cè)诜治稣碣Y料時(shí)所用的格式是不同的: -數(shù)據(jù)集中應(yīng)至少有一個(gè)結(jié)果變量，用

2、于記錄不同處理因素水平下觀察值的大?。?-至少有一個(gè)處理因素變量，用于記錄處理因素的類(lèi)型及其水平數(shù)。例5.1 某職業(yè)病防治院對(duì)31名石棉礦工中的石棉肺患者、可疑患者及非患者進(jìn)行了用力肺活量測(cè)定，請(qǐng)給出數(shù)據(jù)集的結(jié)構(gòu)。解：數(shù)據(jù)集中應(yīng)有兩個(gè)變量，x和group。x記錄肺活量的大??；group取值為1、2或3，分別代表石棉肺患者、可疑患者及非患者。例5.2 某廠醫(yī)務(wù)室測(cè)定了10名氟作業(yè)工人工前、工中及工后4小時(shí)的尿氟濃度，請(qǐng)給出數(shù)據(jù)集的結(jié)構(gòu)。解：數(shù)據(jù)集中應(yīng)有三個(gè)變量，x、group和worker。x記錄尿氟濃度；group取值為1、2或3，分別代表工前、工中及工后；worker取值為1到10，分別代表

3、10名工人。5.1.2 方差分析數(shù)據(jù)集的建立技巧可見(jiàn)方差分析的數(shù)據(jù)集其變量取值有一定的規(guī)律，因此可以利用循環(huán)語(yǔ)句和判斷語(yǔ)句來(lái)簡(jiǎn)化輸入。例5.3 請(qǐng)建立例5.1的數(shù)據(jù)集。例5.4 請(qǐng)建立例5.2的數(shù)據(jù)集。5.2 引 例例5.5 請(qǐng)回答在例5.1中三組石棉礦工的用力肺活量有無(wú)差別？解：對(duì)于單因素的方差分析，在比較之前一般要考慮方差齊性的問(wèn)題。這里可以用GLM過(guò)程中的HOVTEST選項(xiàng)實(shí)現(xiàn)，實(shí)際上這也是在SAS中做方差齊性檢驗(yàn)的唯一途徑。5.3 ANOVA過(guò)程5.3.2 語(yǔ)法說(shuō)明Means語(yǔ)句的選項(xiàng)主要用來(lái)指定兩兩比較的方法和檢驗(yàn)水準(zhǔn)，主要有：兩兩比較方法 BON（Bonferroni檢驗(yàn)）、DUN

4、CAN（新復(fù)極差法）、DUNNETT、GABRIEL、LSD（最小顯著差法 T）、SCHEFFE、SIDAK 、SNK（q檢驗(yàn)）、TUKEY 、WALLER 其它 ALPHA= p值 更改水準(zhǔn)5.3.3 結(jié)果解釋方差分析的統(tǒng)計(jì)結(jié)果部分 5.4 GLM過(guò)程GLM過(guò)程:General Liner Model過(guò)程，它使用最小二乘法對(duì)數(shù)據(jù)擬合廣義線形模型。GLM過(guò)程中可以進(jìn)行回歸分析、方差分析、協(xié)方差分析、劑量反應(yīng)模型分析、多元方差分析和偏相關(guān)分析等等，其功能之強(qiáng)大可見(jiàn)一斑。5.4.1 語(yǔ)法格式這里我們只解釋協(xié)方差分析的做法：由于協(xié)變量并非我們研究的處理因素，因此在CLASS語(yǔ)句中不能有它的位置；但是

5、，另一方面協(xié)變量要影響結(jié)果變量的取值，因此必須在模型中引入，怎么辦？只要在MODEL語(yǔ)句中將其寫(xiě)入即可（寫(xiě)在分析變量的首位）。5.4.2 結(jié)果解釋GLM過(guò)程的結(jié)果和ANOVA過(guò)程是雷同的注意是雷同而不是相同。在處理因素、協(xié)變量及交互作用的檢驗(yàn)結(jié)果部分，GLM過(guò)程會(huì)給出兩種結(jié)果（類(lèi)和類(lèi)結(jié)果）。在我們用到的范圍里，兩種結(jié)果是完全相同的，所以用哪一種都可以。5.4.3 應(yīng)用實(shí)例方差分析基礎(chǔ)常用試驗(yàn)設(shè)計(jì)的ANOVA缺省數(shù)據(jù)的方差分析方差分析基礎(chǔ)單向分組資料的ANOVE 1）樣本容量相等（model y=trt) 2）樣本容量不相等(model y=trt) 3）數(shù)據(jù)轉(zhuǎn)換（model yp=trt) 單一自由度的比較(model y=trt)兩向分組資料的ANOVA(model y=m h)系統(tǒng)分組資料的ANOVA(model y=trt pot(trt)常用試驗(yàn)設(shè)計(jì)的ANOVA完全隨機(jī)化設(shè)計(jì)（單向分組資料）隨機(jī)區(qū)組設(shè)計(jì) 1）單因素（model y=blk trt) 2）兩因素(model y=blk a b a*b) 3）三因素(model y=blk a b a*b c a*c b*c a*b*c) (model y=blk a|b|c)拉丁方設(shè)計(jì)(mode