1 The SAS System 16:03 Monday, May 1, 1995 NOTE: Copyright(c) 1989 by SAS Institute Inc., Cary, NC USA. NOTE: SAS (r) Proprietary Software Release 6.09 TS027 Licensed to UNIVERSITY OF WISCONSIN, Site 0002176014. NOTE: AUTOEXEC processing beginning; file is /usr/sas609/autoexec.sas. NOTE: SAS initialization used: real time 0.30 seconds cpu time 0.12 seconds NOTE: AUTOEXEC processing completed. 1 options nocenter ps=64 ls=79; 2 3 data a; 4 infile 'dhiman.dat' firstobs=2; 5 input cow trt pr carry bcaa1 bcaa2 bcaa3; NOTE: The infile 'dhiman.dat' is: File Name=/usr1/y/yandell/.SAS/dhiman.dat, Owner Name=yandell,Group Name=calsstat, Access Permission=rw-r-----, File Size (bytes)=860 NOTE: 25 records were read from the infile 'dhiman.dat'. The minimum record length was 32. The maximum record length was 32. NOTE: The data set WORK.A has 25 observations and 7 variables. NOTE: DATA statement used: real time 0.29 seconds cpu time 0.07 seconds 6 data dhiman; set a; 7 bcaa = bcaa1; spr = 1; output; 8 bcaa = bcaa2; spr = 2; output; 9 bcaa = bcaa3; spr = 3; output; NOTE: The data set WORK.DHIMAN has 75 observations and 9 variables. NOTE: DATA statement used: real time 0.17 seconds cpu time 0.04 seconds 10 proc sort data=dhiman; by spr; 11 NOTE: The data set WORK.DHIMAN has 75 observations and 9 variables. NOTE: PROCEDURE SORT used: real time 0.18 seconds cpu time 0.02 seconds 12 proc glm; by spr; /* separate analysis by subperiod */ 13 class cow trt pr spr carry; /* check for possible carryover */ 14 model bcaa = cow trt pr carry; /* (no evidence of carryover) */ 15 lsmeans trt pr / stderr pdiff; 2 The SAS System 16:03 Monday, May 1, 1995 16 NOTE: The PROCEDURE GLM printed pages 1-9. NOTE: PROCEDURE GLM used: real time 0.80 seconds cpu time 0.11 seconds 17 proc glm; /* overall split plot analysis */ 18 class cow trt pr spr; 19 model bcaa = cow trt pr cow*trt*pr spr spr*trt spr*pr; 20 random cow cow*trt*pr / test; 21 test h = trt pr e = cow*trt*pr; 22 lsmeans trt pr / e = cow*trt*pr stderr pdiff; 23 lsmeans spr / stderr pdiff; /* not estimable due to proc glm bug */ 24 lsmeans trt*spr / stderr pdiff out=lsmsp; 25 NOTE: TYPE I EMS not available without the E1 option. NOTE: The data set WORK.LSMSP has 15 observations and 5 variables. NOTE: The PROCEDURE GLM printed pages 10-18. NOTE: PROCEDURE GLM used: real time 0.99 seconds cpu time 0.22 seconds 26 proc plot data=lsmsp; 27 plot lsmean*trt=spr; 28 NOTE: The PROCEDURE PLOT printed page 19. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds 29 proc mixed data=dhiman; 30 class cow trt pr spr; 31 model bcaa = trt pr spr spr*trt spr*pr; 32 random cow cow*trt*pr; 33 lsmeans trt*spr; /* LS means done properly */ 34 make 'lsmeans' out=lsmsub; NOTE: 3 observations are not included because of missing values. NOTE: Estimated G matrix is not positive definite. NOTE: The data set WORK.LSMSUB has 15 observations and 6 variables. NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed pages 20-21. NOTE: PROCEDURE MIXED used: real time 0.98 seconds cpu time 0.21 seconds 35 data plotsub; set lsmsub; /* plotting kludge */ 36 trt = substr(level,9,1); 37 spr = substr(level,11,1); NOTE: The data set WORK.PLOTSUB has 15 observations and 8 variables. NOTE: DATA statement used: real time 0.17 seconds 3 The SAS System 16:03 Monday, May 1, 1995 cpu time 0.05 seconds 38 proc plot; 39 plot lsmean*trt=spr; NOTE: The PROCEDURE PLOT printed page 22. NOTE: PROCEDURE PLOT used: real time 0.00 seconds cpu time 0.00 seconds NOTE: The SAS System used: real time 4.61 seconds cpu time 0.91 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414