options ls=85 ps=35; run;
data sasuser.mmex2_40;
input deg_days gas;
lines;
15.6 520
26.8 610
37.8 870
36.4 850
35.5 880
18.6 490
15.3 450
 7.9 250
 0.0 110
40.0   .  /* won't get used in fit, but will ahve predictions calculated */
;
run;
Proc means data=sasuser.mmex2_40;
run;

/* regression, plus graphs via "typewriter" */

Proc reg data=sasuser.mmex2_40;
 model gas = deg_days / CLI clm;
 plot gas*deg_days="*"
      predicted.*deg_days="0" / CLEAR COLLECT OVERLAY;
run;

/* regression, plus high resolution graphs */

Proc reg data=sasuser.mmex2_40 GRAPHICS;
 model gas = deg_days / CLI clm;
 plot gas*deg_days="*"
      predicted.*deg_days="0" / CLEAR COLLECT OVERLAY;
run;