plot figure for mixed model results

[GiaLee]

Hi all,

I have a question on creating a figure using linear mixed model results. The time variable is categorical, and the outcome is a score. My mentor asked me to draw one line in an unadjusted model and then draw one line while adjusting for several covariates. I can easily create the first figure using proc plm (figure 1). However, when I include covariates in the model and use proc plm, it automatically plots the figure and just adjust for effects with the variable after the "time" variable. Could someone assist me with this? Thank you!

unadjusted 

PROC MIXED DATA = a   METHOD = REML COVTEST;
CLASS record_id   time(ref="6")   younger   gender education;
MODEL Score = time   younger gender  education/ SOLUTION OUTpred=PREDDATA;
RANDOM INTERCEPT / SUBJECT = record_id;
repeated time_ / type=ar(1) sub=record_id;
store out=MixedModel;
RUN;
proc plm restore=MixedModel;
effectplot interaction (x=time_);
run;

adjusted 

 

[Ksharp]

Check @Rick_SAS blogs:
https://blogs.sas.com/content/iml/2018/12/19/visualize-mixed-model.html
https://blogs.sas.com/content/iml/2019/05/30/visualize-interaction-effects-regression.html

[GiaLee]

Thanks for the reply. These blogs are really helpful for plotting interaction terms. However, I'm looking for a way to plot the figure between time (category) and score, adjusting for several covariates but not the interaction terms. Thanks.