Documentation Example 1 for PROC GLMPOWER
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: gmpex01 */
/* TITLE: Documentation Example 1 for PROC GLMPOWER */
/* (One-Way ANOVA) */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: power */
/* power analysis */
/* sample size */
/* graphs */
/* general linear models */
/* analysis of variance */
/* PROCS: GLMPOWER */
/* DATA: */
/* */
/* SUPPORT: John Castelloe */
/* REF: */
/* MISC: Graphs produced by ODS Graphics in PROC GLMPOWER */
/* 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 GLMPOWER include: STATISTICAL, ANALYSIS, */
/* DEFAULT, LISTING, and HTMLBLUECML. */
/****************************************************************/
data Fluids;
input Fluid $ LacticAcid1 LacticAcid2 CellWgt;
datalines;
Water 35.6 35.6 2
EZD1 33.7 33.7 1
EZD2 30.2 30.2 1
LZ1 29 28 1
LZ2 25.9 25.9 1
;
proc glmpower data=Fluids;
class Fluid;
model LacticAcid1 LacticAcid2 = Fluid;
weight CellWgt;
contrast "Water vs. others" Fluid -1 -1 -1 -1 4;
contrast "EZD vs. LZ" Fluid 1 1 -1 -1 0;
contrast "EZD1 vs. EZD2" Fluid 1 -1 0 0 0;
contrast "LZ1 vs. LZ2" Fluid 0 0 1 -1 0;
power
stddev = 3.75
alpha = 0.025
ntotal = .
power = 0.9;
run;
*
ods listing style=htmlbluecml;
ods graphics on;
proc glmpower data=Fluids plotonly;
class Fluid;
model LacticAcid1 LacticAcid2 = Fluid;
weight CellWgt;
contrast "Water vs. others" Fluid -1 -1 -1 -1 4;
contrast "EZD vs. LZ" Fluid 1 1 -1 -1 0;
contrast "EZD1 vs. EZD2" Fluid 1 -1 0 0 0;
contrast "LZ1 vs. LZ2" Fluid 0 0 1 -1 0;
power
stddev = 3.75
alpha = 0.025
ntotal = .
power = 0.9;
plot x=power min=.5 max=.95;
run;
proc glmpower data=Fluids plotonly;
class Fluid;
model LacticAcid1 LacticAcid2 = Fluid;
weight CellWgt;
contrast "Water vs. others" Fluid -1 -1 -1 -1 4;
contrast "EZD vs. LZ" Fluid 1 1 -1 -1 0;
contrast "EZD1 vs. EZD2" Fluid 1 -1 0 0 0;
contrast "LZ1 vs. LZ2" Fluid 0 0 1 -1 0;
power
stddev = 3.75
alpha = 0.025
ntotal = 24
power = .;
plot x=n min=24 max=480;
run;
ods graphics off;