/* test macro for continuous variables */
%include "WilcxPowerContinuousNties.sas";
options symbolgen formdlim=']';
filename vlmr '/udd/nhron/macrodevelopment/WilcxRank/data/vlmrccrp.dat';
data vlmrccrp;
infile vlmr firstobs=4;
input trt 1-2 @4 crpchg 7.4;
crpchg1=round(crpchg,0.01);
run;
%WilcxPowerContinuousNties(m=100,n=100,alpha=0.05,theta=, mu=,estimate=Y,dsn=vlmrccrp,var1=trt,var2=crpchg);
%WilcxPowerContinuousNties(m=100,n=100,alpha=0.05,theta=0.574, mu=,estimate=N,dsn=,var1=,var2=);
Output:
The SAS System 10:59 Wednesday, March 12, 2008 2
The FREQ Procedure
Cumulative Cumulative
trt Frequency Percent Frequency Percent
--------------------------------------------------------
1 649 51.18 649 51.18
2 619 48.82 1268 100.00
]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
The SAS System 10:59 Wednesday, March 12, 2008 3
Obs wstat uformula ustat theta
1 382244.5 171319.5 230411.5 0.57355
]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
Power Estimation for Sample1 (100) and Sample2 (100) with theta=0.5735467265 4
10:59 Wednesday, March 12, 2008
m n theta vvar a b z power1 power2 power
100 100 0.57355 .004952032 -0.094759 -2.18503 1.95996 0.46225 0.014443 0.47670
]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
Power Estimation for Sample1 (100) and Sample2 (100) with theta=0.574 5
10:59 Wednesday, March 12, 2008
m n theta vvar a b z power1 power2 power
100 100 0.574 .004945551 -0.088376 -2.19290 1.95996 0.46479 0.014157 0.47895