data cars93a;
 set ma2612.cars93a;
run;
data newcars;
 HIGHMPG=.;
 DISPLACE=2.0;
 HP=110;
 RPM=5000;
run;
proc append base=cars93a data=newcars;
run;
data cars93a;
 set cars93a;
 D2=DISPLACE**2;
 H2=HP**2;
 R2=RPM**2;
run;
/*
And now run the regression with the clm (estimate of mean response)
and cli (prediction) options.
*/
proc reg data=cars93a;
 model highmpg=displace hp rpm d2 h2 r2/clb clm cli alpha=.05;
output out=outpred p=predicted lclm=lclm uclm=uclm lcl=lcl ucl=ucl;
run;
/*
Combine data sets to display regressor and response
along with estimation and confidence limits and 
prediction and prediction limits.
*/
data both;
 merge cars93a outpred;
run;
proc print;run;