【002】SAS 中四格表资料，想求2个group率的差值以及95%怎么处理？

凤鸾

本帖所要讲的如题。
给个例子：
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

膜拜楼主