1                              The SAS System   16:03 Friday, August 5, 1994

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



1          TITLE1 "EXAMPLE2.SAS--Demonstrate LINMOD using MAKESS";
2          
3          * Data are from Danford, Hughes, and McNee,              *;
4          * "On the Analysis of Repeated Measurement Experiments," *;
5          * Biometrics, 16, 1960, pp. 547-565.                     *;
6          
7          * 0. Define raw data file;
8          FILENAME IN01 "C:\PROGRAMS\LINMOD31\EXAMPLES\PAYNE.DAT";
9          DATA PAYNE;
10              INFILE IN01;
11              INPUT @7 GROUP $ INIT SCORE2 SCORE4 SCORE6 SCORE8 SCORE10;
12         * Create indicator variables ;
13         CONTROL =(GROUP="1");
14         LOW     =(GROUP="2");
15         MODERATE=(GROUP="3");
16         HIGH    =(GROUP="4");
17         

NOTE: The infile IN01 is:
      FILENAME=C:\PROGRAMS\LINMOD31\EXAMPLES\PAYNE.DAT,
      RECFM=V,LRECL=132

NOTE: 45 records were read from the infile IN01.
      The minimum record length was 31.
      The maximum record length was 31.
NOTE: The data set WORK.PAYNE has 45 observations and 11 variables.
NOTE: The DATA statement used 2.68 seconds.


18         PROC PRINT DATA=PAYNE;
19         
20         *1. Define directory in which LINMOD source code has been stored;
21         %LET LMDIRECT = C:\PROGRAMS\LINMOD31\SOURCE\    ;
22         
23         *2. Define SAS macro's needed;
24         %INCLUDE "&LMDIRECT.MACROLIB.MAC" / NOSOURCE2 ;
41         
42         * 3. Reduce raw data to a TYPE=CORR file.  See Example 1;
43         *    Not done in this example due to using MAKESS;
44         
45         * 4.1 Start IML with at least minimum space needed;
46         *     Use SHOW SPACE statement and larger WORKSIZE, SYMSIZE 
values
2                              The SAS System   16:03 Friday, August 5, 1994

47         *     to tune performance with bigger problems;

NOTE: The PROCEDURE PRINT printed page 1.
NOTE: The PROCEDURE PRINT used 2.45 seconds.


48         PROC IML WORKSIZE=200 SYMSIZE=350;
Worksize    = 204800
Symbol size = 358400
IML Ready
49         
50         *4.2 Initialize LINMOD code;
51         &LINMOD ;
NOTE: Module _DEFOPT_ defined.
NOTE: Module _EREXIT_ defined.
NOTE: Module _FITMOD_ defined.
NOTE: Module _GETCSS_ defined.
NOTE: Module _INITL_ defined.
NOTE: Module _MAKESS_ defined.
NOTE: Module _NMLIST_ defined.
NOTE: Module _PRPARM_ defined.
NOTE: Module _PVTYPE_ defined.
NOTE: Module _SCAN_ defined.
NOTE: Module _SETOPT_ defined.
NOTE: Module _TSTGLH_ defined.
NOTE: Module _UNIREP_ defined.
NOTE: Module _UNSCAN_ defined.
NOTE: Module _ZPIVOT_ defined.
NOTE: Module FITMODEL defined.
NOTE: Module GETCORSS defined.
NOTE: Module MAKESS defined.
NOTE: Module SETOPT defined.
NOTE: Module TESTGLH defined.
NOTE: Module NAMELIST defined.
NOTE: Module UTREND defined.
NOTE: Module UMEAN defined.
NOTE: Module UPOLY1 defined.
NOTE: Module UPOLY2 defined.
2174       
2175       *5. Change and list options;
2176       OPT_OFF = { "MSH" };
2177       OPT_ON  = { "LISTINFO" "AVAILOPT" };
2178       RUN SETOPT;
2179       
2180       *6.1 Read raw data into IML;
2181       *    Use this approach only for small files;
2182       *    Use PROCSSCP macro otherwise;
2183       USE PAYNE;
2184       READ ALL VAR{"GROUP"} INTO GROUP;
2185       READ ALL VAR{"INIT"} INTO INIT;
3                              The SAS System   16:03 Friday, August 5, 1994

2186       READ ALL VAR{"SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10"} INTO 
Y;
2187       CLOSE PAYNE;
2188       
2189       *6.2 Use IML functions to create indicator variables for design 
matrix;
2190       N=NROW(Y);
2190                  * # observations in sample;
2191       CONSTANT=J(N,1,1);
2191                          *Column of 1's for intercept, etc;
2192       CELLMEAN=DESIGN(GROUP);
2192                                                      *Cell mean coding;
2193       EFFECT  =CONSTANT||DESIGNF(GROUP);
2193                                                      *Effect coding;
2194       REFERENC=CONSTANT||CELLMEAN(|*,2:NCOL(CELLMEAN)|);
2195                                                      *Reference cell 
coding;
2196       *Refer to IML manual for explanation of functions;
2197       *Use PROC PRINT and FREQ on GROUP to determine interpretation of 
BETA;
2198       
2199       *6.3 Assemble all predictors and responses into a single matrix,
2200            which must be called Z (to use MAKESS).  ZNAMES must also 
exist,
2201            with 1 row, with values labeling the columns of Z ;
2202       Z = CELLMEAN || INIT || Y;
2203       ZNAMES = { "CONTROL"  "LOW"  "MODERATE"  "HIGH"  "INIT" }
2204               || { "SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10" } ;
2205       FREE GROUP INIT Y ;
2206       
2207       *7. Create SSCP matrix and associated parameters;
2208       RUN MAKESS;
2209       FREE Z;
2209               *Saves space, but if helpful, use PROCSSCP macro, not 
MAKESS;
2210       
2211       *8. Fit a model;
2212       INDVARS = { "CONTROL" "LOW" "MODERATE" "HIGH" "INIT" };
2213       DEPVARS = { "SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10" };
2214       RUN FITMODEL;
2215       
2216       *9. Conduct any testimation desired;
2217       *C= ---- ;
2218       *U= ---- ;
2219       *RUN TESTGLH;
Exiting IML.
NOTE: The PROCEDURE IML printed pages 2-4.
NOTE: The PROCEDURE IML used 34.0 seconds.


4                              The SAS System   16:03 Friday, August 5, 1994

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