Hi Mrite,   I am doing the same analysis as yours. The hazard analysis shows that the association between Bicarb (LBXSC3SI) and survival is not linear. Thus, I want to plot restricted cubic splines in proportional hazard regression.    I used the following code   PROC PHREG DATA = survey; effect LBXSC3SIs = spline(LBXSC3SI / basis=tpf(noint) NATURALCUBIC details knotmethod=rangefractions(0.05 0.50 0.95) ); model years*dead (0) = LBXSC3SIs / rl; run; ods graphics on;   However, I do not see any plot/graph in the SAS output. Did you find any solution for your code?     ... View more