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计算程序