第七章 方差分析和協(xié)方差分析

1、第七章第七章方差分析和協(xié)方差分析方差分析和協(xié)方差分析（醫(yī)學(xué)統(tǒng)計(jì)之星）上次更新日期：方差分析和協(xié)方差分析在SAS系統(tǒng)中由SASSTAT模塊來(lái)完成，其中我們常用的有ANOVA過(guò)程和GLM過(guò)程。前者運(yùn)算速度較快，但功能較為有限；后者運(yùn)算速度較慢，但功能強(qiáng)大，我們做協(xié)方差分析時(shí)就要用到GLM過(guò)程。本章將首先介紹方差分析所用數(shù)據(jù)集的建立技巧，然后重點(diǎn)介紹這兩個(gè)程序步。其實(shí)，這里的速度快慢只是相對(duì)而言，SAS的處理速度是首屈一指的。舉個(gè)例子，這個(gè)暑

2、假我做了一個(gè)有6600條記錄的，7因素的，交叉設(shè)計(jì)的方差分析（是不是已經(jīng)有人喊頭痛了？），我先是用SPSSFWIN957.5來(lái)做，運(yùn)行了大約10分鐘才出結(jié)果。我又換用SASFWIN956.12來(lái)做，結(jié)果用了――2.47秒！7.1方差分析數(shù)據(jù)集的建立技巧7.1.1方差分析的數(shù)據(jù)集格式方差分析的數(shù)據(jù)集格式統(tǒng)計(jì)分析所用的數(shù)據(jù)格式和我們?cè)诜治稣碣Y料時(shí)所用的格式是不同的。一般來(lái)說(shuō)，數(shù)據(jù)集中應(yīng)至少有一個(gè)結(jié)果變量，用于記錄不同處理因素水平下觀察值的

3、大??；至少有一個(gè)處理因素變量，用于記錄處理因素的類(lèi)型及其水平數(shù)。以單因素方差分析為例，就應(yīng)有一個(gè)結(jié)果變量和一個(gè)處理因素變量；而兩因素的方差分析應(yīng)有一個(gè)結(jié)果變量和兩個(gè)處理因素變量。例7.1某職業(yè)病防治院對(duì)31名石棉礦工中的石棉肺患者、可疑患者及非患者進(jìn)行了用力肺活量測(cè)定，請(qǐng)給出數(shù)據(jù)集的結(jié)構(gòu)（衛(wèi)統(tǒng)p44例5.1）。解：數(shù)據(jù)集中應(yīng)有兩個(gè)變量，x和group。x記錄肺活量的大?。籫roup取值為1、2或3，分別代表石棉肺患者、可疑患者及非患者。

4、解：此例中數(shù)據(jù)較有規(guī)律，各組的例數(shù)均相等，這可正是循環(huán)語(yǔ)句大顯身手的時(shí)候。libnameac:userdataa.wtli5_2dogroup=1to3dowker=1to10這里的SAS語(yǔ)句非常類(lèi)似于Foxbase語(yǔ)句，兩個(gè)循環(huán)分別控制了group和wker的取值，這在輸入數(shù)據(jù)時(shí)是常見(jiàn)的一種技巧。inputvalue@@output兩個(gè)變量只是循環(huán)變量，因此要用output語(yǔ)句寫(xiě)入數(shù)據(jù)集中。endendcards90.5388.434

5、7.37......105.2758.95請(qǐng)注意x變量和處理因素變量的對(duì)應(yīng)關(guān)系，不要改變數(shù)據(jù)塊的次序，否則建立的數(shù)據(jù)集是錯(cuò)的。procprintrun7.2引例例7.5請(qǐng)回答在例7.1中三組石棉礦工的用力肺活量有無(wú)差別？解：由于是三組做比較，因此要用方差分析法。在建立數(shù)據(jù)集后可用ANOVA過(guò)程或GLM過(guò)程進(jìn)行分析，同時(shí)給出選項(xiàng)進(jìn)行各組均數(shù)的兩兩比較（此處用SNK法），兩過(guò)程的輸出格式基本一致。對(duì)于單因素的方差分析，在比較之前一般要考慮方