options nocenter ps=40 ls=72;

data new;
input trt plot sub resp;
cards;
1 1 1 2
1 1 2 3
1 2 1 2
1 2 2 3.5
1 3 1 2.5
1 3 2  3
1 4 1 4.0
1 4 2 2.5
2 1 1 4.5
2 1 2 4.0
2 2 1 3.5
2 2 2 3.5
2 3 1 5
2 3 2 4.5
2 4 1 4.5
2 4 2 4
3 1 1 3.5
3 1 2 1.5
3 2 1 2.0
3 2 2 2.0
3 3 1 1.5
3 3 2 2.0
3 4 1 2.0
3 4 2 2.5
4 1 1 6.0
4 1 2 6.5
4 2 1 5.5
4 2 2 5.5
4 3 1 5.0
4 3 2 5.5
4 4 1 4.5
4 4 2 5.0
;

proc sort;
 by trt plot;

proc means noprint;
 var resp;
 by trt plot;
 output out=new1 mean=respmn;

proc glm;
 class trt;
 model respmn=trt;

proc glm data=new;
 class trt plot sub;
 model resp=trt plot(trt);
 test h=trt e=plot(trt);