topic plotting estimates in proc genomd in Statistical Procedures

plotting estimates in proc genomd

ellen_qi — Sat, 07 Oct 2017 01:22:44 GMT

Hi guys,

I am struggling now to make a plot for the ratios and their confidence intervals from proc genmod.

The code is now like:

proc genmod data=fmdfmd;

class dummy_fmd_bed (ref='0');

model had_procedure (event='1')=age_at_dx dummy_fmd_bed dxtest_renal pmh_dissection pmh_aneurysm pres_hypertension pres_dissection_car_vert meds_no_of_bp_meds/alpha=0.1 dist=bin link=log;

estimate 'age_at_dx' age_at_dx -10/ exp alpha=0.1;
estimate 'dummy_fmd_bed' dummy_fmd_bed 1 -1/exp alpha=0.1;
estimate 'dxtest_renal' dxtest_renal 1 -1/exp alpha=0.1;
estimate 'pmh_dissection' pmh_dissection 1 -1/exp alpha=0.1;
estimate 'pmh_aneurysm' pmh_aneurysm 1 -1/exp alpha=0.1;
estimate 'pres_dissection_car_vert' pres_dissection_car_vert 1 -1/exp alpha=0.1;
estimate 'pres_hypertension' pres_hypertension 1 -1/exp alpha=0.1;
estimate 'meds_no_of_bp_meds' meds_no_of_bp_meds 1 /exp alpha=0.1;

run;

So these exponentiation estimate can give the ratios, but i want a plot of these ratios just like the odds ratio plots. Is there any way to do it?

Re: plotting estimates in proc genomd

Reeza — Sat, 07 Oct 2017 03:56:38 GMT

See the instructions about halfway down the page:

https://blogs.sas.com/content/iml/2015/07/29/or-plots-log-scale.html

@ellen_qi wrote:

Hi guys,

I am struggling now to make a plot for the ratios and their confidence intervals from proc genmod.

The code is now like:

proc genmod data=fmdfmd;

class dummy_fmd_bed (ref='0');

model had_procedure (event='1')=age_at_dx dummy_fmd_bed dxtest_renal pmh_dissection pmh_aneurysm pres_hypertension pres_dissection_car_vert meds_no_of_bp_meds/alpha=0.1 dist=bin link=log;

estimate 'age_at_dx' age_at_dx -10/ exp alpha=0.1;
estimate 'dummy_fmd_bed' dummy_fmd_bed 1 -1/exp alpha=0.1;
estimate 'dxtest_renal' dxtest_renal 1 -1/exp alpha=0.1;
estimate 'pmh_dissection' pmh_dissection 1 -1/exp alpha=0.1;
estimate 'pmh_aneurysm' pmh_aneurysm 1 -1/exp alpha=0.1;
estimate 'pres_dissection_car_vert' pres_dissection_car_vert 1 -1/exp alpha=0.1;
estimate 'pres_hypertension' pres_hypertension 1 -1/exp alpha=0.1;
estimate 'meds_no_of_bp_meds' meds_no_of_bp_meds 1 /exp alpha=0.1;

run;

So these exponentiation estimate can give the ratios, but i want a plot of these ratios just like the odds ratio plots. Is there any way to do it?

Re: plotting estimates in proc genomd

ellen_qi — Sat, 07 Oct 2017 13:20:31 GMT

I know how to get the odds ratio plot under logistic, but mine is under proc genmod, and this literature does not provide methods how to do it.

Re: plotting estimates in proc genomd

Reeza — Sat, 07 Oct 2017 17:10:26 GMT

I don't believe GENMOD supports producing a ODDS RATIO plot by default, therefore the reference to literature that shows you how to create it manually.

EDIT: Note that I've moved this post to the Statistical Procedures forum so someone with more experience can help.

Re: plotting estimates in proc genomd

ellen_qi — Sat, 07 Oct 2017 20:56:49 GMT

Thank you 🙂 First time here for help!