Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- PROC LOGISTIC plots

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

☑ This topic is **solved**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 10-01-2023 05:10 PM
(599 views)

hello all,

I am running this PROC LOGISTIC:

proc logistic data=expi2 plots=effect ( clband=yes SHOWOBS=NO ) ;

model criterion (event= '1') = predictor1 predictor2;

the PLOTS=EFFECT option gives me the plot for predictor1, but not for predictor2.

Is there a way to obtain a plot with both predictors?

I have seen @Rick_SAS syntax for plotting different groups in the same graph, but I have not found a solution for what I need.

The datafile is attached.

Thanks in advance for any suggestions you may have.

Eman

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Producing a 3D plot for two continuous variables is very easy - in fact, it is the default of the EFFECTPLOT statement in that situation. The following produces a contour plot that should be pretty straightforward for any audience to understand.

```
proc logistic data=expi2;
model criterion (event= '1') = predictor1 predictor2 ;
effectplot;
run;
```

12 REPLIES 12

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Of course, I can use the store option and PLM, but that still gives me two separate plots, one for each predictor.

proc logistic data=expi2 plots=effect ( clband=yes SHOWOBS=NO ) ;

model criterion (event= '1') = predictor1 predictor2 ;

store pip;

proc plm restore=pip;

effectplot fit (x=predictor1 );

effectplot fit (x=predictor2 );

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

`proc logistic data=expi2 plots=effect ( clband=yes SHOWOBS=NO x=(predictor1 predictor2)) ;`

By the way, you did not attach the data and in fact, you should NOT attach data. Data should be provided as working SAS data step code (examples and instructions) and not via file attachments or screen captures.

--

Paige Miller

Paige Miller

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

hello.

sorry about the data.

What you suggest produces two separate graphs, one for predictor1 and one for predictor2.

The same results than with proc plm -- but your solution is of course more efficient.

What I would like is to have slopes for both predictors in one graph.

Eman

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I am not aware of any built-in way to have both EFFECT plots on one graph. It doesn't even make sense to do this unless PREDICTOR1 and PREDICTOR2 are on the exact same scale, otherwise you would need two different x-axes.

Of course, you can write your own program to put them on the same graph.

--

Paige Miller

Paige Miller

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

yes, they are on the exact same scale, that is why I would like to put them on the same graph.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

You need to realize that with 2 (or more) x-variables, the only way to do this is to hold X1 constant and then get prediction curve for the X2. And then repeat the process holding the other X2 constant and predict X1. Which may or may not be what you are hoping to get from such a plot. Maybe what you really want is a 3-D plot of probabilities with axes X1 and X2 and predicted probabilities, but I'd have to think of how you could create such a plot. I have done a similar thing using PROC REPORT, showing probabilities in a matrix of X1 and X2, but this can't be called a "plot".

--

Paige Miller

Paige Miller

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I will look into that, although a 3-D plot might be too complicated to read for my audience.

Yes, what you suggest is what I have in mind. I can reorganise my data and obtain the same through other means.

I was just curious if PROC LOGISTIC could handle it itself, but I guess this is not the case. Admittedly, it is not something that is likely to be useful very often, so I am not surprised that it does not exist as a built-in capability.

Thank you all for engaging in the discussion.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Producing a 3D plot for two continuous variables is very easy - in fact, it is the default of the EFFECTPLOT statement in that situation. The following produces a contour plot that should be pretty straightforward for any audience to understand.

```
proc logistic data=expi2;
model criterion (event= '1') = predictor1 predictor2 ;
effectplot;
run;
```

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

That is true, it is not that difficult. I will go with this!

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Sample of a PROC REPORT output for probabilities of two variables, colored (red=high prob, green=low prob) to make it more like a "plot". I have removed the column and row headers for confidentiality. Also, in the real PROC REPORT output, the numbers in the cells are large enough to be readable.

--

Paige Miller

Paige Miller

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

This looks awesome indeed.

I have never used PROC REPORT, but I will look into it.

Thank you!

Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.

**If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website. **

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.