衛(wèi)生部衛(wèi)生經(jīng)濟研究所衛(wèi)生政策和發(fā)展研究室(100191) 王書平

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

原理與方法

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

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

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

舉例及SAS的計算程序

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

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

m，整理表格中的行數(shù);n，整理表格中的列數(shù)。

表2 人工流產(chǎn)史與宮外孕之間關(guān)系的調(diào)查結(jié)果〔6〕

(1)程序的說明。

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

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

討 論

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

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

1.余松林.醫(yī)學現(xiàn)場研究中的統(tǒng)計分析方法(修訂本).武漢:同濟醫(yī)科大學，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配對病例對照研究.中國衛(wèi)生統(tǒng)計，2008，12(25):616-617.

4.金丕煥.醫(yī)用統(tǒng)計方法.上海醫(yī)科大學出版社，2000，9:351-354.

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

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

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