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.25 seconds cpu time 0.25 seconds 1 TITLE1 "EXAMPLEH01.SAS-two within factors, UPOLY2, no between, vary 1 ! SIGMA"; 2 TITLE2 "From Coffey C.S. and Muller K.E. (2003) Properties of "; 3 TITLE3 "internal pilots with the univariate approach to repeated 3 ! measures"; 4 TITLE4 "Statistics in Medicine, 22(15), in press."; 5 6 PROC IML SYMSIZE=1000 WORKSIZE=2000; Worksize = 2048000 Symbol size = 1024000 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. 2 The SAS System 13:47 Monday, August 11, 2003 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. 4951 4952 ALPHA=.04; 4953 OPT_ON={NOPRINT GG HF BOX TOTAL_N UCDF UMETHOD}; 4954 OPT_OFF={WARN ALPHA BETASCAL HLT PBT WLK }; 4955 ROUND=2; 4956 4957 BETASCAL=1; 4958 THETA={.25}#{.5 1 -1 .5}; 4958 ! * =Theta(cr) from 1st sentence *after* 4959 * equation 7, Coffey and Muller (2003); 4960 4961 *Following from Table II in Coffey and Muller (2003); 4962 VARSTARE={.47960 .01000 .01000 .01000}; 4962 ! *epsilon ~ .28; 4963 VARSTARF={.34555 .06123 .05561 .04721}; 4963 ! *epsilon ~ .50; 4964 VARSTARG={.23555 .17123 .05561 .04721}; 4964 ! *epsilon ~ .72; 4965 VARSTARH={.12740 .12740 .12740 .12740}; 4965 ! *epsilon = 1; 4966 VARSTAR=VARSTARE//VARSTARF//VARSTARG//VARSTARH; 4967 4968 SIGSCAL={0.50 1.00 2.00}; 4968 ! * <=> gamma in Coffey and Muller (2003); 4969 4970 *log base 2 spacing Clip (2,4,16) and Region(2,8,32); 4971 *Get orthonormal U matrices; 4972 RUN UPOLY2({1 2 4},"A", {1 3 5},"B", 4973 UA,NMA, UB,NMB, UAB ,NMAB); 4974 U=UAB; 4975 C=1; 4976 4977 ESSENCEX={1}; 4978 REPN={20}; 4979 4980 DO IVAR=1 TO 4 BY 1; 4981 SIGSTAR=DIAG(VARSTAR[IVAR,*]); 4982 4983 SIGMA=U*SIGSTAR*U`; 4983 ! *1st paragraph in section 2.4, Coffey and 3 The SAS System 13:47 Monday, August 11, 2003 4983 ! Muller 2003; 4984 BETA=THETA*U`; 4984 ! *1st paragraph in section 2.4, Coffey and 4984 ! Muller 2003; 4985 4986 DO VERSION=1 TO 2 BY 1; 4986 ! *POWERLIB version; 4987 UCDF=J(4,1,VERSION); 4988 UMETHOD=J(2,1,VERSION); 4989 RUN POWER; 4990 HOLDALL=HOLDALL//_HOLDPOWER; 4991 END; 4992 END; 4993 4994 CREATE ONE VAR _HOLDPOWERLBL; 4995 APPEND FROM HOLDALL; 4996 4997 NOTE: Exiting IML. NOTE: The data set WORK.ONE has 24 observations and 13 variables. NOTE: PROCEDURE IML used: real time 0.76 seconds cpu time 0.72 seconds 4998 PROC SORT DATA=ONE OUT=TWO; 4999 BY UCDF_GG UMETHOD_GG SIGSCAL EPSILON; 5000 NOTE: There were 24 observations read from the data set WORK.ONE. NOTE: The data set WORK.TWO has 24 observations and 13 variables. NOTE: PROCEDURE SORT used: real time 0.03 seconds cpu time 0.03 seconds 5001 PROC PRINT DATA=TWO UNIFORM NOOBS; 5002 BY UCDF_GG UMETHOD_GG UCDF_HF UMETHOD_HF UCDF_BOX TOTAL_N; 5003 PAGEBY UCDF_GG; 5004 TITLE5 "All data in file"; 5005 NOTE: There were 24 observations read from the data set WORK.TWO. NOTE: The PROCEDURE PRINT printed pages 1-2. NOTE: PROCEDURE PRINT used: real time 0.02 seconds cpu time 0.02 seconds 5006 PROC PRINT DATA=TWO(RENAME=(SIGSCAL=GAMMA)) UNIFORM NOOBS; 4 The SAS System 13:47 Monday, August 11, 2003 5007 VAR EPSILON GAMMA POWER_GG POWER_HF ; 5008 BY UCDF_GG UMETHOD_GG UCDF_HF UMETHOD_HF UCDF_BOX TOTAL_N; 5009 PAGEBY UCDF_GG; 5010 TITLE5 "Version 2 far more accurate for Table III, Coffey and Muller 5010 ! (2003)"; NOTE: There were 24 observations read from the data set WORK.TWO. NOTE: The PROCEDURE PRINT printed pages 3-4. NOTE: PROCEDURE PRINT used: real time 0.01 seconds cpu time 0.01 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: The SAS System used: real time 1.41 seconds cpu time 1.05 seconds