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

DATA melanoma; SET melanoma;
    fail=(dc ne 2);
    thick=thick-3;
    age=age/10-5;
    PROC SORT; BY years;
RUN;

proc iml;    
use melanoma;
read all var _num_ into temdat;
dummy=j(205,6,0);
dummy[,1]=temdat[,10]; /* years */
dummy[,2]=temdat[,11]; /* fail */
dummy[,3]=temdat[,6];  /* ulceration */
dummy[,4]=temdat[,7];  /* thickness */
dummy[,5]=temdat[,8];  /* sex */
dummy[,6]=temdat[,9];  /* age */
create mydata from dummy;
append from dummy;
quit;

proc iml;


level=0.05;
unit={"years"};
varlist={"ulceration", "thickness", "sex", "age"};
option={n,n,y,n,n};


****************************************************************************;

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


%additive(mydata,level,unit,varlist,option,dummy1,estimates,dummy2);
quit;

****************************************************************************;

PROC GPLOT DATA=estimates;
PLOT col2*col1 col12*col1 col17*col1/OVERLAY HAXIS=axis1 VAXIS=axis2;
axis1 ORDER=0 TO 10 BY 2
    MINOR=NONE LABEL=('Years');
axis2 ORDER=0 TO 0.4 BY 0.1
    MINOR=NONE LABEL=(A=90 R=0 'Cumulative baseline hazard');
symbol1 V=NONE I=STEPLJ L=1;
symbol2 V=none I=STEPLJ L=33;
symbol3 V=none I=STEPLJ L=33;
RUN;