1 The SAS System 13:47 Monday, August 11, 2003 NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0) Licensed to UNC-CHAPEL HILL-SYSTEMWIDE-T/R, Site 0039466002. NOTE: This session is executing on the WIN_PRO platform. NOTE: SAS initialization used: real time 0.26 seconds cpu time 0.26 seconds 1 *TITLE1 "EXAMPLEA02.SAS--Two sample T-test Plot"; 2 3 FILENAME OUT01 ".\EXAMPLEA02.CGM"; 4 5 *************************************************************; 6 * Performs power calculations for a two sample T test, ; 7 * replicating the plots in "Increasing scientific power with ; 8 * statistical power", by K.E. Muller and V.A. Benignus, ; 9 * Neurotoxicology and Teratology, vol 14, May-June, 1992 ; 10 *************************************************************; 11 12 PROC IML SYMSIZE=1000 WORKSIZE=2000; Worksize = 2048000 Symbol size = 1024000 NOTE: IML Ready 13 %INCLUDE "..\IML\POWERLIB202.IML"/NOSOURCE2; NOTE: Module POWER defined. NOTE: Module _POWER defined. NOTE: Module _RANKSYMM defined. NOTE: Module _SIZECHK defined. NOTE: Module _TYPEMISSCHK defined. NOTE: Module _SETOPT defined. NOTE: Module _HLT defined. NOTE: Module _PBT defined. NOTE: Module _WLK defined. NOTE: Module _SPECIAL defined. NOTE: Module _FIRSTUNI defined. NOTE: Module _HFEXEPS defined. NOTE: Module _GGEXEPS defined. NOTE: Module _LASTUNI defined. NOTE: Module _GLMMPCL defined. NOTE: Module _PROBF defined. NOTE: Module _FINV defined. NOTE: Module _FWARN defined. NOTE: Module _SASDS defined. NOTE: Module _QPROB defined. NOTE: Module _AS defined. NOTE: Module _IROUND defined. NOTE: Module _COUNTR defined. 2 The SAS System 13:47 Monday, August 11, 2003 NOTE: Module _ALOG1 defined. NOTE: Module _EXP1 defined. NOTE: Module _ORDER defined. NOTE: Module _ERRBD defined. NOTE: Module _CTFF defined. NOTE: Module _TRUNCN defined. NOTE: Module _FINDU defined. NOTE: Module _INTEGR defined. NOTE: Module _CFE defined. NOTE: Module NAMELIST defined. NOTE: Module UMEAN defined. NOTE: Module UPOLY1 defined. NOTE: Module UPOLY2 defined. NOTE: Module UPOLY3 defined. 4957 4958 * Define inputs to power program; 4959 4960 SIGMA = {1}; 4961 SIGSCAL= {0.32 1.00 2.05}; 4962 4963 ESSENCEX = I(2); 4964 REPN = { 10 }; 4965 4966 BETA = {0 1}`; 4967 BETASCAL = DO(0, 2.5, 0.05); 4968 C = {1 -1}; 4969 4970 OPT_ON = {DS NOPRINT}; 4971 4972 RUN POWER; NOTE: No data sets in the library WORK. NOTE: The data set WORK._PWRDTMP has 153 observations and 5 variables. 4973 4974 *The DS option creates a SAS file with the power calculation results. 4974 ! ; 4975 *Since no name for the file was specified, WORK.PWRDT1 is created.; 4976 4977 QUIT; NOTE: Exiting IML. NOTE: The PROCEDURE IML printed page 1. NOTE: PROCEDURE IML used: real time 0.76 seconds cpu time 0.73 seconds 4978 4979 PROC CONTENTS DATA=PWRDT1; 4980 RUN; NOTE: PROCEDURE CONTENTS used: 3 The SAS System 13:47 Monday, August 11, 2003 real time 0.02 seconds cpu time 0.02 seconds NOTE: The PROCEDURE CONTENTS printed page 2. 4981 4982 * Create Plots *; 4983 4984 PROC SORT DATA=PWRDT1 4985 OUT=ONE; 4986 BY BETASCAL SIGSCAL; 4987 4988 * create file for power curves of varying VARIANCE *; 4989 NOTE: There were 153 observations read from the data set WORK.PWRDT1. NOTE: The data set WORK.ONE has 153 observations and 5 variables. NOTE: PROCEDURE SORT used: real time 0.02 seconds cpu time 0.02 seconds 4990 PROC TRANSPOSE DATA=ONE OUT=TWO PREFIX=SIGPWR; 4991 VAR POWER; 4992 BY BETASCAL; 4993 4994 * create ANNOTATE dataset and assign symbols for labeling plots *; 4995 NOTE: There were 153 observations read from the data set WORK.ONE. NOTE: The data set WORK.TWO has 51 observations and 5 variables. NOTE: PROCEDURE TRANSPOSE used: real time 0.02 seconds cpu time 0.02 seconds 4996 DATA LABELS (KEEP= X Y XSYS YSYS TEXT STYLE SIZE); 4997 LENGTH TEXT $ 5 STYLE $ 8; 4998 XSYS="2"; YSYS="2"; 4999 5000 X=.26; Y=.95; TEXT="s"; STYLE="cgreek"; *SIZE=1.0; OUTPUT; 5001 X=.31; Y=.97; TEXT="2"; STYLE="TRIPLEX"; SIZE=.75; OUTPUT; 5002 X=.50; Y=.95; TEXT="=0.32"; STYLE="TRIPLEX"; SIZE=1.0; OUTPUT; 5003 X=.76; Y=.70; TEXT="s"; STYLE="cgreek"; SIZE=1.0; OUTPUT; 5004 X=.81; Y=.72; TEXT="2"; STYLE="TRIPLEX"; SIZE=.75; OUTPUT; 5005 X=1.00; Y=.70; TEXT="=1.00"; STYLE="TRIPLEX"; SIZE=1.0; OUTPUT; 5006 X=1.01; Y=.15; TEXT="s"; STYLE="cgreek"; SIZE=1.0; OUTPUT; 5007 X=1.06; Y=.17; TEXT="2"; STYLE="TRIPLEX"; SIZE=.75; OUTPUT; 5008 X=1.25; Y=.15; TEXT="=2.05"; STYLE="TRIPLEX"; SIZE=1.0; OUTPUT; 5009 RUN; 4 The SAS System 13:47 Monday, August 11, 2003 NOTE: The data set WORK.LABELS has 9 observations and 7 variables. NOTE: DATA statement used: real time 0.03 seconds cpu time 0.01 seconds 5010 5011 * The plot will be saved to this file for future inclusion in a 5011 ! document.*; 5012 GOPTIONS GSFNAME=OUT01 DEVICE=CGMOF97P 5013 CBACK=WHITE COLORS=(BLACK) HORIGIN=0IN VORIGIN=0IN 5014 HSIZE=5IN VSIZE=3IN HTEXT=12PT FTEXT=TRIPLEX; 5015 5016 SYMBOL1 I=JOIN V=NONE L=34 W=1.0; 5017 SYMBOL2 I=JOIN V=NONE L=1 W=1.0; 5018 SYMBOL3 I=JOIN V=NONE L=34 W=1.0; 5019 AXIS1 ORDER=(0 TO 1 BY .1) W=1.5 MINOR=NONE MAJOR=(W=1.5) 5020 LABEL=(ANGLE=-90 ROTATE=90); 5021 AXIS2 ORDER=(0 TO 2.5 BY .5) W=1.5 MINOR=NONE MAJOR=(W=1.5); 5022 5023 * The plot overlays power curves for three different variances *; 5024 PROC GPLOT DATA=TWO; 5025 PLOT SIGPWR1*BETASCAL=1 5026 SIGPWR2*BETASCAL=2 5027 SIGPWR3*BETASCAL=3/OVERLAY VAXIS=AXIS1 HAXIS=AXIS2 5027 ! ANNOTATE=LABELS; 5028 LABEL SIGPWR1="Power" SIGPWR2="Power" SIGPWR3="Power" 5029 BETASCAL="Mean Difference"; 5030 RUN; NOTE: 73 RECORDS WRITTEN TO D:\KEM\RESEARCH\SOFTWARE\POWERV2\V202\EXAMPLES\EXAMPLEA02.CGM 5031 QUIT; NOTE: There were 51 observations read from the data set WORK.TWO. NOTE: There were 9 observations read from the data set WORK.LABELS. NOTE: PROCEDURE GPLOT used: real time 0.26 seconds cpu time 0.10 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: The SAS System used: real time 1.92 seconds cpu time 1.22 seconds