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.27 seconds
      cpu time            0.21 seconds
      
1          TITLE1 "EXAMPLEG02.SAS--UPOLY3 for contrasts";
2          TITLE2 "factorial design, repeated measures: A, B, C between, D, E, F
2        !  within";
3          
4          PROC IML SYMSIZE=1000 WORKSIZE=2000;
Worksize    = 2048000
Symbol size = 1024000
NOTE: IML Ready
5          RESET FUZZ NOAUTONAME FW=6 LINESIZE=80;
6          %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.
2 The SAS System                                   13:47 Monday, August 11, 2003

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.
4950       
4951       ALPHA=.05;
4952       
4953       *Choose dimensions of design;
4954       GA=3;
4954     !       *=# groups for between factor A;
4955       GB=3;
4955     !       *=# groups for between factor B;
4956       GC=3;
4956     !       *=# groups for between factor C;
4957       TD=3;
4957     !       *=#Times for within factor D;
4958       TE=3;
4958     !       *=#Times for within factor E;
4959       TF=3;
4959     !       *=#Times for within factor F;
4960       
4961       P=TD#TE#TF;
4962       Q=GA#GB#GC;
4963       ESSENCEX= I(Q);
4964       BETA=J(Q,P,0);
4965       BETA[1,1]=1;
4966       
4967       SIGMA=DIAG(DO(1,P,1));
4967     !                       *Variances are 1,2,3,...p;
4968       
4969       *Get orthonormal U matrices;
4970       CALL UPOLY3 ( (1:TD),"D",  (1:TE),"E",  (1:TF),"F",
4971       		              UD,UDLBL,    UE,UELBL,    UF,UFLBL,
4972                    UDE,UDELBL,  UDF,UDFLBL,  UEF,UEFLBL,  UDEF,UDEFLBL );
4973       
4974       *Get orthonormal C matrices;
4975       CALL UPOLY3 ((1:GA),"A" , (1:GB),"B" , (1:GC),"C",
4976       		             U1,CALBL,   U2,CBLBL,      U3,CCLBL,
4977                  U12,CABLBL, U13,CACLBL, U23,CBCLBL,  U123,CABCLBL);
4978       
4979       BETASCAL={9 18 27};
4980       ROUND=4;
4981       OPT_ON={NOPRINT  GG HF UN  PBT HLT WLK};
4982       OPT_OFF={WARN SIGSCAL ALPHA};
4983       BUG=" ";
4984       
3 The SAS System                                   13:47 Monday, August 11, 2003

4985       C=U1`;
4986       U=UD;
4987         DO REPN=2 TO 12 BY 2;
4988         RUN POWER;
4989         HOLDA=HOLDA//_HOLDPOWER;
4990         END;
4991       PRINT / "AxD";
4992       PRINT HOLDA[COLNAME=_HOLDPOWERLBL ROWNAME=BUG];
4993       
4994       C=U12`;
4995       U=UDE;
4996         DO REPN=2 TO 12 BY 2;
4997         RUN POWER;
4998         HOLDABDE=HOLDABDE//_HOLDPOWER;
4999         END;
5000       PRINT / "AxB x DxE Interaction";
5001       PRINT HOLDABDE[COLNAME=_HOLDPOWERLBL ROWNAME=BUG];
5002       
5003       C=U123`;
5004       U=UDEF;
5005         DO REPN=2 TO 12 BY 2;
5006         RUN POWER;
5007         HABCDEF=HABCDEF//_HOLDPOWER;
5008         END;
5009       PRINT / "AxBxC x DxExF Interaction";
5010       PRINT HABCDEF[COLNAME=_HOLDPOWERLBL ROWNAME=BUG];
NOTE: Exiting IML.
NOTE: The PROCEDURE IML printed pages 1-3.
NOTE: PROCEDURE IML used:
      real time           0.99 seconds
      cpu time            0.97 seconds
      

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