1                              The SAS System  21:19 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 "STORE2.PGM--test POWERLIB IML modules";
2          
3          LIBNAME LOAD01 "\KEM\POWER\STORED\" ;
NOTE: Libref LOAD01 was successfully assigned as follows: 
      Engine:        V608 
      Physical Name: C:\KEM\POWER\STORED
4          
5          PROC IML WORKSIZE=200 SYMSIZE=400;
Worksize    = 204800
Symbol size = 409600
IML Ready
6          
7          RESET STORAGE=LOAD01.POWSTORE;
NOTE: New storage library = LOAD01.POWSTORE
8          LOAD  MODULE = _ALL_           ;
9          
10         
11         ********************************************;
12         
13         OPT_OFF={CASE ALPHA RHOSCAL};
14         
15         *EXAMPLE 3: Multivariate Repeated Measures Model;
16         
17         *Define matrices *;
18         *RHO and BETA from published report or pilot study *;
19         *******************************
20         * BETA ROW1: LOW RISK CONTROL *
21         *      ROW2: LOW RISK TRTMNT  *
22         *      ROW3: HI RISK CONTROL  *
23         *      ROW4: HI RISK TRTMNT   *
24         *   COLUMNS:  18 24 36 MONTHS *
25         *******************************;
26         BETA   = {116.50 108.33 110.83 ,
27                   125.40 115.00 116.40 ,
28                   108.43  77.43  72.14 ,
29                   108.67  95.33  96.67 };
30         
31         RHO = { 1.00 0.79 0.59 ,
32                 0.79 1.00 0.70 ,
33                 0.59 0.70 1.00};
34         
35         ** Variance at 18, 24, and 36 months *;
36         V={256.958 256.958 253.777};
2                              The SAS System  21:19 Tuesday, August 2, 1994

37         D = DIAG(SQRT (V));
38         SIGMA=D*RHO*D;
39         
40         ESSENCEX=I(4);
41         BETASCAL={1};
42         SIGSCAL= {.5 1 };
43         RHOSCAL= {1 .5 };
44         OPT_ON={COLLAPSE};
45         REPN={12 24 48};
46         ROUND=3;
47         
48         * Create an orthnormal U matrix *;
49         TIMESOBS={18 24 36}`;
50         POLYTOBS=ORPOL(TIMESOBS);
51         UTIME=POLYTOBS(|,2:3|);
52         UAVE=POLYTOBS(|,1|);
53         U=UTIME;
54         
55         CMAX  ={1  1 -1 -1,
56                 1 -1  1 -1,
57                 1 -1 -1  1};
58         CTESTS ={1,2,3};
59         
60         *1) Test for any group difference *;
61         C=CMAX;
62         RUN POWER;
63         
64         *2) Tests for individual contrasts *;
65         DO i=1 TO NROW(CTESTS);
66            t= CTESTS(|i,|);
67            C= CMAX(|t,|);
68         
69            RUN POWER;
70         END;
Exiting IML.
NOTE: Storage library LOAD01.POWSTORE closed.
NOTE: The PROCEDURE IML printed pages 1-9.
NOTE: The PROCEDURE IML used 11.32 seconds.


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