data one;
  input university chamber trt $ y;
  cards;
1 1 A 45.8
1 1 B 56.7
1 1 C 15.2
1 2 A 79.8
1 2 B 143.5
1 2 C 71.9
2 1 A 105.3
2 1 B 135.7
2 1 C 111.8
2 2 A 104
2 2 B 132
2 2 C 113.4
3 1 A 114.6
3 1 B 139.4
3 1 C 118.8
3 2 A 92.7
3 2 B 122.2
3 2 C 103
4 1 A 64.2
4 1 B 90.3
4 1 C 74.7
4 2 A 93.4
4 2 B 123
4 2 C 104.7
;

proc glm data=one;
  class university chamber trt;
  model y=university chamber(university) trt trt*university;
  random university chamber(university) trt*university;
  output out=check residual=ehat predicted=yhat;
run;

proc plot;
  plot ehat*yhat=university;
run;