options nodate nonumber ps=55 ls=80;

title 'RPD -- Examples 1.1+ -- Soybean Yield vs. Ozone Level';

data yld_oz;
input ozone yield @@;   
   /* read variable values successively, all the way to the end of each line */
datalines;
0.02 242 0.07 237 0.11 231 0.15 201
;

proc print;
run;


symbol value=x;   /* data points will be plotted as "x" */

proc reg;
model yield = ozone / p r;   
   /* p: output predicted values;  r: perform analysis of residuals */
plot yield*ozone / haxis = 0.0 to 0.16 by 0.02 
                   vaxis = 190 to 260 by 10;
run;