Yes, this blog post was extremely helpful to get me where I'm at now (that's where I got the code I pasted above). Given the multiple lines, it seems to me that the second figure displays the relationship between x (es) and y (gsi) across the different clusters (individuals) -- random effects. But I'm not interested in displaying random effects, just average effects of X on Y by levels of moderator. I'm looking for something very similar to the first figure I posted (i.e. predictions of y [GSI] by x [ES] grouped by arbitrary levels of the moderator [TIS] -- like mean-1SD and mean+1SD), but would like to edit axis labels, legends, line colors. I read in other posts to the community that it is a bit harder to do that in PROC PLM than in PROC SGPLOT, but I haven't been able to send output of my models to PROC SGPLOT. I've been struggling to get at that with the code and examples in the blog post you pasted. I even thought about creating line graphs in Excel based on coefficients in the output, but was unsure if in mixed-effects models the coefficient of X would be (b+d) as in OLS regression (Y = a + bX + cM + dX*M) given the random intercept in my model. Plus I am really interested in learning to do this in SAS given that I have to do similar graphs often in my research.
... View more