options nocenter ls=75;

data new;
infile "coating.dat";
input site batch tablet resp;

proc glm;
class site batch;
model resp=site batch(site);
random batch;
test h=site e=batch(site);
output out=new1 p=pred r=res;

symbol1 v=circle;
proc gplot;
plot res*pred;

proc glm data=new;
class site batch;
model resp=site batch site*batch;
random batch site*batch;
test h=site e=batch*site;
run;