试验组 | 对照组 | |
有效 | 168(98.82%) | 55(96.49%) |
无效 | 2(1.18%) | 2(3.51%) |
合计 | 170 | 57 |
- 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]复制代码