Here is a worked out example based on @StatDave advice, including how to flip the graph
proc logistic data=sashelp.heart;
where BP_status ne "Optimal";
effect spl = spline(weight / details naturalcubic basis=tpf(noint)
knotmethod=percentiles(5) );
model BP_status(event="High") = spl / selection=none alpha=0.1;
units weight=10;
oddsratio weight / at(weight=110 to 210 by 2) cl=pl;
ods output ORPlot=orp;
quit;
data orpGraph;
set orp;
Weight = input(scan(DisplayLabel, 2, "="), best.);
label weight="Weight (lbs)";
run;
proc sgplot data=orpGraph noautolegend;
band x=weight upper=upperCLdisplay lower=lowerCLdisplay;
series x=weight y=OddsRatioEstDisplay;
yaxis label="Odds Ratio of High BP for 10 lbs Weight Increase";
run;
... View more
