data example10_1;
alpha=0.05;                                /*指定一类错误α值*/
beta=0.2;                                   /*指定二类错误β值*/
rt=2;
rho=0.6;
dif=0.4;
s=0.6; 
nc=ceil((2*(probit(1-alpha/2)+probit(1-beta))**2*(s)**2*(1+(rt-1)*rho))/(rt*(dif)**2));
label alpha="α";
label beta="β";
label rt="重复次数";
label rho="相关系数";
label dif="预期两组差值";
label s="合并标准差";
label nc="每组例数";
proc print label;
var rt rho dif s nc;
run;
