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;
     

/* block analysis using proc glm */         
proc glm;
  class block trt;
  model height =  block trt;
  lsmeans trt /stderr;
  estimate 'diff' trt 1 -1;
  output out = resids r = resid p = yhat;
  
  title 'Fixed blocks';

run;

proc plot;
  plot resid*yhat;
  title 'Residual plot';
  
run;