GENERATING PREDICTIONS AND HAZARD FUNCTIONS

After fitting a model with PROC LIFEREG, it's sometimes desirable to generate predicted survival times for the observations in the data set. If you want a single point estimate for each individual, the predicted median survival time is probably the best. You can get this easily with the OUTPUT statement, as shown in the following example:

PROC LIFEREG DATA=recid;
   MODEL week*arrest(0)=fin age race wexp mar paro prio
         / D=WEIBULL;
   OUTPUT OUT=a P=median STD=s;
RUN;
PROC PRINT DATA=a;
   VAR week arrest _prob_ medians;
RUN;