1、单组配对定量资料假设检验POWER
proc power;
onesamplemeans test=t
mean=25
stddev=22
ntotal=12
sides=1
power=.;
run;
quit;
2、多个样本均数比较的POWER
proc power ;
onewayanova
groupmeans =288.9 | 657.9 | 762.9
stddev = 173.6 154.7 127.1
groupweights = (1 1 1)
alpha = 0.05
ntotal =30
power =. ;
plot x =n min=6 max=45 key =bycurve(pos=inset numbers=off);
run;
quit;
3、定性资料比较两组优效性POWER
%let n_t=120;
%let n_r=120;
%let p_t=0.9667;
%let p_r=0.9250;
%let delta_U=0.03;
%let alpha=0.05;
data DA8_8;
power=probnorm(((&p_t-&p_r)-&delta_U)/sqrt(&p_t*( 1-&p_t)/&n_t+&p_r*(1-&p_r)/&n_r)-probit(1-&alpha));
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;

4、定性资料比较两组非劣效性POWER
%let n_t=150;
%let n_r=150;
%let p_t=0.58;
%let p_r=0.46;
%let delta_L=-0.20;
%let delta_U=0.20;
%let alpha=0.05;
data DA8_4;
power=probnorm((&delta_U-(&p_t-&p_r))/sqrt(&p_t*( 1-&p_t)/&n_t+&p_r*(1-&p_r)/&n_r)-probit(1-&alpha/ 2))+probnorm((-&delta_L+(&p_t-&p_r))/sqrt(&p_t*(1-&p_t)/&n_t+&p_r*( 1-&p_r)/&n_r)-probit(1-&alpha/2))- 1;
if power< 0 then power=0;
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;
5、定性资料比较两组等效性
%let n_t=150;
%let n_r=150;
%let p_t=0.58;
%let p_r=0.46;
%let delta_L=-0.20;
%let delta_U=0.20;
%let alpha=0.05;
data DA8_4;
power=probnorm((&delta_U-(&p_t-&p_r))/sqrt(&p_t*( 1-&p_t)/&n_t+&p_r*(1-&p_r)/&n_r)-probit(1-&alpha/ 2))+probnorm((-&delta_L+(&p_t-&p_r))/sqrt(&p_t*(1-&p_t)/&n_t+&p_r*( 1-&p_r)/&n_r)-probit(1-&alpha/2))- 1;
if power< 0 then power=0;
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;

6、定量资料比较两组优效性POWER
%let n_t=39;
%let n_r=43;
%let x_bar_t=1.25;
%let x_bar_r=0.64;
%let s_t=1.13;
%let s_r=0.51;
%let delta_U=0.28;
%let alpha=0.05;
data DA8_6;
power=probt(((&x_bar_t-&x_bar_r)-&delta_U)/sqrt(( 1/&n_t+1/&n_r)*((&n_t-1)*&s_t** 2+(&n_r-1)*&s_r**2)/(&n_t+&n_r- 2))-tinv(1-&alpha,&n_t+&n_r-2),&n_t+&n_r- 2);
if power< 0 then power=0;
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;
7、定量资料比较两组非劣效性POWER
%let n_t=339;
%let n_r=113;
%let x_bar_t=8.58;
%let x_bar_r=9.12;
%let s_t=3.97;
%let s_r=3.61;
%let delta_L=-1.77;
%let alpha=0.05;
data DA8_10;
power=probt(((&x_bar_t-&x_bar_r)-&delta_L)/sqrt(( 1/&n_t+1/&n_r)*((&n_t-1)*&s_t** 2+(&n_r-1)*&s_r**2)/(&n_t+&n_r- 2))-tinv(1-&alpha,&n_t+&n_r-2),&n_t+&n_r- 2);
if power< 0 then power=0;
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;
8、成组资料比较两组等效性POWER
%let n_t=160;
%let n_r=160;
%let x_bar_t=13.29;
%let x_bar_r=14.87;
%let s_t=6.10;
%let s_r=5.84;
%let delta_L=-5;
%let delta_U=5;
%let alpha=0.05;
data DA8_2;
power=probt((-&delta_L+(&x_bar_t-&x_bar_r))/sqrt(( 1/&n_t+1/&n_r)*((&n_t-1)*&s_t** 2+(&n_r-1)*&s_r**2)/(&n_t+&n_r- 2))-tinv(1-&alpha/2,&n_t+&n_r- 2),&n_t+&n_r-2)+probt((&delta_U-(&x_bar_t-&x_bar_r))/sqrt(( 1/&n_t+1/&n_r)*((&n_t-1)*&s_t** 2+(&n_r-1)*&s_r**2)/(&n_t+&n_r- 2))-tinv(1-&alpha/2,&n_t+&n_r- 2),&n_t+&n_r-2)-1;
if power< 0 then power=0;
run;
ods html ;
proc print noobs ;
var power;
run;
ods html close ;

9、成组设计两样本率的POWER
Proc power ;
TwoSampleFreq test=PChi
GroupProportions=( 0.700 0.844 )
Groupns=50| 45
Sides=2
Alpha=0.05
Power=.;
Run;
Quit;
10、成组设计两样本均数的POWER
proc power ;
twosamplemeans test =diff
groupmeans = 1.50 | 1.20
stddev =0.83
groupns = (30 30)
power = . ;
run;
quit;