topic Re: use proc surveylogistic to draw RCS in Statistical Procedures

use proc surveylogistic to draw RCS

qiqiqi — Tue, 02 Nov 2021 13:24:12 GMT

I want to use proc surveylogistic to draw a restrictive cubic spline curve (Y-axis is the odds ratio and x-axis is the exposure), but I don't know how to use the command.

At present, I can only draw a restrictive cubic spline curve with Y-axis as the possibility. Would like to Can someone tell me how to solve it? Thank you！

Re: use proc surveylogistic to draw RCS

SteveDenham — Wed, 03 Nov 2021 12:08:01 GMT

Not sure if I understand the question correctly, but have you considered using the EFFECT statement to spline your X axis values?

SteveDenham

Re: use proc surveylogistic to draw RCS

qiqiqi — Wed, 03 Nov 2021 12:18:00 GMT

Thank you very much for your answer!

I used the EFFECT statement. However, I don't know how to further export the OR and 95CI% for each exposure.

Here are the related commands I used:

proc surveylogistic data=analysis total=totals;

strata shengfen;

cluster site_ no;

weight weight_ quanzhong;

effect dmdurationsp=spline(dmduration/naturalcubic basis=tpf(noint) knotmethod=percentiles(5) details );

class DR2fen gender/param=ref;
model DR2fen (event="1")= gender dmdurationsp age /expb;
store DR2fenmodelsp;

output out=predicted_ values pred=pred lower=lower upper=upper;

run;

Re: use proc surveylogistic to draw RCS

SteveDenham — Fri, 05 Nov 2021 12:51:33 GMT

The LSMEANS statement with the ODDSRATIO option should give what you want.

SteveDenham

Re: use proc surveylogistic to draw RCS

qiqiqi — Fri, 05 Nov 2021 13:28:23 GMT

Thank you very much again! I have read the LSMEANS statement , but there are somethings that I still didn't quite understand how to use LSMEANS statement to get these parameters. If I want to use the model information, and draw the RCS graph with X-axis as the dmduration(the above variable) and y-axis as the corresponding OR and 95CI%(DR2fen, the above variable).