SAS-卡方检验

1、独立样本卡方检验：

独立样本卡方检验
expected：输出理论频率（期望频率）
exact：计算精确概率

数据集如下：

检验结果如下：

2、配对样本卡方检验：

大样本配对卡方检验

数据集如下：

检验结果如下：

3、配对卡方检验（小样本）：

配对四个表卡方检验（小样本）
SAS中McNemar检验不考虑校正问题，当b+c<40时，需校正。需要自己写公式计算。

数据集如下：

代码

data chitest1;
do r = 1 to 2;
do c = 1 to 2;
input f@@;
output;
end;
end;
cards;
64 21
51 33
;
run;

proc freq data=chitest1;
weight f;
tables r*c/chisq expected exact ; /** chisq：卡方检验 **/
run;

data chitest2;
input r c f;
cards;
1 1 80
1 2 10
2 1 31
2 2 11
;
run;

proc freq data=chitest2;
weight f;
tables r*c/agree; /** agree：配对卡方检验（McNemar） **/
run;

data chitest3;
input f11 f12 f21 f22;
if f12+f21<40 then chisq = (abs(f12-f21)-1)**2/(f12+f21) ; /**b+c<40,需校正 **/
else chisq = abs(f12-f21)**2/(f12+f21);
p=1-probchi(chisq,1);
cards;
80 10 31 11
;
run;