data example10_2;
alpha=0.05;                                /*指定一类错误α值*/
beta=0.1;                                  /*指定二类错误β值*/
rt=4;                                       /*指定重复测量的次数*/
rho=0.4;                                   /*指定多次测量间的相关系数*/
p1=0.6;                                    /*指定第1组的阳性率*/
p2=0.5;                                    /*指定第2组的阳性率*/
q1=1-p1;
q2=1-p2;
p=(p1+p2)/2;
q=(q1+q2)/2;
nc=ceil((probit(1-alpha/2)*sqrt(2*p*q)+probit(1-beta)*sqrt(p1*q1+p2*q2))**2*(1+(rt-1)*rho)/(rt*(p1-p2)**2));
label alpha="α";
label beta="β";
label rt="重复次数";
label rho="相关系数";
label p1="第1组阳性率";
label p2="第2组阳性率";
label nc="每组例数";
proc print label;
var rt rho p1 p2 nc;
run;
