1                              The SAS System  21:18 Tuesday, August 2, 1994

NOTE: Copyright(c) 1989 by SAS Institute Inc., Cary, NC USA. 
NOTE: SAS (r) Proprietary Software Release 6.08  TS405
      Licensed to UNIVERSITY OF NORTH CAROLINA AT CHAPEL HILL, Site 
      0003944001.



1          TITLE1 "STORE3.PGM--run without storing for timing comparison";
2          
3          PROC IML WORKSIZE=200 SYMSIZE=400;
Worksize    = 204800
Symbol size = 409600
IML Ready
4          
5          %INCLUDE "\KEM\POWER\DISTRIB\POWERLIB.IML";
NOTE: Module _INPTCHK defined.
NOTE: Module _OPTCHK defined.
NOTE: Module _SNGRCHK defined.
NOTE: Module _SIZECHK defined.
NOTE: Module _TYPECHK defined.
NOTE: Module _SETOPT defined.
NOTE: Module _PROBF defined.
NOTE: Module _HLTMOD defined.
NOTE: Module _PBTMOD defined.
NOTE: Module _WLKMOD defined.
NOTE: Module _SPECIAL defined.
NOTE: Module _UNIMOD defined.
NOTE: Module _GGEXEPS defined.
NOTE: Module _HFEXEPS defined.
NOTE: Module _LASTUNI defined.
NOTE: Module _SASDS defined.
NOTE: Module _POWER defined.
NOTE: Module POWER defined.
1728       
1729       ********************************************;
1730       
1731       OPT_OFF={CASE ALPHA RHOSCAL};
1732       
1733       *EXAMPLE 3: Multivariate Repeated Measures Model;
1734       
1735       *Define matrices *;
1736       *RHO and BETA from published report or pilot study *;
1737       *******************************
1738       * BETA ROW1: LOW RISK CONTROL *
1739       *      ROW2: LOW RISK TRTMNT  *
1740       *      ROW3: HI RISK CONTROL  *
1741       *      ROW4: HI RISK TRTMNT   *
1742       *   COLUMNS:  18 24 36 MONTHS *
1743       *******************************;
1744       BETA   = {116.50 108.33 110.83 ,
2                              The SAS System  21:18 Tuesday, August 2, 1994

1745                 125.40 115.00 116.40 ,
1746                 108.43  77.43  72.14 ,
1747                 108.67  95.33  96.67 };
1748       
1749       RHO = { 1.00 0.79 0.59 ,
1750               0.79 1.00 0.70 ,
1751               0.59 0.70 1.00};
1752       
1753       ** Variance at 18, 24, and 36 months *;
1754       V={256.958 256.958 253.777};
1755       D = DIAG(SQRT (V));
1756       SIGMA=D*RHO*D;
1757       
1758       ESSENCEX=I(4);
1759       BETASCAL={1};
1760       SIGSCAL= {.5 1 };
1761       RHOSCAL= {1 .5 };
1762       OPT_ON={COLLAPSE};
1763       REPN={12 24 48};
1764       ROUND=3;
1765       
1766       * Create an orthnormal U matrix *;
1767       TIMESOBS={18 24 36}`;
1768       POLYTOBS=ORPOL(TIMESOBS);
1769       UTIME=POLYTOBS(|,2:3|);
1770       UAVE=POLYTOBS(|,1|);
1771       U=UTIME;
1772       
1773       CMAX  ={1  1 -1 -1,
1774               1 -1  1 -1,
1775               1 -1 -1  1};
1776       CTESTS ={1,2,3};
1777       
1778       *1) Test for any group difference *;
1779       C=CMAX;
1780       RUN POWER;
1781       
1782       *2) Tests for individual contrasts *;
1783       DO i=1 TO NROW(CTESTS);
1784          t= CTESTS(|i,|);
1785          C= CMAX(|t,|);
1786       
1787          RUN POWER;
1788       END;
Exiting IML.
NOTE: The PROCEDURE IML printed pages 1-9.
NOTE: The PROCEDURE IML used 18.71 seconds.


NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414