SAS宏入门讲解.pptVIP

%macro pw (file,a,b,c,po); data file; mean1=a; mean2=b; v=c; means=abs(mean1-mean2); do n=10 to 500 by 10; n1=n/2; n2=n/2; df=n-2; sn=1/n1+1/n2; t=means/v*sqrt(sn); ualpha=probit(0.975); ubeta=t-ualpha; po=probnorm(ubeta); beta=1-po; output; end; run; %mend pw; %pw (p1,105,115,12,power1); proc print noobs; run; %macro po (x,p,y,z); data x (keep=n p); pai1=y; pai2=z; meanpai=(pai1+pai2)/2; pai=abs(pai1-pai2); mse=meanpai*(1-meanpai); do n=10 to 1000 by 20; ns=sqrt(n);r=2*sqrt(mse); ubeta=ns*pai/r-1.96; p=probnorm(ubeta); output; end; run; %mend po; %macro sorting (file); proc sort data=file; by n; run; %mend sorting; /* Invoking macro */ %po (p1,power1,0.5,0.1); %po (p2,power2,0.5,0.2); %po (p3,power3,0.5,0.3); %po (p4,power4,0.5,0.35); %po (p5,power5,0.5,0.4); %po (p6,power6,0.5,0.45); /*sorting for data */ %sorting (p1); %sorting (p2); %sorting (p3); %sorting (p4); %sorting (p5); %sorting (p6); data total; merge p1 p2 p3 p4 p5 p6; by n; run; quit; goptions device=WIN ctext=blue graphrc interpol=join; axis1 color=blue width=2.0 length=40 ; axis2 color=blue width=2.0 length=15 ; proc gplot data=WORK.TOTAL; plot (POWER1 POWER2 POWER3 POWER4 POWER5 POWER6) * N / overlay haxis=axis1 vaxis=axis2 ; run; quit; /*************************************************/ /* 两样本均数比较的POWER计算程序 */ /*************************************************/ %macro pw; data file; mean1=a; mean2=b; v=c; means=abs(mean1-mean2); do n=10 to 500 by 10; n1=n/2; n2=n/2; df=n-2; sn=1/n1+1/n2; t=means/(v*sqrt(sn)); ualpha=probit(0.975); ubeta=t-ualpha; po=probnorm(ubeta); beta=1-po; output; end; run; %mend pw; %macro sorting; proc sort data=f; by n; run; %mend sorting; %macro analysit (file,a,b,c,po,f); %pw;