1                                                          The SAS System                                  11:51 Sunday, May 7, 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.10 seconds


NOTE: AUTOEXEC processing completed.

1          options ls=79 ps=100;
2          
3             /* B=BARN M=MIXED A=ALFALFA */
4          /*
5          data a;
6             infile 'depies.dat' firstobs=2;
7             input cow trt $ pas barn pad parity $ wk kmy kfat pfat kpro ppro 
scc fcm;
8          proc sort; by parity trt pad wk cow;
9          
10         proc glm;
11            class barn pas pad parity wk cow;
12            model kmy = barn pas(barn) pad(barn*pas)
13               parity parity*barn parity*pas(barn) cow(parity*pad*barn*pas)
14               wk wk*barn wk*pas(barn) wk*parity wk*barn*parity 
wk*parity*pas(barn);
15            random pad(barn*pas) cow(parity*pad*barn*pas);
16            test h = barn pas(barn) e = pad(barn*pas);
17            test h = parity parity*barn parity*pas(barn) e = 
cow(parity*pad*barn*pas);
18            lsmeans barn pas(barn) / e = pad(barn*pas) pdiff stderr;
19         
20         proc mixed;
21            class barn pas pad parity wk cow;
22            model kmy = barn|pas(barn)|parity|wk;
23            random pad(barn*pas) cow(parity*pad*barn*pas);
24            lsmeans pas(barn) barn;
25         
26         proc mixed;
27            class barn pas pad parity wk cow;
28            model kmy = barn|pas(barn)|parity|wk;
29            random pad(barn*pas) cow(parity*pad*barn*pas);
30            lsmeans pas(barn) barn;
31         
32         proc mixed;
33            class barn pas pad parity wk cow;
34            model kmy = barn|pas(barn)|parity|wk;
35            random pad(barn*pas);
36            repeated / type=cs subject = cow(parity*pad);
37            lsmeans barn pas(barn);
38         */
39         
40         data b;
41            infile 'depies3.dat' firstobs=2;
42            input cow trt $ pad parity $
43               kmy1 kmy2 kmy3 kmy4 kmy5 kmy6 kmy7 kmy8 kmy9 kmy10 kmy11 
kmy12 kmy13;
44            /* trt: M=MIXED A=ALFALFA */
45         

NOTE: The infile 'depies3.dat' is:
      File Name=/usr1/y/yandell/.SAS/depies3.dat,
      Owner Name=yandell,Group Name=calsstat,
      Access Permission=rw-r--r--,
      File Size (bytes)=3047

NOTE: 30 records were read from the infile 'depies3.dat'.
      The minimum record length was 0.
      The maximum record length was 101.
NOTE: SAS went to a new line when INPUT statement reached past the end of a 
      line.
NOTE: The data set WORK.B has 29 observations and 17 variables.
NOTE: DATA statement used:
      real time           0.37 seconds
      cpu time            0.09 seconds


46         proc glm;
47            class trt pad parity cow;
48            model kmy1-kmy13 = trt pad(trt) parity parity*trt;
49            test h = trt e = pad(trt);
50            lsmeans parity*trt / out=lsm;
51         

NOTE: The data set WORK.LSM has 52 observations and 5 variables.
NOTE: The PROCEDURE GLM printed pages 1-15.
NOTE: PROCEDURE GLM used:
      real time           0.77 seconds
      cpu time            0.12 seconds


2                               The SAS System        11:51 Sunday, May 7, 1995

52         data c; set lsm;
53            week = substr(_NAME_,4,2);

NOTE: The data set WORK.C has 52 observations and 6 variables.
NOTE: DATA statement used:
      real time           0.15 seconds
      cpu time            0.03 seconds


54         proc sort; by trt;

NOTE: The data set WORK.C has 52 observations and 6 variables.
NOTE: PROCEDURE SORT used:
      real time           0.13 seconds
      cpu time            0.01 seconds


55         proc plot; by trt;
56            plot lsmean*week=parity;
57         

NOTE: The PROCEDURE PLOT printed pages 16-17.
NOTE: PROCEDURE PLOT used:
      real time           0.01 seconds
      cpu time            0.02 seconds


58         proc glm data=b;
59            class trt pad parity cow;
60            model kmy1-kmy13 = trt pad(trt) parity parity*trt / nouni;
61            repeated week 13 (1 2 3 4 5 6 7 8 9 10 11 12 13) polynomial / 
summary;
NOTE: The PROCEDURE GLM printed pages 18-23.
NOTE: PROCEDURE GLM used:
      real time           0.67 seconds
      cpu time            0.09 seconds


NOTE: The SAS System used:
      real time           2.44 seconds
      cpu time            0.51 seconds

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