options nodate nonumber ps=55 ls=80;

title 'RPD -- Example 1.1 -- Y=Yield, X=Ozone';

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

proc print;
run;

proc plot;
plot yield*ozone;
run;

proc reg;
model yield=ozone / p r clm cli;
run;