turn on suggestions

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

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- SAS/GRAPH and ODS Graphics
- /
- Graph probabilities in proc GLIMMIX

Topic Options

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

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

2 weeks ago - last edited 2 weeks ago

Hi,

I am trying to graph the probabilities of an event happening in proc GLIMMIX. Does anyone know if the GLIMMIX procedure generates this type of graph?

Basically, I want to plot an S curve, as one does in logistic regressions.

Regards,

B

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

2 weeks ago

I don't think GLIMMIX generates a FitPlot directly, but you can use the OUTPUT stmt to create an output data set that has the predicted values, then use PROC SGPLOT to create the graph. The trick is that you need to specify what sort of predicted values you want from (BLUP or NOBLUP) and (ILINK or NOILINK). See the GLIMMIX doc for what these terms.

Here is code for an "S-like curve" that you can study and modify:

```
proc glimmix data=sashelp.cars(where=(origin^='Europe'));
model Origin = mpg_city / dist=binary;
output out=modelout pred(noblup ilink)=prob
lcl(noblup ilink)=lcl
ucl(noblup ilink)=ucl;
run;
proc sort data=modelout;
by mpg_city;
run;
proc sgplot data=modelout;
band x=mpg_city lower=lcl upper=ucl;
series x=mpg_city y=prob;
yaxis min=0 max=1;
refline 0 1 / axis=y;
run;
```

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

2 weeks ago

Thank you for the answer, very much appreciated!

I will try that.

Regards,

B