Is there a way to save results from ESTIMATE statement to a file? I have a product of 2 multi-category variables and would like to plot the estimates on a graph. Attaching a code for sample data.
data Neuralgia;
input Treatment $ Sex $ Age Duration Pain $ @@;
datalines;
P F 68 1 No B M 74 16 No P F 67 30 No
P M 66 26 Yes B F 67 28 No B F 77 16 No
A F 71 12 No B F 72 50 No B F 76 9 Yes
A M 71 17 Yes A F 63 27 No A F 69 18 Yes
B F 66 12 No A M 62 42 No P F 64 1 Yes
A F 64 17 No P M 74 4 No A F 72 25 No
P M 70 1 Yes B M 66 19 No B M 59 29 No
A F 64 30 No A M 70 28 No A M 69 1 No
B F 78 1 No P M 83 1 Yes B F 69 42 No
B M 75 30 Yes P M 77 29 Yes P F 79 20 Yes
A M 70 12 No A F 69 12 No B F 65 14 No
B M 70 1 No B M 67 23 No A M 76 25 Yes
P M 78 12 Yes B M 77 1 Yes B F 69 24 No
P M 66 4 Yes P F 65 29 No P M 60 26 Yes
A M 78 15 Yes B M 75 21 Yes A F 67 11 No
P F 72 27 No P F 70 13 Yes A M 75 6 Yes
B F 65 7 No P F 68 27 Yes P M 68 11 Yes
P M 67 17 Yes B M 70 22 No A M 65 15 No
P F 67 1 Yes A M 67 10 No P F 72 11 Yes
A F 74 1 No B M 80 21 Yes A F 69 3 No
;
proc logistic data=Neuralgia;
class Treatment Sex /param=ref;
model Pain= Treatment Sex Age ;
estimate 'Pairwise B vs P' Treatment 0 1 / exp CL;
estimate 'Pairwise A vs B' Treatment 1 -1 / exp CL;
estimate 'Female vs Male' Sex 1 / exp CL;
run;
how about adding ODS output statement:
ods output Estimates=work.E;
proc logistic data=...
...
run;
before proc logistic, and then to use data from WORK.E dataset?
Bart
I can't edit my original post here are more details.
I can't use effectplot. I am doing a conditonal logit and have interaction effect between a time-invariant and time varying covariate. Hence I need to use the ESTIMATE to calculate the effect sizes.
PROC LOGISTIC;
STRATA=id,
MODEL outcome(EVENT='1')=bmi bmi*race;
RUN;
I missed the CLASS statement.
PROC LOGISTIC;
STRATA=id,
CLASS race;
MODEL outcome(EVENT='1')=bmi bmi*race;
RUN;
how about adding ODS output statement:
ods output Estimates=work.E;
proc logistic data=...
...
run;
before proc logistic, and then to use data from WORK.E dataset?
Bart
Perfect! This works beautifully. I could not find any documentation on this. I guess I was not looking in the right place. Thank you so so much!!!
Join us for SAS Innovate 2025, our biggest and most exciting global event of the year, in Orlando, FL, from May 6-9. Sign up by March 14 for just $795.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
Ready to level-up your skills? Choose your own adventure.