options ls=130 ps=600;

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

proc means noprint;
   by parity trt pad wk cow;
   var kmy kfat pfat kpro ppro scc fcm;
   output out = newdhi1 mean = ;
data a; set newdhi1;
   barn = 0;
   if trt = 'B' then barn = 1;
   pas = 0;
   if trt = 'M' then pas =1;
   if trt = 'A' then pas = -1;
   if cow = 3612 then cow = 3647;
proc sort; by wk pas trt pad parity;
proc print;
   var cow trt pas barn pad parity wk kmy kfat pfat kpro ppro scc fcm;