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

5801

收藏 2018-08-28

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

结果：率差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.

扫码加我拉你入群

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

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

被偷走的地下铁

2018-9-12 09:37:21

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）

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群