请选择 进入手机版 | 继续访问电脑版
4669 4

[程序分享] sas9.2临床试验中有效率率差的95%可信区间 [推广有奖]

  • 11关注
  • 0粉丝

本科生

27%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
4.2000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
766 点
帖子
16
精华
0
在线时间
161 小时
注册时间
2018-8-17
最后登录
2024-3-28

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币


试验组对照组
有效168(98.82%)55(96.49%)
无效2(1.18%)2(3.51%)
合计17057
结果:率差2.33(-2.71,7.38)在网上找了《两组率同为100%或0%时率差置信区间》,就试了试,但是结果算出来的和上面不同。求指导,大恩不言谢,拜托了。
  • data _null_;
  •    call symput("z",-probit(0.025));
  • run;
  • data unpaired;
  • a1=168;a2=2;a3=55;a4=2;
  • n1=a1+a2; p1=a1/n1;
  • n2=a3+a4; p2=a3/n2;
  • l1=(2*a1+&z.**2-&z.*sqrt(&z.**2+4*a1*a2/n1))/(2*(n1+&z.**2));*单样本下限;
  • u1=(2*a1+&z.**2+&z.*sqrt(&z.**2+4*a1*a2/n1))/(2*(n1+&z.**2));*单样本上限;
  • l2=(2*a3+&z.**2-&z.*sqrt(&z.**2+4*a3*a4/n2))/(2*(n2+&z.**2));*单样本下限;
  • u2=(2*a3+&z.**2+&z.*sqrt(&z.**2+4*a3*a4/n2))/(2*(n2+&z.**2));*单样本上限;
  • d=p1-p2;
  • l_unpaired=d-sqrt((p1-l1)**2+(u2-p2)**2);*率差下限;
  • u_unpaired=d+sqrt((p2-l2)**2+(u1-p1)**2);*率差上限;
  • proc print;
  • run;

[color=rgb(51, 102, 153) !important]复制代码


二维码

扫码加我 拉你入群

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

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

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

suzhzh 发表于 2018-9-12 09:31:22 |显示全部楼层 |坛友微信交流群
I guess this is is because when the percentage is 0 or 100% the method here is not appropriate.

使用道具

suzhzh 发表于 2018-9-12 09:31
I guess this is is because when the percentage is 0 or 100% the method here is not appropriate.
thanks for your answer.

使用道具

myzhang1982 在职认证  发表于 2023-7-18 16:32:47 |显示全部楼层 |坛友微信交流群
上述SAScode的算法是基于newcombe方法,结果是-0.0161 0.1079 ,-2.71,7.38是基于未经过连续性校正的Wald公式计算

使用道具

hexh 发表于 2023-7-23 15:41:47 |显示全部楼层 |坛友微信交流群
data propd;
   input trt $ res $ f;
datalines;
TR Y 168  
TR N  2
CK Y 55
CK N  2
;
proc freq order=data;
   tables trt*res / riskdiff(column=1 equal var=sample);
   weight f;
run;
quit;

trt * res”表的统计量

列 1 风险估计值
  风险 渐近标准误差 (渐近的) 95%
置信限 (精确)95%
置信限
第 1 行 0.9882 0.0083 0.9720 1.0000 0.9581 0.9986
第 2 行 0.9649 0.0244 0.9171 1.0000 0.8789 0.9957
合计 0.9824 0.0087 0.9653 0.9995 0.9555 0.9952
差值 0.0233 0.0257 -0.0271 0.0738     
差值为(行 1 - 行 2)

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 03:45