Fit plot with observations GLIMMIX

kivester — Thu, 13 Apr 2017 20:20:56 GMT

I am using PROC GLIMMIX (SAS v9.3) to fit a mixed model in which I am most interested in the relationship between the outcome and a continuous covariate. I can get various graphical displays from PROC PLM, but I am having trouble finding anyway to plot the original observations on the effectplots. This is so easy in PROC GENMOD, I can't believe there isn't a way to acheive this with a model fit in GLIMMIX.

Thank you for any help!

Katy

Re: Fit plot with observations GLIMMIX

Rick_SAS — Thu, 13 Apr 2017 23:47:59 GMT

PROC PLM doesn't have accesss to the original data, which is why it behaves differently. PLM only has access to the model.

In fact, the documentation for the OBS option in the EFFECTPLOT statement states "OBS<(obs-options)>
displays observations on the effect plots. An input data set is required; hence the OBS option is not available with PROC PLM."

Just use an OUTPUT statement to create a SAS data set that contains the predicted values. There are several predictions and scales (BLUP, ILINK,...) to choose from. Then you can sort by the covariate and use PROC SGPLOT to overlay the SERIES plot of the predictors on a SCATTER plot of the observations.

topic Re: Fit plot with observations GLIMMIX in Statistical Procedures

Fit plot with observations GLIMMIX

Re: Fit plot with observations GLIMMIX