1:M配對病例-對照研究中兩種可信區間的SAS分析方法

衛生部衛生經濟研究所衛生政策和發展研究室(100191) 王書平

1:M配比的病例對照研究方法資料單因素常用Mantel-Haenzel卡方檢驗〔1-3〕進行分析，但是目前通用的統計軟件大多沒有提供計算程序，而手工計算又十分復雜，本文就如何應用SAS軟件對1:M配比資料的Mantel-Haenzel卡方檢驗的計算方法及對不同的方法計算OR可信區間進行討論和比較。

原理與方法

當1:M配對的病例-對照研究僅存在一個配比因素時可以采用分層分析，即 Mantel-Haenzel卡方檢驗，1:M配對的病例-對照研究資料可按照表1的格式整理，并按照計算公式計算χ2M-H和 ORMH。

表1 1:M病例對照研究資料整理表

ORMH的可信區間方法一的公式〔4〕:首先計算ln(OR)的標準誤SlnOR，然后可信區間。

舉例及SAS的計算程序

1.使用Mantel-Haenzel卡方檢驗來計算χ2M-H，方法一的公式計算ORMH可信區間。

(1)程序說明及程序見附錄。

m，整理表格中的行數;n，整理表格中的列數。

表2 人工流產史與宮外孕之間關系的調查結果〔6〕

(1)程序的說明。

計算χ2和校正χ2及OR與方法一程序相同，只有計算OR的可信區間的程序有所不同:程序見附錄中按照方法二的公式計算所示。

(2)結果的輸出如下:Chisq、CChisq、OR值、Conf95l、Conf95u的含義與方法一解釋相同，Chisq、CChisq、OR值三者結果與方法一的結果相同，Conf95l=5.95和Conf95u=183.06是按照可信區間的公式二計算得出，與方法一的結果不同。

討 論

本文主要論述了配對的病例-對照研究的分析方法，并且對于每種情況的原理和SAS的程序進行論述，同時加以實例來說明問題。就1:1配對的情況來講，使用SPSS中的MeNemar只能給出具體的P值而不能給出卡方值，并且其給出的OR值及OR的可信區間均是按照成組的病例-對照的公式給出的，如果使用SAS編程中的Table后Agree All可以計算MeNemar的卡方值，但其均是未經過校正的卡方值，所給出的OR及可信區間也是按照成組計算。本文通過1:M配對病例-對照的計算的原始公式進行編程，按照不同的情況分別的計算卡方值和OR及可信區間，只需要對程序中的宏定義中的n改為2，按照表1錄入數據即可以計算出相應的卡方值、校正卡方值及OR和OR的可信區間。

對于1:M的配對的病例對照研究Mantel-Haenszel卡方檢驗計算OR的可信區間，本文分別按照計算OR可信區間的方法兩個公式分別編程計算，如果使用方法一的公式得出的OR的95%的可信區間是2.91～373.95，而使用方法二的公式得出的可信區間是5.95～183.06。通過比較，可以得出方法二求出的可信區間的范圍比較窄，因此其計算可信區間的精確度比較高。從方法和原理可以看出方法一計算OR的可信區間主要使用表格中的原始數據，而方法二計算OR的可信區間直接使用未經過校正的卡方值，和卡方值的聯系比較大。

1.余松林.醫學現場研究中的統計分析方法(修訂本).武漢:同濟醫科大學，1985:171-172，194-196.

2.Schlesselman.Case-control studies:design conduct and analysis.New York Oxford University Press，1982:213-220.

3.王旭善，吳德林，張曉峰，等.贛榆縣食管癌流行危險因素1:1配對病例對照研究.中國衛生統計，2008，12(25):616-617.

4.金丕煥.醫用統計方法.上海醫科大學出版社，2000，9:351-354.

5.李立明.流行病學.人民衛生出版社，2004，9:99-100.

6.時景璞，董衛.臨床醫學研究方法——設計、測量與評價.中國醫科大學教材出版社，2006，7:70-71.

附錄 Mantel-Haenzel卡方檢驗和OR及可信區間的SAS計算程序