/* %last test = 113/.99 = what given last month
 cow milk(lb) %last fat% prot% scc(1000) name dim lact farm trt
*/
options ls=79 ps=100;

data shaver;
   infile 'shaver.dat' missover firstobs=2;
   input cow milk lstpct fpct ppct scc id $ dim lact farm $ test ;
   mlm = milk + ((dim-150)*.0029*milk);
   fcm = .4*milk + 15*(milk*(fpct/100));
   mlfcm = fcm + ((dim-150)*.0029*fcm);
   if lact=1 then parity='first'; else parity='mature';
proc sort; by farm test;

/* looks like 4-5 treatments per named cow (cow and name same, id may
   change).  what are the key questions?
 */