I am using complex survey data and want to use restricted cubic splines to plot the change in the Odds Ratio over a Splined Predictor using PROC SURVEYLOGISTIC. I have tried using PROC PLM and the EFFECTPLOT statement, however, the Y axis displays "predicted probabilities", rather than displaying odds ratios. How do I use PROC SURVEYLOGISTIC to display splines that can easily be interpreted as odds ratios? Code that I have tried: TITLE "Splines of effect of age on opioid prescribing"; PROC SURVEYLOGISTIC DATA=work.ed2014_2021_v4; EFFECT spl_age = spline(age / DETAILS NATURALCUBIC BASIS=TPF KNOTMETHOD=PERCENTILELIST(5 25 50 75 95)); MODEL opioids1ormore (EVENT='1') = spl_age/NOINT CLPARM; CLUSTER cpsum; STRATA cstratm; WEIGHT patwt1k; STORE splineModel; OUTPUT OUT=SplineOut PREDICTED=fit; RUN; PROC PLM RESTORE=splineModel PLOTS=ALL NOCLPRINT; EFFECTPLOT FIT(X=age); ODS OUTPUT fitplot=fitplot; SCORE DATA=work.ed2014_2021_v4 OUT=plm_out PREDICTED=_predicted LCLM=lower UCLM=upper; RUN;
... View more