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

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

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

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

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

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