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 "EXAMPLED01.SAS--Power Confidence Limit due to varhat";
2          *TITLE2 "Figure 1 Taylor and Muller, 1995, Amer Statistician, 49,
2        ! p43-47";
3          
4          FILENAME OUT01 ".\EXAMPLED01.CGM";
5          
6          PROC IML WORKSIZE=1000 SYMSIZE=2000;
Worksize    = 1024000
Symbol size = 2048000
NOTE: IML Ready
7          %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.
NOTE: Module _ALOG1 defined.
NOTE: Module _EXP1 defined.
NOTE: Module _ORDER defined.
NOTE: Module _ERRBD defined.
NOTE: Module _CTFF defined.
2 The SAS System                                   13:47 Monday, August 11, 2003

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.
4951       
4952       OPT_OFF = {ALPHA};
4952     !                         *Turn options off;
4953       OPT_ON  = {NOPRINT DS};
4953     !                         *Turn options on;
4954       
4955       *Model is Y = XB + E;
4956       *Specify essence of X matrix, the basic design;
4957       ESSENCEX=I(2);
4957     !                 *Balanced two group t test, cell mean coding;
4958       REPN={12};
4959       
4960       *Hypothesis, as usual in ANOVA, is grand mean invariant so subtract;
4961       BETA = {0 1}`;
4962       BETASCAL=DO(0,.75,.01);
4963       *Gives power for Beta={0  0}` , {0  .01}` , . . . {0  .75}`;
4964       
4965       *Hypothesis is Ho: CBU=Theta0 ;
4966       ALPHA={.01};
4966     !               *Test size for hypothesis test;
4967       C    ={1 -1};
4967     !               * U=I and Theta0=0 are defaults; *This is t-test;
4968       
4969       SIGMA  = {.068};
4969     !                  *Multivariate notation=Variance in univariate;
4970       SIGSCAL= {1}   ;
4970     !                  *Set of covariance matrix multipliers;
4971       
4972       *Statements to create confidence limits;
4973       CLTYPE=1;
4974       N_EST=21;
4974     !                  *# Obs for variance estimate;
4975       RANK_EST=1;
4975     !                  *# model df for study giving variance estimate;
4976       ALPHA_CL=.025;
4976     !                  *Lower confidence limit tail size;
4977       ALPHA_CU=.025;
4977     !                  *Upper confidence limit tail size;
4978       RUN POWER;
NOTE: No data sets in the library WORK.
NOTE: The data set WORK._PWRDTMP has 76 observations and 9 variables.
3 The SAS System                                   13:47 Monday, August 11, 2003

4978     !            *IML and power code end here.  Power program creates file
4978     ! PWRDT1;
4979       **********************************************************;
4980       
4981       *Remainder of code creates plot;
4982       GOPTIONS GSFNAME=OUT01 DEVICE=CGMOF97P
4983       CBACK=WHITE COLORS=(BLACK) HORIGIN=0IN VORIGIN=0IN
4984       HSIZE=5IN VSIZE=3IN HTEXT=12PT FTEXT=TRIPLEX;
4985       
4986       SYMBOL1 I=SPLINE V=NONE L=34 W=2;
4987       SYMBOL2 I=SPLINE V=NONE L= 1 W=2;
4988       SYMBOL3 I=SPLINE V=NONE L=34 W=2;
4989       AXIS1 ORDER=(0 TO 1 BY .2)    W=3 MINOR=NONE MAJOR=(W=2)
4990             LABEL=(ANGLE=-90 ROTATE=90);
4991       AXIS2 ORDER=(0 TO .75 BY .25) W=3 MINOR=NONE MAJOR=(W=2)
4992             LABEL=("Mean Difference, 1/Cr (dL/mg)");
4993       
NOTE: Exiting IML.
NOTE: The PROCEDURE IML printed page 1.
NOTE: PROCEDURE IML used:
      real time           0.92 seconds
      cpu time            0.74 seconds
      

4994       PROC GPLOT DATA=PWRDT1;
4995       PLOT (POWER_L POWER POWER_U)*BETASCAL/ OVERLAY NOFRAME HREF=.5
4996              VZERO VAXIS=AXIS1 HZERO HAXIS=AXIS2 NOLEGEND;
4997       LABEL POWER_L="Power"  POWER  ="Power"  POWER_U="Power"
4998             SIGSCAL="Variance Multiplier";
NOTE:  134 RECORDS WRITTEN TO 
      D:\KEM\RESEARCH\SOFTWARE\POWERV2\V202\EXAMPLES\EXAMPLED01.CGM
NOTE: There were 76 observations read from the data set WORK.PWRDT1.
NOTE: PROCEDURE GPLOT used:
      real time           0.31 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.65 seconds
      cpu time            1.11 seconds