应用SAS软件处理1∶R配比χ2检验
作者：龚洁
单位：武汉市卫生防疫站学校卫生科(430022)
关键词：
中国卫生统计990427 1∶R配比的χ2检验的计算一直是一个比较复杂的问题，在目前几种常用的统计软件包中没有现存的计算机程序。本文介绍应用分层分析的χ2检验的SAS程序处理1∶R配比的χ2检验。
原理与方法
对于分层分析，其第g(g=1,2,…，m)层的资料表格形式为〔1〕：
表1 第g层的资料表格形式
危险因素
, 百拇医药
合计
有
无
发病人数
ag
bg
m1g
未发病人数
cg
dg
m0g
合 计
, 百拇医药
n1g
n0g
ng
其Mantel-Haenszel χ2为：
(1)
统计量χ2MHS服从自由度为1的χ2分布。而对于1∶R配比的χ2检验来说，其第g(g=1,2,…，m)配对组的资料表格形式为〔1，2，3〕：
表2 第g配比组的资料表格式
, http://www.100md.com
接触危险因素
不接触危险因素
合计
病例
yg
1-yg
1
对照
xg
R-xg
R
合计
, http://www.100md.com
xg+yg
1+R-(xg+yg)
1+R
其中，yg是一个指示变量，yg为1时，表示病例接触危险因素，yg为0时表示病例不接触危险因素。将表2的数据代入公式(1)，有：
此公式与1∶R配比的χ2检验的计算公式完全一致。
举例及其SAS计算程序
欲研究儿童意外损伤与某行为方式之间的关系，选择与病例同性别、同年龄的正常儿童作为对照，进行1∶4配比的病例对照研究。共选择了20对，资料见表3。
, 百拇医药
表3 儿童意外损伤与某行为方式关系资料 匹配号
病例
对照
匹配号
病例
对照
1
2
3
4
1
2
3
, http://www.100md.com
4
1
1
1
0
0
0
11
1
0
0
1
0
2
, http://www.100md.com
1
0
1
0
0
12
0
0
0
0
0
3
1
0
, http://www.100md.com
0
0
0
13
1
1
0
1
1
4
0
0
0
0
, 百拇医药
0
14
1
0
0
0
0
5
0
1
0
0
0
15
, http://www.100md.com
1
0
0
0
0
6
1
0
0
0
0
16
1
0
, 百拇医药
1
0
0
7
1
0
1
0
0
17
0
0
0
0
, 百拇医药
0
8
0
0
0
0
0
18
1
1
1
0
1
9
, http://www.100md.com
1
1
1
0
1
19
1
1
0
0
1
10
1
0
, 百拇医药
0
0
0
20
1
0
1
0
0
*：0表示无某种行为方式，1表示有某种行为方式 1.程序说明
g——配比组号
y——为1表示病例，为0表示对照
x——为1表示接触危险因素，为0表示不接触危险因素
, 百拇医药
2.程序
data pd;
input g y x@@;
cards;
1 1 1 1 0 1 1 0 0 1 0 0 1 0 0
2 1 1 2 0 0 2 0 1 2 0 0 2 0 0……；
proc freq;
tables g*y*x/all noprint;
run;
3.结果输出：
利用此程序，算得Mantel Haeszel χ2值为15.362，OR值为72.75。值得一提的是，此方法只在二分类资料(即某危险因素的有无)中验证过，是否适用于其他资料，还有待进一步证实。
参考文献
1.余松林.医学现场研究中的统计分析方法(修订本).武汉：同济医科大学，1985，171～172，194～196.
2.方积乾，等.医学统计学与电脑实验.上海：上海科学技术出版社，1997，218～212.
3.Schlesselman J J.Case-control Studies:Design Conduct,Analysis.New York:Oxford University Press,1982:213～220., 百拇医药
单位：武汉市卫生防疫站学校卫生科(430022)
关键词：
中国卫生统计990427 1∶R配比的χ2检验的计算一直是一个比较复杂的问题，在目前几种常用的统计软件包中没有现存的计算机程序。本文介绍应用分层分析的χ2检验的SAS程序处理1∶R配比的χ2检验。
原理与方法
对于分层分析，其第g(g=1,2,…，m)层的资料表格形式为〔1〕：
表1 第g层的资料表格形式
危险因素
, 百拇医药
合计
有
无
发病人数
ag
bg
m1g
未发病人数
cg
dg
m0g
合 计
, 百拇医药
n1g
n0g
ng
其Mantel-Haenszel χ2为：
(1)
统计量χ2MHS服从自由度为1的χ2分布。而对于1∶R配比的χ2检验来说，其第g(g=1,2,…，m)配对组的资料表格形式为〔1，2，3〕：
表2 第g配比组的资料表格式
, http://www.100md.com
接触危险因素
不接触危险因素
合计
病例
yg
1-yg
1
对照
xg
R-xg
R
合计
, http://www.100md.com
xg+yg
1+R-(xg+yg)
1+R
其中，yg是一个指示变量，yg为1时，表示病例接触危险因素，yg为0时表示病例不接触危险因素。将表2的数据代入公式(1)，有：
此公式与1∶R配比的χ2检验的计算公式完全一致。
举例及其SAS计算程序
欲研究儿童意外损伤与某行为方式之间的关系，选择与病例同性别、同年龄的正常儿童作为对照，进行1∶4配比的病例对照研究。共选择了20对，资料见表3。
, 百拇医药
表3 儿童意外损伤与某行为方式关系资料 匹配号
病例
对照
匹配号
病例
对照
1
2
3
4
1
2
3
, http://www.100md.com
4
1
1
1
0
0
0
11
1
0
0
1
0
2
, http://www.100md.com
1
0
1
0
0
12
0
0
0
0
0
3
1
0
, http://www.100md.com
0
0
0
13
1
1
0
1
1
4
0
0
0
0
, 百拇医药
0
14
1
0
0
0
0
5
0
1
0
0
0
15
, http://www.100md.com
1
0
0
0
0
6
1
0
0
0
0
16
1
0
, 百拇医药
1
0
0
7
1
0
1
0
0
17
0
0
0
0
, 百拇医药
0
8
0
0
0
0
0
18
1
1
1
0
1
9
, http://www.100md.com
1
1
1
0
1
19
1
1
0
0
1
10
1
0
, 百拇医药
0
0
0
20
1
0
1
0
0
*：0表示无某种行为方式，1表示有某种行为方式 1.程序说明
g——配比组号
y——为1表示病例，为0表示对照
x——为1表示接触危险因素，为0表示不接触危险因素
, 百拇医药
2.程序
data pd;
input g y x@@;
cards;
1 1 1 1 0 1 1 0 0 1 0 0 1 0 0
2 1 1 2 0 0 2 0 1 2 0 0 2 0 0……；
proc freq;
tables g*y*x/all noprint;
run;
3.结果输出：
利用此程序，算得Mantel Haeszel χ2值为15.362，OR值为72.75。值得一提的是，此方法只在二分类资料(即某危险因素的有无)中验证过，是否适用于其他资料，还有待进一步证实。
参考文献
1.余松林.医学现场研究中的统计分析方法(修订本).武汉：同济医科大学，1985，171～172，194～196.
2.方积乾，等.医学统计学与电脑实验.上海：上海科学技术出版社，1997，218～212.
3.Schlesselman J J.Case-control Studies:Design Conduct,Analysis.New York:Oxford University Press,1982:213～220., 百拇医药