options ls=75 nodate;

data one;   
   infile "bashor.dat" firstobs=3;
   input day site y;
   ly=log(y);
   sy=sqrt(y);
   location=site;
run;


proc glm;
   class day location;
   model y ly=location|day;
   output out=two r=residual residual_log p=predicted predicted_log;
run;
options ls=64 ps=45;


symbol1 value=dot color=black;
symbol2 value=square color=black;
symbol3 value=triangle color=black;
symbol4 value=diamond color=black;

axis1 offset=(1,1) label=(height=3);
axis2 offset=(1,1) label=(height=3 angle=90);

legend1 label=(height=2);
/*
proc gplot;
   plot y*location=day/haxis=axis1 vaxis=axis2; 
run;
*/
/*
goptions  dev=psepsf;

proc gplot data=two; 
   title "residuals versus predicted - log transform";
   *plot residual_log*predicted_log=location/haxis=axis1 vaxis=axis2 legend=legend1; 
   plot residual*predicted=location/haxis=axis1 vaxis=axis2 legend=legend1; 
run;

proc gplot; 
   title "residuals versus predicted - log transform";
   plot residual_log*predicted_log=location/haxis=axis1 vaxis=axis2 legend=legend1; 
run;
*/
proc mixed method=type3 cl;
   class day location;
   model ly=location/ddfm=satterth outp=predz;
   random day day*location;
   lsmeans location/adj=tukey;
run;
goptions  dev=psepsf;
symbol i=join r=4;
proc gplot;
   plot pred*location=day/haxis=axis1 vaxis=axis2 legend=legend1;
run;
endsas;
proc mixed cl;
   class day location;
   model ly=location/ddfm=satterth;
   random day day*location;
   lsmeans location/adj=bon;
run;