Hi, I would appreciate a help with creating a caterpillar plot from the generalized linear mixed (binomial distribution) models with random intercepts output. I came across with an article on creating a caterpillar plot from proc mcmc output using the cater macro:
ods graphics on;
%CATER(data=postout, var=delta: );
ods graphics off;
I want to do a similar plot but with the proc glimmix output (SolutionR=random effect estimates
and plotting the ‘wiskers’ as +/- 1.96 x the standard errors).
Great. Glad you solved your problem. If you have other questions, please post them. Otherwise, please mark the thread as solved.
If you post your code, we can be more helpful. What statistic are you trying the plot? The parameter estimates and their standard errors?
If so, see this article about regression coefficient plots.
As shown in the article, the basic idea is to output the information you want to plot (estimates and stderrs) to a SAS data set. then use a DATA step to create two variables: Lower=Estimate-1.96*StdErr and Upper=Estimate-+.96*StdErr.
Use the SCATTER statement in PROC SGPLOT and the XERRRORLOW= and XERRORUPPER= options to display the 1.96*StdErr intervals.
Other articles that cerate similar plots are
"Coverage probability of confidence intervals"
and
Thank you Rick! Your response is extremely helpful.
Great. Glad you solved your problem. If you have other questions, please post them. Otherwise, please mark the thread as solved.
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
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.