options ls=75 nodate;

data one;
   infile "milk.dat" firstobs=4;
   input sample lab y;
   ly=log(y);
run;

*goptions colors=(black) device=pslepsf;

proc gplot;
   plot y*sample=lab;
   plot ly*sample=lab;
run;

/*proc print;run;*/

proc glm;
   title "log-transformed data";
   class lab sample;
   *model y=sample|lab;
   model y=sample|lab;
   random sample lab sample*lab;
   test h=lab sample e=sample*lab;
   output out=two p=p r=r;
   lsmeans sample*lab;
run;
proc gplot;
   plot r*p;
run;
proc varcomp;
   title "log-transformed data";
   class sample lab;
   model ly=sample|lab;
run;

proc mixed cl method=type3; 
   title "log-transformed data";
   class sample lab;
   model ly=/s ddfm=satterth cl;
   random sample lab sample*lab;
run;
   
/*proc gplot;
   title "log-transformed data";
   plot r*p;
run;

proc univariate normal;
   title "log-transformed data";
   hist r;
   var r;
run;

proc glm;
   title "raw data";
   class lab sample;
   model y=sample|lab;
   random sample lab sample*lab;
   test h=lab sample e=sample*lab;
   output out=two p=p r=r;
   lsmeans sample*lab;
run;
proc varcomp;
   title "raw data";
   class sample lab;
   model y=sample|lab;
run;*/