1                                The SAS System    12:21 Monday, August 24, 1998

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



1          TITLE1 "EXAMPLE4.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 "..\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\LINMOD32\EXAMPLES\PAYNE.DAT,
      RECFM=V,LRECL=256

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 1.93 seconds.


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


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


49         PROC IML WORKSIZE=2000 SYMSIZE=4000;
Worksize    = 2048000
Symbol size = 4096000
IML Ready
50         
51         *4.2 Initialize LINMOD code;
52         &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.
2286       
2287       %INCLUDE "\KEM\TESTPROG\UNIREPX.IML";
NOTE: Redefining existing module _UNIREP_.
NOTE: Module _UNIREP_ defined.
2398       
2399       
2400       *5. Change and list options;
2401       OPT_ON  = { NOPRINT };
2402       RUN SETOPT;
2403       
2404       *6.1 Read raw data into IML;
2405       *    Use this approach only for small files;
2406       *    Use PROCSSCP macro otherwise;
3                                The SAS System    12:21 Monday, August 24, 1998

2407       USE PAYNE;
2408       READ ALL VAR{"GROUP"} INTO GROUP;
2409       READ ALL VAR{"INIT"} INTO INIT;
2410       READ ALL VAR{"SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10"} INTO Y;
2411       CLOSE PAYNE;
2412       
2413       *6.2 Use IML functions to create indicator variables for design 
matrix;
2414       N=NROW(Y);
2414                  * # observations in sample;
2415       CONSTANT=J(N,1,1);
2415                          *Column of 1's for intercept, etc;
2416       CELLMEAN=DESIGN(GROUP);
2416                                                      *Cell mean coding;
2417       EFFECT  =CONSTANT||DESIGNF(GROUP);
2417                                                      *Effect coding;
2418       REFERENC=CONSTANT||CELLMEAN(|*,2:NCOL(CELLMEAN)|);
2419                                                      *Reference cell 
coding;
2420       *Refer to IML manual for explanation of functions;
2421       *Use PROC PRINT and FREQ on GROUP to determine interpretation of 
BETA;
2422       
2423       *6.3 Assemble all predictors and responses into a single matrix,
2424            which must be called Z (to use MAKESS).  ZNAMES must also exist,
2425            with 1 row, with values labeling the columns of Z ;
2426       Z = CELLMEAN || INIT || Y;
2427       ZNAMES = { "CONTROL"  "LOW"  "MODERATE"  "HIGH"  "INIT" }
2428               || { "SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10" } ;
2429       FREE GROUP INIT Y ;
2430       
2431       *7. Create SSCP matrix and associated parameters;
2432       RUN MAKESS;
2433       FREE Z;
2433               *Saves space, but if helpful, use PROCSSCP macro, not MAKESS;
2434       
2435       *8. Fit a model;
2436       INDVARS = { "CONTROL" "LOW" "MODERATE" "HIGH" "INIT" };
2437       DEPVARS = { "SCORE2" "SCORE4" "SCORE6" "SCORE8" "SCORE10" };
2438       RUN FITMODEL;
2439       
2440       *9. Conduct any testimation desired;
2441       C= J(4,1,-1)||I(4);
2442       U=I(5)  ;
2443       RESET PRINT;
2444       RUN TESTGLH;
2445       PRINT "Following print from user control, not LINMOD:  " ,
2446         _STMAT1_  , _URESUL_;
ERROR: Matrix _URESUL_ has not been set to a value.

4                                The SAS System    12:21 Monday, August 24, 1998

 statement : PRINT           at line  2445 column   1
Exiting IML.
NOTE: The PROCEDURE IML printed pages 2-30.
NOTE: The PROCEDURE IML used 4.28 seconds.


ERROR: Errors printed on page 3.

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