1 The SAS System 13:07 Sunday, June 16, 1996 NOTE: Copyright (c) 1989-1996 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 6.11 TS040 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. The command line option "-nonews" will prevent this display. NOTE: AUTOEXEC processing beginning; file is /afs/.cs.wisc.edu/s/sas-6.11/alpha_osf20/sas611/autoexec.sas. NOTE: SAS initialization used: real time 1.76 seconds cpu time 0.60 seconds NOTE: AUTOEXEC processing completed. 1 options nocenter ps=64 ls=80; 2 3 data phoffman; 4 infile 'season.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 */ 7 lscc = log10(scc); NOTE: The infile 'season.dat' is: File Name=/afs/cs.wisc.edu/p/stat/src/statistics/pda/data/Season/season.dat, Owner Name=yandell,Group Name=yandell, Access Permission=rw-r--r--, File Size (bytes)=27030 NOTE: 442 records were read from the infile 'season.dat'. The minimum record length was 60. The maximum record length was 61. NOTE: Missing values were generated as a result of performing an operation on missing values. Each place is given by: (Number of times) at (Line):(Column). 2 at 7:11 NOTE: The data set WORK.PHOFFMAN has 442 observations and 12 variables. NOTE: DATA statement used: real time 1.70 seconds cpu time 0.48 seconds 8 proc sort; by age; 9 NOTE: The data set WORK.PHOFFMAN has 442 observations and 12 variables. NOTE: PROCEDURE SORT used: real time 0.75 seconds cpu time 0.13 seconds 10 proc glm; 2 The SAS System 13:07 Sunday, June 16, 1996 11 classes herd age trt cow season; 12 model lscc milk = herd age trt age*trt cow(herd*age*trt) 13 season season*trt season*age season*age*trt / ss3; 14 test h=herd age trt age*trt e=cow(herd*age*trt); 15 means age*trt*season; 16 /* 17 lsmeans trt / e=cow(herd*trt) pdiff stderr; 18 lsmeans trt*season / pdiff stderr out=lsm; 19 output out=resid r=remilk p=prmilk; 20 21 proc plot data=resid; by age; 22 plot remilk*prmilk=trt; 23 plot remilk*stage=trt; 24 25 proc plot data=lsm; by age; 26 plot lsmean*season=trt; 27 */ 28 NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. NOTE: The PROCEDURE GLM printed pages 1-5. NOTE: PROCEDURE GLM used: real time 12.23 seconds cpu time 8.08 seconds 29 proc mixed data=phoffman; 30 class herd age trt cow season; 31 model milk = age|trt|season; 32 random herd cow(herd*age*trt); 33 * lsmeans trt*season; 34 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 6-7. NOTE: PROCEDURE MIXED used: real time 4.45 seconds cpu time 2.36 seconds 35 proc mixed data=phoffman; 36 class herd age trt cow season; 37 model milk = age|trt|season; 38 random herd; 39 repeated / type=cs subject=cow(herd*age*trt); 40 * lsmeans trt*season; 41 42 /* 43 proc glm data=phoffman; by age; 44 classes cow; 45 model milk = cow; 46 output out=rmilk r=rmilk; 47 proc glm; by age; 48 classes trt cow season; 49 model rmilk = cow(trt) trt stage season trt*season / ss3; 50 lsmeans season trt*season; 51 means season trt*season; 52 output out=resids r=rrmilk p=pmilk; 3 The SAS System 13:07 Sunday, June 16, 1996 53 54 proc plot; by age; 55 plot rrmilk*pmilk=trt; 56 plot rrmilk*stage=trt; 57 */ 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 8-9. NOTE: PROCEDURE MIXED used: real time 4.86 seconds cpu time 3.13 seconds NOTE: The SAS System used: real time 26.45 seconds cpu time 15.26 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414