sas9.2求临床试验中有效率率差（试验组）的95%可信区间

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

	试验组	对照组
有效	168（98.82%）	55（96.49%）
无效	2（1.18%）	2（3.51%）
合计	170	57

结果：率差2.33（-2.71,7.38）在网上找了《两组率同为100%或0%时率差置信区间》，就试了试，但是结果算出来的和上面不同。求指导，大恩不言谢，拜托了。

1. data _null_;
   
2.    call symput("z",-probit(0.025));
   
3. run;
   
4. 
   
5. data unpaired;
   
6. a1=168;a2=2;a3=55;a4=2;
   
7. n1=a1+a2; p1=a1/n1;
   
8. n2=a3+a4; p2=a3/n2;
   
9. l1=(2*a1+&z.**2-&z.*sqrt(&z.**2+4*a1*a2/n1))/(2*(n1+&z.**2));*单样本下限;
   
10. u1=(2*a1+&z.**2+&z.*sqrt(&z.**2+4*a1*a2/n1))/(2*(n1+&z.**2));*单样本上限;
    
11. l2=(2*a3+&z.**2-&z.*sqrt(&z.**2+4*a3*a4/n2))/(2*(n2+&z.**2));*单样本下限;
    
12. u2=(2*a3+&z.**2+&z.*sqrt(&z.**2+4*a3*a4/n2))/(2*(n2+&z.**2));*单样本上限;
    
13. d=p1-p2;
    
14. l_unpaired=d-sqrt((p1-l1)**2+(u2-p2)**2);*率差下限;
    
15. u_unpaired=d+sqrt((p2-l2)**2+(u1-p1)**2);*率差上限;
    
16. proc print;
    
17. run;

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：置信区间 单样本 有效率 对照组 试验组 sas

这是stata版
不是sas版

你让版主给你移到sas版询问吧

蓝色 发表于 2018-8-28 14:39
这是stata版本
不是sas版本

谢谢，我刚发现了，已经挪过去了