data a;
input water rsp temp rz @@;
cards;
20        35        95        353
20        35        125        183
20        45        95        284
20        45        125        152
24        35        95        471
24        35        125        424
24        45        95        396
24        45        125        373
18        40        110        186
26        40        110        462
22        30        110        368
22        50        110        305
22        40        85        377
22        40        135        365
22        40        110        390
22        40        110        390
22        40        110        390
22        40        110        390
22        40        110        390
22        40        110        390
;
proc rsreg data=a;
  model rz =water rsp temp/lackfit;
run;
data grid;
do;
  rz= . ;
  rsp=40;
  do  water=18 to 28 by 2;
    do temp = 70 to 140 by 5;
             output;
    end;
  end;
end;
data grid;
set a grid;
proc rsreg data=grid out =predict noprint;
    model rz=water rsp temp/predict;
run;
data plot;
set predict；
if rsp=40;
proc g3d data=plot;
plot water*temp=rz / rotate=30 tilt=75 yticknum=6 xticknum=6 zticknum=9
zmin=500 zmax=100 ctop=red cbottom=blue caxis=black;
run;