options nodate nonumber ls=76 ps=55;

title 'Kuehl -- Example 6.5 -- Shrimp Culture';

data shrimp;
do temp=25 to 35 by 10;
do density=80 to 160 by 80;
do saline=10 to 40 by 15;
do rep=1 to 3;
input wtgain @@;
output;
end; end; end; end;
datalines;
86 52 73 544 371 482 390 290 397 53 73 86 393 398 208 249 265 243
439 436 349 249 245 330 247 277 205 324 305 364 352 267 316 188 223 281
;

proc glm;
class temp density saline;
model wtgain = temp density saline temp*density temp*saline
   density*saline temp*density*saline;
run;



quit;