topic Re: Data from proc phreg with frailty models in Statistical Procedures

Data from proc phreg with frailty models

Anton64 — Sat, 25 Jul 2020 14:50:37 GMT

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

Re: Data from proc phreg with frailty models

SteveDenham — Mon, 27 Jul 2020 11:38:44 GMT

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

SteveDenham

Re: Data from proc phreg with frailty models

Anton64 — Wed, 29 Jul 2020 08:37:58 GMT

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.

Re: Data from proc phreg with frailty models

SteveDenham — Mon, 03 Aug 2020 12:46:14 GMT

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.