@sseb  What you're trying to estimate doesn't really make sense unless you have an interaction between sex and BMI (without an interaction the HR is the same for both males and females). If you include an interaction I think it would be something like:   estimate "&value. &year. sex level 1" bmis [-1, &ref.] [1, &value.] sex*bmis [-1, &ref. 1] [1, &value. 1]/ exp cl cov; estimate "&value. &year. sex level 2" bmis [-1, &ref.] [1, &value.] sex*bmis [-1, &ref. 2] [1, &value. 2]/ exp cl cov where sex is a class variable. This should be relatively easy to check since the interaction only matters for the non-reference level, so for the reference level of sex the two lines estimate "&value. &year. sex level 2 (ref)" bmis [-1, &ref.] [1, &value.] sex*bmis [-1, &ref. 2] [1, &value. 2]/ exp cl cov estimate "&value. &year." bmis [-1, &ref.] [1, &value.] / exp cl cov should produce identical results, so you can just change which level is the reference and use the simpler code. ... View more

Hi, I think something like this should do it: proc phreg data=temp;     effect bmiS = spline(bmi / basis=tpf(noint) NATURALCUBIC details knotmethod=list(18.0175 20.5289 22.4059 27.4693) );     model time*event(0)=bmiS  / rl=wald ties=EFRON ;     store sasuser.coxr; run; %macro est(ref=20, start=15, end=35, by=0.1); %Do i = 1 %To %eval(%SysFunc( Ceil( %SysEvalF( ( &End - &Start ) / &By ) ) ) +1) ;    %Let value=%SysEvalF( ( &Start - &By ) + ( &By * &I ) ) ;     estimate "&value." bmis [-1, &ref] [1, &value] / exp cl;     %end; %mend est; ods html select none; ods rtf select none; ods dataset Estimates=Estimates; proc plm restore=sasuser.coxr;     %est(ref=20, start=16, end=35, by=0.01); run; ods html select all; ods rtf select all; data estimates;     set estimates;     BMI=label*1; run; proc sgplot data=estimates NOAUTOLEGEND Noborder;     Series y=ExpEstimate x=BMI / ;     Series y=LowerExp x=BMI / LINEATTRS=(pattern=ShortDash color=Black THICKNESS=1);     Series y=UpperExp x=BMI / LINEATTRS=(pattern=ShortDash color=Black THICKNESS=1);     REFLINE 1 / axis=y;     REFLINE 20 / Axis=X LINEATTRS=(pattern=ThinDot color=Black THICKNESS=1);;     yaxis Values=(0.8 1 1.5 2 3 4 6 8 10) Label="Hazard ratio" Type=LOG LABELATTRS=(weight=BOLD);     xaxis min=16 VALUES=(16 to 35 by 1) LABELATTRS=(weight=BOLD); run; ... View more