SAS-卡方检验

c573489167

已于 2023-08-17 10:43:20 修改

阅读量1.2w

点赞数 1

分类专栏： SAS 文章标签： sas

于 2021-06-25 14:38:25 首次发布

本文链接：https://blog.csdn.net/c573489167/article/details/118220383

版权

SAS 专栏收录该内容

17 篇文章

订阅专栏

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;