* Plot the reference limits;



%macro QPlot(dataset,label1,label2);

axis1 minor=none label=(angle=90 "QTc interval (msec)") 

    order=(420 to 480 by 20);

axis2 minor=none label=("Age (years)") order=(10 to 60 by 10);

symbol1 value=none color=black i=join line=1;

data temp;

    set &dataset;

    format estimate 3.0 x 3.0;

proc sort data=temp; 

    by x;

data annotate;

    length text $50;

    xsys="1"; ysys="1"; hsys="4"; x=50; y=16; position="5";  

    size=0.8; text="&label1"; function="label"; output;

    xsys="1"; ysys="1"; hsys="4"; x=50; y=8; position="5";  

    size=0.8; text="&label2"; function="label"; output;

proc gplot data=temp anno=annotate;

    plot estimate*x/frame haxis=axis2 vaxis=axis1;

    run;

    quit;

%mend QPlot;