options nodate nonumber ps=55 ls=80;

title 'RPD -- Example 1.12 -- Rel. Risk-1 vs. Dust Exposure';

data rr_dust;
input dust rr @@;
rr_1 = rr-1;
datalines;
75 1.10 100 1.05 150 0.97 350 1.90 600 1.83 900 2.45 1300 3.70 1650 3.52
2250 4.16
;

proc print;
run;


symbol value=dot;

proc reg;
model rr_1 = dust / noint p r;   /* noint: remove intercept */
plot rr_1*dust / conf
                 haxis = 0 to 2400 by 400
                 vaxis = -0.5 to 4 by 0.5;
run;


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

proc gplot;
plot rr_1*dust / haxis = 0 to 2400 by 400
                 vaxis = -0.5 to 4 by 0.5;
run;