topic Re: How do I output the Influence statistics in proc logistic in a dataset in Statistical Procedures

How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 14:19:41 GMT

Re: How do I output the Influence statistics in proc logistic in a dataset

StatDave — Tue, 16 Aug 2022 14:37:12 GMT

All of those statistics can be requested in the OUTPUT statement which will put them into a data set. For example, this will save the the DFBeta values in data set MyInf:

output data=MyInf dfbeta=_all_;

See the description of the OUTPUT statement in the LOGISTIC documentation for details.

Or you can save all of them into a data set with this statement:

ods output influence=MyInf;

See this note on saving results from any procedure in a data set.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 14:41:23 GMT

Thanks for your reply but that only gives me output for the variables in my model. What I am looking for is the Influence statistics generated from the influence option in proc logistic. I want the output shown on page 5938 Chapter 76: The LOGISTIC Procedure 15.1.

Re: How do I output the Influence statistics in proc logistic in a dataset

StatDave — Tue, 16 Aug 2022 14:46:59 GMT

The ODS OUTPUT statement I showed saves all of the influence statistics like how the INFLUENCE option displays them.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 15:27:35 GMT

Thanks so much for your help thus far. I am getting output but not like the output shown on page 5938 Chapter 76: The LOGISTIC Procedure 15.1. That's the output I need for my analysis.

Re: How do I output the Influence statistics in proc logistic in a dataset

StatDave — Tue, 16 Aug 2022 15:33:16 GMT

I'm afraid a page number does not help me locate what you are looking at. Please refer to the title of the example or section of the LOGISTIC documentation you are referring to.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 15:43:25 GMT

Hi,
The example is in the Logistic Procedure 15.1. The output example of the Influence statistics in on page # 5938, the example shown is Output 76.6.2 Regression Diagnostics from the INFLUENCE Option.

Re: How do I output the Influence statistics in proc logistic in a dataset

StatDave — Tue, 16 Aug 2022 15:54:12 GMT

That is the example titled "Logistic Regression Diagnostics" which is what I assumed you were looking at and is why I suggested using the ODS statement: ods output influence=MyInf; . The resulting data set (MyInf) has all of the same columns of values as the displayed output that you refer to. So, I don't understand how that doesn't provide what you want.

Re: How do I output the Influence statistics in proc logistic in a dataset

Reeza — Tue, 16 Aug 2022 16:10:12 GMT

Show the code used please.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 16:58:16 GMT

Hi StatDave,

I really appreciate you going the extra mile to help. For some reason I am not getting that output in the output dataset as shown in the Logistic Procedure 15.1 on page 5938 Output 76.6.2 Regression Diagnostics from the INFLUENCE Option. Once again I appreciate your help.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 17:04:09 GMT

Hi,

Here is the code example.

ods output Influence = myinf; 'Not getting same output in an output dataset as shown in Logistic Procedure 15.1 page 5938';
ods graphics on;
proc logistic data = xxx.xxx ;
model purch_ind (event = '1') = 'Not Real Variables'
a
b
c
d
e
f
g
/ influence aggregate scale = none;
/*output out = MyInf dfbeta =_all_;*/
run;
ods graphics off;

Re: How do I output the Influence statistics in proc logistic in a dataset

StatDave — Tue, 16 Aug 2022 17:07:25 GMT

Here is the full example in the documentation with the ODS OUTPUT statement added and a PROC PRINT step which shows the contents of the created data set. Try that and you should get the same basic display from PROC PRINT as from PROC LOGISTIC.

data vaso;
   length Response $12;
   input Volume Rate Response @@;
   LogVolume=log(Volume);
   LogRate=log(Rate);
   datalines;
3.70  0.825  constrict       3.50  1.09   constrict
1.25  2.50   constrict       0.75  1.50   constrict
0.80  3.20   constrict       0.70  3.50   constrict
0.60  0.75   no_constrict    1.10  1.70   no_constrict
0.90  0.75   no_constrict    0.90  0.45   no_constrict
0.80  0.57   no_constrict    0.55  2.75   no_constrict
0.60  3.00   no_constrict    1.40  2.33   constrict
0.75  3.75   constrict       2.30  1.64   constrict
3.20  1.60   constrict       0.85  1.415  constrict
1.70  1.06   no_constrict    1.80  1.80   constrict
0.40  2.00   no_constrict    0.95  1.36   no_constrict
1.35  1.35   no_constrict    1.50  1.36   no_constrict
1.60  1.78   constrict       0.60  1.50   no_constrict
1.80  1.50   constrict       0.95  1.90   no_constrict
1.90  0.95   constrict       1.60  0.40   no_constrict
2.70  0.75   constrict       2.35  0.03   no_constrict
1.10  1.83   no_constrict    1.10  2.20   constrict
1.20  2.00   constrict       0.80  3.33   constrict
0.95  1.90   no_constrict    0.75  1.90   no_constrict
1.30  1.625  constrict
;
proc logistic data=vaso;
   model Response=LogRate LogVolume/influence;
   ods output influence=in;
run;
proc print label; run;

Re: How do I output the Influence statistics in proc logistic in a dataset

Reeza — Tue, 16 Aug 2022 17:13:25 GMT

Your reference numbers (15.1) and page numbers don't align as we don't have whatever book you're using. Can you link to the reference in the online documentation.

https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_logistic_toc.htm

Or whatever version of SAS/STAT you have.

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 18:09:05 GMT

Hi Reeza,

Thanks for your willingness to help, but I think I got it now resolved based on StatDave's post. Regarding your reply it is the SAS/STAT 15.1 User's Guide The Logistic Procedure downloaded from the SAS procedures. The page is 5938, the example on that page is titled "Output 76.6.2 Regression Diagnostics from the INFLUENCE Option".

Re: How do I output the Influence statistics in proc logistic in a dataset

timothymbanks — Tue, 16 Aug 2022 18:11:47 GMT

Hi StatDave,
Thanks once again for your help. It worked finally!!!. Don't know why it was not outputting before but I do get the results now. Thanks a million,

Timothy