样本量计算SAS程序大全

=(1/&q1+1/&q2)*((abs(probit(&alpha))+probit((1-&beta)))*&std/&delta)**2;Р N=int(n)+1;Р n1=int(N*&q1)+1;Р n2=int(N*&q2)+1;Р proc print data=d1_4;Р run;Р %mend;Р %twosample_mean1(alpha=0.05,beta=0.10,delta=5,std=12,q1=0.333,q2=0.667);Р /*---配对二分类资料分析时的样本含量计算---*/Р %macro parisample_ratio(p1=,p2=,alpha=,beta=,pk=);Р data d1_5;Р %let pc=(&p1+&p2-2*&Pk)/2;Р n=(((abs(probit(&alpha/2)))*(sqrt(2*&pc))+(probit(1-Р &beta))*sqrt(2*(&p1-&pk)*(&p2-&pk)/(&pc)))/(&p1-&p2))**2;Р N=int(n)+1;Р proc print data=d1_5;Р run;Р %mend;Р %parisample_ratio(p1=0.48,p2=0.30,alpha=0.05,beta=0.10,pk=0.25);Р /*---两样本率分析时样本含量计算---*/Р %macro twosample_ratio(p1=,p2=,alpha=,beta=,q1=,q2=);Р data d1_6;Р %let pc=&p1*&q1+&p2*&q2;