options ls=79 ps=100;

data a;
   infile 'depies.dat' firstobs=2;
   input cow trt $ pas barn pad parity $ wk kmy kfat pfat kpro ppro scc fcm;
   /* B=BARN M=MIXED A=ALFALFA */
proc sort; by parity trt pad wk cow;

/* model kmy kfat pfat kpro ppro scc fcm = */
proc glm;
   class barn pas pad parity wk cow;
   model kmy = barn pas(barn) pad(barn*pas)
      parity parity*barn parity*pas(barn) cow(parity*pad*barn*pas)
      wk wk*barn wk*pas(barn) wk*parity wk*barn*parity wk*parity*pas(barn);
   random pad(barn*pas) cow(parity*pad*barn*pas);
   test h = barn pas(barn) e = pad(barn*pas);
   test h = parity parity*barn parity*pas(barn) e = cow(parity*pad*barn*pas);
   lsmeans barn pas(barn) / e = pad(barn*pas) pdiff stderr;

proc mixed;
   class barn pas pad parity wk cow;
   model kmy = barn|pas(barn)|parity|wk;
   random pad(barn*pas) cow(parity*pad*barn*pas);
   lsmeans pas(barn) barn;

proc mixed;
   class barn pas pad parity wk cow;
   model kmy = barn|pas(barn)|parity|wk;
   random pad(barn*pas) cow(parity*pad*barn*pas) / type=ar(2);
   lsmeans pas(barn) barn;

proc mixed;
   class barn pas pad parity wk cow;
   model kmy = barn|pas(barn)|parity|wk;
   random pad(barn*pas);
   repeated / type=cs subject = cow(parity*pad);
   lsmeans barn pas(barn);