多组样本平均值比较.pdf

中國醫藥大學 生物統計中心                                                  2009/7  第三章  多組樣本平均值比較  邱顯財  統計分析師    在前一章中，我們介紹了如何比較兩組資料的集中趨勢。當資料分佈符合常態時，可 使用獨立樣本t 檢定(Independent t‐test) ，而當資料分佈違反常態時，則應使用魏克森等級 和檢定(Wilcoxon rank sum test) 。常見的情況為研究者有三組或以上的資料需要作比較，如 果單純地進行兩兩比較，會增加整體的型一誤差(type I error) ，這時需要提供一個整體檢 定的統計方法進行多組資料的集中趨勢檢定。  本章將分三個部分進行：(1) 如何使用百分位數分組；(2) 如何使用統計圖來呈現資 料、檢定分佈的形狀(shape) ；(3) 如何進行多組平均值或中位數的統計檢定。本章將以一 個研究問題來討論這三個部份。    1‐0. 使用百分位數分組      我們想知道血鈣濃度對洗腎病人的存活時間有無影響。因為鈣濃度為連續變項，我們 希望知道低血鈣、中血鈣與高血鈣的洗腎病人，其存活時間有無差異。在對連續變項分組 時，若沒有臨床切點，使用百分位數分組是常見的做法，例如三等分法、四等分法、五等 分法。以下將介紹如何使用PROC UNIVARIATE 產生百分位數。    1‐1. 以PROC UNIVARIATE 產生百分位數       由於我們希望將血鈣濃度分成三組，在這裡我們採用三分位數(tertiles)作為分組的依 據，也就是以33 與66 百分位數作為切點，SAS 語法如【程式01‐1】。      /* 以百分位數(percentile)切分資料  */     /* 33, 66 percentiles */    PROC UNIVARIATE DATA=esrd NOPRINT;    VAR ca;    OUTPUT OUT=b PCTLPRE=Ter PCTLPTS=33,66;    【程式01‐1】、PROC UNIVARIATE 語法  圖01‐1 、【程式01-1】 RUN;    PROC PRINT; RUN;  的輸出結果    【程式01‐1】、PROC UNIVARIATE 語法配合               output 敘述輸出百分位數於資料檔b    在圖01‐1 中我們看到血鈣濃度33 與66 百分位數分別是8.2 與8.9 。        3‐1  SAS教戰手冊 3‐1    中國醫藥大學 生物統計中心                                                  2009/7  第1‐1 節重要指令說明：  1. OUTPUT : 在 UNIVARIATE 程序中，這個敘述可以將統計分析結果輸出於資料檔中。  2. OUT=資料檔名稱：將指定統計值輸出到特定檔案。這裡指定百分位數。  3. PCTLPRE=變項名稱：指定百分位數變項的字首。這個範例裡使用ter 這個字首。  4. PCTLPTS=數值1,數值2,…   (1 至100 皆可) ：欲觀察之百分位數。例如希望知道20 百分 位數為何，就在此輸入20  。    1‐2. 使用「IF…THEN / ELSE… 」敘述將資料分組      我們已經得知血鈣濃度的33 與66 百分位數分別是8.2 與8.9 ，接下來我們將進行分 組的動作。這裡我們在資料步驟(data step)中，執行「IF…THEN / ELSE… 」敘述，產生將資 料分為三組的新變項(ca_1) ：血鈣濃度(ca)8.2 、8.2