1:M配对病例-对照研究中两种可信区间的SAS分析方法
2011-07-27卫生部卫生经济研究所卫生政策和发展研究室100191王书平
卫生部卫生经济研究所卫生政策和发展研究室(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值而不能给出卡方值,并且其给出……
