Data from proc phreg with frailty models

[Anton64]

Hi,

I have to generate a graph showing the cumulative hazard function coming from the following model:

proc phreg data=out.DS ;

     class ARM PARTICIP;

     model end_et*endptfl(0) = ARM BL_NH3 / rl ;

               ID PARTICIP ;

               RANDOM PARTICIP / DIST= GAMMA ;

     hazardratio 'Frailty Model Analysis' ARM;

     run ;

Having specified the RANDOM statement, any ASSESS, BASELINE, and OUTPUT statements are ignored, so I don’t know how I can obtain the data to generate the graph.

thanks for any help

Alessandro Antonellini

[SteveDenham]

Make sure ODS graphics are enabled and add PLOTS=CUMHAZ to the PROC PHREG statement.

SteveDenham

[Anton64]

Hi Steve,

if I apply your suggestions, the graph is still not generated and the log is returning the note

The PLOTS= option (PROC statement) is ignored for a frailty model analysis.

[SteveDenham]

Well, that shoots that down @Anton64 , Does your code produce a table with the hazard ratios?  If so, you could use an ODS OUTPUT statement to get those into a dataset, and then use PROC SGPLOT for the graph.

SteveDenham.