第七章方差分析和協(xié)方差分析（精品）

《第七章方差分析和協(xié)方差分析（精品）》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

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

2、ORWIN957.5來(lái)做，運(yùn)行了大約10分鐘才出結(jié)果。我又換用SASFORWIN956.12來(lái)做結(jié)果用了一一2.47秒!§7.1方差分析數(shù)據(jù)集的建立技巧7.1.1方差分析的數(shù)據(jù)集格式統(tǒng)計(jì)分析所用的數(shù)據(jù)格式和我們?cè)诜治稣碣Y料時(shí)所用的格式是不同的。般來(lái)說(shuō)，數(shù)據(jù)集中應(yīng)至少有一個(gè)結(jié)果變量，用于記錄不同處理因素水平下觀察值的大小；至少有一個(gè)處理因素變量，用于記錄處理因素的類型及其水平數(shù)。以單因素方差分析為例，就應(yīng)有一個(gè)結(jié)果變量和一個(gè)處理因索變量；而兩因索的方差分析應(yīng)有一個(gè)結(jié)果變量和兩個(gè)處理因素變量。例7.1某職業(yè)病防治院對(duì)31名石棉

3、礦工中的石棉肺患者、可疑患者及非患者進(jìn)行了用力肺活量測(cè)定，請(qǐng)給出數(shù)據(jù)集的結(jié)構(gòu)（衛(wèi)統(tǒng)P44例5?1）。解：數(shù)據(jù)集中應(yīng)有兩個(gè)變量，x和groupox記錄肺活量的大小；group取值為1、2或3,分別代表石棉肺患者、可疑患者及非患者。例7.2某廠醫(yī)務(wù)室測(cè)定了10名氟作業(yè)工人工前、工中及工后4小時(shí)的尿氟濃度,請(qǐng)給出數(shù)據(jù)集的結(jié)構(gòu)（衛(wèi)統(tǒng)p46例5.2）。解：數(shù)據(jù)集中應(yīng)有三個(gè)變量，x、group和*70廠1<8廠。x記錄尿氟濃度；group取值為1、2或3,分別代表工前、工中及工后；worker取值為1到10,分別代表10名人。7.1.

4、2方差分析數(shù)據(jù)集的建立技巧可見方弄分析的數(shù)據(jù)集其變量取值有一定的規(guī)律，因此可以利用循環(huán)語(yǔ)句和判斷語(yǔ)句來(lái)簡(jiǎn)化輸入。例7.3請(qǐng)建立例6.1的數(shù)據(jù)集。libnamea'c:user';dataa.wtli5_l;group=l;if_n_>llthenif_n_>20thengroup=3;elsegroup=2;解：我們可以一個(gè)一個(gè)的輸入變量group的值，但這里給出一種更酷的建立方法。此處的_n_為SAS系統(tǒng)變量，表明讀入的數(shù)據(jù)為第幾條記錄，由于石棉肺患者只有X人，因此第12條記錄只能是可疑患者的記錄，同理，從第20條記錄開

5、始的只能是非患者的記錄。這里利用判斷語(yǔ)句來(lái)完成從_n_到group的轉(zhuǎn)換。inputvaluecards;1.81.41.52.13.33.5procprint;run;例7.4請(qǐng)建立例7.2的數(shù)據(jù)集。解：此例中數(shù)據(jù)較有規(guī)律，各組的例數(shù)均相等，這可正是循環(huán)語(yǔ)句大顯身手的時(shí)候。libnameaBc:user1;dataa.wtli5_2;dogroup=lto3;這里的SAS語(yǔ)句非常類似于Foxbase語(yǔ)句.兩個(gè)循環(huán)分別控doworker=lto10;制了group和worker的取值，這在輸入數(shù)據(jù)時(shí)是常見的一種技巧0inp

6、utvalueoutput;兩個(gè)變量只是循環(huán)變量，因此要用output語(yǔ)句寫入數(shù)據(jù)集中.end;end;cards;90.5388.4347.37???請(qǐng)注意X變童和處理因素變量的對(duì)應(yīng)關(guān)系，不要改變數(shù)據(jù)塊的次序，???105.2758.95■否則建立的數(shù)據(jù)集是錯(cuò)的.procprint;run;§7.2引例例7?5請(qǐng)回答在例7.1中三組石棉礦工的用力肺活量有無(wú)差別？解：由于是三組做比較，因此要用方差分析法。在建立數(shù)據(jù)集后可用ANOVAi±程或GLM過程進(jìn)行分析，同時(shí)給出選項(xiàng)進(jìn)行各組均數(shù)的兩兩比較（此處用SNK法），兩過程的輸出

7、格式基本一致。對(duì)于單因素的方差分析，在比較之前一般要考慮方差齊性的問題。這里可以用GLM過程中的HOVTEST選項(xiàng)實(shí)現(xiàn)，實(shí)際上這也是在SAS中做方差齊性檢驗(yàn)的唯一途徑。①設(shè)定數(shù)據(jù)庫(kù)環(huán)境:LIBNAMEA‘C:USEFT;②數(shù)據(jù)步，建立數(shù)據(jù)集（略）：③ANOVA或GUV1過程，進(jìn)行方差分析并同時(shí)做兩兩比較（SNK法）：PROCANOVADATA=A.WTLI5_1;PROCGLMDATA=A.WTLI5_1;CLASSGROUP;CLASSGROUP;MODELVALUE=GROUP;MODELVALUE=GROUP/SS1

8、;MEANGROUP/SNK;MEANGROUP/SNKHOVTEST：RUN;RUN;3這里的GLM過程中有兩個(gè)選項(xiàng)，SS1要求只輸出方差分析SS1模型的結(jié)果（默認(rèn)輸出SS1和SS3兩種模型的結(jié)果），HOVTEST則要求進(jìn)行方差齊性檢驗(yàn)（該選項(xiàng)只能在SAS6.12以上的版本中實(shí)現(xiàn)）。§7