Proc plm Scatterplot and line graph?

[TamaraTavares]

Hello,

I used proc mixed command to conduct linear mixed models to examine whether genetic status (mutation carrier vs. non-mutation carrier) is associated with increased ventricular volume (continous variable). I found a main effect of genetic status and another predictor variable (years to onset) and I am interested in plotting my data & the model in order to visualize it.  

I have ran the following to create a graph:

proc plm source=graph6;
title plot Genetic Status*Years From Symptom Onset; effectplot slicefit (x=YrsfromAV_AAO_Time1 sliceby=GS_Time1);
run;

I was wondering if anyone knows of a way to create this line graph that encompasses the model, in addition to a scatter plot indicating each data point overlaid on it.

Thank you in advance,

Tamara

[Rick_SAS]

The "item store" that PROC PLM reads does not contain observation information, just a summary of the model. So you can't overlay the fit on a scatter plot with PROC PLM. However, you can use the OUTPUT statement to save out the predictions and then use PROC SGPLOT to visualize the model. See the example code in 

https://communities.sas.com/t5/SAS-Statistical-Procedures/How-to-plot-regression-line-from-mixed-mod...