统计211

标题: 【002】SAS 中四格表资料，想求2个group率的差值以及95%怎么处理？ [打印本页]

作者: 凤鸾 时间: 2012-1-6 13:41
标题: 【002】SAS 中四格表资料，想求2个group率的差值以及95%怎么处理？
本帖最后由凤鸾于 2012-1-6 13:44 编辑

本帖所要讲的如题。
给个例子：
Treatment    治疗有效治疗无效 Total    有效 Proportion
--------------    -----------    -------------- ------    -------------------------
Placebo             16          48                64       0.250
Test                   40          20                60       0.667
-------------    ------------ ------------- -------    --------------------------
Total                   56          68                124    0.452
虽然我们可以用x2检验出组别不同他们的治疗效果是不相同的。
但如果我们更感兴趣的是test 与placebo的治疗率差以及它的95%的ci，那我们要怎么操作呢？
也就是我想看下test比placebo它的治疗效果好多少，怎么办？

在SAS中我们可以用riskdiff来解决。下面就是程序的code.
data respire2;
input treat $ outcome $ count @@;
datalines;
test f 40 test u 20
placebo f 16 placebo u 48
;
proc freq order=data; /*注：order=data指的是按照我们data的排序进行分析，这里test排练在placebo之前，那么我们的riskdiff 也就是率差计算的就是test-placebo*/
weight count;
tables treat*outcome / riskdiff measures;/*riskdiff 的结果会出现Column 1 和Column 2的 Risk Estimates，Column1  Risk Estimates即为我们所要研究的治疗为有效的差的相关结果*/
run;

我把riskdiff的SAS OUTPUT贴出来进行讲解:
Statistics for Table of treat by outcome
Column 1 Risk Estimates
                                                            (Asymptotic) 95%          (Exact) 95%
            Risk                ASE                   Confidence Limits          Confidence Limits
-----------------------------------------------------------------------------------------------------------
Row 1    0.6667          0.0609                0.5474       0.7859       0.5331 0.7831
Row 2    0.2500          0.0541                0.1439       0.3561       0.1502 0.3740
Total    0.4516          0.0447                0.3640       0.5392       0.3621 0.5435
Difference 0.4167       0.0814                0.2570       0.5763

Difference is (Row 1 - Row 2)

可以看出 Difference =(Row 1 - Row 2)=0.4167即test 与placebo的治疗率的差值是=0.4167，test比placebo治疗效果好41.67%。率的差值的95%的ci为[0.2570 0.5763]
值得注意的是：当样本量较小时，计算出来的CI要比用计算公式计算出来的CI范围窄了一些。

作者: wangming 时间: 2012-1-29 20:35
膜拜楼主

欢迎光临统计211 (http://www.tj211.com/)