簡易流行病學與生物統計學 二十一、作卡方檢定

簡易流行病學與生物統計學

二十一、作卡方檢定

作者：吳聰賢醫師

SAS（Statistical analysis system）作卡方檢定

SAS（Statistical analysis system）起源於1976年，為一套常被使用的統計分析系統，目前尚未有中文版，在功能表列的Window有三種可以互相切換的視窗Program editor、Log、Output：

   Program editor（編輯視窗）：提供程式的編修、及程式反白框起來執行

   Log（運算紀錄視窗）：顯示程式執行過程、及程式語法錯誤訊息；會因多次程式的執行，累積很多運算紀錄，可用開新檔案的方式消除

   Output（運算結果視窗）：顯示程式執行後的結果；也會因多次程式的執行，累積很多運算紀錄，可用開新檔案的方式消除

以上述喝可樂與中毒的例子書寫運算程式：

               Program editor（編輯視窗）

data intoxi; input cola poison count; cards; 0 0 107 0 1 53 1 0 2 1 1 9 ; proc freq; weight count; tables cola*poison/chisq; run;

  data是資料檔案，intoxi代表檔案檔名，是intoxication（中毒）的縮寫，因為檔名不可以超過8個字元；input代表輸入變項，cola、和poison均是變項名稱，變項命名有三個禁忌：

1.  不可超過8個字元，如intoxication需改成intoxi

2.  不可以用阿拉伯數字開頭，如5yrdeath需改成death5yr

3.  不可以用【-】此符號，需改成其他符號

 count是變項累積的數目。

【0， 0， 107】第一個數字代表喝可樂，0是有喝可樂，1是沒喝可樂；第二個數字代表中毒，0是有中毒，1是沒中毒；後面數字是累積人數，運算結果如下：

      Output（運算結果視窗）

Statistics for table of cola by poison Statistic                    DF     Value         Prob Chi-square                    1      10.559        0.001 Likelihood ratio Chi-square   1      10.320        0.001 Continuity adj. Chi-square    1       8.557        0.003 Mantel-Haenszel Chu-square    1      10.497        0.001 Fisher’s exact test（left）                        1.000               （right）                          2.00E-03               （2-tail）                         2.00E-03

      DF是degree of freedom（自由度）

      Value是x² value（卡方值）

      Prob是probability value（P值）

      2.00E-03代表2.00 ×10^-3            

  上述喝可樂與中毒的例子，卡方檢定結果卡方值10.559，P值0.001達統計上意義。注意！著不到勝算比（odds ratio）、和勝算比95%信賴區間，為什麼？