topic Re: Save results from ESTIMATE statement to a file in SAS Programming

Save results from ESTIMATE statement to a file

PamG — Thu, 14 Mar 2024 20:48:46 GMT

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;

Re: Save results from ESTIMATE statement to a file

PamG — Thu, 14 Mar 2024 20:52:47 GMT

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;

Re: Save results from ESTIMATE statement to a file

PamG — Thu, 14 Mar 2024 20:54:28 GMT

I missed the CLASS statement.

PROC LOGISTIC;

STRATA=id,

CLASS race;

MODEL outcome(EVENT='1')=bmi bmi*race;

RUN;

Re: Save results from ESTIMATE statement to a file

yabwon — Thu, 14 Mar 2024 21:00:41 GMT

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

Re: Save results from ESTIMATE statement to a file

PamG — Thu, 14 Mar 2024 21:36:15 GMT

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!!!