/*  SAS code for analyzing data from the gasoline
    additive data from problem 4 on assignment 8  */

data set1;
  infile 'c:\courses\st511\sas\gas.additive.dat';
  input driver auto additive $ y;
  run;


/*  Print the data  */

proc print data=set1;
run;



/*  UNIX users should delete the next line and replace
    it with the lines in the next comment  */

goptions device=WIN target=WINPRTC rotate=landscape;



  proc mixed data=set1;
    class driver auto additive;
    model y = additive / solution
                e1 e3 outpm=set1 outp=set2;;
    random driver auto / solution;
    ods output solutionr=setr;
	ods output covparms=setv;
    lsmeans additive / pdiff tdiff;
    lsmeans additive / pdiff adjust=Tukey; 
    run;

	proc print data=setr; run;

	proc print data=setv; run;

data driver; set setr;
  if(Effect="driver"); run;

data auto; set setr;
  if(Effect="auto"); run;

proc rank data=driver normal=blom out=driver ; 
  var driver; ranks q;  run;

axis1 label=(f=swiss h=2.5)
      value=(f=swiss h=2.0)  w=3.0  length= 6.5 in;

axis2 label=(f=swiss h=2.0 r=90 a=270)
      value=(f=swiss h=2.0) w= 3.0  length = 5.0 in;

 SYMBOL1 V=circle H=2.0 w=3 l=1 i=none ;
      

 proc gplot data=driver;
   plot driver*q / vaxis=axis2 haxis=axis1;
   title  H=3.0 F=swiss "Normal Plot: Driver Effects";
       label q='Standard Normal Quantiles ';
       label driver = 'Effects ';
 run;

proc rank data=auto normal=blom out=auto; 
  var auto; ranks q;  run;

proc gplot data=auto;
   plot auto*q / vaxis=axis2 haxis=axis1;
   title  H=3.0 F=swiss "Normal Plot: Auto Effects";
       label q='Standard Normal Quantiles ';
       label auto = 'Effects ';
 run;




proc rank data=set1 normal=blom out=set1 ; 
  var resid; ranks q;  run;

axis1 label=(f=swiss h=2.5)
      value=(f=swiss h=2.0)  w=3.0  length= 6.0 in;

axis2 label=(f=swiss h=2.0 r=90 a=270)
      value=(f=swiss h=2.0) w= 3.0  length = 4.5 in;

 SYMBOL1 V=circle H=2.0 w=3 l=1 i=none ;
      

 proc gplot data=set1;
   plot resid*q / vaxis=axis2 haxis=axis1;
   title  H=3.0 F=swiss "Normal Plot: Residuals";
       label q='Standard Normal Quantiles ';
       label resid = 'Residuals';
 run;

 
 proc gplot data=set1;
   plot resid*pred / vaxis=axis2 haxis=axis1;
   title  H=3.0 F=swiss "Residual Plot";
       label pred='Estimated Means';
       label resid = 'Residuals';
 run;