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.22 seconds
      
1          TITLE1 "EXAMPLEC01.SAS--power for Time by Treatment Interaction";
2          TITLE2 "Compares power of tests when SIGMA most favorable to UN";
3          
4          PROC IML SYMSIZE=1000 WORKSIZE=2000;
Worksize    = 2048000
Symbol size = 1024000
NOTE: IML Ready
5          %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.
NOTE: Module _TRUNCN defined.
NOTE: Module _FINDU defined.
NOTE: Module _INTEGR defined.
2 The SAS System                                   13:47 Monday, August 11, 2003

NOTE: Module _CFE defined.
NOTE: Module NAMELIST defined.
NOTE: Module UMEAN defined.
NOTE: Module UPOLY1 defined.
NOTE: Module UPOLY2 defined.
NOTE: Module UPOLY3 defined.
4949       
4950       OPT_OFF={ALPHA};
4951       OPT_ON = {ORTHU UN HF GG BOX  HLT PBT WLK MMETHOD  UMETHOD MMETHOD};
4952       * Specifying the option ORTHU in OPT_ON allows the program to provide
4952     ! ;
4953       * an orthonormal U matrix if one is not given by the user;
4954       * This is the case for the following code;
4955       
4956       P=3;
4957       Q=4;
4958       C=J(Q-1,1,1)||(-I(Q-1));
4959       U=( J(P-1,1,1)||(-I(P-1)) )`;
4960       ALPHA=.01;
4961       
4962       VARIANCE=1;
4962     !              RHO=0.4;
4963       SIGMA=VARIANCE#(I(P)#(1-RHO) + J(P,P,RHO));
4963     !                                             *Compound symmetry;
4964       SIGSCAL={1, 2};
4965       
4966       ESSENCEX=I(Q);
4966     !                REPN={5,10};
4967       BETA=J(Q,P,0);
4968       BETA[1,1]=1;
4969       BETASCAL=DO(0, 2.0 , 0.50);
4970       
4971       MMETHOD={4,4,4};
4971     !                    *two moment null approximations + OS noncen
4971     ! multiplier ON;
4972       
4973       UCDF = {4,2,2,4};
4973     !                    * UN and Box (4) exact via Davies' algorithm
4973     ! (1980), as in;
4974                          * Muller Edwards Taylor (2003). If exact fails
4974     ! then ;
4975                          * switch to approximation 2, MET 2003;
4976                          * HF and GG: (2) Muller, Edwards, and Taylor
4976     ! (2004) approx;
4977       
4978       RUN POWER;
4979       QUIT;
NOTE: Exiting IML.
NOTE: The PROCEDURE IML printed pages 1-5.
NOTE: PROCEDURE IML used:
3 The SAS System                                   13:47 Monday, August 11, 2003

      real time           0.72 seconds
      cpu time            0.68 seconds
      

NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414
NOTE: The SAS System used:
      real time           1.18 seconds
      cpu time            0.91 seconds