【原创】临床试验计算样本量的SAS程序---以文献为例

如题，在研究《样本量估计及其在nQuery和SAS软件上的实现4——率的比较(三)》这篇文章时，对例子2-9进行的SAS编程实践，最终得出和文献一样的结果，自己写个程序不容易，以前有人发过相关SAS程序，但是根本无法运行，还是自己动手比较好，现在也把正确的SAS程序文档上传附件，供大家使用~~~

SAS9.2软件程序如下（都是我手动输入的，只需替换不同的a,s,p1,p2,power值即可）：
%macro PTT0(a,s,p1,p2,power);
data PTT0;
p1=&p1;
p2=&p2;
s=&s;
a=&a;
power=&power;
q1=1-p1;
q2=1-p2;
p=1/2*(p1+p2);
q=1-p;
za=probit(1-a/s);
zb=probit(power);
n=(za*sqrt(2*p*q)+zb*sqrt(p1*q1+p2*q2))**2/(p1-p2)**2;
n=ceil(n);
run;
proc printdata=PTT0 label;
var a s p1 p2power n;
label a='Test significance level'
s='1 or 2 sided test'
p1='proportion1'
p2='proportion2'
power='power'
n='n'n;quit;
%mend PTT0;
%PTT0(a=0.05,s=2,p1=0.45,p2=0.65,power=0.95);
输出的结果是：