06-16-2016 01:19 PM
I want to display estimates from an HLM model for two different groups (variables: out= my outcome, ID = case identifyer, age = measure of time, diagnosis = grouping variable).
class = ID;
model out = age diagnosis / solution outp=pred;
random= intercept /SUBJECT = ID TYP = UN GCORR;
proc sgplot data=pred;
series y=pred x=age /group = diagnosis;
Resulting plot has lines for individual cases and I'm trying to create a graph with one line per group. Added bonus if you can help with overlaying raw data within same graph and/or adding confidence intervals. Thanks in advance!
06-16-2016 01:39 PM
proc sort data=pred; by diagnosis age; run;
before the call to proc sgplot.
scatter x=age y=out / group=diagnosis;
before run; to get the raw data points.