data example7_5;
input temp water time y@;
cards;
210	2.0	30	209
210	4.0	40	208
220	2.0	40	215
220	4.0	30	220
230	2.0	40	252
230	4.0	30	248
240	2.0	30	237
240	4.0	40	231
;
proc glm;
class temp water time;
model y=temp water time;
lsmeans temp/adjust=bon tdiff;
run;
