options ls=75; /* Sensitivity for a one-sample (paired) or two-sample t-test. Input parameters: p=alpha level (.05 is default) power=power(.80 is default) m=Sample Size sample=1 or 2 for one vs. 2 sample test side=1 or 2 (One sided vs. 2 sided) Output parameters: Sample size, actual power, cutoff for significance, Effect Size in Standard Deviations, and Percentile of the Median (alternate Hypothesis). */ data t;input p power m sample side; if p<0 or p>1 then p=.05; if power<0 or power>1 then power=.8; if side=. then side=2; if sample ne 1 then sample=2; if abs(side-1)>.1 then side=2; za=probit(1-(p/side));zb=probit(1-power); delinit=(abs(za)+abs(zb))*sqrt(sample/m)-.001; d=.001;i=0; t=za-.001; test=delinit;delinit=test; aa:t=t+d; i=i+1; y=side*probt(-t,sample*(m-1)); if y-p<0 or i>3000 then goto bb; goto aa; bb:del=delinit-.001; count2=0; cc:;count2=count2+1; del=del+.001; nc=sqrt(m/sample)*del; ypow=1-probt(t,sample*(m-1),nc); if ypow