Documentation Example 8 for PROC POWER
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: pwrex08 */
/* TITLE: Documentation Example 8 for PROC POWER */
/* (Customizing Plots) */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: power */
/* sample size */
/* power analysis */
/* plot statement */
/* graphs */
/* PROCS: POWER */
/* DATA: */
/* */
/* SUPPORT: John Castelloe */
/* REF: */
/* MISC: Graphs produced by ODS Graphics in PROC POWER */
/* require a style that distinguishes curves based */
/* on a combination of color, line style, and symbol */
/* marker. Styles that are well-suited for use in */
/* PROC POWER include: STATISTICAL, ANALYSIS, DEFAULT, */
/* LISTING, and HTMLBLUECML. */
/****************************************************************/
proc power;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = 0.9
ntotal = .;
run;
*
ods listing style=htmlbluecml;
ods graphics on;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = 0.9
ntotal = .;
plot x=power min=0.5 max=0.95;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = 0.9
ntotal = .;
plot y=power min=0.5 max=0.95;
run;
proc power;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = 0.5 0.95
ntotal = .;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 200;
plot x=n min=20 max=500;
run;
proc power plotonly;
twosamplemeans test=diff
meandiff = 3 6
stddev = 7 9
power = 0.9
ntotal = .;
plot x=effect min=3 max=6;
run;
proc power plotonly;
ods output plotcontent=PlotData;
twosamplemeans test=diff
groupmeans = 12 | 18
stddev = 7
groupweights = 2 | 1
power = .
ntotal = 20;
plot x=n min=20 max=50 npoints=20;
run;
proc print data=PlotData;
var NominalNTotal NTotal;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 18
stddev = 7
groupweights = 2 | 1
power = .
ntotal = 20;
plot x=n min=20 max=50 npoints=5;
run;
proc power plotonly;
twosamplemeans test=diff
nfractional
groupmeans = 12 | 18
stddev = 7
groupweights = 2 | 1
power = .
ntotal = 20;
plot x=n min=20 max=50 npoints=20;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 100;
plot x=n min=20 max=500
yopts=(ref=0.8 0.9);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 100;
plot x=n min=20 max=500
yopts=(ref=0.8 0.9 crossref=yes);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 100;
plot x=n min=20 max=500
xopts=(ref=100 crossref=yes);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
alpha = 0.01 0.025 0.1
power = .
ntotal = 100;
plot x=n min=20 max=500;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
alpha = 0.01 0.025 0.1
power = .
ntotal = 100;
plot x=n min=20 max=500
vary (linestyle, symbol, color);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 18 15
stddev = 7 9
alpha = 0.01 0.025 0.1
sides = 1 2
power = .
ntotal = 100;
plot x=n min=20 max=500
vary (linestyle by stddev,
symbol by alpha sides,
panel by groupmeans);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 200;
plot x=n min=20 max=500
key = byfeature(pos=inset);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 200;
plot x=n min=20 max=500
key = bycurve;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 200;
plot x=n min=20 max=500
key = bycurve(numbers=off pos=inset);
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 200;
plot x=n min=20 max=500
key = oncurves;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 232 382 60 98;
plot x=n min=20 max=500
markers=analysis;
run;
proc power plotonly;
twosamplemeans test=diff
groupmeans = 12 | 15 18
stddev = 7 9
power = .
ntotal = 232 382 60 98;
plot x=n min=20 max=500
markers=nice;
run;
ods graphics off;