options formdlim = '-';

data mead;
     
  infile 'mead53.txt' firstobs = 2;

/* for pc, replace with something like:
  infile 'c:/Documents and Settings/pdixon.IASTATE/Desktop/mead53a.txt';
    or (for file in my documents):
  infile 'c:/Documents and Settings/pdixon.IASTATE/My Documents/mead53a.txt';
  */
  
  input block trt $ height;
  run;
     

/* blocks as random effect using proc mixed */         
proc mixed method = type3;
  class block trt;
  model height =  trt /outp = resids;
  random block;  
  lsmeans trt;
  estimate 'diff' trt 1 -1;
  title 'Random blocks';
run;

proc print data = resids;
  title "Residual data set";
run;
     
proc plot;
  plot resid*pred;
run;