Wanting to plot restricted cubic splines in proportional hazards regression

[giants619]

I would like to create a cubic spline plot for proportional hazards regression. It should look something like this:

I have tried the following code, which did not produce any graphics:

ods graphics on;
proc icphreg data=updated2;
model time*mortstat(0) = RDW / basehaz=splines;
hazardratio RDW;

run;

ods graphics on;
PROC PHREG DATA = updated2;
effect RDWs = spline(RDW / basis=tpf(noint) NATURALCUBIC details knotmethod=rangefractions(0.05 0.50 0.95) );
model TIME*MORTSTAT(0) = RDWs / rl;
run;

I saw that several other individuals asked the same question on these forums years ago but did not get any response.

[sbxkoenk]

Hello @giants619 ,

Both the PROC ICPHREG and PROC PHREG have a PLOTS= option on the PROC statement.

If you do not use it, you cannot get any plots.

PLOTS= Controls the plots that are produced through ODS Graphics.

Check the documentation 

https://go.documentation.sas.com/doc/en/helpcenterwlcm/1.0/home.htm

of both procedures and add the PLOTS= option.

Cheers,

Koen