FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;

/*
Dataset containing points to be plotted
 */

DATA points;
    DO rr = 1 TO 10 BY .01; /* risk ratios */     
        ARRAY p {*} p1-p10; /* p-values */
        ARRAY or {*} or1-or10; /* odds ratios for each p-value */
        DO j = 1 TO 10 BY 1;
            p(j) = 0.02*j;
            or(j) = (1/p(j) - 1)/(1/p(j) - rr)*rr;
            OUTPUT;
            END;
        END;
RUN;

/*
Dataset for annotation of diagonal
*/
    
DATA diagonal;
    FUNCTION='MOVE'; X=0; Y=0; XSYS='2'; YSYS='2';
    OUTPUT;
    FUNCTION='DRAW'; X=10; Y=10; XSYS='2'; YSYS='2';
    LINE=2;
    OUTPUT;
RUN;

/*
The plot
*/

PROC GPLOT DATA = points ANNOTATE = diagonal; 
    PLOT (or1-or10)*rr=1
        /OVERLAY HAXIS = AXIS1 VAXIS = AXIS2; 
    AXIS1
        MINOR = NONE
        ORDER = (0 TO 10 BY 1)
        LABEL = ('Risk ratio');
    AXIS2
        MINOR = NONE
        ORDER = (0 TO 10 BY 1)
        LABEL = (A=90 R=0 'Odds ratio');
    SYMBOL1 VALUE=NONE I=SPLINE LINE=1 COLOR=BLACK;
RUN;
QUIT;