1 The SAS System 15:19 Sunday, April 30, 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.24 seconds cpu time 0.11 seconds NOTE: AUTOEXEC processing completed. 1 options nocenter ps=64 ls=80; 2 3 data phoffman; 4 infile 'phoffman.dat' firstobs=2; 5 input age herd trt cow season milk fat prot scc lscore stage; 6 /* season: 1=spring 2=summer 3=fall 4=winter */ NOTE: The infile 'phoffman.dat' is: File Name=/usr1/y/yandell/.SAS/phoffman.dat, Owner Name=yandell,Group Name=calsstat, Access Permission=rw-r--r--, File Size (bytes)=27030 NOTE: 442 records were read from the infile 'phoffman.dat'. The minimum record length was 60. The maximum record length was 61. NOTE: The data set WORK.PHOFFMAN has 442 observations and 11 variables. NOTE: DATA statement used: real time 0.33 seconds cpu time 0.11 seconds 7 proc sort; by age; 8 NOTE: The data set WORK.PHOFFMAN has 442 observations and 11 variables. NOTE: PROCEDURE SORT used: real time 0.21 seconds cpu time 0.04 seconds 9 proc glm; by age; 10 classes herd trt cow season; 11 model milk = herd trt cow(herd*trt) season trt*season / ss3; 12 test h=herd trt e=cow(herd*trt); 13 means season trt*season; 14 lsmeans trt / e=cow(herd*trt) pdiff stderr; 15 lsmeans trt*season / pdiff stderr out=lsm; 16 output out=resid r=remilk p=prmilk; 17 NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: The above message was for the following by-group: AGE=1 2 The SAS System 15:19 Sunday, April 30, 1995 NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: The above message was for the following by-group: AGE=2 NOTE: The data set WORK.RESID has 442 observations and 13 variables. NOTE: The data set WORK.LSM has 16 observations and 6 variables. NOTE: The PROCEDURE GLM printed pages 1-10. NOTE: PROCEDURE GLM used: real time 1.68 seconds cpu time 0.86 seconds 18 proc plot data=resid; by age; /* diagnostic plots */ 19 plot remilk*prmilk=trt; 20 plot remilk*stage=trt; 21 NOTE: The PROCEDURE PLOT printed pages 11-14. NOTE: PROCEDURE PLOT used: real time 0.04 seconds cpu time 0.04 seconds 22 proc plot data=lsm; by age; /* interaction plots */ 23 plot lsmean*season=trt; 24 NOTE: The PROCEDURE PLOT printed pages 15-16. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds 25 proc mixed data=phoffman; by age; 26 class herd trt cow season; 27 model milk = trt|season; 28 random herd cow(herd*trt); 29 lsmeans trt*season; 30 NOTE: Estimated G matrix is not positive definite. NOTE: The above message was for the following by-group: AGE=2 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 17-20. NOTE: PROCEDURE MIXED used: real time 0.69 seconds cpu time 0.51 seconds 31 proc mixed data=phoffman; by age; 32 class herd trt cow season; 33 model milk = trt|season; 34 random herd; 35 repeated / type=cs subject=cow(herd*trt); 36 lsmeans trt*season; 37 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 21-24. 3 The SAS System 15:19 Sunday, April 30, 1995 NOTE: PROCEDURE MIXED used: real time 1.12 seconds cpu time 0.94 seconds 38 proc glm data=phoffman; by age; 39 classes cow; 40 model milk = cow; 41 output out=rmilk r=rmilk; NOTE: The data set WORK.RMILK has 442 observations and 12 variables. NOTE: The PROCEDURE GLM printed pages 25-28. NOTE: PROCEDURE GLM used: real time 1.11 seconds cpu time 0.36 seconds 42 proc glm; by age; 43 classes trt cow season; 44 model rmilk = cow(trt) trt stage season trt*season / ss3; 45 lsmeans season trt*season; 46 means season trt*season; 47 output out=resids r=rrmilk p=pmilk; 48 NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: The above message was for the following by-group: AGE=1 NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: The above message was for the following by-group: AGE=2 NOTE: The data set WORK.RESIDS has 442 observations and 14 variables. NOTE: The PROCEDURE GLM printed pages 29-36. NOTE: PROCEDURE GLM used: real time 1.34 seconds cpu time 0.58 seconds 49 proc plot; by age; 50 plot rrmilk*pmilk=trt; 51 plot rrmilk*stage=trt; 52 53 NOTE: The PROCEDURE PLOT printed pages 37-40. NOTE: PROCEDURE PLOT used: real time 0.04 seconds cpu time 0.04 seconds NOTE: The SAS System used: real time 7.25 seconds cpu time 3.73 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414