sas统计学习笔记（四）

有时候列联表的情况很负责，比如某个因素的每个水平上都有一个2×2列联表，一种用于分析多层2×2表的方法是Mantel-Haenszel卡方统计量，这个技术有时也叫元分析。用sas软件计算分层表的卡方值非常简单，只需要在PROC FREQ中输入三向表，并加上选项ALL。

下面有一个例子，假设我们有两个2×2表，一个针对男性，一个针对女性，每个表给出了睡眠时间和体能测试达标之间的关系。

  男性检验结果   女性检验结果   失败通过  失败通过睡眠低20100睡眠低30100 高15150 高25200

***Program to compute a Mantel-Haenszel Chi-squre Test;DATA ABILITY;  INPUT GENDER $ RESULTS $ SLEEP $ COUNT;DATALINES;BOYS FALL 1-LOW 20BOYS FALL 2-HIGH 15BOYS PASS 1-LOW 100BOYS PASS 2-HIGH 150GIRLS FALL 1-LOW 30GIRLS FALL 2-HIGH 25GIRLS PASS 1-LOW 100GIRLS PASS 2-HIGH 200;PROC FREQ DATA=ABILITY;  TITLE "Mantel-Haenszel Chi-squre Test";  TABLES GENDER*SLEEP*RESULTS/ALL;  *第一个为分情况讨论的变量，gender在第一个则分为男女两个表格输出sleep*results；  WEIGHT COUNT;RUN;

（部分结果输出图）

可以看到，男性卡方值为3.7013（p=0.0544），相对风险为1.8333，女性卡方值为9.0106（p=0.0027），相对风险为2.0769.

两个表的综合结果如下：

检验优势比的同质性Breslow-Day统计量不显著（p=0.698），表示我们可以综合这两个表。