BookmarkSubscribeRSS Feed
boban
Calcite | Level 5

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

 

2 REPLIES 2
Rick_SAS
SAS Super FREQ

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;
boban
Calcite | Level 5

Thank you for the answer, very much appreciated!

 

I will try that.

 

Regards,

B

SAS INNOVATE 2024

Innovate_SAS_Blue.png

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. 

Register now!

How to Concatenate Values

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.

Get the $99 certification deal.jpg

 

 

Back in the Classroom!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 2 replies
  • 1681 views
  • 0 likes
  • 2 in conversation