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. ... View more