sas统计学习笔记(四)
2017-08-18 17:14
211 查看
有时候列联表的情况很负责,比如某个因素的每个水平上都有一个2×2列联表,一种用于分析多层2×2表的方法是Mantel-Haenszel卡方统计量,这个技术有时也叫元分析。用sas软件计算分层表的卡方值非常简单,只需要在PROC FREQ中输入三向表,并加上选项ALL。
下面有一个例子,假设我们有两个2×2表,一个针对男性,一个针对女性,每个表给出了睡眠时间和体能测试达标之间的关系。
(部分结果输出图)
可以看到,男性卡方值为3.7013(p=0.0544),相对风险为1.8333,女性卡方值为9.0106(p=0.0027),相对风险为2.0769.
两个表的综合结果如下:
检验优势比的同质性Breslow-Day统计量不显著(p=0.698),表示我们可以综合这两个表。
下面有一个例子,假设我们有两个2×2表,一个针对男性,一个针对女性,每个表给出了睡眠时间和体能测试达标之间的关系。
男性检验结果 | 女性检验结果 | ||||||
失败 | 通过 | 失败 | 通过 | ||||
睡眠 | 低 | 20 | 100 | 睡眠 | 低 | 30 | 100 |
高 | 15 | 150 | 高 | 25 | 200 |
***Program to compute a Mantel-Haenszel Chi-squre Test; DATA ABILITY; INPUT GENDER $ RESULTS $ SLEEP $ COUNT; DATALINES; BOYS FALL 1-LOW 20 BOYS FALL 2-HIGH 15 BOYS PASS 1-LOW 100 BOYS PASS 2-HIGH 150 GIRLS FALL 1-LOW 30 GIRLS FALL 2-HIGH 25 GIRLS PASS 1-LOW 100 GIRLS 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),表示我们可以综合这两个表。
相关文章推荐
- SAS统计分析学习笔记(八)——T检验和非参数比较
- 统计学习方法学习笔记
- 机器学习-李航-统计学习方法学习笔记之感知机(1)
- 算法导论 学习笔记 第九章 中值和顺序统计
- SAS学习笔记:读取混合格式的数据
- linux运用shell脚本统计列出目录文件个数(学习笔记)
- 统计学习笔记--感知机模型
- ≪统计学习精要(The Elements of Statistical Learning)≫课堂笔记(十三)
- 统计学习笔记之感知机
- Hadoop学习笔记—20.网站日志分析项目案例(三)统计分析
- ≪统计学习精要(The Elements of Statistical Learning)≫课堂笔记(三)
- BZOJ 1036: [ZJOI2008]树的统计Count [树链剖分]【学习笔记】
- 统计学习方法学习笔记:概论
- C++学习笔记27——泛型算法之统计(accumate count)
- SAS学习笔记之输出到Excel
- [统计学习方法学习笔记]1.1 统计学习
- R语言学习笔记3--描述统计
- 统计学习笔记——k近邻法与kd树
- 统计学习笔记(4) 线性回归(1)
- 统计学习方法学习笔记一