1 The SAS System 16:06 Sunday, April 7, 1996 NOTE: Copyright (c) 1989-1993 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 6.10 TS018 Licensed to UNIVERSITY OF WISCONSIN, Site 0002176032. This message is contained in the SAS news file, and is presented upon initialization. Edit the files "news" in the "misc/base" directory to display site-specific news and information in the program log. All products for this release are production except for the following: SAS/CALC - EXPERIMENTAL SAS/CONNECT - BETA SAS/EIS - EXPERMENTAL SAS/OR - BETA SAS/QC - BETA SAS/SHARE - BETA NOTE: AUTOEXEC processing beginning; file is /afs/cs.wisc.edu/s/sas-6.10/alpha_osf20/sas610/autoexec.sas. NOTE: SAS initialization used: real time 0.68 seconds cpu time 0.55 seconds NOTE: AUTOEXEC processing completed. 1 options nocenter ps=64 ls=80; /* Two-Way Random Effects Model */ 2 data twoway; 3 infile 'random.dat'; /* MJ 18/19 */ 4 input row col y; 5 NOTE: The infile 'random.dat' is: File Name=/afs/cs.wisc.edu/p/stat/Data/MJ/random.dat, Owner Name=yandell,Group Name=yandell, Access Permission=rw-r-----, File Size (bytes)=98 NOTE: 14 records were read from the infile 'random.dat'. The minimum record length was 6. The maximum record length was 6. NOTE: The data set WORK.TWOWAY has 14 observations and 3 variables. NOTE: DATA statement used: real time 0.80 seconds cpu time 0.28 seconds 6 proc glm; /* THE OLD PROC GLM */ 7 class row col; 8 model y = row col row*col / ss1 ss2 ss3 solution; 9 random row col row*col / test; 10 estimate 'grand mean' intercept 1; 11 estimate 'row 1 col 1' intercept 1 row 1 col 1 row*col 1; 12 estimate 'row 1 col 2' intercept 1 row 1 col 0 1 row*col 0 1; 13 estimate 'row 1 col 3' intercept 1 row 1 col 0 0 1 row*col 0 0 1; 14 estimate 'row 2 col 1' intercept 1 row 0 1 col 1 row*col 0 0 0 1; 15 estimate 'row 2 col 2' intercept 1 row 0 1 col 0 1 row*col 0 0 0 0 1; 2 The SAS System 16:06 Sunday, April 7, 1996 16 estimate 'row 2 col 3' intercept 1 row 0 1 col 0 0 1 row*col 0 0 0 0 0 1; 17 NOTE: TYPE I EMS not available without the E1 option. NOTE: TYPE II EMS not available without the E2 option. NOTE: The PROCEDURE GLM printed pages 1-5. NOTE: PROCEDURE GLM used: real time 2.53 seconds cpu time 1.06 seconds 18 proc varcomp; /* THE BETTER PROC VARCOMP */ 19 class row col; 20 model y = row col row*col; NOTE: The PROCEDURE VARCOMP printed pages 6-7. NOTE: PROCEDURE VARCOMP used: real time 0.63 seconds cpu time 0.16 seconds 21 proc varcomp method=type1; 22 class row col; 23 model y = row col row*col; NOTE: The PROCEDURE VARCOMP printed pages 8-9. NOTE: PROCEDURE VARCOMP used: real time 0.40 seconds cpu time 0.06 seconds 24 proc varcomp method=ml; 25 class row col; 26 model y = row col row*col; NOTE: The PROCEDURE VARCOMP printed pages 10-11. NOTE: PROCEDURE VARCOMP used: real time 0.56 seconds cpu time 0.08 seconds 27 proc varcomp method=reml; 28 class row col; 29 model y = row col row*col; 30 NOTE: The PROCEDURE VARCOMP printed pages 12-13. NOTE: PROCEDURE VARCOMP used: real time 0.68 seconds cpu time 0.11 seconds 31 proc mixed; /* THE NEW PROC MIXED */ 32 class row col; 33 model y = ; 34 random row col row*col / solution; 35 estimate 'grand mean' intercept 1; 36 estimate 'row 1 col 1' intercept 1 | row 1 col 1 row*col 1; 37 estimate 'row 1 col 2' intercept 1 | row 1 col 0 1 row*col 0 1; 38 estimate 'row 1 col 3' intercept 1 | row 1 col 0 0 1 row*col 0 0 1; 3 The SAS System 16:06 Sunday, April 7, 1996 39 estimate 'row 2 col 1' intercept 1 | row 0 1 col 1 row*col 0 0 0 1; 40 estimate 'row 2 col 2' intercept 1 | row 0 1 col 0 1 row*col 0 0 0 0 1; 41 estimate 'row 2 col 3' intercept 1 | row 0 1 col 0 0 1 row*col 0 0 0 0 0 1; NOTE: Estimated G matrix is not positive definite. NOTE: The PROCEDURE MIXED printed pages 14-15. NOTE: PROCEDURE MIXED used: real time 0.85 seconds cpu time 0.50 seconds NOTE: The SAS System used: real time 9.45 seconds cpu time 3.38 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414