options nodate nonumber ps=55 ls=80;

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

data yld_oz;
input ozone yield @@;
datalines;
0.02 242 0.07 237 0.11 231 0.15 201
;


symbol value=x;

proc reg;
model yield = ozone / p r clb clm cli;
   /* clb: CIs for betas;  clm: CIs for mean responses;  cli: PIs */
plot yield*ozone / conf      /* add confidence limits for mean */
                   haxis = 0.0 to 0.16 by 0.02
                   vaxis = 170 to 290 by 10;
run;


/* To produce a graph with nicer confidence limits ... */

symbol interpol=rlclm95;   /* regression, linear, ... */
                           /* confidence limits (for) mean, ... */
                           /* 95% (limits) */

proc gplot;
plot yield*ozone / haxis = 0.0 to 0.16 by 0.02
                   vaxis = 170 to 290 by 10;
run;